431 research outputs found
Probing Extreme-Density Matter with Gravitational Wave Observations of Binary Neutron Star Merger Remnants
We present a proof-of-concept study, based on numerical-relativity
simulations, of how gravitational waves (GWs) from neutron star merger remnants
can probe the nature of matter at extreme densities. Phase transitions and
extra degrees of freedom can emerge at densities beyond those reached during
the inspiral, and typically result in a softening of the equation of state
(EOS). We show that such physical effects change the qualitative dynamics of
the remnant evolution, but they are not identifiable as a signature in the GW
frequency, with the exception of possible black-hole formation effects. The EOS
softening is, instead, encoded in the GW luminosity and phase and is in
principle detectable up to distances of the order of several Mpcs with advanced
detectors and up to hundreds of Mpcs with third generation detectors. Probing
extreme-density matter will require going beyond the current paradigm and
developing a more holistic strategy for modeling and analyzing postmerger GW
signals.Comment: 5 pages, 3 figures. Matches version accepted on ApJ
Dynamical Mass Ejection from Binary Neutron Star Mergers
We present fully general-relativistic simulations of binary neutron star
mergers with a temperature and composition dependent nuclear equation of state.
We study the dynamical mass ejection from both quasi-circular and
dynamical-capture eccentric mergers. We systematically vary the level of our
treatment of the microphysics to isolate the effects of neutrino cooling and
heating and we compute the nucleosynthetic yields of the ejecta. We find that
eccentric binaries can eject significantly more material than quasi-circular
binaries and generate bright infrared and radio emission. In all our
simulations the outflow is composed of a combination of tidally- and
shock-driven ejecta, mostly distributed over a broad angle from
the orbital plane, and, to a lesser extent, by thermally driven winds at high
latitudes. Ejecta from eccentric mergers are typically more neutron rich than
those of quasi-circular mergers. We find neutrino cooling and heating to
affect, quantitatively and qualitatively, composition, morphology, and total
mass of the outflows. This is also reflected in the infrared and radio
signatures of the binary. The final nucleosynthetic yields of the ejecta are
robust and insensitive to input physics or merger type in the regions of the
second and third r-process peaks. The yields for elements on the first peak
vary between our simulations, but none of our models is able to explain the
Solar abundances of first-peak elements without invoking additional first-peak
contributions from either neutrino and viscously-driven winds operating on
longer timescales after the mergers, or from core-collapse supernovae.Comment: 19 pages, 10 figures. We corrected a problem in the formulation of
the neutrino heating scheme and re-ran all of the affected models. The main
conclusions are unchanged. This version also contains one more figure and a
number of improvements on the tex
Signatures of hypermassive neutron star lifetimes on r-process nucleosynthesis in the disk ejecta from neutron star mergers
We investigate the nucleosynthesis of heavy elements in the winds ejected by
accretion disks formed in neutron star mergers. We compute the element
formation in disk outflows from hypermassive neutron star (HMNS) remnants of
variable lifetime, including the effect of angular momentum transport in the
disk evolution. We employ long-term axisymmetric hydrodynamic disk simulations
to model the ejecta, and compute r-process nucleosynthesis with tracer
particles using a nuclear reaction network containing species. We
find that the previously known strong correlation between HMNS lifetime,
ejected mass, and average electron fraction in the outflow is directly related
to the amount of neutrino irradiation on the disk, which dominates mass
ejection at early times in the form of a neutrino-driven wind. Production of
lanthanides and actinides saturates at short HMNS lifetimes ( ms),
with additional ejecta contributing to a blue optical kilonova component for
longer-lived HMNSs. We find good agreement between the abundances from the disk
outflow alone and the solar r-process distribution only for short HMNS
lifetimes ( ms). For longer lifetimes, the rare-earth and third
r-process peaks are significantly under-produced compared to the solar pattern,
requiring additional contributions from the dynamical ejecta. The
nucleosynthesis signature from a spinning black hole (BH) can only overlap with
that from a HMNS of moderate lifetime ( ms). Finally, we show that
angular momentum transport not only contributes with a late-time outflow
component, but that it also enhances the neutrino-driven component by moving
material to shallower regions of the gravitational potential, in addition to
providing additional heating.Comment: 18 pages, 11 figures, published version with small change
General Relativistic Three-Dimensional Multi-Group Neutrino Radiation-Hydrodynamics Simulations of Core-Collapse Supernovae
We report on a set of long-term general-relativistic three-dimensional (3D)
multi-group (energy-dependent) neutrino-radiation hydrodynamics simulations of
core-collapse supernovae. We employ a full 3D two-moment scheme with the local
M1 closure, three neutrino species, and 12 energy groups per species. With
this, we follow the post-core-bounce evolution of the core of a nonrotating
- progenitor in full unconstrained 3D and in octant symmetry for
. We find the development of an asymmetric runaway
explosion in our unconstrained simulation. We test the resolution dependence of
our results and, in agreement with previous work, find that low resolution
artificially aids explosion and leads to an earlier runaway expansion of the
shock. At low resolution, the octant and full 3D dynamics are qualitatively
very similar, but at high resolution, only the full 3D simulation exhibits the
onset of explosion.Comment: Accepted to Ap
R-process Nucleosynthesis from Three-Dimensional Magnetorotational Core-Collapse Supernovae
We investigate r-process nucleosynthesis in three-dimensional (3D)
general-relativistic magnetohydrodynamic simulations of rapidly rotating
strongly magnetized core collapse. The simulations include a microphysical
finite-temperature equation of state and a leakage scheme that captures the
overall energetics and lepton number exchange due to postbounce neutrino
emission and absorption. We track the composition of the ejected material using
the nuclear reaction network SkyNet. Our results show that the 3D dynamics of
magnetorotational core-collapse supernovae (CCSN) are important for their
nucleosynthetic signature. We find that production of r-process material beyond
the second peak is reduced by a factor of 100 when the magnetorotational jets
produced by the rapidly rotating core undergo a kink instability. Our results
indicate that 3D magnetorotationally powered CCSNe are a robust r-process
source only if they are obtained by the collapse of cores with unrealistically
large precollapse magnetic fields of order G. Additionally, a
comparison simulation that we restrict to axisymmetry, results in overly
optimistic r-process production for lower magnetic field strengths.Comment: 10 pages, 9 figures, 2 tables. submitted to Ap
Low mass binary neutron star mergers : gravitational waves and neutrino emission
Neutron star mergers are among the most promising sources of gravitational
waves for advanced ground-based detectors. These mergers are also expected to
power bright electromagnetic signals, in the form of short gamma-ray bursts,
infrared/optical transients, and radio emission. Simulations of these mergers
with fully general relativistic codes are critical to understand the merger and
post-merger gravitational wave signals and their neutrinos and electromagnetic
counterparts. In this paper, we employ the SpEC code to simulate the merger of
low-mass neutron star binaries (two neutron stars) for a set of
three nuclear-theory based, finite temperature equations of state. We show that
the frequency peaks of the post-merger gravitational wave signal are in good
agreement with predictions obtained from simulations using a simpler treatment
of gravity. We find, however, that only the fundamental mode of the remnant is
excited for long periods of time: emission at the secondary peaks is damped on
a millisecond timescale in the simulated binaries. For such low-mass systems,
the remnant is a massive neutron star which, depending on the equation of
state, is either permanently stable or long-lived. We observe strong
excitations of l=2, m=2 modes, both in the massive neutron star and in the form
of hot, shocked tidal arms in the surrounding accretion torus. We estimate the
neutrino emission of the remnant using a neutrino leakage scheme and, in one
case, compare these results with a gray two-moment neutrino transport scheme.
We confirm the complex geometry of the neutrino emission, also observed in
previous simulations with neutrino leakage, and show explicitly the presence of
important differences in the neutrino luminosity, disk composition, and outflow
properties between the neutrino leakage and transport schemes.Comment: Accepted by PRD; 23 pages; 24 figures; 4 table
Neutrino-heated winds from millisecond protomagnetars as sources of the weak r-process
We explore heavy element nucleosynthesis in neutrino-driven winds from rapidly rotating,
stronglymagnetized protoneutron stars (‘millisecond protomagnetars’) forwhich themagnetic
dipole is aligned with the rotation axis, and the field is assumed to be a static force-free configuration.
We process the protomagnetar wind trajectories calculated by Vlasov, Metzger &
Thompson through the r-process nuclear reaction network SkyNet using contemporary models
for the evolution of the wind electron fraction during the protoneutron star cooling phase.
Although we do not find a successful second or third-peak r-process for any rotation period
P, we show that protomagnetars with P ∼ 1–5 ms produce heavy element abundance distributions
that extend to higher nuclear mass number than from otherwise equivalent spherical
winds (with the mass fractions of some elements enhanced by factors of �100–1000). The
heaviest elements are synthesized by outflows emerging along flux tubes that graze the closed
zone and pass near the equatorial plane outside the light cylinder. Due to dependence of the
nucleosynthesis pattern on the magnetic field strength and rotation rate of the protoneutron
star, natural variations in these quantities between core collapse events could contribute to the
observed diversity of the abundances of weak r-process nuclei in metal-poor stars. Further
diversity, including possibly even a successful third-peak r-process, could be achieved for
misaligned rotators with non-zero magnetic inclination with respect to the rotation axis. If
protomagnetars are central engines for GRBs, their relativistic jets should contain a high-mass
fraction of heavy nuclei of characteristic mass number ¯A ≈ 100, providing a possible source
for ultrahigh energy cosmic rays comprised of heavy nuclei with an energy spectrum that
extends beyond the nominal Grezin–Zatsepin–Kuzmin cut-off for protons or iron nuclei
Post-merger evolution of a neutron star-black hole binary with neutrino transport
We present a first simulation of the post-merger evolution of a black
hole-neutron star binary in full general relativity using an energy-integrated
general relativistic truncated moment formalism for neutrino transport. We
describe our implementation of the moment formalism and important tests of our
code, before studying the formation phase of a disk after a black hole-neutron
star merger. We use as initial data an existing general relativistic simulation
of the merger of a neutron star of 1.4 solar mass with a black hole of 7 solar
mass and dimensionless spin a/M=0.8. Comparing with a simpler leakage scheme
for the treatment of the neutrinos, we find noticeable differences in the
neutron to proton ratio in and around the disk, and in the neutrino luminosity.
We find that the electron neutrino luminosity is much lower in the transport
simulations, and that the remnant is less neutron-rich. The spatial
distribution of the neutrinos is significantly affected by relativistic
effects. Over the short timescale evolved, we do not observe purely
neutrino-driven outflows. However, a small amount of material (3e-4Msun) is
ejected in the polar region during the circularization of the disk. Most of
that material is ejected early in the formation of the disk, and is fairly
neutron rich. Through r-process nucleosynthesis, that material should produce
high-opacity lanthanides in the polar region, and could thus affect the
lightcurve of radioactively powered electromagnetic transients. We also show
that by the end of the simulation, while the bulk of the disk is neutron-rich,
its outer layers have a higher electron fraction. As that material would be the
first to be unbound by disk outflows on longer timescales, the changes in Ye
experienced during the formation of the disk could have an impact on the
nucleosynthesis outputs from neutrino-driven and viscously-driven outflows.
[Abridged]Comment: 29 pages, 25 figure
- …