Hypermassive Neutron Star Research Articles

Coalescence of binary neutron stars in a scalar-tensor theory of gravity

We carry out numerical-relativity simulations of coalescing binary neutron stars in a scalar-tensor theory that admits spontaneous scalarization. We model neutron stars with realistic equations of state. We choose the free parameters of the theory taking into account the constraints imposed by the latest observations of neutron-star-- white-dwarf binaries with pulsar timing. We show that even within those severe constraints, scalarization can still affect the evolution of the binary neutron stars not only during the late inspiral, but also during the merger stage. We also confirm that even when both neutron stars have quite small scalar charge at large separations, they can be strongly scalarized dynamically during the final stages of the inspiral. In particular, we identify the binary parameters for which scalarization occurs either during the late inspiral or only after the onset of the merger when a remnant, supramassive or hypermassive neutron star is formed. We also discuss how those results can impact the extraction of physical information on gravitational waves once they are detected.

MAGNETICALLY DRIVEN WINDS FROM DIFFERENTIALLY ROTATING NEUTRON STARS AND X-RAY AFTERGLOWS OF SHORT GAMMA-RAY BURSTS

Besides being among the most promising sources of gravitational waves, merging neutron star binaries also represent a leading scenario to explain the phenomenology of short gamma-ray bursts (SGRBs). Recent observations have revealed a large subclass of SGRBs with roughly constant luminosity in their X-ray afterglows, lasting $10\!-\!10^4$ s. These features are generally taken as evidence of a long-lived central engine powered by the magnetic spin-down of a uniformly rotating, magnetized object. We propose a different scenario in which the central engine powering the X-ray emission is a differentially rotating hypermassive neutron star (HMNS) that launches a quasi-isotropic and baryon-loaded wind driven by the magnetic field, which is built-up through differential rotation. Our model is supported by long-term, three-dimensional, general-relativistic, and ideal magnetohydrodynamic simulations, showing that this isotropic emission is a very robust feature. For a given HMNS, the presence of a collimated component depends sensitively on the initial magnetic field geometry, while the stationary electromagnetic luminosity depends only on the magnetic energy initially stored in the system. We show that our model is compatible with the observed timescales and luminosities and express the latter in terms of a simple scaling relation.

Linking electromagnetic and gravitational radiation in coalescing binary neutron stars

We expand on our study of the gravitational and electromagnetic emissions from the late stage of an inspiraling neutron star binary as presented in Ref. \cite{Palenzuela:2013hu}. Interactions between the stellar magnetospheres, driven by the extreme dynamics of the merger, can yield considerable outflows. We study the gravitational and electromagnetic waves produced during the inspiral and merger of a binary neutron star system using a full relativistic, resistive MHD evolution code. We show that the interaction between the stellar magnetospheres extracts kinetic energy from the system and powers radiative Poynting flux and heat dissipation. These features depend strongly on the configuration of the initial stellar magnetic moments. Our results indicate that this power can strongly outshine pulsars in binaries and have a distinctive angular and time-dependent pattern. Our discussion provides more detail than Ref. \cite{Palenzuela:2013hu}, showing clear evidence of the different effects taking place during the inspiral. Our simulations include a few milliseconds after the actual merger and study the dynamics of the magnetic fields during the formation of the hypermassive neutron star. We also briefly discuss the possibility of observing such emissions.

Remnant massive neutron stars of binary neutron star mergers: Evolution process and gravitational waveform

Massive (hypermassive and supramassive) neutron stars are likely to be often formed after the merger of binary neutron stars. We explore the evolution process of the remnant massive neutron stars and gravitational waves emitted by them, based on numerical-relativity simulations for binary neutron star mergers employing a variety of equations of state and choosing a plausible range of the neutron star mass of binaries. We show that the lifetime of remnant hypermassive neutron stars depends strongly on the total binary mass and also on the equations of state. Gravitational waves emitted by the remnant massive neutron stars universally have a quasiperiodic nature of an approximately constant frequency although the frequency varies with time. We also show that the frequency and time-variation feature of gravitational waves depend strongly on the equations of state. We derive a fitting formula for the quasiperiodic gravitational waveforms, which may be used for the data analysis of a gravitational-wave signal.

Magnetorotational instability in relativistic hypermassive neutron stars

