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.