Svoboda | Graniru | BBC Russia | Golosameriki | Facebook
Next Article in Journal
A Novel Scalable Quantum Protocol for the Dining Cryptographers Problem
Previous Article in Journal
Development and Validation of a Compressible Reacting Gas-Dynamic Flow Solver for Supersonic Combustion
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Dynamics
Volume 4
Issue 1
10.3390/dynamics4010009
dynamics-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phase Diagram of Nuclear Pastas in Neutron Star Crusts

by
Jorge A. Muñoz
* and
Jorge A. López
*
Department of Physics, The University of Texas at El Paso, El Paso, TX 79968, USA
*
Authors to whom correspondence should be addressed.
Dynamics 2024, 4(1), 157-169; https://doi.org/10.3390/dynamics4010009
Submission received: 23 October 2023 / Revised: 23 January 2024 / Accepted: 25 January 2024 / Published: 22 February 2024
Browse Figures
Versions Notes

Abstract

:
Two neural networks were trained to predict, respectively, the Euler characteristic and the curvature of nuclear pastas in neutron star crust conditions generated by molecular dynamics simulations of neutron star matter with 0.1 < x < 0.5, 0.040 fm−3 < ρ < 0.085 fm−3 (0.68 × 1014 g/cm3 < ρ < 1.43 × 1014 g/cm3), and 0.2 MeV < T < 4.0 MeV, where x is proton content, the density is ρ , and the temperature is T. The predictions of the two networks were combined to determine the nuclear pasta phase that is thermodynamically stable at a given x, ρ , and T, and a three-dimensional phase diagram that extrapolated slightly the regions of existing molecular dynamics data was computed. The jungle gym and anti-jungle gym structures are prevalent at high temperature and low density, while the anti-jungle gym and anti-gnocchi structures dominate at high temperature and high density. A diversity of structures exist at low temperatures and intermediate density and proton content. The trained models used in this work are open access and available at a public repository to promote comparison to pastas obtained with other models.
Keywords:
neutron star matter; nuclear pasta; Minkowski functionals; molecular dynamics; nuclear phases; neural networks
PACS:
26.60.+c; 26.60.−c; 02.70.Ns; 21.65.+f

1. Introduction

The structure of neutron star crusts is important for understanding neutrino flow [1] and the formation of nuclei in neutron star mergers [2]. Recent studies have shown that “pasta” structures can be formed both in nuclear matter (NM) [3] and neutron star matter (NSM) [4]. Those studies found that such structures are formed through dynamical processes, and vary depending on the proton content, density, and temperature. For a recent review see Ref. [5]. Here we continue a previous study [6,7] to relate the shape of the pasta structures to thermodynamic state variables. The hundreds of pastas analyzed were obtained through classical molecular dynamics (CMD) simulations in [4]. We use machine learning (ML) to learn the relationship between the Minkowski functionals of the different structure morphologies to the thermodynamics state variables of the CMD simulations that generated them.

1.1. Astrophysical Caveats

As this is not an article about astrophysics, but about a technique to categorize the nuclear pastas (that may exist in astrophysical objects), we believe it is convenient to state the range of applicability of the present work.
The composition of neutron star crusts is a complex and fascinating topic in astrophysics. The crust of a neutron star is a thin layer of matter that is about 1 kilometer thick, made up of a variety of elements, including iron, nickel, calcium, etc. forming a crystal lattice. The inner crust starts with Fe 56 , which has a proton fraction of Z/A = 26/56 = 0.46, and ends with Kr 118 , with a proton fraction of Z/A = 36/118 = 0.3 [8]. The drop of proton fraction with larger densities is due to the increase of electron capture by protons producing a large neutron excess and a correspondingly low proton fraction.
Thus, as it will be explained below and in Appendix A, the classical molecular dynamics model uses protons and neutrons and, as is, cannot be used to study the crust composed of a crystal lattice made of iron, nickel, calcium, etc.
At greater depths, that is, between the stellar core and the outer crust, there exists a transition region in which the density increases up to saturation density, and reaches higher values at the core. In this region, when the density is about half of the nuclear density, the nuclear system begins to undergo cluster formation, forming regions with nuclear density (clusters) and low-density regions (dilute neutron gas); such a system is known as the pasta. Since most field theories are unable to study clustering and other critical phenomena, the exact density of the crust-to-core transition has not been determined; Thomas-Fermi calculations place it at 0.067 fm−3 (or about 1.13 × 10 14 g/cm3) [9], while molecular dynamics simulations put it at about 0.08 fm−3 (or about 1.34 × 10 14 g/cm3) [5].
Similarly, the value of the proton fraction in this pasta region is also not well established; However, it is believed that the fraction of neutrons in the vapor and in the pastas is dictated by the nuclear symmetry energy [10]. Unfortunately, the symmetry energy is not firmly constrained by nuclear mass models (see the seriousness of this problem by a comparison of hundreds of models in Refs. [11,12]). The problem in this region is even more complicated for clustered systems, as the mean field results predict E s y m → 0 at low densities, whereas the experimental results [13] (and CMD calculations [5]) indicate E s y m → 5 to 10 MeV.
Thus, models that misestimate the value of E s y m m at low densities predict proton ratios between 1 % and 10 % (e.g., Monte Carlo simulations [10]), or between 45 % at very low densities (e.g., Brueckner-Thomas-Fermi model within a energy density functional). Whereas CMD, which yields correct values for E s y m m at low densities, predicts an unbound nuclear system at values of Z / A ≲ 20 % , with the possibility of finding pastas at higher values of Z / A [5].
In a different astrophysical scenario, the conditions inside a supernova are very similar to some of the conditions inside a neutron star, and nuclear pasta may also form inside supernovae. In this scenario, the pasta has to be warm and neutrinos that are trapped could play a role in setting the proton fraction [14]. Supernovae can be expected to have proton fractions of about 0.33 at high temperatures and at relatively low densities [15].
Another handicap of the CDM model is the impossibility of producing new particles through reactions, such as mesons and hyperons. Furthermore, when neutrons are spin-paired in neutron stars, they experience a decrease in energy due to the formation of Cooper pairs which leads to superfluidity; this phenomenon is not taken into account by CMD. However, it is expected that neutron pairing will not affect the pasta structures, as the energy change that occurs when neutrons pair is of the order of 1 MeV per neutron, much smaller than the binding energies of ≈ 10 MeV per nucleon in the pasta structures.
In conclusion, in this work we study pastas with densities between 0.02 fm−3 ≤ ρ ≤ 0.12 fm−3 (or 0.34 × 1014 g/cm3 ≤ ρ ≤ 2.02 × 10 14 g/cm3), temperatures 0.2 MeV ≤ 4 MeV, and Z / A between 0.5 and 0.3. [By extrapolation these ranges will be extended to 0.015 fm−3 ≤ ρ ≤ 0.12 fm−3, 0.0 MeV ≤ T ≤ 5 MeV, and 0 ≤ Z / A ≤ 0.5 ]. This hopefully covers the ranges of interest for both neutron star crusts and supernovae.

