The observation of a compact object with a mass of $2.50-2.67M_{\odot}$ on August 14, 2019, by the LIGO Scientific and Virgo collaborations (LVC) has the potential to improve our understanding of the supranuclear equation of state. While the gravitational-wave analysis of the LVC suggests that GW190814 likely was a binary black hole system, the secondary component could also have been the heaviest neutron star observed to date. We use our previously derived nuclear-physics--multi-messenger astrophysics framework to address the nature of this object. Based on our findings, we determine GW190814 to be a binary black hole merger with a probability of $>99.9\%$. Even if we weaken previously employed constraints on the maximum mass of neutron stars, the probability of a binary black hole origin is still $\sim 86\%$. Furthermore, we study the impact that this observation has on our understanding of the nuclear equation of state by analyzing the allowed region in the mass-radius diagram of neutron stars for both a binary black hole or neutron star - black hole scenario. We find that the unlikely scenario in which the secondary object was a neutron star requires rather stiff equations of state with a maximum speed of sound $c_s\geq \sqrt{0.6}$ times the speed of light, while the binary black hole scenario does not offer any new insight.

The multi-messenger analysis of neutron-star mergers is not only of astrophysical interest, it also allows studying the behavior of supranuclear dense matter, probing the principles of general relativity, and measuring the expansion rate of the Universe. We perform a multi-messenger analysis of the gravitational-wave signal GW170817 and its electromagnetic counterparts AT2017gfo and GRB170817A, as well as the gravitational-wave signal GW190425. By incorporating information from the NICER observation of PSR J0030+0451, the radio observations of PSR J0740+6620, PSR J0348+4032, PSR J1614-2230, and nuclear-theory computations using chiral effective field theory, we provide new constraints on the neutron-star equation of state. Within our analysis, we determine the radius of a $1.4M_\odot$ NS to be $R_{1.4M_\odot}=11.74^{+0.98}_{-0.79}\ \rm km$ at $90\%$ confidence and estimate the Hubble constant to be $H_0=66.2^{+4.4}_{-4.2}\ \rm km /Mpc/s$ at $1\sigma$ uncertainty. The same multi-messenger-analysis techniques will be applicable for future observations, improving our understanding of the most extreme objects in the Universe.

The detection of GW170817 is revolutionizing many areas of astrophysics with the joint observation of gravitational waves and electromagnetic emissions. These multi-messenger events provide a new approach to determine the Hubble constant, thus, they are a promising candidate for mitigating the tension between measurements of Type Ia supernovae via the local distance ladder and the Cosmic Microwave Background. In addition to the "standard siren" provided by the gravitational-wave measurement, the kilonova itself has characteristics that allow to improve existing measurements or to perform yet another, independent measurement of the Hubble constant without gravitational-wave information. Here, we employ standardization techniques borrowed from the type-Ia community and apply them to kilonovae, not using any information from the gravitational-wave signal. We use two versions of this technique, one derived from direct observables measured from the lightcurve, and the other based on inferred ejecta parameters, e.g., mass, velocity, and composition, for two different models. These lead to Hubble Constant measurements of $H_0 = 109^{+49}_{-35}$\,km $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$ for the measured analysis, and $H_0 = 85^{+22}_{-17}$\,km $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$ and $H_0 = 79^{+23}_{-15}$\,km $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$ for the inferred analyses. This measurement has error bars within ~$\sim$\,2 to the gravitational-wave measurements ($H_0=74^{+16}_{-8}$\,km $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$), showing its promise as an independent constraint on $H_0$.