Next Article in Journal
De Sitter-Invariant Black Holes
Next Article in Special Issue
Estimates of the Surface Magnetic Field Strength of Radio Pulsars
Previous Article in Journal
Characterizing Low-Energy Charged Particles in the Magnetosphere with the LEM CubeSat Spectrometer Project: Detector Concept and Hardware Characterisation
Previous Article in Special Issue
A Short Survey of Matter-Antimatter Evolution in the Primordial Universe
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Neutron Star Binaries Produced by Binary-Driven Hypernovae, Their Mergers, and the Link between Long and Short GRBs

by
Laura M. Becerra
1,2,*,†,
Chris Fryer
3,†,
Jose F. Rodriguez
1,2,*,†,
Jorge A. Rueda
2,4,5,6,7,*,† and
Remo. Ruffini
2,4,8,†
1
GIRG, Escuela de Física, Universidad Industrial de Santander, Bucaramanga 680002, Colombia
2
ICRANet, Piazza della Repubblica 10, 65122 Pescara, Italy
3
CCS-2, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
4
ICRA, Dip. di Fisica, Sapienza Università di Roma, P.le Aldo Moro 5, 00185 Rome, Italy
5
ICRANet-Ferrara, Dip. di Fisica e Scienze Della Terra, Università Degli Studi di Ferrara, Via Saragat 1, 44122 Ferrara, Italy
6
Dip. di Fisica e Scienze della Terra, Università degli Studi di Ferrara, Via Saragat 1, 44122 Ferrara, Italy
7
INAF, Istituto de Astrofisica e Planetologia Spaziali, Via Fosso del Cavaliere 100, 00133 Rome, Italy
8
INAF, Viale del Parco Mellini 84, 00136 Rome, Italy
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Universe 2023, 9(7), 332; https://doi.org/10.3390/universe9070332
Submission received: 15 May 2023 / Revised: 8 July 2023 / Accepted: 9 July 2023 / Published: 12 July 2023
(This article belongs to the Special Issue Remo Ruffini Festschrift)

Abstract

:
The binary-driven hypernova (BdHN) model explains long gamma-ray bursts (GRBs) associated with supernovae (SNe) Ic through physical episodes that occur in a binary composed of a carbon-oxygen (CO) star and a neutron star (NS) companion in close orbit. The CO core collapse triggers the cataclysmic event, originating the SN and a newborn NS (hereafter ν NS) at its center. The ν NS and the NS accrete SN matter. BdHNe are classified based on the NS companion fate and the GRB energetics, mainly determined by the orbital period. In BdHNe I, the orbital period is of a few minutes, so the accretion causes the NS to collapse into a Kerr black hole (BH), explaining GRBs of energies > 10 52 erg. BdHN II, with longer periods of tens of minutes, yields a more massive but stable NS, accounting for GRBs of 10 50 10 52 erg. BdHNe III have still longer orbital periods (e.g., hours), so the NS companion has a negligible role, which explains GRBs with a lower energy release of < 10 50 erg. BdHN I and II might remain bound after the SN, so they could form NS-BH and binary NS (BNS), respectively. In BdHN III, the SN likely disrupts the system. We perform numerical simulations of BdHN II to compute the characteristic parameters of the BNS left by them, their mergers, and the associated short GRBs. We obtain the mass of the central remnant, whether it is likely to be a massive NS or a BH, the conditions for disk formation and its mass, and the event’s energy release. The role of the NS nuclear equation of state is outlined.

1. Introduction