1.2. Terminology

“Nuclear matter” is composed of protons and neutrons. Likewise, “neutron star matter” is nuclear matter as it would be found on the crust of a neutron star, i.e., embedded in a gas of electrons. The “neutron star crust” is the outermost layer of a neutron star, with a thickness of about 1 km, where the density varies from zero to normal nuclear density. “Saturation density” is the density at which nuclear matter is at its most stable state, at zero temperature, ρ 0 = 0.16 fm−3 (2.74 × 1014 g/cm3). The “neutrino flow” in a neutron star is important for cooling and for transporting energy from the core to the crust. The “pasta” is a variety of exotic structures that can form in NM and NSM, and are called as such after their resemblance to pasta shapes, such as spaghetti, gnocchi, and lasagna. The “Minkowski functionals” of each structure are their respective values of the Euler characteristic, the volume and area occupied by the pastas, and the mean integral curvature. Finally, the “proton content” x = Z / A is the fraction of protons Z to the total number of nucleons A, which includes both neutrons and protons.

1.3. Nuclear Pasta

Neutron stars are the dense cadavers of large stars that have undergone a supernova explosion. Due to the intense gravitational pull of their compact mass, neutron stars overpower the electron degeneracy pressure, forcing electrons to fuse with protons, forming neutrons. This results in a tightly packed sphere of nearly pure neutron matter with protons and electrons scattered throughout, occupying a space several thousand times smaller than the original star.
At the neutron star surface, the pressure is relatively low, allowing nuclei from helium to iron to exist without being crushed together. However, in the crust, as the density increases, nuclear attraction and Coulomb repulsion forces become comparable in strength, and the interplay between them, and the pressure in the crust, leads to the formation of what has been dubbed “nuclear pasta” due to their resemblance to various types of pasta. If nuclear pasta does indeed exist, it would be the universe’s strongest material [16].
Formally, nuclear pasta is a theorized arrangement of protons and neutrons that may exist in the crusts of neutron stars [17]. The study of structures of nuclear pasta at zero temperature began in the 1980s through static methods like the liquid drop model [17,18], Thomas-Fermi models [19] and mean field theories [20]. At higher temperatures, the dynamics of nuclear pastas have been studied through quantum molecular dynamics simulations [21,22,23] and classical potential models [24,25]. Recently, classical molecular dynamics (CMD), which can find local minima by cooling, have studied non-traditional phases [26].
Various properties of clustered media, such as nuclear pastas, and their phase transitions into different pasta structures [26], have been calculated using CMD. In particular, recent studies have investigated the pasta-like structures and phase transitions of NM and NSM at sub-saturation densities, low temperatures, and proton content from 0.1 to 0.5 [1,2,3,4,5,6].
As an illustration, Figure 1 shows characteristic nuclear matter pastas obtained with CMD for 2000 nucleons at different values of density and temperatures, and for an equal number of protons and neutrons. The colors used represent protons (orange) and neutrons (purple), except at T = 0.25 MeV, where some protons appear red and some neutrons blue.
The examples at 0.02 fm−3 (0.34 ×   10 14 g/cm3) vary from “spaghetti” (tubes) at T = 1.00 MeV to “gnocchi” (bundles) at T = 0.01 MeV. At 0.04 fm−3 (0.68 ×   10 14 g/cm3), the pastas go from “anti-gnocchi” (i.e., holes) at T = 1.00 MeV to “jungle gym” (cross-linked structures similar to a lattice of rods) at T = 0.01 MeV. These structures vary with the proton content and with the addition of the electron gas for NSM.

1.4. The Importance of This Study

The present work, although it is not an astrophysics study, it will develop a method to compare the structures obtained with various models at different temperatures, densities and proton content, that may apply to different astrophysical scenarios. The study will demonstrate the utility of the Minkowski functionals to characterize the nuclear pastas and determine phase changes among them. A public code based on machine learning technology will be available in a public repository.

2. Materials and Methods

2.1. Research Design

The goal of this project is to connect the structure of nuclear pastas to measurable quantities such as the Euler characteristic. This was achieved through the following steps:
(i)
Generate pastas.
(ii)
Compute the Minkowski functionals.
(iii)
Interpolate and extrapolate the values of the Minkowski functionals using neural networks.
(iv)
Plot the Minkowski functionals for the different types of pastas.
In this work we use the pastas obtained with CMD simulations in Refs. [6,7], compute the Minkowski functionals using cubic voxels, and use two neural networks to interpolate and extrapolate the value of the Euler characteristic and mean curvature for arbitrary values of proton content, temperature and density. Here we provide some details of the methods employed, more complete information can be found in Appendix A and Appendix B.

2.2. Classical Molecular Dynamics

