No Arabic abstract
Jets can become collimated as they propagate through dense environments and understanding such interactions is crucial for linking physical models of the environments to observations. In this work, we use 3D special-relativistic simulations to study how jets propagate through the environment created around a neutron star merger remnant by neutrino-driven winds. We simulate four jets with two different initial structures, top-hat and Gaussian, and two luminosities. After jet breakout, we study the angular jet structures and the resulting afterglow light curves. We find that the initial angular structures are efficiently washed out during the propagation, despite the small wind mass of only $sim 10^{-3}$ M$_odot$. The final structures depend, however, on the jet luminosity, as less energetic jets are more strongly collimated. Although entrainment of baryons leads to only moderate outflow Lorentz factors ($approx 40$), all simulated jets can well reproduce the afterglow observed in the aftermath of GW170817. The inferred physical parameters (e.g. inclination angle, ambient particle number density), however, vary substantially between the fits and appear to be sensitive to smaller details of the angular jet shape, indicating that observationally inferred parameters may depend sensitively on the employed jet models.
Very long base interferometry (VLBI) radio images recently proved to be essential in breaking the degeneracy in the ejecta model for the neutron star merger event GW170817. We discuss the properties of synthetic radio images of merger jet afterglow by using semi-analytic models of laterally spreading or non-spreading jets. The image centroid initially moves away from the explosion point in the sky with an apparent superlumianal velocity. After reaching a maximum displacement its motion is reversed. This behavior is in line with that found in full hydrodynamics simulations. Since the evolution of the centroid shift and jet image size are significantly different in the two jet models, observations of these characteristics for very bright events might be able to confirm or constrain the lateral expansion law of merger jets. We explicitly demonstrate how $theta_{rm obs}$ is obtained by the centroid shift of radio images or its apparent velocity provided the ratio of the jet core size $theta_{c}$ and the viewing angle $theta_{rm obs}$ is determined by afterglow light curves. We show that a simple method based on a point-source approximation provides reasonable angular estimates ($10-20%$ errors at most). By taking a sample of structured Gaussian jet results, we find that the model with $theta_{rm obs} sim 0.32$ rad can explain the main features of the GW170817 afterglow light curves and the radio images.
The first neutron star-neutron star (NS-NS) merger was discovered on August 17, 2017 through gravitational waves (GW170817) and followed with electromagnetic observations. This merger was detected in an old elliptical galaxy with no recent star formation. We perform a suite of numerical calculations to understand the formation mechanism of this merger. We probe three leading formation mechanisms of double compact objects: classical isolated binary star evolution, dynamical evolution in globular clusters and nuclear cluster formation to test whether they are likely to produce NS-NS mergers in old host galaxies. Our simulations with optimistic assumptions show current NS-NS merger rates at the level of 10^-2 yr^-1 from binary stars, 5 x 10^-5 yr^-1 from globular clusters and 10^-5 yr^-1 from nuclear clusters for all local elliptical galaxies (within 100 Mpc^3). These models are thus in tension with the detection of GW170817 with an observed rate 1.5 yr^-1 (per 100 Mpc^3; LIGO/Virgo estimate). Our results imply that either (i) the detection of GW170817 by LIGO/Virgo at their current sensitivity in an elliptical galaxy is a statistical coincidence; or that (ii) physics in at least one of our three models is incomplete in the context of the evolution of stars that can form NS-NS mergers; or that (iii) another very efficient (unknown) formation channel with a long delay time between star formation and merger is at play.
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.
VLBI and JVLA observations revealed that GW170817 involved a narrow jet ($ theta_j approx 4^circ $) that dominated the afterglow peak at our viewing angle, $ theta_{rm obs} approx 20^circ $. This implies that at the time of the afterglow peak, the observed signal behaved like an afterglow of a top-hat jet seen at $ theta_{rm obs} gg theta_j $, and it can be modeled by analytic expressions that describe such jets. We use a set of numerical simulations to calibrate these analytic relations and obtain generic equations for the peak time and flux of such an afterglow as seen from various observing angles. Using the calibrated equations and the estimated parameters of GW170817, we estimate the detectability of afterglows from future double neutron star mergers during the Advanced LIGO/Virgo observation run O3. GW170817 took place at a relatively low-density environment. Afterglows of similar events will be detectable only at small viewing angles, $ theta_{rm obs} lesssim 20^circ $, and only $sim 20% $ of the GW detections of these events will be accompanied by a detectable afterglow. At higher densities, more typical to sGRB sites, up to $ 70% $ of the GW detections are expected to be followed by a detectable afterglow, typically at $ theta_{rm obs} sim 30^circ $. We also provide the latest time one should expect an afterglow detection. We find that for typical parameters, if the jet emission had not been detected within about a year after the merger, it is unlikely to be ever detected.
Although the main features of the evolution of binary neutron star systems are now well established, many details are still subject to debate, especially regarding the post-merger phase. In particular, the lifetime of the hyper-massive neutron stars formed after the merger is very hard to predict. In this work, we provide a simple analytic relation for the lifetime of the merger remnant as function of the initial mass of the neutron stars. This relation results from a joint fit of data from observational evidence and from various numerical simulations. In this way, a large range of collapse times, physical effects and equation of states is covered. Finally, we apply the relation to the gravitational wave event GW170817 to constrain the equation of state of dense matter.