ABSTRACT
We present a new moment-based energy-integrated neutrino transport code for neutron star merger simulations in general relativity. In the merger context, ours is the first code to include ...Doppler effects at all orders in υ/c, retaining all non-linear neutrino–matter coupling terms. The code is validated with a stringent series of tests. We show that the inclusion of full neutrino–matter coupling terms is necessary to correctly capture the trapping of neutrinos in relativistically moving media, such as in differentially rotating merger remnants. We perform preliminary simulations proving the robustness of the scheme in simulating ab-initio mergers to black hole collapse and long-term neutron star remnants up to ${\sim }70\,$ ms. The latter is the longest dynamical space-time, 3D, general relativistic simulations with full neutrino transport to date. We compare results obtained at different resolutions and using two different closures for the moment scheme. We do not find evidences of significant out-of-thermodynamic equilibrium effects, such as bulk viscosity, on the post-merger dynamics or gravitational wave emission. Neutrino luminosities and average energies are in good agreement with theory expectations and previous simulations by other groups using similar schemes. We compare dynamical and early wind ejecta properties obtained with M1 and with our older neutrino treatment. We find that the M1 results have systematically larger proton fractions. However, the differences in the nucleosynthesis yields are modest. This work sets the basis for future detailed studies spanning a wider set of neutrino reactions, binaries, and equations of state.
We investigate r-process nucleosynthesis in 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 robust r-process sources only if they are obtained by the collapse of cores with unrealistically large precollapse magnetic fields of the order of 1013 G. Additionally, a comparison simulation that we restrict to axisymmetry results in overly optimistic r-process production for lower magnetic field strengths.
We investigate the influence of magnetic fields on the evolution of binary neutron star (BNS) merger remnants via three-dimensional (3D) dynamical-spacetime general-relativistic magnetohydrodynamic ...(MHD) simulations. We evolve a post-merger remnant with an initial poloidal magnetic field, resolve the magnetoturbulence driven by shear flows, and include a microphysical finite-temperature equation of state. A neutrino leakage scheme that captures the overall energetics and lepton number exchange is also included. We find that turbulence induced by the magnetorotational instability in the hypermassive neutron star (HMNS) amplifies magnetic field to beyond magnetar strength (1015 G). The ultra-strong toroidal field is able to launch a relativistic jet from the HMNS. We also find a magnetized wind that ejects neutron-rich material with a rate of . The total ejecta mass in our simulation is 5 × 10−3 M . This makes the ejecta from the HMNS an important component in BNS mergers and a promising source of r-process elements that can power a kilonova. The jet from the HMNS reaches a terminal Lorentz factor of ∼5 in our highest-resolution simulation. The formation of this jet is aided by neutrino cooling preventing the accretion disk from protruding into the polar region. As neutrino pair-annihilation and radiative processes in the jet (which were not included in the simulations) will boost the Lorentz factor in the jet further, our simulations demonstrate that magnetars formed in BNS mergers are a viable engine for short gamma-ray bursts.
As gravitational-wave detectors become more sensitive and broaden their frequency bandwidth, we will access a greater variety of signals emitted by compact binary systems, shedding light on their ...astrophysical origin and environment. A key physical effect that can distinguish among different formation scenarios is the misalignment of the spins with the orbital angular momentum, causing the spins and the binary's orbital plane to precess. To accurately model such precessing signals, especially when masses and spins vary in the wide astrophysical range, it is crucial to include multipoles beyond the dominant quadrupole. Here, we develop the first multipolar precessing waveform model in the effective-one-body (EOB) formalism for the entire coalescence stage (i.e., inspiral, merger and ringdown) of binary black holes: SEOBNRv4PHM. In the nonprecessing limit, the model reduces to SEOBNRv4HM, which was calibrated to numerical-relativity (NR) simulations, and waveforms from black-hole perturbation theory. We validate SEOBNRv4PHM by comparing it to the public catalog of 1405 precessing NR waveforms of the Simulating eXtreme Spacetimes (SXS) collaboration, and also to 118 SXS precessing NR waveforms, produced as part of this project, which span mass ratios 1-4 and (dimensionless) black-hole's spins up to 0.9. We stress that SEOBNRv4PHM is not calibrated to NR simulations in the precessing sector. We compute the unfaithfulness against the 1523 SXS precessing NR waveforms, and find that, for 94% (57%) of the cases, the maximum value, in the total mass range 20 − 200 M⊙, is below 3% (1%). Those numbers change to 83% (20%) when using the inspiral-merger-ringdown, multipolar, precessing phenomenological model IMRPhenomPv3HM. We investigate the impact of such unfaithfulness values with two Bayesian, parameter-estimation studies on synthetic signals. We also compute the unfaithfulness between those waveform models as a function of the mass and spin parameters to identify in which part of the parameter space they differ the most. We validate them also against the multipolar, precessing NR surrogate model NRSur7dq4, and find that the SEOBNRv4PHM model outperforms IMRPhenomPv3HM.
Neutron star mergers are among the most promising sources of gravitational waves for advanced ground-based detectors. Simulations of these mergers with fully general relativistic codes are critical ...to understand the merger and postmerger gravitational wave signals and their neutrinos and electromagnetic counterparts. In this paper, we employ the Spectral Einstein Code to simulate the merger of low mass neutron star binaries for a set of three nuclear-theory-based, finite temperature equations of state. We show that the frequency peaks of the postmerger gravitational wave signal are in good agreement with predictions obtained from recent simulations using a simpler treatment of gravity. 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.
ABSTRACT 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 380 ms. 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.
We present ENIGMA, a time domain, inspiral-merger-ringdown waveform model that describes nonspinning binary black holes systems that evolve on moderately eccentric orbits. The inspiral evolution is ...described using a consistent combination of post-Newtonian theory, self-force and black hole perturbation theory. Assuming eccentric binaries that circularize prior to coalescence, we smoothly match the eccentric inspiral with a stand-alone, quasicircular merger, which is constructed using machine learning algorithms that are trained with quasicircular numerical relativity waveforms. We show that ENIGMA reproduces with excellent accuracy the dynamics of quasicircular compact binaries. We validate ENIGMA using a set of Einstein Toolkit eccentric numerical relativity waveforms, which describe eccentric binary black hole mergers with mass-ratios between 1≤q≤5.5, and eccentricities e0≲0.2 ten orbits before merger. We use this model to explore in detail the physics that can be extracted with moderately eccentric, nonspinning binary black hole mergers. In particular, we use ENIGMA to show that the gravitational wave transients GW150914, GW151226, GW170104, GW170814 and GW170608 can be effectively recovered with spinning, quasicircular templates if the eccentricity of these events at a gravitational wave frequency of 10 Hz satisfies e0≤{0.175,0.125,0.175,0.175,0.125}, respectively. We show that if these systems have eccentricities e0∼0.1 at a gravitational wave frequency of 10 Hz, they can be misclassified as quasicircular binaries due to parameter space degeneracies between eccentricity and spin corrections. Using our catalog of eccentric numerical relativity simulations, we discuss the importance of including higher-order waveform multipoles in gravitational wave searches of eccentric binary black hole mergers.
ABSTRACT We present results from high-resolution semiglobal simulations of neutrino-driven convection in core-collapse supernovae. We employ an idealized setup with parameterized neutrino ...heating/cooling and nuclear dissociation at the shock front. We study the internal dynamics of neutrino-driven convection and its role in redistributing energy and momentum through the gain region. We find that even if buoyant plumes are able to locally transfer heat up to the shock, convection is not able to create a net positive energy flux and overcome the downward transport of energy from the accretion flow. Turbulent convection does, however, provide a significant effective pressure support to the accretion flow as it favors the accumulation of energy, mass, and momentum in the gain region. We derive an approximate equation that is able to explain and predict the shock evolution in terms of integrals of quantities such as the turbulent pressure in the gain region or the effects of nonradial motion of the fluid. We use this relation as a way to quantify the role of turbulence in the dynamics of the accretion shock. Finally, we investigate the effects of grid resolution, which we change by a factor of 20 between the lowest and highest resolution. Our results show that the shallow slopes of the turbulent kinetic energy spectra reported in previous studies are a numerical artifact. Kolmogorov scaling is progressively recovered as the resolution is increased.
We present results from a controlled numerical experiment investigating the effect of stellar density gas on the coalescence of binary black holes (BBHs) and the resulting gravitational waves (GWs). ...This investigation is motivated by the proposed stellar core fragmentation scenario for BBH formation and the associated possibility of an electromagnetic counterpart to a BBH GW event. We employ full numerical relativity coupled with general-relativistic hydrodynamics and set up a 30+30 M_{⊙} BBH (motivated by GW150914) inside gas with realistic stellar densities. Our results show that at densities ρ≳10^{6}-10^{7} g cm^{-3} dynamical friction between the BHs and gas changes the coalescence dynamics and the GW signal in an unmistakable way. We show that for GW150914, LIGO observations appear to rule out BBH coalescence inside stellar gas of ρ≳10^{7} g cm^{-3}. Typical densities in the collapsing cores of massive stars are in excess of this density. This excludes the fragmentation scenario for the formation of GW150914.