Neutron Star Mergers Research Articles

Thermal Conductivity and Thermal Hall Effect in Dense Electron-Ion Plasma

In this study, we examine thermal conductivity and the thermal Hall effect in electron-ion plasmas relevant to hot neutron stars, white dwarfs, and binary neutron star mergers, focusing on densities found in the outer crusts of neutron stars and the interiors of white dwarfs. We consider plasma consisting of single species of ions, which could be either iron [56]Fe or carbon [12]C nuclei. The temperature range explored is from the melting temperature of the solid T∼109 K up to 1011 K. This covers both degenerate and non-degenerate electron regimes. We find that thermal conductivity increases with density and temperature for which we provide analytical scaling relations valid in different regimes. The impact of magnetic fields on thermal conductivity is also analyzed, showing anisotropy in low-density regions and the presence of the thermal Hall effect characterized by the Righi–Leduc coefficient. The transition froma degenerate to non-degenerate regime is characterized by a minimum ratio of thermal conductivity to temperature, which is analogous to the minimum observed already in the case of electrical conductivity. We provide also formulas fit to our numerical results, which can be used in dissipative magneto-hydrodynamics simulations of warm compact stars.

Binary neutron star mergers in massive scalar-tensor theory: Properties of postmerger remnants

We investigate the properties of postmerger remnants of binary neutron star mergers in the framework of Damour–Esposito-Farese-type scalar-tensor theory of gravity with a massive scalar field by numerical relativity simulation. It is found that the threshold mass for prompt collapse is raised in the presence of the excited scalar field. Our simulation results also suggest the existence of a long-lived ϕ mode in hypermassive neutron stars due to the presence of the massive scalar field that enhances the quasiradial oscillation in the remnant. We investigate the descalarization condition in hypermassive neutron stars and discover a distinctive signature in postmerger gravitational waves. Published by the American Physical Society 2024

Incidence of afterglow plateaus in gamma-ray bursts associated with binary neutron star mergers

One of the most surprising gamma-ray burst (GRB) features discovered with the Swift X-ray telescope (XRT) is a plateau phase in the early X-ray afterglow light curves. These plateaus are observed in the majority of long GRBs, while their incidence in short GRBs (SGRBs) is still uncertain due to their fainter X-ray afterglow luminosity with respect to long GRBs. An accurate estimate of the fraction of SGRBs with plateaus is of utmost relevance given the implications that the plateau may have for our understanding of the jet structure and possibly of the nature of the binary neutron star (BNS) merger remnant. This work presents the results of an extensive data analysis of the largest and most up-to-date sample of SGRBs observed with the XRT, and for which the redshift has been measured. We find a plateau incidence of 18-37<!PCT!> in SGRBs, which is a significantly lower fraction than that measured in long GRBs (>50<!PCT!>). Although still debated, the plateau phase could be explained as energy injection from the spin-down power of a newly born magnetized neutron star (NS; magnetar). We show that this scenario can nicely reproduce the observed short GRB (SGRBs) plateaus, while at the same time providing a natural explanation for the different plateau fractions between short and long GRBs. In particular, our findings may imply that only a minority of BNS mergers generating SGRBs leave behind a sufficiently stable or long-lived NS to form a plateau. From the probability distribution of the BNS remnant mass, a fraction 18-37<!PCT!> of short GRB plateaus implies a maximum NS mass in the range $ odot

Simulating Short Gamma-Ray Burst Jets in Realistic Late Binary Neutron Star Merger Environments

The electromagnetic emission and the afterglow observations of the binary neutron star merger event GW170817A confirmed the association of the merger with a short gamma-ray burst (GRB) harboring a narrow (5°–10°) and powerful (1049–1050 erg) jet. Using the 1 s long neutrino-radiation general relativistic MHD simulation of coalescing neutron stars of K. Kiuchi et al., and following the semi-analytical estimates of M. Pais et al., we inject a narrow, powerful, unmagnetized jet into the post-merger phase. We explore different opening angles, luminosities, central engine durations, and times after the merger. We explore early (0.1 s following the merger) and late (1 s) jet launches; the latter is consistent with the time delay of ≈1.74 s observed between GW170817 and GRB 170817A. We demonstrate that the semi-analytical estimates correctly predict the jets’ breakout and collimation conditions. When comparing our synthetic afterglow light curves to the observed radio data of GW170807, we find a good agreement for a 3 × 1049 erg jet launched late with an opening angle in the range ≃5°–7°.

Magnetic dissipation in short gamma-ray-burst jets

Aims. Short gamma-ray bursts originate when relativistic jets emerge from the remnants of binary neutron star (BNS) mergers, as observed in the first multi-messenger event GW170817–GRB 170817A, which coincided with a gravitational wave signal. Both the jet and the remnant are believed to be magnetized, and the presence of magnetic fields is known to influence the jet propagation across the surrounding post-merger environment. In the magnetic interplay between the jet and the environment itself, effects due to a finite plasma conductivity may be important, especially in the first phases of jet propagation. We aim to investigate these effects, from jet launching to its final breakout from the post-merger environment. Methods. We used the PLUTO numerical code to perform 2D axisymmetric and full 3D resistive relativistic magnetohydrodynamic (MHD) simulations, employing spherical coordinates with spatial radial stretching. We considered and compared different models for physical resistivity, which must be small but still dominating over the intrinsic numerical dissipation (which yields unwanted smearing of structures in any ideal MHD code). Stiff terms in the current density are treated with IMplicit-EXplicit Runge Kutta algorithms for time-stepping. A Synge-like gas (Taub equation of state) is also considered. All simulations were performed using an axisymmetric analytical model for both the jet propagation environment and the jet injection; we leave the case of jet propagation in a realistic environment (i.e., imported from an actual BNS merger simulation) to a future study. Results. As expected, no qualitative differences are detected due to the effect of a finite conductivity, but significant quantitative differences in the jet structure and induced turbulence are clearly seen in 2D axisymmetric simulations, and we also compare different resistivity models. We see the formation of regions with a resistive electric field parallel to the magnetic field, and nonthermal particle acceleration may be enhanced there. The level of dissipated Ohmic power is also dependent on the various recipes for resistivity. Most of the differences arise before the breakout from the inner environment, whereas once the jet enters the external weakly magnetized environment region, these differences are preserved during further propagation despite the lower grid refinement. Finally, we show and discuss the 3D evolution of the jet within the same environment in order to highlight the emergence of nonaxisymmetric features.

Swiftly Chasing Gravitational Waves across the Sky in Real Time

We introduce a new capability of the Neil Gehrels Swift Observatory, dubbed “continuous commanding,” that achieves 10 s latency response time on orbit to unscheduled target-of-opportunity requests received on the ground. We show that this will allow Swift to respond to premerger (early-warning) gravitational-wave (GW) detections, rapidly slewing the Burst Alert Telescope (BAT) across the sky to place the GW origin in the BAT field of view at or before merger time. This will dramatically increase the GW/gamma-ray burst (GRB) codetection rate and enable prompt arcminute localization of a neutron star merger. We simulate the full Swift response to a GW early-warning alert, including input sky maps produced at different early-warning times, a complete model of the Swift attitude control system, and a full accounting of the latency between the GW detectors and the spacecraft. 60 s of early warning can double the rate of a prompt GRB detection with arcminute localization, and 140 s guarantees observation anywhere on the unocculted sky, even with localization areas ≫1000 deg2. While 140 s is beyond current GW detector sensitivities, 30–70 s is achievable today. We show that the detection yield is now limited by the latency of LIGO/Virgo cyberinfrastructure and motivate a focus on its reduction. Continuous commanding has been integrated as a general capability of Swift, significantly increasing its versatility in response to the growing demands of time-domain astrophysics. We demonstrate this potential on an externally triggered fast radio burst (FRB), slewing 81° across the sky, and collecting X-ray and UV photons from the source position <150 s after the trigger was received from the Canadian Hydrogen Intensity Mapping Experiment, thereby setting the earliest and deepest such constraints on high-energy activity from nonrepeating FRBs. The Swift Team invites the community to consider and propose novel scientific applications of ultra-low-latency UV, X-ray, and gamma-ray observations.

A census of the Sun’s ancestors and their contributions to the Solar System chemical composition

In this work, we compute the rates and numbers of different types of stars and phenomena (supernovae, novae, white dwarfs, merging neutron stars, black holes) that contributed to the chemical composition of the Solar System. During the Big Bang, only light elements formed, while all the heavy ones, from carbon to uranium and beyond, have since been created inside stars. Stars die and release the newly formed elements into the interstellar gas. This process is called ‘chemical evolution’. In particular, we analyse the death rates of stars of all masses, whether they die quiescently or explosively. These rates and total star numbers are computed in the context of a revised version of the two-infall model for the chemical evolution of the Milky Way, which reproduces the observed abundance patterns of several chemical species, the global solar metallicity, and the current gas, stellar, and total surface mass densities relatively well. We also compute the total number of stars ever born and still alive as well as the number of stars born up to the formation of the Solar System with mass and metallicity like those of the Sun. This latter number accounts for all the possible existing Solar systems that can host life in the solar vicinity. We conclude that, among all the stars (from 0.8 to 100 M⊙) that were born and died from the Big Bang up until the Solar System formation epoch and that contributed to its chemical composition, 93.00% were stars that died as single white dwarfs (without interacting significantly with a companion star) and originated in the mass range of 0.8–8 M⊙, while 5.24% were neutron stars and 0.73% were black holes, both originating from core-collapse supernovae (M &gt; 8 M⊙); 0.64% were Type Ia supernovae and 0.40% were nova systems, both originating from the same mass range as the white dwarfs. The number of stars similar to the Sun born from the Big Bang up until the formation of the Solar System, with metallicity in the range 12+log(Fe/H)= 7.50 ± 0.04 dex, is ~31•107, and in particular our Sun is the ~2.61• 107-th star of this kind.

Energy-dependent and Energy-integrated Two-moment General-relativistic Neutrino Transport Simulations of a Hypermassive Neutron Star

We compare two-moment-based energy-dependent and three variants of energy-integrated neutrino transport general-relativistic magnetohydrodynamics simulations of a hypermassive neutron star. To study the impacts due to the choice of the neutrino transport schemes, we perform simulations with the same setups and input neutrino microphysics. We show that the main differences between energy-dependent and energy-integrated neutrino transport are found in the disk and ejecta properties, as well as in the neutrino signals. The properties of the disk surrounding the neutron star and the ejecta in energy-dependent transport are very different from the ones obtained using energy-integrated schemes. Specifically, in the energy-dependent case, the disk is more neutron-rich at early times and becomes geometrically thicker at later times. In addition, the ejecta is more massive and, on average, more neutron-rich in the energy-dependent simulations. Moreover, the average neutrino energies and luminosities are about 30% higher. Energy-dependent neutrino transport is necessary if one wants to better model the neutrino signals and matter outflows from neutron star merger remnants via numerical simulations.

Lower-mass-gap Black Holes in Dense Star Clusters

The existence of compact stellar remnants in the mass range 2–5 M ⊙ has long been debated. This so-called lower-mass gap (LMG) was initially suggested by the lack of low-mass X-ray binary observations with accretors about 2–5 M ⊙, but it has recently been called into question following newer observations, including an LMG candidate with a millisecond pulsar (MSP) companion in the dense globular cluster NGC 1851. Here, we model NGC 1851 with a grid of similar dense star clusters utilizing the state-of-the-art Monte Carlo N-body code Cluster Monte Carlo, and we specifically study the formation of LMG black holes (BHs). We demonstrate that both massive star evolution and dynamical interactions can contribute to forming LMG BHs. In general, the collapse of massive remnants formed through mergers of neutron stars (NSs) or massive white dwarfs produces the largest number of LMG BHs among all formation channels. However, in more massive clusters, supernova core collapse can contribute comparable numbers. Our NGC 1851-like models can reproduce MSP—LMG BH binaries similar to the observed system. Additionally, the LMG BHs can also become components of dynamically assembled binaries, and some will be in merging BH–NS systems similar to the recently detected gravitational wave source GW230529. However, the corresponding merger rate is probably ≲1 Gpc−3 yr−1.

Binary neutron star merger offsets from their host galaxies. GW 170817 as a case study

The locations of binary neutron star (BNS) mergers within their host galaxies encode the systemic kicks that these systems received in the supernova aftermath. We investigate how the galactic potential and the systemic kicks shape the offset distribution of BNS mergers with a case study of GW 170817 and its host NGC 4993. We derived dynamical constraints on the host potential from integral field spectroscopy with Jeans anisotropic modelling. We evolved the trajectories of synthetic BNSs from the BPASS code in the galactic potential, using two different kick prescriptions to investigate how the observed offsets might differentiate between these two possibilities. The simulation was repeated after swapping the host potential with that of a dwarf galaxy, to test the effect of the potential on the offsets. The location of GW 170817 is entirely consistent with our predictions, regardless of large or small kicks, because the strong potential of NGC 4993 is only diagnostic of very large kicks. In galaxies of similar or greater mass, large offsets can constrain large kicks, while small offsets do not provide much information. In an old dwarf galaxy, on the other hand, small offsets can constrain small kicks, while large offsets would prevent host association.

Differential rotation in neutron stars at finite temperatures

IntroductionThis paper investigates the impact of differential rotation on the bulk properties and onset of rotational instabilities in neutron stars at finite temperatures up to 50 MeV.MethodsUtilizing the relativistic Brueckner-Hartree-Fock (RBHF) formalism in full Dirac space, the study constructs equation of state (EOS) models for hot neutron star matter, including conditions relevant for high temperatures. These finite-temperature EOS models are applied to compute the bulk properties of differentially rotating neutron stars with varying structural deformations.ResultsThe findings demonstrate that the stability of these stars against bar-mode deformation, a key rotational instability, is only weakly dependent on temperature. Differential rotation significantly affects the maximum mass and radius of neutron stars, and the threshold for the onset of bar-mode instability shows minimal sensitivity to temperature changes within the examined range.DiscussionThese findings are crucial for interpreting observational data from neutron star mergers and other high-energy astrophysical events. The research underscores the necessity of incorporating differential rotation and finite temperature effects in neutron star models to predict their properties and stability accurately.

Off-Axis Color Characteristics of Binary Neutron Star Merger Events: Applications for Space Multi-Band Variable Object Monitor and James Webb Space Telescope

With advancements in gravitational wave detection technology, an increasing number of binary neutron star (BNS) merger events are expected to be detected. Due to the narrow opening angle of jet cores, many BNS merger events occur off-axis, resulting in numerous gamma-ray bursts (GRBs) going undetected. Models suggest that kilonovae, which can be observed off-axis, offer more opportunities to be detected in the optical/near-infrared band as electromagnetic counterparts of BNS merger events. In this study, we calculate kilonova emission using a three-dimensional semi-analytical code and model the GRB afterglow emission with the open-source Python package afterglowpy at various inclination angles. Our results show that it is possible to identify the kilonova signal from the observed color evolution of BNS merger events. We also deduce the optimal observing window for SVOM/VT and JWST/NIRCam, which depends on the viewing angle, jet opening angle, and circumburst density. These parameters can be cross-checked with the multi-band afterglow fitting. We suggest that kilonovae are more likely to be identified at larger inclination angles, which can also help determine whether the observed signals without accompanying GRBs originate from BNS mergers.

Binary neutron star mergers using a discontinuous Galerkin-finite difference hybrid method

Abstract We present a discontinuous Galerkin-finite difference hybrid scheme that allows high-order shock capturing with the discontinuous Galerkin method for general relativistic magnetohydrodynamics in dynamical spacetimes. We present several optimizations and stability improvements to our algorithm that allow the hybrid method to successfully simulate single, rotating, and binary neutron stars. The hybrid method achieves the efficiency of discontinuous Galerkin methods throughout almost the entire spacetime during the inspiral phase, while being able to robustly capture shocks and resolve the stellar surfaces. We also use Cauchy-Characteristic evolution to compute the first gravitational waveforms at future null infinity from binary neutron star mergers. The simulations presented here are the first successful binary neutron star inspiral and merger simulations using discontinuous Galerkin methods.

Prospect of Detecting Magnetic Fields from Strong-magnetized Binary Neutron Stars

Binary neutron star mergers are unique sources of gravitational waves in multi-messenger astronomy. The inspiral phase of binary neutron stars can emit gravitational waves as chirp signals. The present waveform models of gravitational waves only considered the gravitational interaction. In this paper, we derive the waveform of the gravitational wave signal taking into account the presence of magnetic fields. We found that the electromagnetic interaction and radiation can introduce different frequency-dependent power laws for both the amplitude and frequency of the gravitational wave. We show from the results of the Fisher information matrix that the third-generation observation may detect magnetic dipole moments if the magnetic field is ∼1017 G.

Probing kaon meson condensations through gravitational waves during neutron star inspiral phases

Gravitational waves from binary neutron star merger events have shown unparalleled advantages in deciphering the stellar internal structure, yet it remains underexplored to seek clues about kaon meson condensations inside massive neutron stars through these events. To address this gap, this study examines the binary star inspiral processes with kaon meson condensations, employing various kaon meson potential wells to probe unique imprints in the inspiral gravitational waves. Our analysis reveals that the inspiral process of kaon meson condensation binary stars exhibits lower gravitational wave frequencies and shorter retarded times compared to traditional neutron stars. Furthermore, a deeper kaon meson potential well in these systems further shortens the inspiral duration time, resulting in distinctive gravitational waveforms and notable orbital phase deviations. These findings underscore the pivotal role of kaon meson condensations in modulating the inspiral gravitational waves, suggesting the inspiral gravitational waves could serve as a powerful tool for probing potential signals of kaon mesons within neutron stars.

Stable case BB/BC mass transfer to form GW190425-like massive binary neutron star mergers

On April 25, 2019, the LIGO-Virgo Collaboration discovered a gravitational-wave (GW) signal from a binary neutron star (BNS) merger, that is, GW190425. Due to the inferred large total mass, the origin of GW190425 remains unclear. Assuming GW190425 originated from the standard isolated binary evolution channel, its immediate progenitor is considered to be a close binary system, consisting of a He-rich star and a NS just after the common envelope phase. We aim to study the formation of GW190425 in a solar-like environment by using the detailed binary evolution code MESA We perform detailed stellar structure and binary evolution calculations that take into account mass loss, internal differential rotation, and tidal interactions between a He-rich star and a NS companion. We explore the parameter space of the initial binary properties, including initial NS and He-rich masses and initial orbital period. We find that the immediate post-common-envelope progenitor system, consisting of a primary $ NS and a secondary He-rich star with an initial mass of $ in a close binary with an initial period of $ days days $), that experiences stable Case BB/BC mass transfer (MT) during binary evolution, can reproduce the formation of GW190425-like BNS events. Our studies reveal that the secondary He-rich star of the GW190425's progenitor before its core collapse can be efficiently spun up through tidal interaction, finally remaining as a NS with rotational energy even reaching $ erg $, which is always much higher than the neutrino-driven energy of the supernova (SN) explosion. If the newborn secondary NS is a magnetar, we expect that GW190425 can be the remnant of a magnetar-driven SN, namely a magnetar-driven ultra-stripped SN, a superluminous SN, or a broad-line Type Ic SN. Our results show that GW190425 could be formed through the isolated binary evolution, which involves a stable Case BB/BC MT just after the common envelope phase. On top of that, we show the He-rich star can be tidally spun up, potentially forming a spinning magnetized NS (magnetar) during the second SN explosion.

Measuring neutron star radius with second and third generation gravitational wave detector networks

Abstract The next generation of ground-based interferometric gravitational wave detectors will observe mergers of black holes and neutron stars throughout cosmic time. A large number of the binary neutron star merger events will be observed with extreme high fidelity, and will provide stringent constraints on the equation of state of nuclear matter. In this paper, we investigate the systematic improvement in the measurability of the equation of state with increase in detector sensitivity by combining constraints obtained on the radius of a $1.4 \, \mathrm{M}_{\odot}$ neutron star from a simulated source population. Since the measurability of the equation of state depends on its stiffness, we consider a range of realistic equations of state that span the current observational constraints. We show that a single 40km Cosmic Explorer detector can pin down the neutron star radius for a soft, medium and stiff equation of state to an accuracy of 10m within a decade, whereas the current generation of ground-based detectors like the Advanced LIGO-Virgo network would take $\mathcal{O}(10^5)$ years to do so for a soft equation of state. &amp;#xD;

Robustness of neutron star merger simulations to changes in neutrino transport and neutrino-matter interactions

Robustness of neutron star merger simulations to changes in neutrino transport and neutrino-matter interactions

A simulation study on the sub-threshold joint gravitational wave-electromagnetic wave observation on binary neutron star mergers

Abstract The coalescence of binary neutron stars (BNS) is a prolific source of gravitational waves (GWs) and electromagnetic (EM) radiation, offering a dual observational window into the Universe. Lowering the signal-to-noise ratio (S/N) threshold is a simple and cost-effective way to enhance the detection probability of GWs from BNS mergers. In this study, we introduce a metric of the purity of joint GW and EM detections Pjoint, which is in analogue to Pastro in GW only observations. By simulating BNS merger GWs jointly detected by the HLV network and EM counterparts (kilonovae and short Gamma-ray bursts, sGRBs) with an assumed merger rate density of BNS, we generate catalogs of GW events and EM counterparts. Through this simulation, we analyze joint detection pairs, both correct and misidentified. We find the following: 1. For kilonovae, requiring Pjoint &amp;gt; 95% instead of $P_{\rm astro}&amp;gt;95~{{\%}}$ reduces the S/N from 9.2 to 8.5-8.8, allowing 5-13 additional joint detections per year and increasing the GW detection volume by 9-17%; 2. For sGRBs, requiring Pjoint &amp;gt; 95% instead of Pastro reduces the S/N from 9.2 to 8.1-8.5; 3. Increasing kilonova or sGRB detection capability does not improve Pjoint due to a higher rate of misidentifications. We also show that sub-threshold GW and kilonova detections can reduce the uncertainty in measuring the Hubble constant to 89-92% of its original value, and sub-threshold GW and sGRB observations can enhance the precision of constraining the speed of GWs to 88% of previously established values.

Large eddy simulations of magnetized mergers of black holes and neutron stars

Large eddy simulations of magnetized mergers of black holes and neutron stars