Gravitational waves have been detected from the inspiral of a binary neutron-star, GW170817, which allowed constraints to be placed on the neutron star equation of state. The equation of state can be further constrained if gravitational waves from a post-merger remnant are detected. Post-merger waveforms are currently generated by numerical-relativity simulations, which are computationally expensive. Here we introduce a hierarchical model trained on numerical-relativity simulations, which can generate reliable post-merger spectra in a fraction of a second. Our spectra have mean fitting factors of 0.95, which compares to fitting factors of 0.76 and 0.85 between different numerical-relativity codes that simulate the same physical system. This method is the first step towards generating large template banks of spectra for use in post-merger detection and parameter estimation.

Determining the differential-rotation law of compact stellar objects produced in binary neutron stars mergers or core-collapse supernovae is an old problem in relativistic astrophysics. Addressing this problem is important because it impacts directly on the maximum mass these objects can attain and hence on the threshold to black-hole formation under realistic conditions. Using the results from a large number of numerical simulations in full general relativity of binary neutron star mergers described with various equations of state and masses, we study the rotational properties of the resulting hypermassive neutron stars. We find that the angular-velocity distribution shows only a modest dependence on the equation of state, thus exhibiting the traits of "quasi-universality" found in other aspects of compact stars, both isolated and in binary systems. The distributions are characterized by an almost uniformly rotating core and a "disk". Such a configuration is significantly different from the $j-{\rm constant}$ differential-rotation law that is commonly adopted in equilibrium models of differentially rotating stars. Furthermore, the rest-mass contained in such a disk can be quite large, ranging from $\simeq 0.03\,M_{\odot}$ in the case of high-mass binaries with stiff equations of state, up to $\simeq 0.2\,M_{\odot}$ for low-mass binaries with soft equations of state. We comment on the astrophysical implications of our findings and on the long-term evolutionary scenarios that can be conjectured on the basis of our simulations.

We show how gravitational-wave observations with advanced detectors of tens to several tens of neutron-star binaries can measure the neutron-star radius with an accuracy of several to a few percent, for mass and spatial distributions that are realistic, and with none of the sources located within 100 Mpc. We achieve such an accuracy by combining measurements of the total mass from the inspiral phase with those of the compactness from the postmerger oscillation frequencies. For estimating the measurement errors of these frequencies we utilize analytical fits to postmerger numerical-relativity waveforms in the time domain, obtained here for the first time, for four nuclear-physics equations of state and a couple of values for the mass. We further exploit quasi-universal relations to derive errors in compactness from those frequencies. Measuring the average radius to well within 10% is possible for a sample of 100 binaries distributed uniformly in volume between 100 and 300 Mpc, so long as the equation of state is not too soft or the binaries are not too heavy.

The most popular model for short gamma-ray bursts (sGRBs) involves the coalescence of binary neutron stars. Because the progenitor is actually hidden from view, we must consider under which circumstances such merging systems are capable of producing a successful sGRB. Soon after coalescence, winds are launched from the merger remnant. In this paper, we use realistic wind profiles derived from global merger simulations in order to investigate the interaction of sGRB jets with these winds using numerical simulations. We analyze the conditions for which these axisymmetric winds permit relativistic jets to breakout and produce a sGRB. We find that jets with luminosities comparable to those observed in sGRBs are only successful when their half-opening angles are below ~20{\deg}. This jet collimation mechanism leads to a simple physical interpretation of the luminosities and opening angles inferred for sGRBs. If wide, low luminosity jets are observed, they might be indicative of a different progenitor avenue such as the merger of a neutron star with a black hole. We also use the observed durations of sGRB to place constraints on the lifetime of the wind phase, which is determined by the time it takes the jet to breakout. In all cases we find that the derived limits argue against completely stable remnants for binary neutron star mergers that produce sGRBs.

A number of works have shown that important information on the equation of state of matter at nuclear density can be extracted from the gravitational waves emitted by merging neutron-star binaries. We present a comprehensive analysis of the gravitational-wave signal emitted during the inspiral, merger and post-merger of 56 neutron-star binaries. This sample of binaries, arguably the largest studied to date with realistic equations of state, spans across six different nuclear-physics equations of state and ten masses, allowing us to sharpen a number of results recently obtained on the spectral properties of the gravitational-wave signal. Overall we find that: (i) for binaries with masses differing no more than $20\%$, the frequency at gravitational-wave amplitude's maximum is related quasi-universally with the tidal deformability of the two stars; (ii) the spectral properties vary during the post-merger phase, with a transient phase lasting a few millisecond after the merger and followed by a quasi-stationary phase; (iii) when distinguishing the spectral peaks between these two phases, a number of ambiguities in the identification of the peaks disappear, leaving a simple and robust picture; (iv) using properly identified frequencies, quasi-universal relations are found between the spectral features and the properties of the neutron stars; (v) for the most salient peaks analytic fitting functions can be obtained in terms of the stellar tidal deformability or compactness. Altogether, these results support the idea that the equation of state of nuclear matter can be constrained tightly when a signal in gravitational waves from binary neutron stars is detected.

Extending previous work by a number of authors, we have recently presented a new approach in which the detection of gravitational waves from merging neutron star binaries can be used to determine the equation of state of matter at nuclear density and hence the structure of neutron stars. In particular, after performing a large number of numerical-relativity simulations of binaries with nuclear equations of state, we have found that the post-merger emission is characterized by two distinct and robust spectral features. While the high-frequency peak was already shown to be associated with the oscillations of the hypermassive neutron star produced by the merger and to depend on the equation of state, we have highlighted that the low-frequency peak is related to the merger process and to the total compactness of the stars in the binary. This relation is essentially universal and provides a powerful tool to set tight constraints on the equation of state. We here provide additional information on the extensive analysis performed, illustrating the methods used, the tests considered, as well as the robustness of the results. We also discuss additional relations that can be deduced when exploring the data and how these correlate with various properties of the binary. Finally, we present a simple mechanical toy model that explains the main spectral features of the post-merger signal and can even reproduce analytically the complex waveforms emitted right after the merger.

Determining the equation of state of matter at nuclear density and hence the structure of neutron stars has been a riddle for decades. We show how the imminent detection of gravitational waves from merging neutron star binaries can be used to solve this riddle. Using a large number of accurate numerical-relativity simulations of binaries with nuclear equations of state, we find that the postmerger emission is characterized by two distinct and robust spectral features. While the high-frequency peak has already been associated with the oscillations of the hypermassive neutron star produced by the merger and depends on the equation of state, a new correlation emerges between the low-frequency peak, related to the merger process, and the total compactness of the stars in the binary. More importantly, such a correlation is essentially universal, thus providing a powerful tool to set tight constraints on the equation of state. If the mass of the binary is known from the inspiral signal, the combined use of the two frequency peaks sets four simultaneous constraints to be satisfied. Ideally, even a single detection would be sufficient to select one equation of state over the others. We test our approach with simulated data and verify it works well for all the equations of state considered.

Inspiralling compact binaries as standard sirens will soon become an invaluable tool for cosmology when advanced interferometric gravitational-wave detectors begin their observations in the coming years. However, a degeneracy in the information carried by gravitational waves between the total rest-frame mass $M$ and the redshift $z$ of the source implies that neither can be directly extracted from the signal, but only the combination $M(1+z)$, the redshifted mass. Recent work has shown that for binary neutron star systems, a tidal correction to the gravitational-wave phase in the late-inspiral signal that depends on the rest-frame source mass could be used to break the mass-redshift degeneracy. We propose here to use the signature encoded in the post-merger signal to deduce the redshift to the source. This will allow an accurate extraction of the intrinsic rest-frame mass of the source, in turn permitting the determination of source redshift and luminosity distance solely from gravitational-wave observations. This will herald a new era in precision cosmography and astrophysics. Using numerical simulations of binary neutron star mergers of very slightly different mass, we model gravitational-wave signals at different redshifts and use Bayesian parameter estimation to determine the accuracy with which the redshift can be extracted for a source of known mass. We find that the Einstein Telescope can determine the source redshift to $sim 10$--$20%$ at redshifts of $z<0.04$.

Determining the conditions under which a black hole can be produced is a long-standing and fundamental problem in general relativity. We use numerical simulations of colliding selfgravitating fluid objects to study the conditions of black-hole formation when the objects are boosted to ultrarelativistic speeds. Expanding on previous work, we show that the collision is characterized by a type-I critical behaviour, with a black hole being produced for masses above a critical value, M_c, and a partially bound object for masses below the critical one. More importantly, we show for the first time that the critical mass varies with the initial effective Lorentz factor <\gamma> following a simple scaling of the type M_c ~ K <\gamma>^{-1.0}, thus indicating that a black hole of infinitesimal mass is produced in the limit of a diverging Lorentz factor. Furthermore, because a scaling is present also in terms of the initial stellar compactness, we provide a condition for black-hole formation in the spirit of the hoop conjecture.

The stability properties of relativistic stars against gravitational collapse to black hole is a classical problem in general relativity. A sufficient criterion for secular instability was established by Friedman, Ipser and Sorkin (1988), who proved that a sequence of uniformly rotating barotropic stars is secularly unstable on one side of a turning point and then argued that a stronger result should hold: that the sequence should be stable on the opposite side, with the turning point marking the onset of secular instability. We show here that this expectation is not met. By computing in full general relativity the $F$-mode frequency for a large number of rotating stars, we show that the neutral-stability point, i.e., where the frequency becomes zero, differs from the turning point for rotating stars. Using numerical simulations we validate that the new criterion can be used to assess the dynamical stability of relativistic rotating stars.

We examine the propagation of collisionless particles emitted from a spherical shell to infinity. The number distribution at infinity, calculated as a function of the polar angle, exhibits a small deviation from uniformity. The number of particles moving from the polar region toward the equatorial plane is slightly larger than that of particles in the opposite direction, for an emission radius $ > 4.5M$ in extreme Kerr space-time. This means that the black hole spin exerts an anti-collimation effect on the particles stream propagating along the rotation axis. We also confirm this property in the weak field limit. The quadrupole moment of the central object produces a force toward the equatorial plane. For a smaller emission radius $r<4.5M$, the absorption of particles into the black hole, the non-uniformity and/or the anisotropy of the emission distribution become much more important.

We present a systematic study of the steep decay emission from gamma-ray bursts (GRBs) observed by the Swift X-Ray Telescope (XRT). In contrast to the analysis described in recent literature, we produce composite Burst Alert Telescope (BAT) and XRT light curves by extrapolating the XRT data (2-10 keV) into the BAT energy range (15-25 keV) rather than extrapolating the BAT data into the XRT energy band (0.3-10 keV). Based on the fits to the composite light curves, we have confirmed the existence of an exponential decay component which smoothly connects the BAT prompt data to the XRT steep decay for several GRBs. We also find that the XRT steep decay for some of the bursts can be well fit by a combination of a power-law with an exponential decay model. We discuss this exponential component within the frame work of both the internal and the external shock model.

We show that the jet structure of gamma-ray bursts (GRBs) can be investigated with the tail emission of the prompt GRB. The tail emission which we consider is identified as a steep-decay component of the early X-ray afterglow observed by the X-ray Telescope onboard Swift. Using a Monte Carlo method, we derive, for the first time, the distribution of the decay index of the GRB tail emission for various jet models. The new definitions of the zero of time and the time interval of a fitting region are proposed. These definitions for fitting the light curve lead us an unique definition of the decay index, which is useful to investigate the structure of the GRB jet. We find that if the GRB jet has a core-envelope structure, the predicted distribution of the decay index of the tail has a wide scatter and has multiple peaks, which cannot be seen for the case of the uniform and the Gaussian jet. Therefore, the decay index distribution tells us the information on the jet structure. Especially, if we observe events whose decay index is less than about 2, both the uniform and the Gaussian jet models will be disfavored according to our simulation study.

The observational diversity of optical emission, which coincides with prompt gamma-ray bursts (GRBs), has been discovered in the recent Swift era. We show that on the assumption of the synchrotron radiation for the observed energy range below the X-ray band, the observed diversity can be explained using the internal shock model by taking into account a high-latitude emission and the spectral change due to the synchrotron self-absorption. It may even be possible in our model to include bright optical flashes found, e.g., in GRB~990123. The prediction of our model is that the spectral index in the optical band is dependent on whether the optical light curve correlates with those in the X-rays and/or gamma-rays or not, which will be tested in the near future.