In recent years, new astrophysical observations have provided a wealth of exciting input for nuclear physics. For example, the observations of two-solar-mass neutron stars put strong constraints on ...possible phase transitions to exotic phases in strongly interacting matter at high densities. Furthermore, the recent observation of a neutron-star merger in both the electromagnetic spectrum and gravitational waves has provided compelling evidence that neutron-star mergers are an important site for the production of extremely neutron-rich nuclei within the r-process. In the coming years, an abundance of new observations is expected, which will continue to provide crucial constraints on the nuclear physics of these events. To reliably analyze such astrophysical observations and extract information on nuclear physics, it is very important that a consistent approach to nuclear systems is used. Such an approach consists of a precise and accurate method to solve the nuclear many-body problem in nuclei and nuclear matter, combined with modern nuclear Hamiltonians that allow to estimate the theoretical uncertainties. Quantum Monte Carlo methods are ideally suited for such an approach and have been successfully used to describe atomic nuclei and nuclear matter. In this contribution, I will present a detailed description of Quantum Monte Carlo methods focusing on the application of these methods to astrophysical problems. In particular, I will discuss how to use Quantum Monte Carlo methods to describe nuclear matter of relevance to the physics of neutron stars.
Recent observations of neutron stars with gravitational waves and X-ray timing provide unprecedented access to the equation of state (EoS) of cold dense matter at densities difficult to realize in ...terrestrial experiments. At the same time, predictions for the EoS equipped with reliable uncertainty estimates from chiral effective field theory (xEFT) allow us to bound our theoretical ignorance. In this work, we analyze astrophysical data using a nonparametric representation of the neutron-star EoS conditioned on xEFT to directly constrain the underlying physical properties of the compact objects without introducing modeling systematics. We discuss how the data alone constrain the EoS at high densities when we condition on xEFT at low densities. Here, we also demonstrate how to exploit astrophysical data to directly test the predictions of xEFT for the EoS up to twice nuclear saturation density, in order to estimate the density at which these predictions might break down. We nd that the existence of massive pulsars, gravitational waves from GW170817, and NICER observations of PSR J0030+0451 favor xEFT predictions for the EoS up to nuclear saturation density over a more agnostic analysis by as much as a factor of 7 for the quantum Monte Carlo (QMC) calculations used in this work. While xEFT predictions using QMC are fully consistent with gravitational-wave data up to twice nuclear saturation density, NICER observations suggest that the EoS stiffens relative to these predictions at or slightly above nuclear saturation density. Additionally, for these QMC calculations, we marginalize over the uncertainty in the density at which xEFT begins to break down, constraining the radius of a 1.4M$\bigodot$ neutron star to R1.4 = $11.40^{+1.38}_{–1.04}$ ($12.54^{+0.71}_{–0.63}$) km and the pressure at twice nuclear saturation density to p(2nsat) = $14.2^{+18.1}_{–8.4}$ ($28.7^{+15.3}_{–15.0}$) MeV=fm3 with massive pulsar and gravitational-wave (and NICER) data.
We propose the existence of a lower bound on the energy of pure neutron matter (PNM) on the basis of unitary-gas considerations. We discuss its justification from experimental studies of cold atoms ...as well as from theoretical studies of neutron matter. We demonstrate that this bound results in limits to the density-dependent symmetry energy, which is the difference between the energies of symmetric nuclear matter and PNM. In particular, this bound leads to a lower limit to the volume symmetry energy parameter S0. In addition, for assumed values of S0 above this minimum, this bound implies both upper and lower limits to the symmetry energy slope parameter L,which describes the lowest-order density dependence of the symmetry energy. A lower bound on neutron-matter incompressibility is also obtained. These bounds are found to be consistent with both recent calculations of the energies of PNM and constraints from nuclear experiments. Our results are significant because several equations of state that are currently used in astrophysical simulations of supernovae and neutron star mergers, as well as in nuclear physics simulations of heavy-ion collisions, have symmetry energy parameters that violate these bounds. Furthermore, below the nuclear saturation density, the bound on neutron-matter energies leads to a lower limit to the density-dependent symmetry energy, which leads to upper limits to the nuclear surface symmetry parameter and the neutron-star crust-core boundary. We also obtain a lower limit to the neutron-skin thicknesses of neutron-rich nuclei. Above the nuclear saturation density, the bound on neutron-matter energies also leads to an upper limit to the symmetry energy, with implications for neutron-star cooling via the direct Urca process.
Observations of neutron-star mergers with distinct messengers, including gravitational waves and electromagnetic signals, can be used to study the behavior of matter denser than an atomic nucleus and ...to measure the expansion rate of the Universe as quantified by the Hubble constant. We performed a joint analysis of the gravitational-wave event GW170817 with its electromagnetic counterparts AT2017gfo and GRB170817A, and the gravitational-wave event GW190425, both originating from neutron-star mergers. We combined these with previous measurements of pulsars using x-ray and radio observations, and nuclear-theory computations using chiral effective field theory, to constrain the neutron-star equation of state. We found that the radius of a 1.4-solar mass neutron star is Formula: see text km at 90% confidence and the Hubble constant is Formula: see text at 1σ uncertainty.
To obtain an understanding of the structure and reactions of nuclear systems from first principles has been a long-standing goal of nuclear physics. In this respect, few- and many-body systems ...provide a unique laboratory for studying nuclear interactions. During the past decades, the development of accurate representations of the nuclear force has undergone substantial progress. Particular emphasis has been devoted to chiral effective field theory (EFT), a low-energy effective representation of quantum chromodynamics (QCD). Within chiral EFT, many studies have been carried out dealing with the construction of both the nucleon-nucleon (NN ) and three-nucleon (3N ) interactions. The aim of the present article is to give a detailed overview of the chiral interaction models that are local in configuration space, and show recent results for nuclear systems obtained by employing these local chiral forces.
Neutron stars are astrophysical objects of extremes, reaching the highest densities we can observe in the cosmos, and probing matter under conditions that cannot be recreated in terrestrial ...experiments. In August 2017, the first neutron-star merger has been observed, which provided compelling evidence that these events are an important site for r-process nucleosynthesis. Furthermore, the gravitational-wave signal of such events might shed light upon the nature of strongly interacting matter in the neutron-star core. To understand these remarkable events, reliable nuclear physics input is essential. In this contribution, I explain how to use chiral effective field theory and advanced many-body methods to provide a consistent and systematic approach to strongly inter- acting systems from nuclei to neutron stars with controlled theoretical uncertainties. I will discuss recent results for the equation of state relevant for the nuclear astrophysics of neutron stars and neutron-star mergers.
Abstract
In the past few years, new observations of neutron stars (NSs) and NS mergers have provided a wealth of data that allow one to constrain the equation of state (EOS) of nuclear matter at ...densities above nuclear saturation density. However, most observations were based on NSs with masses of about 1.4
M
⊙
, probing densities up to ∼three to four times the nuclear saturation density. Even higher densities are probed inside massive NSs 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 multimessenger astrophysics constraints and derive updated constraints on the EOS describing the NS 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 AT 2017gfo, and the gamma-ray burst GRB 170817A, we are able to estimate the radius of a typical 1.4
M
⊙
NS to be
11.94
−
0.87
+
0.76
km
at 90% confidence. Our analysis allows us to revisit the upper bound on the maximum mass of NSs and disfavors the presence of a strong first-order phase transition from nuclear matter to exotic forms of matter, such as quark matter, inside NSs.
We perform a joint Bayesian inference of neutron-star mass and radius constraints based on GW170817, observations of quiescent low-mass x-ray binaries (QLMXBs), photospheric radius expansion x-ray ...bursting sources, and x-ray timing observations of J0030+0451. With this dataset, the form of the prior distribution still has an impact on the posterior mass-radius curves and equation of state (EOS), but this impact is smaller than recently obtained when considering QLMXBs alone. We analyze the consistency of the electromagnetic data by including an "intrinsic scattering" contribution to the uncertainties, and find only a slight broadening of the posteriors. This suggests that the gravitational-wave and electromagnetic observations of neutron-star structure are providing a consistent picture of the neutron-star mass-radius curve and the EOS.
Abstract
The observation of a compact object with a mass of 2.50–2.67
M
⊙
on 2019 August 14, 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-multimessenger 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 ∼81%. 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
times the speed of light, while the binary black hole scenario does not offer any new insight.