Gamma-ray bursts (GRBs) are classified using the time (in the observer’s frame) T 90 , in which 90 % of the observed isotropic energy ( E iso ) in the gamma-rays is released. Long GRBs have T 90 > 2 s and short GRBs, T 90 < 2 s [1,2,3,4,5]. The two types of sources, short and long GRBs, are thought to be related to phenomena occurring in gravitationally collapsed objects, e.g., stellar-mass black holes (BHs) and neutron stars (NSs).
For short GRBs, mergers of binary NSs (BNSs) and/or NS-BH were soon proposed as progenitors [6,7,8,9]). For long bursts, the core-collapse of a single massive star leading to a BH (or a magnetar), a collapsar [10], surrounded by a massive accretion disk has been the traditional progenitor (see, e.g., [11,12], for reviews). The alternative binary-driven hypernova (BdHN) model exploits the increasing evidence for the relevance of a binary progenitor for long GRBs, e.g., their association with Ic-type supernovae (SNe) [13,14,15,16], proposing a binary system composed of a carbon-oxygen star (CO) and an NS companion for long GRBs. We refer the reader to [17,18,19,20,21,22,23] for theoretical details on the model.
In this article, we are interested in the direct relationship between long and short GRBs predicted by the BdHN scenario. The CO undergoes core collapse, ejecting matter in a supernova (SN) explosion and forming a newborn NS ( ν NS) at its center. The NS companion attracts part of the ejected material leading to an accretion process with high infalling rates. Also, the ν NS gains mass via a fallback accretion process. The orbital period is the most relevant parameter for the CO-NS system’s fate. In BdHN of type I, the NS reaches the critical mass, gravitationally collapsing into a Kerr BH. It occurs for short orbital periods (usually a few minutes) and explains GRBs with energies above 10 52 erg. In BdHN II, the orbital period is larger, up to a few tens of minutes, so the accretion rate decreases, and the NS becomes more massive but remains stable. These systems explain GRBs with energies 10 50 10 52 erg. In BdHN III, the orbital separation is still larger; the NS companion does not play any role, and the energy release is lower than 10 50 erg. If the binary is not disrupted by the mass loss in the SN explosion (see [20] for details), a BdHN I produces a BH-NS, whereas a BdHN II produces a BNS. In BdHN III, the SN is expected to disrupt the system. Therefore, in due time, the mergers of NS-BHs left by BdHNe I and of BNS left by BdHNe II are expected to lead to short GRBs.
Short GRBs from BNS mergers have been classified into short gamma-ray flashes (S-GRFs) and authentic short GRBs (S-GRBs), depending on whether the central remnant is an NS or a BH, respectively [24]. Two different subclasses of short GRBs from BNS mergers have been electromagnetically proposed [20,24,25]:
(1) Authentic short GRBs (S-GRBs): short bursts with isotropic energy E iso 10 52 erg and peak energy E p , i 2 MeV. They occur when a BH is formed in the merger, which is revealed by the onset of a GeV emission (see [25,26,27]). Their electromagnetically inferred isotropic occurrence rate is ρ S GRB 1.9 1.1 + 1.8 × 10 3 Gpc 3 year 1 [24]. The distinct signature of the formation of the BH, namely the observation of the 0.1 –100 GeV emission by the Fermi-LAT, needs the presence of baryonic matter interacting with the newly-formed BH, e.g., via an accretion process (see, e.g., [26,28]).
(2) Short gamma-ray flashes (S-GRFs): short bursts with E iso 10 52 erg and E p , i 2 MeV. They occur when no BH is formed in the merger, i.e., when it leads to a massive NS. Their U-GRB electromagnetically inferred isotropic occurrence rate is ρ S GRF 3.6 1.0 + 1.4 Gpc 3 year 1 [24].
(3) Ultrashort gamma-ray flashes (U-GRFs): in [20], it has been advanced a new class short bursts, the ultrashort GRBs (U-GRBs) produced by NS-BH binaries when the merger leaves the central BH with very little or completely without surrounding matter. An analogous system could be produced in BNS mergers. We shall call these systems ultrashort GRFs, for short U-GRFs. Their gamma-ray emission is expected to occur in a prompt short radiation phase. The post-merger radiation is drastically reduced, given the absence of baryonic matter to power an extended emission. A kilonova can still be observed days after the merger, in the infrared, optical, and ultraviolet wavelengths, produced by the radioactive decay of r-process yields [29,30,31,32]. Kilonova models used a dynamical ejecta composed of matter expelled by tides prior or during the merger, and a disk-wind ejecta by matter expelled from post-merger outflows in accretion disks [33], so U-GRFs are expected to have only the dynamical ejecta kilonova emission.
We focus on the BNSs left by BdHNe II and discuss how their properties impact the subsequent merger process and the associated short GRB emission, including their GW radiation. Since an accretion disk around the central remnant of a BNS merger, i.e., a newborn NS or a BH, is an important ingredient in models of short GRBs (see, e.g., [34] and references therein), we give some emphasis to the conditions and consequences for the merger leaving a disk. We study BNSs formed through binary evolution channels. Specifically, we expect these systems to form following a binary evolution channel similar to that of two massive stars leading to stripped-envelope binaries, described in previous studies (e.g., [35,36]). In this process, the CO star undergoes mass loss in multiple mass-transfer and common-envelope phases through interactions with the NS companion (see, e.g., [37,38,39]). This leads to removing the H/He layers of the secondary star, which ends up as a CO star. Recently, significant progress has been made in the study of alternative evolution channels for the progenitor of BNSs, such as hierarchical systems involving triple and quadrupole configurations [40,41], which are motivated by the presence of massive stars in multiple systems [42]. These systems are out of the scope of this study.
The article is organized as follows. In Section 2, we discuss the numerical simulations of BdHNe and specialize in an example of a BNS led by a BdHN II. Section 3 introduces a theoretical framework to analyze the BNS merger outcome configuration properties based on the conservation laws of baryon number, angular momentum, and mass-energy. We present in Section 4 a specific example analyzing a BNS merger using the above-mentioned theoretical framework, including estimates of the energy and angular momentum release. We include the radiation in gravitational waves (GWs) and estimate its detection by current facilities. Section 5 presents a summary and the conclusions of this work.

2. A BNS Left by a BdHN II

Figure 1 shows a snapshot of the mass density with the vector velocity field at the binary’s equatorial plane some minutes after the CO collapse and the expansion of the SN ejecta. The system’s evolution was simulated with an SPH code, where the NS companion and the ν NS are point particles that interact gravitationally with the SPH particles of the SN ejecta. For details of these numerical simulations, we refer to [23,43]. In these simulations, the influence of the star’s magnetic field has be disregarded, as the magnetic pressure remains significantly lower than the random pressure exerted on the infalling material. The simulation of Figure 1 corresponds to a CO-NS for a CO star evolved from a zero-age main-sequence (ZAMS) star of M zams = 15 M . The CO mass is about 3.06 M , whose core collapse leaves a 1.4 M ν NS and ejects 1.66 M . The NS companion’s initial mass is 1.4 M , and the initial binary period of the system is about 4.5 min.
From the accretion rate on the NSs, we have calculated the evolution of the mass and angular momentum of the binary components (see [43], for details). Table 1 summarizes the final parameters of the ν NS and the NS, including the gravitational mass, m, dimensionless angular momentum, j, angular velocity, Ω , equatorial radius, R eq , and moment of inertia, I. These structure parameters have been calculated with the RNS code [44] and using the GM1 [45,46] and TM1 [47] EOS (see Table 2 for details of the EOS). The BNS left by the BdHN II event has a period P orb = 14.97 min, orbital separation a orb 2 × 10 10 cm, and eccentricity e = 0.45 .

3. Inferences from Conservation Laws

We analyze the properties of the central remnant NS formed after the merger. We use the conservation laws of baryon number, energy, and angular momentum for this aim.

3.1. Baryon Number Conservation

The total baryonic mass of the system must be conserved, so the binary baryonic mass, M b , will redistribute among that of the postmerger’s central remnant, m b , c ; the ejecta’s mass, m ej , which is unbound to the system; and the matter kept bound to the system, e.g., in the form of a disk of mass m d . Therefore, we have the constraint
M b = m b , c + m ej + m d , M b = m b , 1 + m b , 2 .
For a uniformly rotating NS, the relation among its baryonic mass, m b , i , gravitational mass, m i , and angular momentum J i , is well represented by the simple function
m b , i M m i M + 13 200 m i M 2 1 1 130 j i 1.7 , i = 1 , 2 , c ,
where j i c J i / ( G M 2 ) , which fits numerical integration solutions of the axisymmetric Einstein equations for various nuclear EOS, with a maximum error of 2 % [48]. Thus, Equation (2) is a nearly universal, i.e., EOS-independent, formula. Equation (2) applies to the merging components ( i = 1 , 2 ) as well as to the central remnant ( i = c ).

3.2. Angular Momentum Conservation

We can make more inferences about the merger’s fate from the conservation of angular momentum. The angular momentum of the binary during the inspiral phase is given by
J = μ r 2 Ω + J 1 + J 2 , J i = 2 5 κ i m i R i 2 Ω i , i = 1 , 2 ,
where r is the orbital separation, μ = m 1 m 2 / M is the reduced mass, M = m 1 + m 2 is the total binary mass, and Ω = G M / r 3 is the orbital angular velocity. The gravitational mass and stellar radius of the i-th stellar component are, respectively, m i and R i ; J i is its angular momentum, Ω i its angular velocity, and κ i is the ratio between its moment of inertia to that of a homogeneous sphere. We adopt the convention m 2 m 1 . After the merger, the angular momentum is given by the sum of the angular momentum of the central remnant, the disk, and the ejecta. Angular momenta conservation implies that the angular momenta at merger, J merger , equals that of the final configuration plus losses:
J merger = J c + J d + Δ J ,
where J c and J d are, respectively, the angular momenta of the central remnant and the eventual surrounding disk, Δ J accounts for angular momentum losses, e.g., via gravitational waves, and we have neglected the angular momentum carried out by the ejecta since it is expected to have small mass ∼ 10 4 10 2 M . Simulations suggest that this ejecta comes from interface of the merger, where matter is squeezed and ejected perpendicular to the orbital plane, see, e.g., [49,50]. The definition of the merger point will be discussed below.
The angular momentum of the binary at the merger point is larger than the maximum value a uniformly rotating NS can attain, i.e., the angular momentum at the Keplerian/mass-shedding limit, J K . Thus, the remnant NS should evolve first through a short-lived phase that radiates the extra angular momentum over that limit and enters the rigidly rotating stability phase from the mass-shedding limit. Thus, we assume the remnant NS after that transition phase starts its evolution with angular momentum
J c = J K 0.7 G m c 2 c .
Equation (5) fits the angular momentum of the Keplerian sequence from full numerical integration of the Einstein equations and is nearly independent of the nuclear EOS (see, e.g., [48] and references therein). Therefore, the initial dimensionless angular momentum of the central remnant is
j c = c J c G M 2 0.7 m c M 2 .
We model the disk’s angular momentum as a ring at the remnant’s inner-most stable circular orbit (ISCO). Thus, we use the formula derived in Cipolletta et al. [51], which fits, with a maximum error of 0.3 % , the numerical results of the angular momentum per unit mass of a test particle circular orbit in the general relativistic axisymmetric field of a rotating NS. Within this assumption, the disk’s angular momentum is given by
J d = J ISCO G c m c m d 2 3 0.37 j c m c / M 0.85 .
Notice that Equation (7) reduces to the known result for the Schwarzschild metric for vanishing angular momentum, as it must. However, it differs from the result for the Kerr metric, which tells us that the Kerr metric does not describe the exterior spacetime of a rotating NS (see [51] for a detailed discussion).
The estimate of J merger requires the knowledge of the merger point, which depends on whether or not the binary secondary becomes noticeably deformed by the tidal forces. When the binary mass ratio q m 2 / m 1 is close or equal to 1, the stars are only deformed before the point of contact [52]. Therefore, for q 1 , we can assume the point of the merger as the point of contact
r merger r cont = ( C 2 + q C 1 ) ( 1 + q ) C 1 C 2 G M c 2 ,
where C 1 , 2 G m 1 , 2 / ( c 2 R 1 , 2 ) is the compactness of the BNS components.
When the masses are different, if we model the stars as Newtonian incompressible spheroids, there is a minimal orbital separation r ms , below which no equilibrium configuration is attainable, i.e., one star begins to shed mass to the companion due to the tidal forces. In this approximation, r ms 2.2 q 1 / 3 R 2 [53]. Numerical relativity simulations of BH-NS quasi-equilibrium states suggest that the mass-shedding occurs at a distance (see [54] and references therein) of
r ms ( 0.270 ) 2 / 3 q 1 / 3 R 2 .
Our analysis adopts the mass-shedding distance of Equation (9). For a system with q = 0.7 (similar mass ratio of the one in Table 1), we have found that the less-compact star begins to shed mass before the point of contact, independently of the EOS, which agrees with numerical relativity simulations. Consequently, for non-symmetric binaries q < 1 , we define the merging at the point as the onset of mass-shedding, r merger r ms .
Based on the above two definitions of merger point, Equations (8) and (9), the angular momentum at the merger is given by
J merger = ν C 2 + q C 1 ( 1 + q ) C 1 C 2 G M 2 c , q 1 , ν q 1 / 3 [ ( 1 + q ) C 2 ] 1 / 2 G M 2 c , q < 1 ,
where we have introduced the so-called symmetric mass-ratio parameter, ν q / ( 1 + q ) 2 .

3.3. Mass-Energy Conservation

The conservation of mass-energy before and after the merger implies the energy released equals the mass defect of the system, i.e.,
E GW + E other = Δ M c 2 = [ M ( m c + m ej + m d ) ] c 2 ,
where Δ M is the system’s mass defect. We have also defined E GW = E GW insp + E GW pm the total energy emitted in GWs in the inspiral regime, E GW insp , and in the merger and post-merger phases, E GW pm . The energy E other is radiated in channels different from the GW emission, e.g., electromagnetic (photons) and neutrinos.

4. A Specific Example of BNS Merger

We analyze the merger of the 1.505 + 1.404 M BNS in Table 1. For these component masses, the inferred orbital separation of a orb 2 × 10 10 cm and eccentricity e = 0.45 , the merger is expected to be driven by GW radiation on a timescale [55] of
τ GW = c 5 G 3 5 256 a orb 4 μ M 2 F ( e ) 73.15 kyr , F ( e ) = 48 19 1 g ( e ) 4 0 e g ( e ) 4 ( 1 e 2 ) 5 / 2 e ( 1 + 121 304 e 2 ) d e 0.44 ,
where g ( e ) = e 12 / 19 ( 1 e 2 ) 1 ( 1 + 121 e 2 / 304 ) 870 / 2299 .
From Equations (2), (7) and (10), and the conservation Equations (1), (4) and (11), we can obtain the remnant and disk’s mass as a function of the angular momentum losses, Δ J , as well as an estimate of the energy and angular momentum released in the cataclysmic event. We use the NS structure parameters obtained for the GM1 EOS and the TM1 EOS. The total gravitational mass of the system is M = m 1 + m 2 = 2.909 M , so using Equation (2), we obtain the total baryonic mass of the binary, M b = m b , 1 + m b , 2 3.184 M . The binary’s mass fraction is q = 0.933 , so we assume the merger starts at the contact point. With this, the angular momentum at the merger, as given by Equation (10), for the GM1 and TM1 EOS is, respectively, J merger 5.65 G M 2 / c and J merger 5.73 G M 2 / c .
Figure 2 shows that the disk’s mass versus the central remnant’s mass for selected values of the angular momentum loss for the two EOS. The figure shows the system’s final parameters lie between two limiting cases: zero angular momentum loss leading to maximal disk mass and maximal angular momentum loss leading to zero disk mass.

4.1. Maximal Disk Mass

We obtain the configuration corresponding to the maximum disk mass switching off angular momentum losses. Let us specialize in the GM1 EOS. By setting Δ J = 0 , the solution of the system of equations formed by the baryon number and angular conservation equations leads to the central remnant’s mass, m c = 2.697 M , and disk’s mass, m d = 0.073 M . This limiting case switches off the GW emission, so it also sets an upper limit to the energy released in mechanisms different than GWs. Thus, Equation (11) implies that E other = Δ M c 2 = [ M ( m c + m ej + m d ) ] c 2 ( M m c m d ) c 2 0.139 M c 2 2.484 × 10 53 erg of energy are carried out to infinity by a mechanism different than GWs and not accompanied by angular momentum losses.

4.2. Zero Disk Mass

The other limiting case corresponds when the angular momentum loss and the remnant mass are maximized, i.e., when no disk is formed (see Figure 2). By setting m d = 0 , the solution of the conservation equations leads to the maximum angular momentum loss, Δ J = 0.331 G M 2 / c , and the maximum remnant’s mass, m c = 2.756 M .
Thus, the upper limit to the angular momentum carried out by GWs is given by the maximum amount of angular momentum losses, i.e., Δ J GW 0.331 G M 2 / c . In the inspiral phase of the merger, the system releases
E GW insp G m 1 m 2 2 r cont = q C 1 C 2 M c 2 2 ( 1 + q ) ( C 2 + q C 1 ) + 1 2 j 1 | Ω 1 | + j 2 | Ω 2 | G M 2 c .
For the binary we are analyzing, E GW insp 0.0194 M c 2 0.0563 M c 2 1.0073 × 10 53 erg. The transitional non-axisymmetric object (e.g., triaxial ellipsoid) formed immediately after the merger mainly generates these GWs, and their emission ends when the stable remnant NS is finally formed. We can model such a rotating object as a compressible ellipsoid with a polytropic EOS of index n = 0.5 –1 [56]. The object will spin up by angular momentum loss to typical frequencies of 1.4 2.0 kHz. The energy emitted in GWs is E GW pm 0.0079 M c 2 1.404 × 10 52 erg. Therefore, the energy released in GWs is, E GW = E GW insp + E GW pm 0.0642 M c 2 1.147 × 10 53 erg. If no disk is formed, i.e., for a U-GRF, the mass-energy defect is Δ M c 2 = [ M ( m c + m ej ) ] c 2 ( M m c ) c 2 0.153 M c 2 2.734 × 10 53 erg. This implies that E other = Δ M c 2 E GW 0.089 M c 2 1.591 × 10 53 erg are released in forms of energy different than GW radiation.
Therefore, combining the above two results, we conclude that for the present merger, assuming the GM1 EOS, the merger releases 0 < E GW 1.147 × 10 53 erg in GWs and 1.591 × 10 53 E other < 2.484 × 10 53 erg are released in other energy forms. The energy observed in short GRBs and further theoretical analysis, including numerical simulations of the physical processes occurring during the merger, will clarify the efficiency of converting E other into observable radiation. Since no BH is formed (in this GM1 EOS analysis), the assumption that the merger leads to an S-GRF suggests an efficiency lower than 10 % .
We now estimate the detection efficiency of the GW radiation released by the system in the post-merger phase when angular momentum losses are maximized, i.e., in the absence of a surrounding disk. We find the root-sum-squared strain of the signal, i.e.,
h rss = 2 | h ˜ + | 2 + | h ˜ × | 2 d f 1 π d f ¯ G E GW pm c 3 ,
where h ˜ + and h ˜ × are the Fourier transforms of the GW polarizations, d is the distance to the source, f ¯ is the mean GW frequency in the postmerger phase. These signals are expected to be detected with a 50 % of efficiency by the LIGO/Virgo pipelines [57] when h rss 10 22 Hz 1 / 2 [58]. For the energy release in the post-merger phase, we have f ¯ = 1671.77 Hz, so these signals could be detected up to a distance of d 10 Mpc.

5. Discussion and Conclusions

As some BdHN I and II systems remain bound after the GRB-SN event, the corresponding NS-BH and BNS systems, driven by GW radiation, will merge and lead to short GRBs. For a few minutes binary, the merger time is of the order of 10 4 year. This implies that the binaries will still be close to the long GRB site by the merger time, which implies a direct link between long and short GRBs [20].
The occurrence rate of long and short bursts, however, should differ as the SN explosion likely disrupts the binaries with long orbital periods. We are updating our previous analysis on this interesting topic reported in [59]. We refer the reader to Bianco et al. [60] for a preliminary discussion.
As a proof of concept, this article examined this unique connection between long and short GRBs predicted by the BdHN scenario, emphasizing the case of mergers of BNS left by BdHNe II. For this particular case, the simulations predict that the outcome system will be a NSs binary with the star spins anti-aligned. The application of the present theoretical framework to the analysis of other merging binaries, such as the BH-NS binaries produced by BdHN I (see [20] for a general discussion), will be addressed in a separate work.
We have carried out a numerical SPH simulation of a BdHN II occurring in a CO-NS of orbital period 4.5 min. The mass of the CO is 3.06 M and that of the NS companion, 1.4 M . The CO is the pre-SN star obtained from a ZAMS star of M zams = 15 M simulated from MESA code. The SPH simulation follows [23,43]. It computes the accretion rate onto the ν NS (left by the CO core collapse) and the NS companion while the ejecta expands within the binary. For the event that left a ν NS-NS eccentric binary of 1.505 + 1.404 M , orbital separation 2 × 10 10 cm, orbital period of ≈15 min and eccentricity e = 0.45 . The SN ejecta matter forms a disk around both stars with opposite spins, so we expect that the ν -NS binary will also have anti-aligned spins as well. The above parameters suggest the BNS merger leading to a short GRB occurs in ≈73 kyear after the BdHN II event.
Whether or not the central remnant of the BNS merger will be a Kerr BH or a massive, fast-rotating NS depends on the nuclear EOS. For instance, we have shown the GM1 EOS leads to the latter while the TM1 EOS leads to the former. As an example of the theoretical framework presented in this article, we quantify the properties of the merger using the GM1 EOS. We infer the mass of the NS central remnant and the surrounding disk as a function of the angular momentum losses. We then emphasize the merger features in the limiting cases of maximum and zero angular momentum loss, corresponding to a surrounding disk’s absence or maximum mass. We estimated the maximum energy and angular momentum losses in GWs. We showed that the post-merger phase could release up to ≈ 10 52 erg in ≈ 1.7 kHz GWs, and LIGO/Virgo could, in principle, detect such emissions for sources up to ≈10 Mpc. We assessed that up to a few 10 53 erg of energy could be released in other forms of energy, so a ≲ 10 % of efficiency of its conversion into observable electromagnetic radiation would lead to an S-GRF.
The direct link between long and short GRB progenitors predicted by the BdHN model opens the way to exciting astrophysical developments. For instance, the relative rate of BdHNe I and II and S-GRBs and S-GRFs might give crucial information on the nuclear EOS of NSs and the CO-NS parameters. At the same time, this information provides clues for the stellar evolution path of the binary progenitors leading to the CO-NS binaries of the BdHN scenario. Although challenging because of their expected ultrashort duration, observing a U-GRF would also be relevant for constraining the EOS of NS matter. An extended analysis is encouraged, including additional BNS parameters obtained from SPH simulations of BdHNe for various CO-NS systems and nuclear EOS.

Author Contributions

Conceptualization, L.M.B., C.F., J.F.R., J.A.R. and R.R.; methodology, L.M.B., J.F.R. and J.A.R.; formal analysis, L.M.B., J.F.R. and J.A.R.; investigation, L.M.B., J.F.R. and J.A.R.; writing—original draft, L.M.B., C.F., J.F.R., J.A.R. and R.R.; writing—review and editing, L.M.B., J.F.R. and J.A.R. All authors have read and agreed to the published version of the manuscript.

Funding

L.M.B. is supported by the Vicerrectoría de Investigación y Extensión—Universidad Industrial de Santander Postdoctoral Fellowship Program No. 2023000359. J.F.R. has received funding/support from the Patrimonio Autónomo—Fondo Nacional de Financiamiento para la Ciencia, la Tecnología y la Innovación Francisco José de Caldas (MINCIENCIAS–COLOMBIA) Grant No. 110685269447 RC-80740-465-2020, Project No. 69553.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BdHNBinary-driven hypernova
BHBlack hole
BNSBinary neutron star
COCarbon-oxygen
EOSEquation of state
GRBGamma-ray burst
GWGravitational wave
ISCOInnermost stable circular orbit
NSNeutron star
ν NSNewborn neutron star
S-GRBShort gamma-ray burst
S-GRFShort gamma-ray flash
SNSupernova
U-GRBUltrashort gamma-ray burst
U-GRFUltrashort gamma-ray flash
ZAMSZero-age main-sequence

References

  1. Mazets, E.P.; Golenetskii, S.V.; Ilinskii, V.N.; Panov, V.N.; Aptekar, R.L.; Gurian, I.A.; Proskura, M.P.; Sokolov, I.A.; Sokolova, Z.I.; Kharitonova, T.V. Catalog of cosmic gamma-ray bursts from the KONUS experiment data. I. Astrophys. Space Sci. 1981, 80, 3–83. [Google Scholar] [CrossRef]
  2. Klebesadel, R.W. The durations of gamma-ray bursts. In Gamma-Ray Bursts—Observations, Analyses and Theories; Ho, C., Epstein, R.I., Fenimore, E.E., Eds.; Cambridge University Press: Cambridge, UK, 1992; pp. 161–168. [Google Scholar]
  3. Dezalay, J.P.; Barat, C.; Talon, R.; Syunyaev, R.; Terekhov, O.; Kuznetsov, A. Short cosmic events—A subset of classical GRBs? In American Institute of Physics Conference Series; Paciesas, W.S., Fishman, G.J., Eds.; AIP Publishing: Melville, NY, USA, 1992; Volume 265, pp. 304–309. [Google Scholar]
  4. Kouveliotou, C.; Meegan, C.A.; Fishman, G.J.; Bhat, N.P.; Briggs, M.S.; Koshut, T.M.; Paciesas, W.S.; Pendleton, G.N. Identification of two classes of gamma-ray bursts. Astrophys. J. 1993, 413, L101–L104. [Google Scholar] [CrossRef]
  5. Tavani, M. Euclidean versus Non-Euclidean Gamma-Ray Bursts. Astrophys. J. 1998, 497, L21–L24. [Google Scholar] [CrossRef] [Green Version]
  6. Goodman, J. Are gamma-ray bursts optically thick? Astrophys. J. 1986, 308, L47. [Google Scholar] [CrossRef]
  7. Paczynski, B. Gamma-ray bursters at cosmological distances. Astrophys. J. 1986, 308, L43–L46. [Google Scholar] [CrossRef]
  8. Eichler, D.; Livio, M.; Piran, T.; Schramm, D.N. Nucleosynthesis, neutrino bursts and gamma-rays from coalescing neutron stars. Nature 1989, 340, 126–128. [Google Scholar] [CrossRef]
  9. Narayan, R.; Piran, T.; Shemi, A. Neutron star and black hole binaries in the Galaxy. Astrophys. J. 1991, 379, L17–L20. [Google Scholar] [CrossRef]
  10. Woosley, S.E. Gamma-ray bursts from stellar mass accretion disks around black holes. Astrophys. J. 1993, 405, 273–277. [Google Scholar] [CrossRef]
  11. Mészáros, P. Theories of Gamma-Ray Bursts. Annu. Rev. Astron. Astrophys. 2002, 40, 137. [Google Scholar] [CrossRef] [Green Version]
  12. Piran, T. The physics of gamma-ray bursts. Rev. Mod. Phys. 2004, 76, 1143–1210. [Google Scholar] [CrossRef] [Green Version]
  13. Galama, T.J.; Vreeswijk, P.M.; van Paradijs, J.; Kouveliotou, C.; Augusteijn, T.; Bohnhardt, H.; Brewer, J.P.; Doublier, V.; Gonzalez, J.-F.; Leibundgut, B.; et al. An unusual supernova in the error box of the γ-ray burst of 25 April 1998. Nature 1998, 395, 670–672. [Google Scholar] [CrossRef]
  14. Woosley, S.E.; Bloom, J.S. The Supernova Gamma-Ray Burst Connection. Annu. Rev. Astron. Astrophys. 2006, 44, 507–556. [Google Scholar] [CrossRef] [Green Version]
  15. Della Valle, M. Supernovae and Gamma-Ray Bursts: A Decade of Observations. Int. J. Mod. Phys. D 2011, 20, 1745–1754. [Google Scholar] [CrossRef]
  16. Hjorth, J.; Bloom, J.S. The Gamma-Ray Burst—Supernova Connection. In Gamma-Ray Bursts; Cambridge Astrophysics Series; Kouveliotou, C., Wijers, R.A.M.J., Woosley, S., Eds.; Cambridge University Press: Cambridge, UK, 2012; Volume 51, Chapter 9; pp. 169–190. [Google Scholar]
  17. Rueda, J.A.; Ruffini, R. On the Induced Gravitational Collapse of a Neutron Star to a Black Hole by a Type Ib/c Supernova. Astrophys. J. 2012, 758, L7. [Google Scholar] [CrossRef] [Green Version]
  18. Izzo, L.; Rueda, J.A.; Ruffini, R. GRB 090618: A candidate for a neutron star gravitational collapse onto a black hole induced by a type Ib/c supernova. Astron. Astrophys. 2012, 548, L5. [Google Scholar] [CrossRef] [Green Version]
  19. Fryer, C.L.; Rueda, J.A.; Ruffini, R. Hypercritical Accretion, Induced Gravitational Collapse, and Binary-Driven Hypernovae. Astrophys. J. 2014, 793, L36. [Google Scholar] [CrossRef] [Green Version]
  20. Fryer, C.L.; Oliveira, F.G.; Rueda, J.A.; Ruffini, R. Neutron-Star-Black-Hole Binaries Produced by Binary-Driven Hypernovae. Phys. Rev. Lett. 2015, 115, 231102. [Google Scholar] [CrossRef] [Green Version]
  21. Becerra, L.; Cipolletta, F.; Fryer, C.L.; Rueda, J.A.; Ruffini, R. Angular Momentum Role in the Hypercritical Accretion of Binary-driven Hypernovae. Astrophys. J. 2015, 812, 100. [Google Scholar] [CrossRef] [Green Version]
  22. Becerra, L.; Bianco, C.L.; Fryer, C.L.; Rueda, J.A.; Ruffini, R. On the Induced Gravitational Collapse Scenario of Gamma-ray Bursts Associated with Supernovae. Astrophys. J. 2016, 833, 107. [Google Scholar] [CrossRef] [Green Version]
  23. Becerra, L.; Ellinger, C.L.; Fryer, C.L.; Rueda, J.A.; Ruffini, R. SPH Simulations of the Induced Gravitational Collapse Scenario of Long Gamma-Ray Bursts Associated with Supernovae. Astrophys. J. 2019, 871, 14. [Google Scholar] [CrossRef] [Green Version]
  24. Ruffini, R.; Rueda, J.A.; Muccino, M.; Aimuratov, Y.; Becerra, L.M.; Bianco, C.L.; Kovacevic, M.; Moradi, R.; Oliveira, F.G.; Pisani, G.B.; et al. On the Classification of GRBs and Their Occurrence Rates. Astrophys. J. 2016, 832, 136. [Google Scholar] [CrossRef]
  25. Ruffini, R.; Muccino, M.; Kovacevic, M.; Oliveira, F.G.; Rueda, J.A.; Bianco, C.L.; Enderli, M.; Penacchioni, A.V.; Pisani, G.B.; Wang, Y.; et al. GRB 140619B: A short GRB from a binary neutron star merger leading to black hole formation. Astrophys. J. 2015, 808, 190. [Google Scholar] [CrossRef] [Green Version]
  26. Ruffini, R.; Muccino, M.; Aimuratov, Y.; Bianco, C.L.; Cherubini, C.; Enderli, M.; Kovacevic, M.; Moradi, R.; Penacchioni, A.V.; Pisani, G.B.; et al. GRB 090510: A Genuine Short GRB from a Binary Neutron Star Coalescing into a Kerr-Newman Black Hole. Astrophys. J. 2016, 831, 178. [Google Scholar] [CrossRef] [Green Version]
  27. Ruffini, R.; Muccino, M.; Aimuratov, Y.; Amiri, M.; Bianco, C.L.; Chen, Y.C.; Eslam Panah, B.; Mathews, G.J.; Moradi, R.; Pisani, G.B.; et al. On the universal GeV emission in short GRBs. arXiv 2018, arXiv:1802.07552. [Google Scholar]
  28. Aimuratov, Y.; Ruffini, R.; Muccino, M.; Bianco, C.L.; Penacchioni, A.V.; Pisani, G.B.; Primorac, D.; Rueda, J.A.; Wang, Y. GRB 081024B and GRB 140402A: Two Additional Short GRBs from Binary Neutron Star Mergers. Astrophys. J. 2017, 844, 83. [Google Scholar] [CrossRef] [Green Version]
  29. Li, L.X.; Paczyński, B. Transient Events from Neutron Star Mergers. Astrophys. J. 1998, 507, L59–L62. [Google Scholar] [CrossRef] [Green Version]
  30. Metzger, B.D.; Martínez-Pinedo, G.; Darbha, S.; Quataert, E.; Arcones, A.; Kasen, D.; Thomas, R.; Nugent, P.; Panov, I.V.; Zinner, N.T. Electromagnetic counterparts of compact object mergers powered by the radioactive decay of r-process nuclei. Mon. Not. R. Astron. Soc. 2010, 406, 2650–2662. [Google Scholar] [CrossRef]
  31. Tanvir, N.R.; Levan, A.J.; Fruchter, A.S.; Hjorth, J.; Hounsell, R.A.; Wiersema, K.; Tunnicliffe, R.L. A ‘kilonova’ associated with the short-duration γ-ray burst GRB 130603B. Nature 2013, 500, 547–549. [Google Scholar] [CrossRef] [Green Version]
  32. Berger, E.; Fong, W.; Chornock, R. An r-process Kilonova Associated with the Short-hard GRB 130603B. Astrophys. J. 2013, 774, L23. [Google Scholar] [CrossRef] [Green Version]
  33. Metzger, B.D. Kilonovae. Living Rev. Relativ. 2017, 20, 3. [Google Scholar] [CrossRef] [Green Version]
  34. Berger, E. Short-Duration Gamma-Ray Bursts. Annu. Rev. Astron. Astrophys. 2014, 52, 43–105. [Google Scholar] [CrossRef] [Green Version]
  35. Tauris, T.M.; Langer, N.; Podsiadlowski, P. Ultra-stripped supernovae: Progenitors and fate. Mon. Not. R. Astron. Soc. 2015, 451, 2123–2144. [Google Scholar] [CrossRef] [Green Version]
  36. Tauris, T.M.; Kramer, M.; Freire, P.C.C.; Wex, N.; Janka, H.T.; Langer, N.; Podsiadlowski, P.; Bozzo, E.; Chaty, S.; Kruckow, M.U.; et al. Formation of Double Neutron Star Systems. Astrophys. J. 2017, 846, 170. [Google Scholar] [CrossRef] [Green Version]
  37. Nomoto, K.; aki Hashimoto, M. Presupernova evolution of massive stars. Phys. Rep. 1988, 163, 13–36. [Google Scholar] [CrossRef]
  38. Kim, H.J.; Yoon, S.C.; Koo, B.C. Observational Properties of Type Ib/c Supernova Progenitors in Binary Systems. Astrophys. J. 2015, 809, 131. [Google Scholar] [CrossRef] [Green Version]
  39. Yoon, S.C. Evolutionary Models for Type Ib/c Supernova Progenitors. Publ. Astron. Soc. Aust. 2015, 32, e015. [Google Scholar] [CrossRef] [Green Version]
  40. Hamers, A.S.; Fragione, G.; Neunteufel, P.; Kocsis, B. First- and second-generation black hole and neutron star mergers in 2+2 quadruples: Population statistics. Mon. Not. R. Astron. Soc. 2021, 506, 5345–5360. [Google Scholar] [CrossRef]
  41. Vynatheya, P.; Hamers, A.S. How Important Is Secular Evolution for Black Hole and Neutron Star Mergers in 2+2 and 3+1 Quadruple-star Systems? Astrophys. J. 2022, 926, 195. [Google Scholar] [CrossRef]
  42. Sana, H.; de Mink, S.E.; de Koter, A.; Langer, N.; Evans, C.J.; Gieles, M.; Gosset, E.; Izzard, R.G.; Le Bouquin, J.B.; Schneider, F.R.N. Binary Interaction Dominates the Evolution of Massive Stars. Science 2012, 337, 444. [Google Scholar] [CrossRef] [Green Version]
  43. Becerra, L.M.; Moradi, R.; Rueda, J.A.; Ruffini, R.; Wang, Y. First minutes of a binary-driven hypernova. Phys. Rev. D 2022, 106, 083002. [Google Scholar] [CrossRef]
  44. Stergioulas, N.; Friedman, J.L. Comparing Models of Rapidly Rotating Relativistic Stars Constructed by Two Numerical Methods. Astrophys. J. 1995, 444, 306. [Google Scholar] [CrossRef]
  45. Glendenning, N.K.; Moszkowski, S.A. Reconciliation of neutron-star masses and binding of the Lambda in hypernuclei. Phys. Rev. Lett. 1991, 67, 2414–2417. [Google Scholar] [CrossRef] [PubMed]
  46. Pal, S.; Bandyopadhyay, D.; Greiner, W. Antikaon condensation in neutron stars. Nucl. Phys. A 2000, 674, 553–577. [Google Scholar] [CrossRef] [Green Version]
  47. Sugahara, Y.; Toki, H. Relativistic mean-field theory for unstable nuclei with non-linear σ and ω terms. Nucl. Phys. A 1994, 579, 557–572. [Google Scholar] [CrossRef]
  48. Cipolletta, F.; Cherubini, C.; Filippi, S.; Rueda, J.A.; Ruffini, R. Fast rotating neutron stars with realistic nuclear matter equation of state. Phys. Rev. D 2015, 92, 023007. [Google Scholar] [CrossRef] [Green Version]
  49. Oechslin, R.; Janka, H.T.; Marek, A. Relativistic neutron star merger simulations with non-zero temperature equations of state. I. Variation of binary parameters and equation of state. Astron. Astrophys. 2007, 467, 395–409. [Google Scholar] [CrossRef] [Green Version]
  50. Bauswein, A.; Goriely, S.; Janka, H.T. Systematics of Dynamical Mass Ejection, Nucleosynthesis, and Radioactively Powered Electromagnetic Signals from Neutron-star Mergers. Astrophys. J. 2013, 773, 78. [Google Scholar] [CrossRef] [Green Version]
  51. Cipolletta, F.; Cherubini, C.; Filippi, S.; Rueda, J.A.; Ruffini, R. Last stable orbit around rapidly rotating neutron stars. Phys. Rev. D 2017, 96, 024046. [Google Scholar] [CrossRef] [Green Version]
  52. Bernuzzi, S. Neutron star merger remnants. Gen. Relativ. Gravit. 2020, 52, 108. [Google Scholar] [CrossRef]
  53. Chandrasekhar, S. Ellipsoidal Figures of Equilibrium; Dover: New York, NY, USA, 1969. [Google Scholar]
  54. Shibata, M.; Taniguchi, K. Coalescence of Black Hole-Neutron Star Binaries. Living Rev. Relat. 2011, 14, 6. [Google Scholar] [CrossRef] [Green Version]
  55. Maggiore, M. Gravitational Waves: Volume 1: Theory and Experiments; Oxford University Press: Oxford, UK, 2007. [Google Scholar]
  56. Lai, D.; Shapiro, S.L. Gravitational radiation from rapidly rotating nascent neutron stars. Astrophys. J. 1995, 442, 259–272. [Google Scholar] [CrossRef] [Green Version]
  57. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Abernathy, M.R.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; et al. All-sky search for long-duration gravitational wave transients in the first Advanced LIGO observing run. Class. Quantum Gravity 2018, 35, 065009. [Google Scholar] [CrossRef] [Green Version]
  58. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; Adya, V.B.; et al. Search for Post-merger Gravitational Waves from the Remnant of the Binary Neutron Star Merger GW170817. Astrophys. J. 2017, 851, L16. [Google Scholar] [CrossRef] [Green Version]
  59. Ruffini, R.; Rodriguez, J.; Muccino, M.; Rueda, J.A.; Aimuratov, Y.; Barres de Almeida, U.; Becerra, L.; Bianco, C.L.; Cherubini, C.; Filippi, S.; et al. On the Rate and on the Gravitational Wave Emission of Short and Long GRBs. Astrophys. J. 2018, 859, 30. [Google Scholar] [CrossRef] [Green Version]
  60. Bianco, C.L.; Mirtorabi, M.T.; Moradi, R.; Rastegarnia, F.; Rueda, J.A.; Ruffini, R.; Wang, Y.; Della Valle, M.; Li, L.; Zhang, S.R. Probing electromagnetic-gravitational wave emission coincidence in type I binary-driven hypernova family of long GRBs at very-high redshift. arXiv 2023, arXiv:2306.05855. [Google Scholar] [CrossRef]