CMD is a mature technique that has been applied to determine the pasta structures, both in NM and NSM [1,25,26]. In CMD, nucleons are treated as classical particles that interact through pair potentials. The dynamics of the nucleons are determined by their equations of motion. Other publications [25,26] describe the details and main advantages of CMD; the model optimized for NSM was introduced in Ref. [4]. The CMD model treats neutron star matter as a mixture of protons p, neutrons n, and electrons, with different potentials for the interactions between each pair of particles, p p , n n and n p . The electron gas is mimicked by a screening potential. Here, we use CMD to investigate pasta in NSM; the CMD model is presented in Appendix A.
To simulate neutron star matter, we used the LAMMPS code [27] with the potentials described in Appendix A. A = 4000 nucleons were prepared in a cubic cell under periodic boundary conditions with proton ratios of x = Z / A = 0.1 , 0.2, 0.3, 0.4, and 0.5, densities ranging between 0.04 fm−3 < ρ < 0.085 fm−3 (0.68 ×   10 14 g/cm3 < ρ < 1.43   ×   10 14 g/cm3), and temperatures going from T = 4.0 MeV down to 0.2 MeV. Systems were prepared with nucleons placed at random with a minimum distance of 0.01 fm, and with speeds distributed according to a Maxwell-Boltzmann distribution at a given temperature. The evolution of the system at a given density was studied by solving the equations of motion. The position, momenta, and energy of the nuclei were stored at fixed time-steps. The pressure, temperature, and density were also recorded. At the end, the resulting pasta structures at the various combinations of values of x, ρ and T were obtained.

2.3. Minkowski Functionals

The pasta, as any other geometrical figures, can be characterized by means of the Minkowski functionals [28]. The volume (V), surface area (S), Euler characteristic ( χ ), and integral mean curvature (B) can quantify the shape and connectivity of the nuclear structures. Previous studies [25,29,30,31] have used the Minkowski functionals to investigate the pasta. The Euler characteristic is a topological invariant that describes the shape of a geometric body with a topological invariant number, χ , that will not change if the body is stretched or bent. χ can be used to distinguish between different topological spaces. For example, a sphere has χ = 1 while a torus has χ = 0, so they cannot be topologically equivalent. A more thorough description can be found in Appendix B and in Refs. [3,25].
After the work of [25], it is known that the pasta structures can be classified in terms of the curvature and the Euler characteristic, as illustrated in Table 1. Notice that the term “anti” refers to the situation in which occupied regions are inverted into empty ones and vice versa. Likewise, the term “jungle gym” denotes a 3D trellis; the signs of these functionals are correlated with the different phases of the pasta [32]. In this work, pasta structures produced by CMD are “voxelized”, i.e., represented by “voxels”, or regular units of volume of d = 2.35 fm per side. The V, S, and B are then computed, and χ is obtained according to Equation (A4). See the appendix of Ref. [3] for complete details.

2.4. Neural Network Interpolation

To interpolate the value of the two Minkowski functionals, the integral mean curvature and Euler characteristic, for arbitrary values of temperature and density, we trained fully-connected feed-forward neural networks of each functional using the package scikit-learn [33]. We used a network architecture of 4 hidden layers with 50 neurons per layer and with a rectified linear unit (ReLU) activation function. The weights of the network were determined by means of the method Adam of stochastic gradient-based optimization [34]. The data were partitioned 80% (randomly selected) for training the network, and 20% used to test the performance of the model. The normalized mean absolute error (NMAE) and the root-mean-square error (NRMSE) for the selected models on the test and the full data set are shown in Table 2. In the CMD data, the range of the integral mean curvature was − 1894 ≤ B ≤ 1556 and the range of the Euler characteristics was − 1093 ≤ χ ≤ 396 .

3. Results

In the study ( i ) we use the 9600 pastas obtained with CMD simulations of the cooling of 4000 nucleons from an initial temperature of 4.0 MeV to different final temperatures, at different proton content and densities [6,7], ( i i ) compute the Minkowski functionals, and ( i i i ) interpolate the value of the Minkowski functionals using neural networks for arbitrary values of T, x and ρ . The end product of steps ( i ) and ( i i ) is a dataset of values for the two Minkowski functionals in a dense data grid; these values are then used for training neural networks. This procedure is used to identify the structure of any pasta based on its proton content, density and temperature, without having to resort to the cumbersome voxeling. The methodology is independent of the method used to produce the pasta.

Phase Diagram

Sixteen constant-density cuts of the three-dimensional T − ρ − x phase diagram of nuclear pastas in neutron star crusts predicted by our trained neural networks are shown in Figure 2. To generate the cuts, the predictions of each neural network were aligned with matplotlib [35] color maps: Greens, which is a sequential color map that starts with the color white and linearly moves towards green, was used for the mean curvature; seismic, which is a divergent color map with blue and red at the two extremes and white in the middle, was used for the Euler characteristic. Consequentially, blue is associated with the jungle-gym and anti-jungle gym structures (negative χ ); red with the gnocchi and anti-gnocchi structures (positive χ ); green with gnocchi, spaghetti, and jungle gym structures (positive B); and lighter colors are associated with anti-gnocchi, anti-spaghetti, and anti-jungle gym structures (negative B). Each cut was computed by making individual predictions at 500 points along the proton proportion axis, and 100 points along the temperature axis. Black squares show points at which either mean curvature or Euler characteristic change sign, and these changes were detected by scanning the data along both the proton proportion rows and temperature columns. The following predictions in the regions are extrapolations: ρ < 0.04 fm−3 (0.68 ×   10 14 g/cm3), ρ > 0.085 fm−3, (1.43 ×   10 14 g/cm3), x < 0.1 , T < 0.2 MeV, and T > 4 MeV; are extrapolations.

4. Discussion

It is easy to see that at T > 2.0 MeV most structures are jungle gym, anti-jungle gym or anti-gnocchi. The rest of the other pasta structures exist only at lower temperatures, T < 2.0 MeV, and higher proton proportions, 0.2 < x < 0.5 . The lasagna is the least prominent of the structures, but it corresponds to the region near the boundary between any anti-structure and the corresponding structure; also notice that the lasagna is equivalent to its anti-lasagna. It is reassuring that the phase changes occur in the same temperature-density region where the caloric curve changes slopes, cf. Figures 48–50 in Ref. [5]; the change of slope could indicate the pasta structure changes are of first order. It must be mentioned that our results go beyond previous studies of the pasta diagram performed with static Hartree-Fock at fixed proton content x = 0.3 and non-zero temperature [32], and for T = 0 and x = 0.5 [36], and with quantum molecular dynamics for T = 0 and x = 0.3 [37] for cold matter.
In Figure 3, we present a comparison of our (extrapolated) results for T = 0.01 MeV with those of Ref. [36] for zero temperature obtained with static Hartree-Fock calculations. Using the labels and colors of Ref. [36] and assuming that χ ∼ 0 ( B ∼ 0 ) when their values are within ± 12 percent ( ± 6 percent) of their maximum values in our training set, we denote gnocchi as SPH (for sphere), spaghetti (cylinder-like) as ROD1, planar meshes of orthogonal rods as ROD2, jungle gym (three-dimensional grids of rods) as ROD3, lasagna as SLAB, anti-spaghetti (rod-bubble) as ROD1B, anti-planar meshes of orthogonal rods (ROD2 bubble) as ROD2B, and anti-gnocchi (sphere bubble) as SPHB. Notice that the slab and jungle gym are symmetric under the exchange of matter and voids. The uniform phase, i.e., without the pasta structure, is denoted as UNI. The difference between the results obtained by the two models underscores the advantage that CMD has in finding meta-stable minima by dynamical means. Other recent studies [38] show the same effect, obtaining large differences in pasta structures as a function of the generalized Skyrme effective interactions used. It is worth mentioning that the variety of pasta phases exists, say, in the density range of up to 0.08 fm−3 (1.43 × 1014 g/cm3), in agreement with [9].
In passing, we remark that other classical dynamic models, c.f. [24,30,31] use different potentials, but produce the same type of pastas, although not necessarily in the same order in density and temperature. Likewise, it is worth mentioning that superfluidity, although relevant at T ≤ 1 MeV and ρ ≲ ρ 0 / 3 [39] has not been shown to have any effect on the structure of the pastas.
To promote the comparison of the structures obtained with various models at different temperatures, densities and proton content, the code used and ML models used in this work are available in a public repository [40] under the name “Dangelo”, in honor of the Italian painter Sergio Dangelo, founder of the arte nucleare movement, and who could very well have painted the phase diagrams of Figure 2. The program produces phase diagrams of nuclear pastas using the neural network models trained on the classical molecular dynamics data presented in this work. The models are open access as well.
These findings are of interest to study the structure of the crust of neutron stars as a function of depth. In future work, we will combine these results with a Lane–Emden polytropic equation of state appropriate for a neutron star to determine the structure, i.e., the type of pasta, that exists at different depths, i.e., at different densities, temperatures and proton content, from the surface of the neutron star.

5. Conclusions

In this study we go beyond previous investigations of the structures attained by neutron star matter at varying conditions of T, ρ and x [6,7]. Using molecular dynamics simulations, 9600 pastas were produced, analyzed with respect to their Minkowski functionals, and the Minkowski functionals used to train neural network models, which were then used to explore the pasta configuration of the 3D space T − ρ − x . We found that the phase diagram of neutron star matter corresponds, in qualitative terms, to phase changes detected through the caloric curve in previous studies.
We conclude that machine learning techniques help estimate the Minkowski functionals for values of T, x, and ρ . The evolution of the pasta structures as a function of x, T and ρ is linked to the 3D representation of the functionals, that is, as χ goes from positive to negative, the structure goes from anti-pasta to pasta, etc. The linkage of the pasta structures to these functionals will be studied in a forthcoming article for neutron star matter and nuclear matter.
We encourage other investigators to compare the structures obtained with other models with those obtained with CMD by using the code “Dangelo” [40].

Author Contributions

Conceptualization, J.A.L. and J.A.M.; methodology, J.A.L. and J.A.M.; software, J.A.M.; validation, J.A.L. and J.A.M.; formal analysis, J.A.L. and J.A.M.; investigation, J.A.L. and J.A.M.; resources, J.A.L. and J.A.M.; data curation, J.A.L. and J.A.M.; writing—original draft preparation, J.A.L.; writing—review and editing, J.A.L. and J.A.M.; visualization, J.A.L. and J.A.M.; supervision, J.A.L. and J.A.M.; project administration, J.A.L. and J.A.M.; funding acquisition, J.A.L. and J.A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This material is based upon work supported by the U.S. Department of Energy, Office of Science, Reaching a New Energy Sciences Workforce (RENEW) program. J.A.M. was supported by RENEW Award Number DE-SC-0021994 and J.A.L. by RENEW Award Number CW46261 via Oak Ridge National Laboratory.

Data Availability Statement

The data presented in this study are available on request from the corresponding authors.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this article:
CMDClassical molecular dynamics
MDMolecular dynamics
MLMachine learning
NMNuclear matter
NSMNeutron star matter
NMAENormalized mean absolute error
ReLURectified linear unit
NRMSENormalized root mean square error

Appendix A. Molecular Dynamics

Appendix A.1. Why Use Classical Molecular Dynamics

In nuclear physics theoretical models are indispensable tools for comprehending heavy-ion collisions at intermediate energies. The majority of researchers employ models where nucleons (or “test particles”) interact through an average mean field that allows them to incorporate some quantum and semi-quantum features. Unfortunately, such an approach sacrifices very important features of the reactions, such as clusterization, fragment formation, phase transitions, and other crucial phenomena.
Conversely, classical dynamics models, while failing to capture quantum effects, excel at preserving nucleon-nucleon correlations, statistical fluctuations, phase transitions, and other critical phenomena essential to understanding the later stages of fragmentation reactions.
In a paradoxical interplay, the quantum aspects of mean field models facilitate a realistic energy distribution in the early stages of the reaction but falter in fragmenting the system in the later stages. On the other hand, classical methods do not represent accurately the reaction dynamics in the early part, but excel in the latter phase of the collision.
Furthermore, after six decades of evolution of the independent mean field models, the community of practitioners has not agreed on a single model to use. Indeed, calculations from the two fundamental families, Quantum Molecular Dynamics (QMD) and Boltzmann-Uhling-Uhlenbeck (BUU), still exhibit grave inconsistencies in their predictions. A comparison of different models revealed variations in flow observables, such as the effectiveness of Pauli blocking in nucleon-nucleon collisions, the stability of initialized nuclei, and others, primarily attributed to differences in the initialization of the system, collision integral treatment, and other intertwined effects [41].
The latter discrepancies likely stem from the operator’s selection of “test particles” and the fluctuations manually introduced in both QMD and the BUU Langevin framework. While fragmentation was not examined in this comparison, it remains a critical challenge for all such models [42,43,44,45,46]. To compare model predictions to final yields, the resulting matter density distributions must be fed into “afterburner” computer programs such as the Statistical Multifragmentation Model [47] and similar tools.
In contrast, despite not respecting quantum principles, classical dynamics models (such as CMD) can capture both the equilibrium and non-equilibrium phases of a collision. CMD models can indeed describe hydrodynamic flow, non-equilibrium dynamics, and phase transitions with non-adjustable parameters. They can also simulate phase transitions (critical phenomena), neck fragmentation [48], obtain isoscaling and caloric curves in reactions, as well as nuclear pasta formation in infinite systems [4,5,49].
Thus, despite its failures in accounting for quantum effects (the limitations of which have been assessed in [50]), we believe we err less in studying nuclear systems using classical molecular dynamics than other models that do not include the many-body correlations necessary to study phase changes and a plethora of critical phenomena. Furthermore, it is clear that this study could have been performed with pastas produced by other methods; unfortunately, no other technique can generate the thousands of pastas needed to train a neural network at such a large range of T, x and ρ , and without a large consumption of computer time. We understand that CMD is not the choice of most colleagues, and we ask them to have an open mind, and accept that other models can be used to explore nuclear phenomena.

Appendix A.2. The Classical Molecular Dynamics Model

In classical molecular dynamics (CMD) nucleons interact through a modified version of the Pandharipande Medium potential. The Pauli repulsion can be simulated in classical mechanics via a neutron-proton, attractive at large separations, and repulsive at small ones [51]. The potentials are phenomenological, constructed with combinations of Yukawa interactions, with parameters adjusted as to approximate the behavior of the Vlasov-Nordheim in the classical limit [52], and yield at saturation density, ρ 0 = 0.16 fm−3 (2.7 × 10 14 g/cm3), a binding energy E ( ρ 0 ) = − 16 MeV/nucleon, and a compressibility of about 250 MeV.
The corresponding mathematical expressions are
V n p ( r ) = V r r e − μ r r − V r r c e − μ r r c − V a r e − μ a r + V a r c e − μ a r c V n n ( r ) = V 0 r e − μ 0 r − V 0 r c e − μ 0 r c
where, after the cutoff radius r c the potentials are set to zero. Notice that in this phenomenological approach, to have nuclear matter bound but purely neutron matter unbound, V n p is attractive, and V n n is purely repulsive. Furthermore, in studies of infinite systems, the proton-proton Coulomb interaction is ignored and V n n = V p p ; the Coulomb interaction is only used in reactions.
The values of μ r , μ a , μ 0 and V r , V a , V 0 were first determined by the Urbana-Champaign group of V.J. Pandharipande for cold nuclear matter [51,52,53], and were later improved to reproduce cold NM binding energies more accurately [2]; the values of the parameters are listed in Table A1.
The potentials developed by Pandharipande [51,52,53] were aimed at reproducing the experimental cross sections in nucleon–nucleon collisions. By design, these potentials did not contain the Coulomb potential of the proton–proton interactions, which is assumed to be embedded in the potentials. Furthermore, there is no momentum dependence in the potentials, and thus, one cannot extract the momentum dependence of the compressibility, saturation density, symmetry energy, etc.
In the neutron star crust, electrons fill the voids between protons and neutrons. This means that the effect of the Coulomb potential of the electron gas must be included in any model of the crust. The principal effect of the electron cloud is to screen the Coulomb potential of the protons. This screening effect reduces the strength of the Coulomb forces between the protons, which is important for understanding the structure and dynamics of the crust. The long range of the Coulomb interaction makes it difficult to simulate infinite systems. Some approximation is therefore needed when modeling the neutron star crust. One common approximation is to use the Thomas-Fermi approximation, which assumes that the electron gas is in a state of uniform density.
The Coulomb effect of the gas of electrons is implemented through the Thomas-Fermi screening potential [21,24], which takes the gas of electrons as a Fermi gas with the same density of the protons. The electron gas becomes polarized screening the protons. The corresponding Poisson equation yields the potential
V T F = q 2 r e − r / λ ,
with a cutoff distance for V T F of 20 fm, and a screening length of λ = 10 fm, to have the properties of the resulting pasta essentially constant (see Ref. [54]).
The pasta configurations were obtained by cooling a system of 4000 nucleons of a given percentage of protons and neutrons and at a certain density, from T = 4.0 MeV down to a final temperature of 0.2 MeV. The proportion of protons to neutrons used was of 0.1 ≤ X = Z / A ≤ 0.5 , and the densities were 0.04 fm−3 ≤ ρ ≤ 0.085 fm−3 (0.68 × 1014 g/cm3 < ρ < 1.43 × 10 14 g/cm3). The size of the cell is determined by the number of particles and the density; however, using it with periodic boundary conditions, it is large enough as to avoid finite size effects of smaller lengths than the cell size [54].
The pastas were extracted from a Markov chain produced by CMD. To avoid correlations between the positions of the nucleons, each of the pastas was separated from the previous one by 10,000 time steps. The total number of structures obtained was 9600.
Table A1. Parameter set for the CMD computations.
Table A1. Parameter set for the CMD computations.
ParameterValueParameterValue
V r 3097.0 MeV μ r 1.648 fm−1
V a 2696.0 MeV μ a 1.528 fm−1
V 0 379.5 MeV μ 0 1.628 fm−1
r c 5.4/20 fm

Appendix B. Minkowski Functionals

Geometries of complex structures can be decomposed into simple structures that are connected side by side. The simple structures used are the voxels. The connectivity across voxels can be quantified by the number of vertices, ( n v ), edges ( n e ), faces ( n f ) and volume ( n c ) of the external shape of the structure formed with voxels.
The Minkowski functionals, are defined in terms of n v , n e , n f and n c as follows [28]:
V = n c , S = − 6 n c + 2 n f , 2 B = 3 n c − 2 n f + n e , χ = − n c + n f − n e + n v
V is the volume of the intersection of the voxels, and S is the external surface area of the intersections of the voxels; these are two important properties that determine its shape. The mean breadth (B) functional is the mean external curvature of the collection of voxels, with positive values of B corresponding to smooth, round surfaces, such as spheres and cylinders, and negative values of B corresponding to bodies with spiky or irregular surfaces, such as bodies with inner bubbles or hollow cylinders [32,37]. A value of B = 0 (with χ = 0 ) corresponds to bodies with flat surfaces like slabs.
The Euler characteristic χ is the Euler number of the external shape of the voxels, which can be written as
χ = isolated regions + cavities − tunnels
Solid objects, including those with isolated voids, can have positive Euler characteristics. However, structures like sponges or large inner voids, like jungle gym structures, have negative Euler characteristics because the number of voids or tunnels is larger than the number of isolated regions.
The various values of B ,   χ corresponding to the different pastas associated are listed in Table 1 (see Ref. [25] for more information). As it was mentioned previously, slab-like structures (lasagnas) correspond to B ∼ χ ∼ 0 . Spherical-like (Gnocchi) or cylinder-like (spaghetti) structures, correspond to B ≥ 0 and χ ≥ 0 . Likewise, sponge-like (Jungle-gym) structures attain χ < 0 . The structures with matter replaced by voids and vice versa use the prefix anti- in Table 1.
The calculation of the Minkowski functionals require that each nucleon is represented by a 3D structure. Here we bin space using “voxels”, each of which contains a single nucleon. The procedure to do this is explained in Ref. [5].

References

  1. Alcain, P.; Dorso, C. The neutrino opacity of neutron rich matter. Nucl. Phys. A 2017, 961, 183–199. [Google Scholar] [CrossRef]
  2. Alcain, P.N.; Dorso, C.O. Dynamics of fragment formation in neutron-rich matter. Phys. Rev. C 2018, 97, 015803. [Google Scholar] [CrossRef]
  3. Dorso, C.; Frank, G.; López, J. Phase transitions and symmetry energy in nuclear pasta. Nucl. Phys. A 2018, 978, 35–64. [Google Scholar] [CrossRef]
  4. Dorso, C.; Frank, G.; López, J. Symmetry energy in neutron star matter. Nucl. Phys. A 2019, 984, 77–98. [Google Scholar] [CrossRef]
  5. López, J.A.; Dorso, C.O.; Frank, G.A. Properties of nuclear pastas. Front. Phys. 2021, 16, 24301. [Google Scholar] [CrossRef]
  6. López, J.A.; Muñoz, J.A. Analytical expression and neural network study of the symmetry energy. CERN Proc. 2019, 1, 29–34. [Google Scholar]
  7. López, J.; Muñoz, J.; Dorso, C.; Frank, G. Machine learning Minkoswki functionals of neutron star crusts. J. Phys. Conf. Ser. 2020, 1643, 012054. [Google Scholar] [CrossRef]
  8. Utama, R.; Piekarewicz, J.; Prosper, H. Nuclear mass predictions for the crustal composition of neutron stars: A Bayesian neural network approach. Phys. Rev. C 2016, 93, 014311. [Google Scholar] [CrossRef]
  9. Sharma, B.K.; Centelles, M.; Viñas, X.; Baldo, M.; Burgio, G.F. Unified equation of state for neutron stars on a microscopic basis. Astron. Astrophys. 2015, 584, A103. [Google Scholar] [CrossRef]
  10. Piekarewicz, J.; Sánchez, G.T. Proton fraction in the inner neutron-star crust. Phys. Rev. C 2012, 85, 015807. [Google Scholar] [CrossRef]
  11. Dutra, M.; Lourenço, O.; Sá Martins, J.S.; Delfino, A.; Stone, J.R.; Stevenson, P.D. Skyrme interaction and nuclear matter constraints. Phys. Rev. C 2012, 85, 035201. [Google Scholar] [CrossRef]
  12. Dutra, M.; Lourenço, O.; Avancini, S.S.; Carlson, B.V.; Delfino, A.; Menezes, D.P.; Providência, C.; Typel, S.; Stone, J.R. Relativistic mean-field hadronic models under nuclear matter constraints. Phys. Rev. C 2014, 90, 055203. [Google Scholar] [CrossRef]
  13. Hagel, K.; Natowitz, J.B.; Röpke, G. The equation of state and symmetry energy of low-density nuclear matter. Eur. Phys. J. A 2014, 50, 1–16. [Google Scholar] [CrossRef]
  14. Horowitz, M.; Pérez García, M.; Piekarewicz, J. Neutrino-pasta scattering: The opacity of nonuniform neutron-rich matter. Phys. Rev. C 2004, 69, 045804. [Google Scholar] [CrossRef]
  15. Pi, M.; Barranco, M.; Perez-Canyellas, A.; Polls, A. Supernova matter: A semiclassical approach. Astron. Astrophys. Suppl. Ser. 1986, 64, 439–451. [Google Scholar]
  16. Caplan, M.E.; Schneider, A.S.; Horowitz, C.J. Elasticity of Nuclear Pasta. Phys. Rev. Lett. 2018, 121, 132701. [Google Scholar] [CrossRef]
  17. Ravenhall, D.G.; Pethick, C.J.; Wilson, J.R. Structure of Matter below Nuclear Saturation Density. Phys. Rev. Lett. 1983, 50, 2066–2069. [Google Scholar] [CrossRef]
  18. Hashimoto, M.A.; Seki, H.; Yamada, M. Shape of Nuclei in the Crust of Neutron Star. Prog. Theor. Phys. 1984, 71, 320–326. [Google Scholar] [CrossRef]
  19. Williams, R.; Koonin, S. Sub-saturation phases of nuclear matter. Nucl. Phys. A 1985, 435, 844–858. [Google Scholar] [CrossRef]
  20. Page, D.; Lattimer, J.M.; Prakash, M.; Steiner, A.W. Minimal Cooling of Neutron Stars: A New Paradigm. Astrophys. J. Suppl. Ser. 2004, 155, 623. [Google Scholar] [CrossRef]
  21. Maruyama, T.; Niita, K.; Oyamatsu, K.; Maruyama, T.; Chiba, S.; Iwamoto, A. Quantum molecular dynamics approach to the nuclear matter below the saturation density. Phys. Rev. C 1998, 57, 655–665. [Google Scholar] [CrossRef]
  22. Kido, T.; Maruyama, T.; Niita, K.; Chiba, S. MD simulation study for nuclear matter. Nucl. Phys. A 2000, 663, 877c–880c. [Google Scholar] [CrossRef]
  23. Watanabe, G.; Sato, K.; Yasuoka, K.; Ebisuzaki, T. Microscopic study of slablike and rodlike nuclei: Quantum molecular dynamics approach. Phys. Rev. C 2002, 66, 012801. [Google Scholar] [CrossRef]
  24. Horowitz, C.J.; Pérez-García, M.A.; Carriere, J.; Berry, D.K.; Piekarewicz, J. Nonuniform neutron-rich matter and coherent neutrino scattering. Phys. Rev. C 2004, 70, 065806. [Google Scholar] [CrossRef]
  25. Dorso, C.O.; Giménez Molinelli, P.A.; López, J.A. Topological characterization of neutron star crusts. Phys. Rev. C 2012, 86, 055805. [Google Scholar] [CrossRef]
  26. Alcain, P.N.; Giménez Molinelli, P.A.; Dorso, C.O. Beyond nuclear “pasta”: Phase transitions and neutrino opacity of new “pasta” phases. Phys. Rev. C 2014, 90, 065803. [Google Scholar] [CrossRef]
  27. Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comput. Phys. 1995, 117, 1–19. [Google Scholar] [CrossRef]
  28. Michielsen, K.; De Raedt, H. Integral-geometry morphological image analysis. Phys. Rep. 2001, 347, 461–538. [Google Scholar] [CrossRef]
  29. Watanabe, G.; Sato, K.; Yasuoka, K.; Ebisuzaki, T. Structure of cold nuclear matter at subnuclear densities by quantum molecular dynamics. Phys. Rev. C 2003, 68, 035806. [Google Scholar] [CrossRef]
  30. Schneider, A.S.; Horowitz, C.J.; Hughto, J.; Berry, D.K. Nuclear “pasta” formation. Phys. Rev. C 2013, 88, 065807. [Google Scholar] [CrossRef]
  31. Caplan, M.E.; Horowitz, C.J. Colloquium: Astromaterial science and nuclear pasta. Rev. Mod. Phys. 2017, 89, 041002. [Google Scholar] [CrossRef]
  32. Schuetrumpf, B.; Klatt, M.A.; Iida, K.; Maruhn, J.A.; Mecke, K.; Reinhard, P.G. Time-dependent Hartree-Fock approach to nuclear “pasta” at finite temperature. Phys. Rev. C 2013, 87, 055805. [Google Scholar] [CrossRef]
  33. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-learn: Machine Learning in Python. J. Mach. Learn. Res. 2011, 12, 2825–2830. [Google Scholar]
  34. Kingma, D.P.; Ba, J. Adam: A Method for Stochastic Optimization. In Proceedings of the 3rd International Conference on Learning Representations, ICLR 2015, San Diego, CA, USA, 7–9 May 2015. [Google Scholar]
  35. Hunter, J.D. Matplotlib: A 2D graphics environment. Comput. Sci. Eng. 2007, 9, 90–95. [Google Scholar] [CrossRef]
  36. Schuetrumpf, B.; Iida, K.; Maruhn, J.A.; Reinhard, P.G. Nuclear “pasta matter” for different proton fractions. Phys. Rev. C 2014, 90, 055802. [Google Scholar] [CrossRef]
  37. Sonoda, H.; Watanabe, G.; Sato, K.; Yasuoka, K.; Ebisuzaki, T. Phase diagram of nuclear “pasta” and its uncertainties in supernova cores. Phys. Rev. C 2008, 77, 035806. [Google Scholar] [CrossRef]
  38. Shchechilin, N.N.; Pearson, J.; Chamel, N. Nuclear Pasta in Cold Non-Accreting Neutron Stars: Symmetry Energy Effects. Phys. Sci. Forum 2023, 7, 10. [Google Scholar] [CrossRef]
  39. Sedrakian, A.; Clark, J. Superfluidity in nuclear systems and neutron stars. Eur. Phys. J. A 2019, 55, 167. [Google Scholar] [CrossRef]
  40. Muñoz, J.A. Jamunozlab/Dangelo. 2023. Available online: https://zenodo.org/records/10032533 (accessed on 22 October 2023).
  41. Colonna, M.; Zhang, Y.; Wang, Y.; Cozma, D.; Danielewicz, P.; Ko, C.; Zhang, F. Comparison of Heavy-Ion Transport Simulations: Mean-field Dynamics in a Box. Phys. Rev. C 2021, 104, 024603. [Google Scholar] [CrossRef]
  42. Bertsch, G.; Das Gupta, S. A guide to microscopic models for intermediate energy heavy ion collisions. Phys. Rep. 1988, 160, 189–233. [Google Scholar] [CrossRef]
  43. Danielewicz, P. Determination of the mean-field momentum-dependence using elliptic flow. Nucl. Phys. A 2000, 673, 375–410. [Google Scholar] [CrossRef]
  44. Li, B.A. Neutron-Proton Differential Flow as a Probe of Isospin-Dependence of the Nuclear Equation of State. Phys. Rev. Lett. 2000, 85, 4221–4224. [Google Scholar] [CrossRef]
  45. Aichelin, J.; Stöcker, H. Quantum molecular dynamics—A novel approach to N-body correlations in heavy ion collisions. Phys. Lett. B 1986, 176, 14–19. [Google Scholar] [CrossRef]
  46. Ono, A.; Danielewicz, P.; Friedman, W.A.; Lynch, W.G.; Tsang, M.B. Isospin fractionation and isoscaling in dynamical simulations of nuclear collisions. Phys. Rev. C 2003, 68, 051601. [Google Scholar] [CrossRef]
  47. Bondorf, J.; Botvina, A.; Iljinov, A.; Mishustin, I.; Sneppen, K. Statistical multifragmentation of nuclei. Phys. Rep. 1995, 257, 133–221. [Google Scholar] [CrossRef]
  48. Chernomoretz, A.; Gingras, L.; Larochelle, Y.; Beaulieu, L.; Roy, R.; St-Pierre, C.; Dorso, C.O. Quasiclassical model of intermediate velocity particle production in asymmetric heavy ion reactions. Phys. Rev. C 2002, 65, 054613. [Google Scholar] [CrossRef]
  49. Horowitz, C.J.; Pérez-García, M.A.; Berry, D.K.; Piekarewicz, J. Dynamical response of the nuclear “pasta” in neutron star crusts. Phys. Rev. C 2005, 72, 035801. [Google Scholar] [CrossRef]
  50. López, J.A.; Ramírez-Homs, E.; González, R.; Ravelo, R. Isospin-asymmetric nuclear matter. Phys. Rev. C 2014, 89, 024611. [Google Scholar] [CrossRef]
  51. Vicentini, A.; Jacucci, G.; Pandharipande, V.R. Fragmentation of hot classical drops. Phys. Rev. C 1985, 31, 1783–1793. [Google Scholar] [CrossRef]
  52. Lenk, R.J.; Schlagel, T.J.; Pandharipande, V.R. Accuracy of the Vlasov-Nordheim approximation in the classical limit. Phys. Rev. C 1990, 42, 372–385. [Google Scholar] [CrossRef]
  53. Lenk, R.J.; Pandharipande, V.R. Disassembly of hot classical charged drops. Phys. Rev. C 1986, 34, 177–184. [Google Scholar] [CrossRef] [PubMed]
  54. Alcain, P.N.; Giménez Molinelli, P.A.; Nichols, J.I.; Dorso, C.O. Effect of Coulomb screening length on nuclear “pasta” simulations. Phys. Rev. C 2014, 89, 055801. [Google Scholar] [CrossRef]
Dynamics 04 00009 g001
Figure 1. Nuclei arrangements (pastas) stable at a given density and temperature when the proportion of protons is 0.5. See text for details. (color online).
Figure 1. Nuclei arrangements (pastas) stable at a given density and temperature when the proportion of protons is 0.5. See text for details. (color online).
Dynamics 04 00009 g001
Dynamics 04 00009 g002aDynamics 04 00009 g002b
Figure 2. Constant-density temperature vs. proton proportion cuts of the three-dimensional phase diagram of nuclear pasta structures of neutron star matter predicted by the current methodology. See text for details. (color online).
Figure 2. Constant-density temperature vs. proton proportion cuts of the three-dimensional phase diagram of nuclear pasta structures of neutron star matter predicted by the current methodology. See text for details. (color online).
Dynamics 04 00009 g002aDynamics 04 00009 g002b
Dynamics 04 00009 g003
Figure 3. Pasta phases at T ≈ 0 MeV for given densities (in fm−3) and proton proportion predicted by the current methodology (left), and static Hartree-Fock from Ref. [36] (right). (color online).
Figure 3. Pasta phases at T ≈ 0 MeV for given densities (in fm−3) and proton proportion predicted by the current methodology (left), and static Hartree-Fock from Ref. [36] (right). (color online).
Dynamics 04 00009 g003
Table 1. Integral mean curvature (B) and Euler characteristic ( χ ) of the pasta shapes.
Table 1. Integral mean curvature (B) and Euler characteristic ( χ ) of the pasta shapes.
B < 0B∼ 0B > 0
χ > 0 Anti-Gnocchi Gnocchi
χ ∼ 0 Anti-SpaghettiLasagnaSpaghetti
χ < 0 Anti-Jungle Gym Jungle Gym
Table 2. Normalized root-mean-square error (NRMSE) and normalized mean absolute error (NMAE) for the neural network models used in this work.
Table 2. Normalized root-mean-square error (NRMSE) and normalized mean absolute error (NMAE) for the neural network models used in this work.
TestFull
BNRMSE8.8%5.5%
NMAE3.0%4.4%
χ NRMSE11.8%5.4%
NMAE3.9%3.8%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Muñoz, J.A.; López, J.A. Phase Diagram of Nuclear Pastas in Neutron Star Crusts. Dynamics 2024, 4, 157-169. https://doi.org/10.3390/dynamics4010009

AMA Style

Muñoz JA, López JA. Phase Diagram of Nuclear Pastas in Neutron Star Crusts. Dynamics. 2024; 4(1):157-169. https://doi.org/10.3390/dynamics4010009

Chicago/Turabian Style

Muñoz, Jorge A., and Jorge A. López. 2024. "Phase Diagram of Nuclear Pastas in Neutron Star Crusts" Dynamics 4, no. 1: 157-169. https://doi.org/10.3390/dynamics4010009

Article Metrics

Back to TopTop