A differentially rotating hypermassive neutron star (HMNS) is a metastable object which can be formed in the merger of neutron-star binaries. The eventual collapse of the HMNS into a black hole is a key element in generating the physical conditions expected to accompany the launch of a short gamma-ray burst. We investigate the influence of magnetic fields on HMNSs by performing three-dimensional simulations in general-relativistic magnetohydrodynamics. In particular, we provide direct evidence for the occurrence of the magnetorotational instability (MRI) in HMNS interiors. For the first time in simulations of these systems, rapidly-growing and spatially-periodic structures are observed to form with features like those of the channel flows produced by the MRI in other systems. Moreover, the growth time and wavelength of the fastest-growing mode are extracted and compared successfully with analytical predictions. The MRI emerges as an important mechanism to amplify magnetic fields over the lifetime of the HMNS, whose collapse to a black hole is accelerated. The evidence provided here that the MRI can actually develop in HMNSs could have a profound impact on the outcome of the merger of neutron-star binaries and on its connection to short gamma-ray bursts.

COMPACT BINARY PROGENITORS OF SHORT GAMMA-RAY BURSTS

In recent years, detailed observations and accurate numerical simulations have provided support to the idea that mergers of compact binaries containing either two neutron stars (NSs) or an NS and a black hole (BH) may constitute the central engine of short gamma-ray bursts (SGRBs). The merger of such compact binaries is expected to lead to the production of a spinning BH surrounded by an accreting torus. Several mechanisms can extract energy from this system and power the SGRBs. Here we connect observations and numerical simulations of compact binary mergers, and use the current sample of SGRBs with measured energies to constrain the mass of their powering tori. By comparing the masses of the tori with the results of fully general-relativistic simulations, we are able to infer the properties of the binary progenitors which yield SGRBs. By assuming a constant efficiency in converting torus mass into jet energy, epsilon_{jet}=10%, we find that most of the tori have masses smaller than 0.01M_{sun}, favoring "high-mass" binary NSs mergers, i.e., binaries with total masses >~1.5 the maximum mass of an isolated NS. This has important consequences for the gravitational-wave signals that may be detected in association with SGRBs, since "high-mass" systems do not form a long-lived hypermassive NS after the merger. While NS-BH systems cannot be excluded to be the engine of at least some of the SGRBs, the BH would need to have an initial spin of ~0.9, or higher.

Intense electromagnetic outbursts from collapsing hypermassive neutron stars

We study the gravitational collapse of a magnetized neutron star using a novel numerical approach able to capture both the dynamics of the star and the behavior of the surrounding plasma. In this approach, a fully general relativistic magnetohydrodynamics implementation models the collapse of the star and provides appropriate boundary conditions to a force-free model which describes the stellar exterior. We validate this strategy by comparing with known results for the rotating monopole and aligned rotator solutions and then apply it to study both rotating and non-rotating stellar collapse scenarios, and contrast the behavior with what is obtained when employing the electrovacuum approximation outside the star. The non-rotating electrovacuum collapse is shown to agree qualitatively with a Newtonian model of the electromagnetic field outside a collapsing star. We illustrate and discuss a fundamental difference between the force-free and electrovacuum solutions, involving the appearance of large zones of electric-dominated field in the vacuum case. This provides a clear demonstration of how dissipative singularities appear generically in the non-linear time-evolution of force-free fluids. In both the rotating and non-rotating cases, our simulations indicate that the collapse induces a strong electromagnetic transient. In the case of sub-millisecond rotation, the magnetic field experiences strong winding and the transient carries much more energy. This result has important implications for models of gamma-ray bursts.

DYNAMICAL CAPTURE BINARY NEUTRON STAR MERGERS

We study dynamical capture binary neutron star mergers as may arise in dense stellar regions such as globular clusters. Using general-relativistic hydrodynamics, we find that these mergers can result in the prompt collapse to a black hole or in the formation of a hypermassive neutron star, depending not only on the neutron star equation of state but also on impact parameter. We also find that these mergers can produce accretion disks of up to a tenth of a solar mass and unbound ejected material of up to a few percent of a solar mass. We comment on the gravitational radiation and electromagnetic transients that these sources may produce.

Importance of cooling in triggering the collapse of hypermassive neutron stars

The inspiral and merger of a binary neutron star (NSNS) can lead to the formation of a hypermassive neutron star (HMNS). As the HMNS loses thermal pressure due to neutrino cooling and/or centrifugal support due to gravitational wave emission, and/or magnetic breaking of differential rotation, it will collapse to a black hole. To assess the importance of shock-induced thermal pressure and cooling, we adopt an idealized equation of state and perform NSNS simulations in full general relativity through late inspiral, merger, and HMNS formation, accounting for cooling. We show that thermal pressure contributes significantly to the support of the HMNS against collapse and that thermal cooling accelerates its ``delayed'' collapse. Our simulations demonstrate explicitly that cooling can induce the catastrophic collapse of a hot hypermassive neutron star formed following the merger of binary neutron stars. Thus, cooling physics is important to include in NSNS merger calculations to accurately determine the lifetime of the HMNS remnant and to extract information about the neutron star equation of state, cooling mechanisms, bar instabilities and B-fields from the gravitational waves emitted during the transient phase prior to black hole formation.

Gravitational waves, neutrino emissions and effects of hyperons in binary neutron star mergers

Numerical simulations for the merger of binary neutron stars are performed in full general relativity incorporating both nucleonic and hyperonic finite-temperature equations of state (EOS) and neutrino cooling. It is found that for the nucleonic and hyperonic EOS, a hyper-massive neutron star (HMNS) with a long lifetime (tlife ≳ 10 ms) is the outcome for the total mass ≈2.7 M⊙. For the total mass ≈3 M⊙, a long-lived (short-lived with tlife ≈ 3 ms) HMNS is the outcome for the nucleonic (hyperonic) EOS. It is shown that the typical total neutrino luminosity of the HMNS is ∼3–6 × 1053 erg s−1 and the effective amplitude of gravitational waves from the HMNS is 1–4 × 10−22 at f ≈ 2–3.2 kHz for a source of distance of 100 Mpc. During the HMNS phase, characteristic frequencies of gravitational waves shift to a higher frequency for the hyperonic EOS in contrast to the nucleonic EOS in which they remain constant approximately. Our finding suggests that the effects of hyperons are well imprinted in gravitational waves and their detection will give us a potential opportunity to explore the composition of the neutron star matter. We present the neutrino luminosity curve when a black hole is formed as well.

Current status of numerical-relativity simulations in Kyoto

We describe the current status of our numerical simulations for the collapse of a massive stellar core to a black hole (BH) and the merger of binary neutron stars (BNS), performed in the framework of full general relativity incorporating finite-temperature equations of state (EOS) and neutrino cooling. For the stellar core collapse simulation, we present the latest numerical results. We employed a purely nucleonic EOS derived by Shen et al. [Nucl. Phys. A 637, 435 (1998)]. As an initial condition, we adopted a 100 M presupernova model calculated by Umeda and Nomoto [Astrophys. J. 637, 1014 (2008)], which has a massive core (M ≈ 3M) with a high value of entropy per baryon (s ≈ 4kB). Changing the degree of rotation for the initial condition, we clarify the strong dependence of the outcome of the collapse on this. When the rotation is rapid enough, the shock wave formed at the core bounce is deformed to a torus-like shape. Then, the infalling matter accumulates in the central region due to the oblique shock at the torus surface, hitting the proto-neutron star and dissipating the kinetic energy there. As a result, outflows can be launched. The proto-neutron eventually collapses to a BH and an accretion torus is formed around it.We also found that the evolution of the BH and torus depends strongly on the rotation initially given. In the BNS merger simulations, we employ an EOS incorporating a degree of freedom for hyperons derived by Shen et al. [Astrophys. J. Suppl. 197, 20 (2011)], in addition to the purely nucleonic EOS. The numerical simulations show that for the purely nucleonic EOS, a hypermassive neutron star (HMNS) with a long lifetime (10 ms) is the outcome for the total mass M ≤ 3.0M. In contrast, the formed HMNS collapses to a BH in a shorter time scale with the hyperonic EOS for M ≥ 2.7M. It is shown that the typical total neutrino luminosity of the HMNS is ∼(3–10)×1053 ergs/s and the effective amplitude of gravitational waves from the HMNS is (2–6)×10−22 at f ≈ 2–2.5 kHz for a source distance of 100 Mpc.

Effects of Hyperons in Binary Neutron Star Mergers

Numerical simulations for the merger of binary neutron stars are performed in full general relativity incorporating both nucleonic and hyperonic finite-temperature equations of state (EOS) and neutrino cooling. It is found that even for the hyperonic EOS, a hypermassive neutron star is first formed after the merger for the typical total mass ≈2.7M(⊙), and subsequently collapses to a black hole (BH). It is shown that hyperons play a substantial role in the postmerger dynamics, torus formation around the BH, and emission of gravitational waves (GWs). In particular, the existence of hyperons is imprinted in GWs. Therefore, GW observations will provide a potential opportunity to explore the composition of neutron star matter.

Gravitational Waves and Neutrino Emission from the Merger of Binary Neutron Stars

Numerical simulations for the merger of binary neutron stars are performed in full general relativity incorporating a finite-temperature (Shen's) equation of state (EOS) and neutrino cooling for the first time. It is found that for this stiff EOS, a hypermassive neutron star (HMNS) with a long lifetime (≫10 ms) is the outcome for the total mass ≲3.0M(⊙). It is shown that the typical total neutrino luminosity of the HMNS is ∼3-8×10(53) erg/s and the effective amplitude of gravitational waves from the HMNS is 4-6×10(-22) at f=2.1-2.5 kHz for a source distance of 100 Mpc. We also present the neutrino luminosity curve when a black hole is formed for the first time.

Binary neutron star mergers: Dependence on the nuclear equation of state

We perform a numerical-relativity simulation for the merger of binary neutron stars with 6 nuclear-theory-based equations of state (EOSs) described by piecewise polytropes. Our purpose is to explore the dependence of the dynamical behavior of the binary neutron star merger and resulting gravitational waveforms on the EOS of the supernuclear-density matter. The numerical results show that the merger process and the first outcome are classified into three types; (i) a black hole is promptly formed, (ii) a short-lived hypermassive neutron star (HMNS) is formed, (iii) a long-lived HMNS is formed. The type of the merger depends strongly on the EOS and on the total mass of the binaries. For the EOS with which the maximum mass is larger than 2Msun, the lifetime of the HMNS is longer than 10 ms for a total mass m_0=2.7Msun. A recent radio observation suggests that the maximum mass of spherical neutron stars is M_max \geq 1.97\pm 0.04Msun in one \sigma level. This fact and our results support the possible existence of a HMNS soon after the onset of the merger for a typical binary neutron star with m_0=2.7Msun. We also show that the torus mass surrounding the remnant black hole is correlated with the type of the merger process; the torus mass could be large, \geq 0.1Msun, in the case that a long-lived HMNS is formed. We also show that gravitational waves carry information of the merger process, the remnant, and the torus mass surrounding a black hole.

AFTERGLOW OF A BINARY NEUTRON STAR MERGER

The merger of two neutron stars often results in a rapidly and differentially rotating hypermassive neutron star (HMNS). We show by numerical-relativity simulation that the magnetic-field profile around such HMNS is dynamically varied during its subsequent evolution, and as a result, electromagnetic radiation with a large luminosity ∼0.1B2R3Ω is emitted with baryons (B, R, and Ω are poloidal magnetic-field strength at stellar surface, stellar radius, and angular velocity of an HMNS). The predicted luminosity of electromagnetic radiation, which is primarily emitted along the magnetic-dipole direction, is ∼1047(B/1013 G)2(R/10 km)3(Ω/104 rad s−1) erg s−1, which is comparable to the luminosity of quasars.

Formation and evolution of black hole and accretion disk in collapse of massive stellar cores

AbstractWe describe the results of our numerical simulations of the collapse of a massive stellar core to a BH, performed in the framework of full general relativity incorporating finite-temperature equation of state and neutrino cooling. We adopt a 100 M⊙ presupernova model calculated by Umeda & Nomoto (2008), which has a massive core with a high value of entropy per baryon. Changing the degree of rotation for the initial condition, we clarify the dependence of the outcome on this. When the rotation is rapid enough, the shock wave formed at the core bounce is deformed to be a torus-like shape. Then, the infalling matter is accumulated in the central region due to the oblique shock at the torus surface, hitting the hypermassive neutron star (HMNS) and dissipating the kinetic energy there. As a result, outflows can be launched. The HMNS eventually collapses to a BH and an accretion torus is formed around it. We also found that the evolution of the BH and torus depends strongly on the rotation initially given.

Observational constraints on the nature of very short gamma-ray bursts

We discuss a very peculiar subgroup of gamma-ray bursts among the BATSE sources. These bursts are very short ( T 90 ⩽ 0.1 s), hard, and came predominantly from a restricted direction of the sky (close to the Galactic anti-center). We analyze their arrival times and possible correlations, as well as the profiles of individual bursts. We find no peculiarities in the arrival times of Very Short Bursts (VSBs) despite their highly non-uniform spatial distribution. There is no dependence in the burst shapes on location. Bursts coming both from the burst-enhancement Galactic Anticenter region and from all other directions show considerable dispersion in their rise and fall times. Significant fraction of VSBs have multiple peaks despite their extremely short duration. Burst time properties are most likely to be consistent with two origin mechanisms: either with binary NS–NS mergers with low total masses passing through a phase of hypermassive neutron star, or with evaporation of the primordial black holes in the scenario of no photosphere formation.

Accurate evolutions of inspiralling neutron-star binaries: Prompt and delayed collapse to a black hole

Binary neutron-star systems represent primary sources for the gravitational-wave detectors that are presently operating or are close to being operating at the target sensitivities. We present a systematic investigation in full general relativity of the dynamics and gravitational-wave emission from binary neutron stars which inspiral and merge, producing a black hole surrounded by a torus. Our results represent the state of the art from several points of view: (i) We use high-resolution shock-capturing methods for the solution of the hydrodynamics equations and high-order finite-differencing techniques for the solution of the Einstein equations; (ii) We employ adaptive mesh-refinement techniques with ``moving boxes'' that provide high-resolution around the orbiting stars; (iii) We use as initial data accurate solutions of the Einstein equations for a system of binary neutron stars in irrotational quasicircular orbits; (iv) We exploit the isolated-horizon formalism to measure the properties of the black holes produced in the merger; (v) Finally, we use two approaches, based either on gauge-invariant perturbations or on Weyl scalars, to calculate the gravitational waves emitted by the system. Within our idealized treatment of the matter, these techniques allow us to perform accurate evolutions on time scales never reported before (i.e. $\ensuremath{\sim}30\text{ }\text{ }\mathrm{ms}$) and to provide the first complete description of the inspiral and merger of a neutron-star binary leading to the prompt or delayed formation of a black hole and to its ringdown. We consider either a polytropic equation of state or that of an ideal fluid and show that already with this idealized treatment a very interesting phenomenology can be described. In particular, we show that while higher-mass polytropic binaries lead to the prompt formation of a rapidly rotating black hole surrounded by a dense torus, lower-mass binaries give rise to a differentially rotating star, which undergoes large oscillations and emits large amounts of gravitational radiation. Eventually, also the hyper-massive neutron star collapses to a rotating black hole surrounded by a torus. Finally, we also show that the use of a nonisentropic equation of state leads to significantly different evolutions, giving rise to a delayed collapse also with high-mass binaries, as well as to a more intense emission of gravitational waves and to a geometrically thicker torus.

General relativistic simulations of magnetized binary neutron star mergers

Binary neutron stars (NSNS) are expected to be among the leading sources of gravitational waves observable by ground-based laser interferometers and may be the progenitors of short-hard gamma-ray bursts. We present a series of general relativistic NSNS coalescence simulations both for unmagnetized and magnetized stars. We adopt quasiequilibrium initial data for circular, irrotational binaries constructed in the conformal thin-sandwich (CTS) framework. We adopt the BSSN formulation for evolving the metric and a high-resolution shock-capturing scheme to handle the magnetohydrodynamics. Our simulations of unmagnetized binaries agree with the results of Shibata, Taniguchi and Ury\ifmmode \bar{u}\else \={u}\fi{} [M. Shibata, K. Taniguchi, and K. Ury\ifmmode \bar{u}\else \={u}\fi{}, Phys. Rev. D 68, 084020 (2003).]. In cases in which the mergers result in a prompt collapse to a black hole, we are able to use puncture gauge conditions to extend the evolution and determine the mass of the material that forms a disk. We find that the disk mass is less than 2% of the total mass in all cases studied. We then add a small poloidal magnetic field to the initial configurations and study the subsequent evolution. For cases in which the remnant is a hypermassive neutron star, we see measurable differences in both the amplitude and phase of the gravitational waveforms following the merger. For cases in which the remnant is a black hole surrounded by a disk, the disk mass and the gravitational waveforms are about the same as the unmagnetized cases. Magnetic fields substantially affect the long-term, secular evolution of a hypermassive neutron star (driving ``delayed collapse'') and an accretion disk around a nascent black hole.

Collapse of magnetized hypermassive neutron stars in general relativity: Disk evolution and outflows

We study the evolution in axisymmetry of accretion disks formed self-consistently through collapse of magnetized hypermassive neutron stars to black holes. Such stars can arise following the merger of binary neutron stars. They are differentially rotating, dynamically stable, and have rest masses exceeding the mass limit for uniform rotation. However, hypermassive neutron stars are secularly unstable to collapse due to MHD-driven angular momentum transport. The rotating black hole which forms in this process is surrounded by a hot, massive, magnetized torus and a magnetic field collimated along the spin axis. This system is a candidate for the central engine of a short-hard gamma-ray burst (GRB). Our code integrates the coupled Einstein-Maxwell-MHD equations and is used to follow the collapse of magnetized hypermassive neutron star models in full general relativity until the spacetime settles down to a quasistationary state. We then employ the Cowling approximation, in which the spacetime is frozen, to track the subsequent evolution of the disk. This approximation allows us to greatly extend the disk evolutions and study the resulting outflows, which may be relevant to the generation of a GRB. We find that outflows are suppressed when a stiff equation of state is assumed for low-density disk material and are sensitive to the initial magnetic field configuration.