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 observat
ions 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.
At supranuclear densities, explored in the core of neutron stars, a strong phase transition from hadronic matter to more exotic forms of matter might be present. To test this hypothesis, binary neutron-star mergers offer a unique possibility to probe
matter at densities that we can not create in any existing terrestrial experiment. In this work, we show that, if present, strong phase transitions can have a measurable imprint on the binary neutron-star coalescence and the emitted gravitational-wave signal. We construct a new parameterization of the supranuclear equation of state that allows us to test for the existence of a strong phase transition and extract its characteristic properties purely from the gravitational-wave signal of the inspiraling neutron stars. We test our approach using a Bayesian inference study simulating 600 signals with three different equations of state and find that for current gravitational-wave detector networks already twelve events might be sufficient to verify the presence of a strong phase transition. Finally, we use our methodology to analyze GW170817 and GW190425, but do not find any indication that a strong phase transition is present at densities probed during the inspiral.
In recent years, the combination of advanced quantum Monte Carlo (QMC) methods and local interactions derived from chiral effective field theory (EFT) has been shown to provide a versatile and systematic approach to nuclear systems. Calculations at n
ext-to-next-to-leading order in chiral EFT have lead to fascinating results for nuclei and nucleonic matter. On the one hand, ground-state properties of nuclei are well reproduced up to $Aleq16$, even though these potentials have been fit to nucleon-nucleon scattering and few-body observables only. On the other hand, a reasonable description of neutron-matter properties emerges. While regulator functions applied to two- and three-nucleon forces are a necessary ingredient in these many-body calculations, the use of local regulators leads to a substantial residual regulator and cutoff dependence that increases current theoretical uncertainties. In this contribution, we review local chiral interactions, their applications, and QMC results for nuclei and neutron matter. In addition, we address regulator issues for such potentials and present a possible path forward.
