Investigating the crust of neutron stars with neural-network quantum states

Bryce Fore,Jane Kim,Morten Hjorth-Jensen,Alessandro Lovato

2024-07-31

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$^{-3}$ and 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.

Nuclear Theory,High Energy Astrophysical Phenomena,Solar and Stellar Astrophysics

What problem does this paper attempt to address?

The paper aims to address the following key issues:

1. **Accurate Description of Low-Density Nuclear Matter**: To better understand the physical properties of neutron star crusts, precise modeling of low-density nuclear matter is required. Particularly, in the density range of approximately 0.01 fm\(^{-3}\) to 0.1 fm\(^{-3}\), matter transitions from neutron-rich nuclei to various higher-density "pasta-like" structures, eventually reaching a uniform liquid state.

2. **Overcoming Limitations of Traditional Methods**: Existing theoretical methods (such as the compressible liquid drop model, self-consistent mean-field models, etc.) have limitations in describing low-density nuclear matter, especially in capturing deviations caused by quantum mechanical shell effects and strong correlations.

3. **Developing a New Variational Monte Carlo Method**: Researchers propose a variational Monte Carlo method based on neural network quantum states (NQS), which can avoid the fermion sign problem and effectively simulate the transition from uniform phases to neutron-rich nuclei.

4. **Microscopic Simulation of Nuclear Cluster Formation**: Through this method, researchers can microscopically simulate the spontaneous formation and dissolution of nuclear clusters, which is crucial for understanding the complex structure of neutron star crusts.

5. **Improving Energy Calculations**: The work also aims to obtain lower energy values than the existing auxiliary field diffusion Monte Carlo (AFDMC) method, particularly in low-density regions, which helps in more accurately predicting the physical state within neutron star crusts.

6. **Exploring Properties of Nuclear Matter at Different Densities**: By simulating nuclear matter at different densities and proton fractions, the paper investigates symmetric nuclear matter (SNM), pure neutron matter (PNM), and isospin asymmetric cases in between, to reveal the formation mechanisms of nuclear clusters and their dependence on proton fraction.

In summary, the goal of this paper is to improve the theoretical description of low-density nuclear matter within neutron star crusts by introducing a new method based on neural network quantum states, thereby better understanding and predicting the physical behavior of neutron star crusts.