In this work, we study the merger of two neutron stars with a gravitational mass of 1.4Mmiddot in circle each, employing the Shen-Horowitz-Teige equation of state. This equation of state is a corner ...case, allowing the formation of a stable neutron star with the given total baryonic mass of 3.03Mmiddot in circle. We investigate in unprecedented detail the structure of the remnant, in particular the mass distribution, the thermal structure, and the rotation profile. We also compute fluid trajectories both inside the remnant and those relevant for the formation of the disk. We find a peanut-shaped fluid flow inside the remnant following a strong m=2 perturbation. Moreover, the flow is locally compressive, causing the appearance of dynamic hot spots. Further, we introduce new diagnostic measures that are easy to implement in numeric simulations and that allow one to quantify mass and compactness of merger remnants in a well-defined way. As in previous studies of supra- and hypermassive stars, we find a remnant with a slowly rotating core and an outer envelope rotating at nearly Keplerian velocity. We compute a Tolman-Oppenheimer-Volkoff star model which agrees well with that of the remnant in the core, while the latter possesses extensive outer layers rotating close to Kepler velocity. Finally, we extract the gravitational wave signal and discuss the detectability with modern observatories. This study has implications for the interpretation of gravitational wave detections from the postmerger phase and is relevant for short gamma-ray burst models.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
We present general relativistic numerical simulations of binary neutron star (BNS) mergers with different initial spin configurations. We focus on models with stars of mass 1.4 M⊙ each, which employ ...the equation of state (EOS) by Shen, Horowitz, and Teige, and which result in stable NSs as merger remnants. For comparison, we consider two irrotational equal mass (M=1.35 M⊙) and unequal mass (M=1.29, 1.42 M⊙) BNS models using the APR4 EOS, which result in a supramassive merger remnant. We present visualizations of the fluid flow and temperature distribution and find a strong impact of the spin on vortex structure and nonaxisymmetric deformation. We compute the radial mass distribution and the rotation profile in the equatorial plane using recently developed measures independent of spatial gauge, revealing slowly rotating cores that can be well approximated by the cores of spherical stars. We also study the influence of the spin on the inspiral phase and the gravitational wave (GW) signal. Using a newly developed analysis method, we further show that gravitational waveforms from BNS mergers can exhibit one or more phase jumps after merger, which occur together with minima of the strain amplitude. We provide a natural explanation in terms of the remnant’s quadrupole moment, and show that cancellation effects due to phase jumps can have a strong impact on the GW power spectrum. Finally, we discuss the impact of the spin on the amount of ejected matter.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
Recently, general relations among the quadrupole moment (Q), the moment of inertia (I) and the tidal deformability (Love number) of a neutron star were shown to exist. They are nearly independent of ...the nuclear matter equation of state and would be of great aid in extracting parameters from observed gravitational waves and in testing general relativity. These relations, however, do not account for strong magnetic fields. We consider this problem by studying the effect of a strong magnetic field on slowly rotating relativistic neutron stars and show that, for simple magnetic field configurations that are purely poloidal or purely toroidal, the relation between Q and I is again nearly universal. However, different magnetic field geometries lead to different I-Q relations, and, in the case of a more realistic twisted-torus magnetic field configuration, the relation depends significantly on the equation of state, losing its universality. I-Love-Q relations must thus be used with very great care, since universality is lost for stars with long spin periods, i.e. P 10 s, and strong magnetic fields, i.e. B 1012 G.
We present new results of fully general relativistic magnetohydrodynamic simulations of binary neutron star (BNS) mergers performed with the Whisky code. All the models use a piecewise polytropic ...approximation of the APR4 equation of state for cold matter, together with a 'hybrid' part to incorporate thermal effects during the evolution. We consider both equal and unequal-mass models, with total masses such that either a supramassive NS or a black hole is formed after merger. Each model is evolved with and without a magnetic field initially confined to the stellar interior. We present the different gravitational wave (GW) signals as well as a detailed description of the matter dynamics (magnetic field evolution, ejected mass, post-merger remnant/disk properties). Our simulations provide new insights into BNS mergers, the associated GW emission and the possible connection with the engine of short gamma-ray bursts (both in the 'standard' and in the 'time-reversal' scenarios) and other electromagnetic counterparts.
ABSTRACT
Multimessenger astronomy received a great boost following the discovery of kilonova (KN) AT2017gfo, the optical counterpart of the gravitational wave source GW170817 associated with the ...short gamma-ray burst GRB 170817A. AT2017gfo was the first KN that could be extensively monitored in time using both photometry and spectroscopy. Previously, only few candidates have been observed against the glare of short GRB afterglows. In this work, we aim to search the fingerprints of AT2017gfo-like KN emissions in the optical/NIR light curves of 39 short GRBs with known redshift. For the first time, our results allow us to study separately the range of luminosity of the blue and red components of AT2017gfo-like kilonovae in short GRBs. In particular, the red component is similar in luminosity to AT2017gfo, while the blue KN can be more than 10 times brighter. Finally, we exclude a KN as luminous as AT2017gfo in GRBs 050509B and 061201.
The new era of multimessenger astrophysics requires the capability of studying different aspects of the evolution of compact objects. In particular, the merger of neutron star binaries is a strong ...source of gravitational waves and electromagnetic radiation, from radio to γ-rays, as demonstrated by the detection of GW170817 and its electromagnetic counterparts. In order to understand the physical mechanisms involved in such systems, it is necessary to employ fully general relativistic magnetohydrodynamic (GRMHD) simulations able to include the effects of a composition and temperature dependent equation of state describing neutron star matter as well as neutrino emission and reabsorption. Here, we present our new code named Spritz that solves the GRMHD equations in 3D Cartesian coordinates and on a dynamical spacetime. The code can support tabulated equations of state, taking into account finite temperature effects and allowing for the inclusion of neutrino radiation. In this first paper, we present the general features of the code and a series of tests performed in special and general relativity to assess the robustness of the basic GRMHD algorithms implemented. Among these tests, we also present the first comparison between a non-staggered and a staggered formulation of the vector potential evolution, which is used to guarantee the divergence-less character of the magnetic field. With respect to other publicly available GRMHD codes, Spritz combines the robust approach of a staggered formulation of the vector potential together with the use of an equation of state driver (EOS_Omni) that can allow the code to use finite temperature equations of state. A next version of the code will fully test the EOS_Omni driver by coupling it with a neutrino leakage scheme.
The merger of two neutron stars is predicted to give rise to three major detectable phenomena: a short burst of γ-rays, a gravitational-wave signal, and a transient optical-near-infrared source ...powered by the synthesis of large amounts of very heavy elements via rapid neutron capture (the r-process). Such transients, named 'macronovae' or 'kilonovae', are believed to be centres of production of rare elements such as gold and platinum. The most compelling evidence so far for a kilonova was a very faint near-infrared rebrightening in the afterglow of a short γ-ray burst at redshift z = 0.356, although findings indicating bluer events have been reported. Here we report the spectral identification and describe the physical properties of a bright kilonova associated with the gravitational-wave source GW170817 and γ-ray burst GRB 170817A associated with a galaxy at a distance of 40 megaparsecs from Earth. Using a series of spectra from ground-based observatories covering the wavelength range from the ultraviolet to the near-infrared, we find that the kilonova is characterized by rapidly expanding ejecta with spectral features similar to those predicted by current models. The ejecta is optically thick early on, with a velocity of about 0.2 times light speed, and reaches a radius of about 50 astronomical units in only 1.5 days. As the ejecta expands, broad absorption-like lines appear on the spectral continuum, indicating atomic species produced by nucleosynthesis that occurs in the post-merger fast-moving dynamical ejecta and in two slower (0.05 times light speed) wind regions. Comparison with spectral models suggests that the merger ejected 0.03 to 0.05 solar masses of material, including high-opacity lanthanides.
Full text
Available for:
IJS, KISLJ, NUK, SBMB, UL, UM, UPUK
We construct general relativistic models of stationary, strongly magnetized neutron stars. The magnetic field configuration, obtained by solving the relativistic Grad–Shafranov equation, is a ...generalization of the twisted-torus model recently proposed in the literature; the stellar deformations induced by the magnetic field are computed by solving the perturbed Einstein's equations; stellar matter is modelled using realistic equations of state. We find that in these configurations the poloidal field dominates over the toroidal field and that, if the magnetic field is sufficiently strong during the first phases of the stellar life, it can produce large deformations.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, IZUM, KILJ, NLZOH, NUK, OILJ, PILJ, PNG, SAZU, SBCE, SBMB, UL, UM, UPUK
We derive the luminosity function φ(L) and redshift distribution Ψ(z) of short gamma-ray bursts (SGRBs) using all the available observer-frame constraints (i.e. peak flux, fluence, peak energy and ...duration distributions) of the large population of Fermi SGRBs and the rest-frame properties of a complete sample of SGRBs detected by Swift. We show that a steep φ(L) ∝ L− α with α ≥ 2.0 is excluded if the full set of constraints is considered. We implement a Markov chain Monte Carlo method to derive the φ(L) and Ψ(z) functions assuming intrinsic Ep−Liso and Ep−Eiso correlations to hold or, alternatively, that the distributions of intrinsic peak energy, luminosity, and duration are independent. To make our results independent from assumptions on the progenitor (NS−NS binary mergers or other channels) and from uncertainties on the star formation history, we assume a parametric form for the redshift distribution of the population of SGRBs. We find that a relatively flat luminosity function with slope ~0.5 below a characteristic break luminosity ~3 × 1052 erg s-1 and a redshift distribution of SGRBs peaking at z ~ 1.5−2 satisfy all our constraints. These results also hold if no Ep−Liso and Ep−Eiso correlations are assumed and they do not depend on the choice of the minimum luminosity of the SGRB population. We estimate, within ~200 Mpc (i.e. the design aLIGO range for the detection of gravitational waves produced by NS−NS merger events), that there should be 0.007−0.03 SGRBs yr-1 detectable as γ-ray events. Assuming current estimates of NS−NS merger rates and that all NS−NS mergers lead to a SGRB event, we derive a conservative estimate of the average opening angle of SGRBs ⟨ θjet ⟩ ~ 3°−6°. The luminosity function implies a prompt emission average luminosity ⟨L⟩ ~ 1.5 × 1052 erg s-1, higher by nearly two orders of magnitude than previous findings in the literature, which greatly enhances the chance of observing SGRB “orphan” afterglows. Effort should go in the direction of finding and identifying such orphan afterglows as counterparts of GW events.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK