No Arabic abstract
We present a rapid analytic framework for predicting kilonova light curves following neutron star (NS) mergers, where the main input parameters are binary-based properties measurable by gravitational wave detectors (chirp mass and mass ratio, orbital inclination) and properties dependent on the nuclear equation of state (tidal deformability, maximum NS mass). This enables synthesis of a kilonova sample for any NS source population, or determination of the observing depth needed to detect a live kilonova given gravitational wave source parameters in low latency. We validate this code, implemented in the public MOSFiT package, by fitting it to GW170817. A Bayes factor analysis overwhelmingly ($B>10^{10}$) favours the inclusion of an additional luminosity source in addition to lanthanide-poor dynamical ejecta during the first day. This is well fit by a shock-heated cocoon model, though differences in the ejecta structure, opacity or nuclear heating rate cannot be ruled out as alternatives. The emission thereafter is dominated by a lanthanide-rich viscous wind. We find the mass ratio of the binary is $q=0.92pm0.07$ (90% credible interval). We place tight constraints on the maximum stable NS mass, $M_{rm TOV}=2.17^{+0.08}_{-0.11}$ M$_odot$. For a uniform prior in tidal deformability, the radius of a 1.4 M$_odot$ NS is $R_{1.4}sim 10.7$ km. Re-weighting with a prior based on equations of state that support our credible range in $M_{rm TOV}$, we derive a final measurement $R_{1.4}=11.06^{+1.01}_{-0.98}$ km. Applying our code to the second gravitationally-detected neutron star merger, GW190425, we estimate that an associated kilonova would have been fainter (by $sim0.7$ mag at one day post-merger) and declined faster than GW170817, underlining the importance of tuning follow-up strategies individually for each GW-detected NS merger.
In the past few years, new observations of neutron stars and neutron-star mergers have provided a wealth of data that allow one to constrain the equation of state of nuclear matter at densities above nuclear saturation density. However, most observations were based on neutron stars with masses of about 1.4 solar masses, probing densities up to $sim$ 3-4 times the nuclear saturation density. Even higher densities are probed inside massive neutron stars such as PSR J0740+6620. Very recently, new radio observations provided an update to the mass estimate for PSR J0740+6620 and X-ray observations by the NICER and XMM telescopes constrained its radius. Based on these new measurements, we revisit our previous nuclear-physics multi-messenger astrophysics constraints and derive updated constraints on the equation of state describing the neutron-star interior. By combining astrophysical observations of two radio pulsars, two NICER measurements, the two gravitational-wave detections GW170817 and GW190425, detailed modeling of the kilonova AT2017gfo, as well as the gamma-ray burst GRB170817A, we are able to estimate the radius of a typical 1.4-solar mass neutron star to be $11.94^{+0.76}_{-0.87} rm{km}$ at 90% confidence. Our analysis allows us to revisit the upper bound on the maximum mass of neutron stars and disfavours the presence of a strong first-order phase transition from nuclear matter to exotic forms of matter, such as quark matter, inside neutron stars.
The first detection of gravitational waves from a neutron star-neutron star merger, GW170817, has opened up a new avenue for constraining the ultradense-matter equation of state (EOS). The deviation of the observed waveform from a point-particle waveform is a sensitive probe of the EOS controlling the merging neutron stars structure. In this topical review, I discuss the various constraints that have been made on the EOS in the year following the discovery of GW170817. In particular, I review the surprising relationship that has emerged between the effective tidal deformability of the binary system and the neutron star radius. I also report new results that make use of this relationship, finding that the radius inferred from GW170817 lies between 9.8 and 13.2 km at 90% confidence, with distinct likelihood peaks at 10.8 and 12.3 km. I compare these radii, as well as those inferred in the literature, to X-ray measurements of the neutron star radius. I also summarize the various maximum mass constraints, which point towards a maximum mass < 2.3 M_sun, depending on the fate of the remnant, and which can be used to additionally constrain the high-density EOS. I review the constraints on the EOS that have been performed directly, through Bayesian inference schemes. Finally, I comment on the importance of disentangling thermal effects in future EOS constraints from neutron star mergers.
The very first detection of gravitational waves from a neutron star binary merger, GW170817, exceeded all expectations. The event was relatively nearby, which may point to a relatively high merger rate. It was possible to extract finite-size effects from the gravitational-wave signal, which constrains the nuclear equation of state. Also, an electromagnetic counterpart was detected at many wavebands from radio to gamma rays marking the begin of a new multi-messenger era involving gravitational waves. We describe how multi-messenger observations of GW170817 are employed to constrain the nuclear equation of state. Combining the information from the optical emission and the mass measurement through gravitational waves leads to a lower limit on neutron star radii. According to this conservative analysis, which employs a minimum set of assumptions, the radii of neutron stars with typical masses should be larger than about 10.7~km. This implies a lower limit on the tidal deformability of about 210, while much stronger lower bounds are not supported by the data of GW170817. The multi-messenger interpretation of GW170817 rules out very soft nuclear matter and complements the upper bounds on NS radii which are derived from the measurement of finite-size effects during the pre-merger phase. We highlight the future potential of multi-messenger observations and of GW measurements of the postmerger phase for constraining the nuclear equation of state. Finally, we propose an observing strategy to maximize the scientific yield of future multi-messenger observations.
Light axion fields, if they exist, can be sourced by neutron stars due to their coupling to nuclear matter, and play a role in binary neutron star mergers. We report on a search for such axions by analysing the gravitational waves from the binary neutron star inspiral GW170817. We find no evidence of axions in the sampled parameter space. The null result allows us to impose constraints on axions with masses below $10^{-11} {rm eV}$ by excluding the ones with decay constants ranging from $1.6times10^{16} {rm GeV}$ to $10^{18} {rm GeV}$ at $3sigma$ confidence level. Our analysis provides the first constraints on axions from neutron star inspirals, and rules out a large region in parameter space that has not been probed by the existing experiments.
The thermodynamical properties of the equation of state (EoS) of high-density matter (above nuclear saturation density) and the possible existence of exotic states such as phase transitions from nuclear/hadronic matter into quark-gluon plasma, or the appearance of hyperons, may critically influence the stability and dynamics of compact relativistic stars. From a theoretical point of view, establishing the existence of those states requires the analysis of the `convexity of the EoS. We show indications of the existence of regions in the dense-matter EoS where the thermodynamics may be non-convex as a result of a non-monotonic dependence of the sound speed with the rest-mass density. When this happens, non-conventional dynamics may develop. In this paper we investigate the effects of a phenomenological, non-convex EoS on the equilibrium structure of stable compact stars and on the dynamics of unstable neutron stars that collapse gravitationally to black holes, both for spherically symmetric and uniformly-rotating configurations. We show how the dynamics of the collapse with a non-convex EoS departs from the convex case, leaving distinctive imprints on the gravitational waveforms. The astrophysical significance of these results for microphysical EoSs is discussed.