Figure 1. Massdensity snapshots and velocity field on the orbital plane of a BdHN for a CO left by a M zams = 15 M and a 1.4 M NS companion, with an initial orbital period of about 4.5 min. We follow the expansion of the SN ejecta in the presence of the NS companion and the ν NS with a smoothed particle hydrodynamic (SPH) code. It is clear that a disk with opposite spins has formed around both stars.
Figure 1. Massdensity snapshots and velocity field on the orbital plane of a BdHN for a CO left by a M zams = 15 M and a 1.4 M NS companion, with an initial orbital period of about 4.5 min. We follow the expansion of the SN ejecta in the presence of the NS companion and the ν NS with a smoothed particle hydrodynamic (SPH) code. It is clear that a disk with opposite spins has formed around both stars.
Universe 09 00332 g001
Figure 2. Disk mass versus central remnant (NS) mass. Selected values of the angular momentum loss (in units of G M 2 / c ) are shown as points. The initial BNS has a total gravitational mass of 2.909 M and a mass fraction q = 0.933 , so we assume the merger starts at the contact point. The maximum mass along the Keplerian sequence for the GM1 EOS is 2.84 M and for the TM1 EOS it is 2.62 M (see Table 2). Thus, for the former EOS, the central remnant is a massive fast-rotating NS, while the latter suggests a prompt collapse into a Kerr BH.
Figure 2. Disk mass versus central remnant (NS) mass. Selected values of the angular momentum loss (in units of G M 2 / c ) are shown as points. The initial BNS has a total gravitational mass of 2.909 M and a mass fraction q = 0.933 , so we assume the merger starts at the contact point. The maximum mass along the Keplerian sequence for the GM1 EOS is 2.84 M and for the TM1 EOS it is 2.62 M (see Table 2). Thus, for the former EOS, the central remnant is a massive fast-rotating NS, while the latter suggests a prompt collapse into a Kerr BH.
Universe 09 00332 g002
Table 1. BNS produced by a BdHN II originated in a CO-NS with an orbital period of 4.5 min. The CO star mass is 3.06 M , obtained from the stellar evolution of a ZAMS star of M zams = 15 M , and the NS companion has 1.4 M . The numerical smoothed-particle hydrodynamic (SPH) simulation follows the SN produced by the CO core collapse and estimates the accretion rate onto the ν NS and the NS companion. The structure parameters of the NSs are calculated for the GM1 and TM1 EOS. We refer to [43] for additional details.
Table 1. BNS produced by a BdHN II originated in a CO-NS with an orbital period of 4.5 min. The CO star mass is 3.06 M , obtained from the stellar evolution of a ZAMS star of M zams = 15 M , and the NS companion has 1.4 M . The numerical smoothed-particle hydrodynamic (SPH) simulation follows the SN produced by the CO core collapse and estimates the accretion rate onto the ν NS and the NS companion. The structure parameters of the NSs are calculated for the GM1 and TM1 EOS. We refer to [43] for additional details.
mj Ω R eq I Ω R eq I
[ M ] [s 1 ][km][g cm 2 ][s 1 ][km][g cm 2 ]
GM1 EOSTM1 EOS
ν NS 1.505 0.259 1114.6 14.03 2.04 × 10 45 1077.1 14.47 2.11 × 10 45
NS 1.404 0.011 52.14 14.01 1.85 × 10 45 56.6 14.49 1.93 × 10 45
Table 2. Properties of the selected EOS. From left to right: maximum stable mass of non-rotating configurations, uniformly rotating configurations, set by the maximum mass of the Keplerian/mass-shedding sequence and the corresponding angular velocity.
Table 2. Properties of the selected EOS. From left to right: maximum stable mass of non-rotating configurations, uniformly rotating configurations, set by the maximum mass of the Keplerian/mass-shedding sequence and the corresponding angular velocity.
EOS M max j = 0 M max j kep Ω kep max
[ M ] [ M ] [ s 1 ]
GM1 2.38 2.84 1.001 × 10 4
TM1 2.19 2.62 8.83 × 10 3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Becerra, L.M.; Fryer, C.; Rodriguez, J.F.; Rueda, J.A.; Ruffini, R. Neutron Star Binaries Produced by Binary-Driven Hypernovae, Their Mergers, and the Link between Long and Short GRBs. Universe 2023, 9, 332. https://doi.org/10.3390/universe9070332

AMA Style

Becerra LM, Fryer C, Rodriguez JF, Rueda JA, Ruffini R. Neutron Star Binaries Produced by Binary-Driven Hypernovae, Their Mergers, and the Link between Long and Short GRBs. Universe. 2023; 9(7):332. https://doi.org/10.3390/universe9070332

Chicago/Turabian Style

Becerra, Laura M., Chris Fryer, Jose F. Rodriguez, Jorge A. Rueda, and Remo. Ruffini. 2023. "Neutron Star Binaries Produced by Binary-Driven Hypernovae, Their Mergers, and the Link between Long and Short GRBs" Universe 9, no. 7: 332. https://doi.org/10.3390/universe9070332

APA Style

Becerra, L. M., Fryer, C., Rodriguez, J. F., Rueda, J. A., & Ruffini, R. (2023). Neutron Star Binaries Produced by Binary-Driven Hypernovae, Their Mergers, and the Link between Long and Short GRBs. Universe, 9(7), 332. https://doi.org/10.3390/universe9070332

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop