Search for: All records

Abstract Binary neutron star mergers can produce extreme magnetic fields, some of which can lead to strong magnetar-like remnants. While strong magnetic fields have been shown to affect the dynamics of outflows and angular momentum transport in the remnant, they can also crucially alter the properties of nuclear matter probed in the merger. In this work, we provide a first assessment of the latter, determining the strength of the pressure anisotropy caused by Landau-level quantization and the anomalous magnetic moment. To this end, we perform the first numerical relativity simulation with a magnetic polarization tensor and a magnetic-field-dependent equation of state using a new algorithm we present here, which also incorporates a mean-field dynamo model to control the magnetic field strength present in the merger remnant. Our results show that—in the most optimistic case—corrections to the anisotropy can be in excess of 10% and are potentially largest in the outer layers of the remnant. This work paves the way for a systematic investigation of these effects. 

Abstract A significant interest has emerged recently in assessing whether collimated and ultrarelativistic outflows can be produced by a long-lived remnant from a binary neutron star (BNS) merger, with different approaches leading to different outcomes. To clarify some of the aspect of this process, we report the results of long-term (i.e., ∼110 ms) state-of-the-art general relativistic magnetohydrodynamics simulations of the inspiral and merger of a BNS system of magnetized stars. We find that after ∼50 ms from the merger anα–Ω dynamo driven by the magnetorotational instability sets in in the densest regions of the disk and leads to the breakout of the magnetic field lines from the accretion disk around the remnant. The breakout is responsible for generating a collimated, magnetically driven outflow with only mildly relativistic velocities and for a violent eruption of electromagnetic energy. We provide evidence that this outflow is partly collimated via a Blandford–Payne mechanism. Finally, by including or not the radiative transport via neutrinos, we determine the role they play in the launching of the collimated wind. In this way, we conclude that the mechanism of magnetic field breakout we observe is robust and takes place even without neutrinos. Contrary to previous expectations, the inclusion of neutrino absorption and emission leads to a smaller baryon pollution in polar regions and hence accelerates the occurrence of the breakout, yielding a larger electromagnetic luminosity. Given the mildly relativistic nature of these disk-driven breakout outflows, it is difficult to consider them responsible for the jet phenomenology observed in short gamma-ray bursts. 

Abstract Black hole–neutron star binaries are of interest in many ways: they are intrinsically transient, radiate gravitational waves detectable by LIGO, and may produceγ-ray bursts. Although it has long been assumed that their late-stage orbital evolution is driven entirely by gravitational wave emission, we show here that in certain circumstances, mass transfer from the neutron star onto the black hole can both alter the binary's orbital evolution and significantly reduce the neutron star's mass: when the fraction of its mass transferred per orbit is ≳10⁻², the neutron star's mass diminishes by order unity, leading to mergers in which the neutron star mass is exceptionally small. The mass transfer creates a gas disk around the black holebeforemerger that can be comparable in mass to the debris remaining after merger, i.e., ~0.1M_⊙. These processes are most important when the initial neutron star–black hole mass ratioqis in the range ≈0.2–0.8, the orbital semimajor axis is 40 ≲ a₀/r_g ≲ 300 (r_g ≡ GM_BH/c²), and the eccentricity is large ate₀ ≳ 0.8. Systems of this sort may be generated through the dynamical evolution of a triple system, as well as by other means. 

Abstract We perform the first magnetohydrodynamic simulation tracking the magnetosphere of a collapsing magnetar. The collapse is expected for massive rotating magnetars formed in merger events and may occur many hours after the merger. Our simulation suggests a novel mechanism for a gamma-ray burst (GRB), which is uncollimated and forms a delayed high-energy counterpart of the merger gravitational waves. The simulation shows that the collapse launches an outgoing magnetospheric shock, and a hot magnetized outflow forms behind the shock. The outflow is baryon free and uncollimated, and its power peaks on a millisecond timescale. Then, the outflow becomes modulated by the ring-down of the nascent black hole, imprinting its kilohertz quasi-normal modes on the GRB tail. 

Abstract A variety of high-energy astrophysical phenomena are powered by the release—via magnetic reconnection—of the energy stored in oppositely directed fields. Single-fluid resistive magnetohydrodynamic (MHD) simulations with uniform resistivity yield dissipation rates that are much lower (by nearly 1 order of magnitude) than equivalent kinetic calculations. Reconnection-driven phenomena could be accordingly modeled in resistive MHD employing a nonuniform, “effective” resistivity informed by kinetic calculations. In this work, we analyze a suite of fully kinetic particle-in-cell (PIC) simulations of relativistic pair-plasma reconnection—where the magnetic energy is greater than the rest mass energy—for different strengths of the guide field orthogonal to the alternating component. We extract an empirical prescription for the effective resistivity, ${\eta }_{\mathrm{eff}}=\alpha B_0\mid \mathit{J}{\mid }^p/\left(\mid \mathit{J}{\mid }^{p+1}+{\left(e{n}_{t}c\right)}^{p+1}\right)$, whereB₀is the reconnecting magnetic field strength,Jis the current density,n_tis the lab-frame total number density,eis the elementary charge, andcis the speed of light. The guide field dependence is encoded inαandp, which we fit to PIC data. This resistivity formulation—which relies only on single-fluid MHD quantities—successfully reproduces the spatial structure and strength of nonideal electric fields and thus provides a promising strategy for enhancing the reconnection rate in resistive MHD simulations. 

Abstract In nuclear matter in isolated neutron stars, the flavor content (e.g., proton fraction) is subject to weak interactions, establishing flavor (β-)equilibrium. However, there can be deviations from this equilibrium during the merger of two neutron stars. We study the resulting out-of-equilibrium dynamics during the collision by incorporating direct and modified Urca processes (in the neutrino-transparent regime) into general-relativistic hydrodynamics simulations with a simplified neutrino transport scheme. We demonstrate how weak-interaction-driven bulk viscosity in postmerger simulations can emerge and assess the bulk viscous dynamics of the resulting flow. We further place limits on the impact of the postmerger gravitational-wave strain. Our results show that weak-interaction-driven bulk viscosity can potentially lead to a phase shift of the postmerger gravitational-wave spectrum, although the effect is currently on the same level as the numerical errors of our simulation.