An official website of the United States government
Abstract An accurate description of low-density nuclear matter is crucial for explaining the physics of neutron star crusts. In the density range between approximately 0.01 fm−3and 0.1 fm−3, matter transitions from neutron-rich nuclei to various higher-density pasta shapes, before ultimately reaching a uniform liquid. In this work, we introduce a variational Monte Carlo method based on a neural Pfaffian-Jastrow quantum state, which allows us to model the transition from the liquid phase to neutron-rich nuclei microscopically. At low densities, nuclear clusters dynamically emerge from the microscopic interactions among protons and neutrons, which we model based on pionless effective field theory. Our variational Monte Carlo approach represents a significant improvement over the state-of-the-art auxiliary-field diffusion Monte Carlo method, which is severely hindered by the fermion-sign problem in this low-density regime and cannot capture the onset of clusters. In addition to computing the energy per particle of symmetric nuclear matter and pure neutron matter, we analyze an intermediate isospin-asymmetry configuration to elucidate the formation of nuclear clusters. We also provide evidence that the presence of such nuclear clusters influences the amount of protons in the crust compared to protons in beta-equilibrated, neutrino-transparent matter.
more »
« less
Low Density Neutron Star Matter with Quantum Molecular Dynamics: The Role of Isovector Interactions
The effect of isospin-dependent nuclear forces on the inner crust of neutron stars is modeled within the framework of Quantum Molecular Dynamics (QMD). To successfully control the density dependence of the symmetry energy of neutron-star matter below nuclear saturation density, a mixed vector-isovector potential is introduced. This approach is inspired by the baryon density and isospin density-dependent repulsive Skyrme force of asymmetric nuclear matter. In isospin-asymmetric nuclear matter, the system shows nucleation, as nucleons are arranged into shapes resembling nuclear pasta. The dependence of clusterization in the system on the isospin properties is also explored by calculating two-point correlation functions. We show that, as compared to previous results that did not involve such mixed interaction terms, the energy symmetry slope L is successfully controlled by varying the corresponding coupling strength. Nevertheless, the effect of changing the slope of the nuclear symmetry energy L on the crust-core transition density does not seem significant. To the knowledge of the authors, this is the first implementation of such a coupling in a QMD model for isospin asymmetric matter, which is relevant to the inner crust of neutron and proto-neutron stars.
more »
« less
- Award ID(s):
- 1748621
- PAR ID:
- 10418341
- Date Published:
- Journal Name:
- Universe
- Volume:
- 8
- Issue:
- 7
- ISSN:
- 2218-1997
- Page Range / eLocation ID:
- 380
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The formation of clusters at sub-saturation densities, as a result of many-body correlations, constitutes an essential feature for a reliable modelization of the nuclear matter equation of state (EoS). Phenomenological models that make use of energy density functionals (EDFs) offer a convenient approach to account for the presence of these bound states of nucleons when introduced as additional degrees of freedom. However, in these models clusters dissolve, by construction, when the nuclear saturation density is approached from below, revealing inconsistencies with recent findings that evidence the existence of short-range correlations (SRCs) even at larger densities. The idea of this work is to incorporate SRCs in established models for the EoS, in light of the importance of these features for the description of heavy-ion collisions, nuclear structure and in the astrophysical context. Our aim is to describe SRCs at supra-saturation densities by using effective quasi-clusters immersed in dense matter as a surrogate for correlations, in a regime where cluster dissolution is usually predicted in phenomenological models. Within the EDF framework, we explore a novel approach to embed SRCs within a relativistic mean-field model with density dependent couplings through the introduction of suitable in-medium modifications of the cluster properties, in particular their binding energy shifts, which are responsible for describing the cluster dissolution. As a first exploratory step, the example of a quasi-deuteron within the generalized relativistic density functional approach is investigated. The zero temperature case is examined, where the deuteron fraction is given by the density of a boson condensate. For the first time, suitable parameterizations of the cluster mass shift at zero temperature are derived for all baryon densities. They are constrained by experimental results for the effective deuteron fraction in nuclear matter near saturation and by microscopic many-body calculations in the low-density limit. A proper description of well-constrained nuclear matter quantities at saturation is kept through a refit of the nucleon meson coupling strengths. The proposed parameterizations allow to also determine the density dependence of the quasi-deuteron mass fraction at arbitrary isospin asymmetries. The strength of the deuteron-meson couplings is assessed to be of crucial importance. Novel effects on some thermodynamic quantities, such as the matter incompressibility, the symmetry energy and its slope, are finally discerned and discussed. The findings of the present study represent a first step to improve the description of nuclear matter and its EoS at supra-saturation densities in EDFs by considering correlations in an effective way. In a next step, the single-particle momentum distributions in nuclear matter can be explored using proper wave functions of the quasi-deuteron in the medium. The momentum distributions are expected to exhibit a high-momentum tail, as observed in the experimental study of SRCs by nucleon knockout with high-energy electrons. This will be studied in a forthcoming publication with an extensive presentation of the theoretical method and the results.more » « less
-
We present two recent parametrizations of the equation of state (FSU2R and FSU2H models) that reproduce the properties of nuclear matter and finite nuclei, fulfill constraints on high-density matter stemming from heavy-ion collisions, produce 2 M_Sun neutron stars, and generate neutron star radii below 13 km. Making use of these equations of state, cooling simulations for isolated neutron stars are performed. We find that two of the models studied, FSU2R (with nucleons) and, in particular, FSU2H (with nucleons and hyperons), show very good agreement with cooling observations, even without including nucleon pairing. This indicates that cooling observations are compatible with an equation of state that produces a soft nuclear symmetry energy and, thus, generates small neutron star radii. Nevertheless, both schemes produce cold isolated neutron stars with masses above 1.8 M_Sun.more » « less
-
null (Ed.)We report a precision measurement of the parity-violating asymmetry APV in the elastic scattering of longitudinally polarized electrons from 208Pb. We measure APV=550±16(stat)±8(syst) parts per billion, leading to an extraction of the neutral weak form factor FW(Q2=0.00616 GeV2)=0.368±0.013. Combined with our previous measurement, the extracted neutron skin thickness is Rn−Rp=0.283±0.071 fm. The result also yields the first significant direct measurement of the interior weak density of 208Pb: ρ0W=−0.0796±0.0036(exp)±0.0013(theo) fm−3 leading to the interior baryon density ρ0b=0.1480±0.0036(exp)±0.0013(theo) fm−3. The measurement accurately constrains the density dependence of the symmetry energy of nuclear matter near saturation density, with implications for the size and composition of neutron stars.more » « less
-
For densities beyond nuclear saturation, there is still a large uncertainty in the equations of state (EoSs) of dense matter that translate into uncertainties in the internal structure of neutron stars. The MUSES Calculation Engine provides a free and open-source composable workflow management system, which allows users to calculate the EoSs of dense and hot matter that can be used, e.g., to describe neutron stars. For this work, we make use of two MUSES EoS modules, i.e., Crust Density Functional Theory and Chiral Mean Field model, with beta-equilibrium with leptons enforced in the Lepton module, then connected by the Synthesis module using different functions: hyperbolic tangent, generalized Gaussian, bump, and smoothstep. We then calculate stellar structure using the QLIMR module and discuss how the different interpolating functions affect our results.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3390/universe8070380
