Nuclear pasta in hot neutron-star matter and proto-neutron stars

This study investigates nuclear pasta phases in hot neutron-star matter and proto-neutron stars using a compressible liquid-drop model, finding that models with a small symmetry energy slope predict diverse pasta shapes that significantly influence the thermal evolution and structure of the inner crust.

The Gist

Imagine a neutron star not as a smooth, solid ball of rock, but as a cosmic kitchen where the ingredients are being cooked at temperatures hotter than the center of the sun. This paper explores what happens to the "food" inside these stars, specifically focusing on a strange, gooey state of matter called "Nuclear Pasta."

Here is the breakdown of the research in simple terms:

1. The Setting: A Cosmic Pressure Cooker

When a massive star dies and collapses, it becomes a Proto-Neutron Star (PNS). Think of this as a newborn neutron star that is still very hot and cooling down. Inside, the pressure is so immense that atoms are crushed together.

Usually, we think of matter as either a solid (like a rock) or a gas (like steam). But in these stars, the matter is in a weird middle ground. It's a dense liquid mixed with a thin gas, separated by a sharp boundary. Because of the intense competition between the forces holding the nuclei together (surface tension) and the forces pushing them apart (electric repulsion), the matter doesn't just stay as a blob. It shapes itself into weird structures, much like how pasta shapes form when you squeeze dough.

2. The "Pasta" Shapes

As you go deeper into the star (increasing density), the matter changes shape, just like a game of Tetris:

• Droplets: Small, round balls of nuclei (like gnocchi).
• Rods: Long, thin sticks (like spaghetti).
• Slabs: Flat sheets (like lasagna).
• Tubes: Hollow cylinders (like penne).
• Bubbles: Holes in the liquid (like swiss cheese or ravioli).

Scientists call this "Nuclear Pasta." The paper asks: Which shapes actually appear, and does the temperature or the specific type of nuclear physics change the menu?

3. The Two Chefs: TM1 vs. TM1e

The researchers used two different "cookbooks" (mathematical models) to predict what the pasta looks like. Both cookbooks agree on how to handle the basic ingredients (protons and neutrons), but they disagree on one specific spice: Symmetry Energy.

• The "Spicy" Chef (TM1 model): This model assumes a high "symmetry energy slope" (a high value of $L = 110.8$). Think of this as a chef who loves spicy food. This model predicts that the matter stays very rigid. No matter how you cook it, you mostly get droplets (gnocchi). It's stubborn and doesn't want to change shape into rods or sheets.
• The "Mild" Chef (TM1e model): This model uses a lower symmetry energy slope ($L = 40$). This is a more flexible chef. This model predicts a full menu! At low temperatures, the matter happily transforms from droplets into rods, sheets, tubes, and bubbles as you go deeper into the star.

The Big Discovery: The study found that the "Mild" chef (TM1e) is likely the one who actually runs the kitchen in real neutron stars. The "Spicy" chef's prediction (only droplets) seems too limited compared to what we know about the universe.

4. Why Does This Matter? (The Thermal Evolution)

Why should we care if the matter is shaped like spaghetti or lasagna?

Imagine the star is a house. The inner crust (where the pasta lives) is the insulation in the walls.

• If the insulation is made of smooth, uniform blocks (no pasta), heat flows through easily.
• If the insulation is made of complex, wavy pasta shapes, it acts like a maze. It makes it harder for heat and neutrinos (tiny ghost particles) to escape.

The researchers found that this "pasta layer" in the inner crust is about 1.2 kilometers thick. This layer plays a huge role in how fast the star cools down. If the star has a lot of pasta, it might stay hot longer or cool down differently than if it were just smooth matter.

5. The Takeaway

This paper is like a culinary investigation into the most extreme kitchen in the universe.

• The Problem: We didn't know exactly what shapes the matter takes inside hot, newborn neutron stars.
• The Method: They used advanced math to simulate the star, testing two different theories about how nuclear forces work.
• The Result: One theory (TM1e) suggests a rich variety of pasta shapes exists, while the other (TM1) suggests a boring, droplet-only world. The evidence points to the "rich variety" theory being correct.
• The Impact: These shapes act as a thermal blanket, influencing how the star cools and evolves over time.

In short, the universe might be filled with cosmic spaghetti and lasagna, and understanding this "pasta" helps us understand how neutron stars age and cool down.

Technical Summary

1. Problem Statement

The study addresses the complex structure of matter in the inner crust of Proto-Neutron Stars (PNS) and hot neutron-star matter. As a PNS cools following a core-collapse supernova, it transitions from uniform matter to a non-uniform state where heavy nuclei form. Near the crust-core transition density, the competition between surface tension and Coulomb energy leads to the formation of exotic, non-spherical nuclear shapes known as "nuclear pasta" (e.g., rods, slabs, tubes, bubbles).

Key challenges addressed include:

• Determining the specific pasta phases that exist at finite temperatures ($T > 0$) relevant to PNS cooling.
• Understanding how the nuclear symmetry energy and its density dependence (specifically the slope parameter $L$) influence the emergence and stability of these pasta phases.
• Quantifying the impact of these pasta phases on the macroscopic properties of PNSs, such as their mass-radius relations and thermal evolution.

2. Methodology

The authors employ a multi-component theoretical framework combining nuclear interactions with thermodynamic modeling:

• Nuclear Interaction Model: They utilize the Relativistic Mean-Field (RMF) approach. To isolate the effects of symmetry energy, they compare two specific parameterizations:

  • TM1e: An extended model with a small symmetry energy slope ($L = 40$ MeV), consistent with recent astrophysical constraints.
  • TM1: The original model with a large symmetry energy slope ($L = 110.8$ MeV).
  • Both models share identical isoscalar properties, ensuring differences arise solely from the isovector (symmetry energy) sector.

• Thermodynamic Framework: The Compressible Liquid-Drop (CLD) model is used to describe pasta phases.

  • Wigner-Seitz Approximation: Space is divided into charge-neutral cells with specific geometric symmetries (droplet, rod, slab, tube, bubble).
  • Phase Coexistence: Each cell contains a dense liquid phase and a dilute gas phase (containing free nucleons and $\alpha$ particles) separated by a sharp interface.
  • Free Energy Minimization: The equilibrium state is determined by minimizing the total free energy density, which includes bulk, surface, and Coulomb contributions. This leads to equilibrium conditions that differ from standard Gibbs conditions due to finite-size effects.

• Surface Tension Treatment:

  • Surface tension ($\tau$) is calculated self-consistently using the Thomas-Fermi approximation for a 1D nuclear system.
  • A parameterized form of $\tau$ is employed, dependent on temperature ($T$) and the proton fraction of the liquid phase ($Y_p^L$), fitted to Thomas-Fermi results.

• PNS Modeling:

  • The study adopts a constant entropy per baryon ($S$) approach to simulate the neutrino-free cooling stage of a PNS.
  • Isentropic Equations of State (EOS) are constructed for $S = 0.5, 1, 2$ MeV.
  • The Tolman-Oppenheimer-Volkoff (TOV) equations are solved to derive mass-radius relations and internal profiles.

3. Key Contributions

• Systematic Comparison of Symmetry Energy Slopes: The paper provides a rigorous comparison of pasta phases between models with low ($L=40$) and high ($L=110.8$) symmetry energy slopes, isolating the role of $L$ in phase transitions.
• Finite-Temperature Pasta Phase Diagrams: The authors map out the phase diagrams in the baryon density ($n_b$) vs. temperature ($T$) plane, identifying the specific density ranges where different pasta shapes (droplets, rods, slabs, tubes, bubbles) are thermodynamically stable.
• Surface Tension Parameterization: They provide a computationally efficient parameterized surface tension derived from Thomas-Fermi calculations, facilitating the study of pasta phases at finite temperatures without the prohibitive cost of full self-consistent TF calculations at every step.
• Macroscopic Impact on PNS: The study quantifies how the inclusion of pasta phases in the inner crust alters the global structure of a cooling proto-neutron star.

4. Key Results

• Dependence on Symmetry Energy Slope ($L$):

  • TM1e ($L=40$ MeV): Predicts a rich variety of pasta shapes (droplets $\to$ rods $\to$ slabs $\to$ tubes $\to$ bubbles) at low temperatures. The non-uniform matter region is extensive.
  • TM1 ($L=110.8$ MeV): Predicts only the droplet configuration up to the crust-core transition density. The region of non-uniform matter is significantly smaller, and the critical temperature for the disappearance of pasta phases is lower.
  • Mechanism: A smaller $L$ favors more asymmetric matter, leading to larger proton fractions in the liquid phase and a broader density range for non-uniform structures.

• Thermodynamic Properties:

  • Proton Fraction ($Y_p$): In the TM1e model, the average proton fraction in non-uniform matter is significantly higher than in uniform matter and higher than in the TM1 model.
  • Free Energy: The formation of pasta structures lowers the free energy compared to uniform matter at subsaturation densities. However, bulk equilibrium calculations (neglecting surface/Coulomb terms) underestimate the free energy compared to the CLD results.
  • Temperature Effects: As temperature increases, the density range of non-uniform matter shrinks. At the critical temperature $T_c$, the system transitions to uniform matter.

• Proto-Neutron Star Structure:

  • Inner Crust Thickness: Using the TM1e model with $S=1$, nuclear pasta phases are found to exist in the inner crust with a thickness of approximately 1.2 km (extending from $\sim 12.1$ km to $\sim 13.3$ km in a $1.4 M_\odot$ star).
  • Mass-Radius Relations: The presence of pasta phases results in larger radii for PNSs compared to models using only uniform matter. This effect is more pronounced for lower-mass stars.
  • Cooling Evolution: As the entropy $S$ decreases (cooling), the PNS contracts, and the central density increases. The pasta phases play a crucial role in the thermal evolution by potentially altering thermal conductivity and neutrino transport.

5. Significance

• Astrophysical Constraints: The findings highlight that the symmetry energy slope $L$ is a critical parameter for determining the internal structure of neutron stars. Models with large $L$ (like TM1) may fail to predict complex pasta structures, potentially leading to inaccurate EOS tables for supernova simulations and binary merger modeling.
• Thermal Evolution: The existence of a thick pasta layer ($\sim 1.2$ km) in the inner crust suggests that pasta phases could significantly influence the thermal conductivity and neutrino opacity of the star, thereby affecting the cooling timescale of proto-neutron stars.
• Methodological Advancement: The successful application of the CLD method with parameterized surface tension offers a robust and efficient tool for generating EOS tables that include pasta phases for use in large-scale hydrodynamic simulations of core-collapse supernovae and neutron star mergers.

In conclusion, the paper demonstrates that nuclear pasta phases are a robust feature of hot neutron-star matter when the symmetry energy slope is small ($L \approx 40$ MeV), significantly altering the inner crust structure and the macroscopic properties of proto-neutron stars.