From supernovae to neutron stars: crust formation time

This paper presents a simple analytic estimate for the onset time of neutron star crust formation during the late post-convective cooling phase, deriving closed-form expressions that predict the first solid phase typically appears between 100 and 500 seconds after birth, with specific dependence on the protoneutron star's mass, radius, and composition.

The Gist

Imagine a newborn neutron star not as a cold, dead rock, but as a scorching, chaotic ball of soup. This is a Protoneutron Star (PNS). When it is first born in a supernova explosion, it is incredibly hot and filled with a "soup" of particles, including heavy atomic nuclei floating freely like islands in a hot ocean. At this stage, the star is entirely fluid; it cannot hold its shape against stress because the heat is too strong.

This paper asks a simple question: How long does it take for this hot, fluid soup to cool down enough to turn into a solid crust?

Think of it like a pot of hot soup cooling on a stove. Eventually, the surface gets cold enough that the ingredients stop swirling and start locking together into a solid layer. For a neutron star, this "solid layer" is called a crust, and its formation is a major milestone in the star's life.

Here is how the authors figured out the timing, using simple analogies:

1. The Cooling Process (The Leaky Bucket)

The star cools down by shooting out invisible particles called neutrinos. Imagine the star as a hot bucket with a leaky bottom. The faster the water (heat) leaks out, the faster the bucket cools.

• The authors used a mathematical "leak rate" based on how heavy the star is and how big it is.
• They calculated that as time passes, the "soup" inside gets less chaotic (lower entropy) and the temperature drops.

2. The "Freezing Point" (The Crystal Lattice)

In a normal freezer, water turns to ice when it hits 0°C. In a neutron star, the "freezing point" is different. It depends on how strongly the heavy atomic nuclei (the "islands" in the soup) attract each other.

• If the nuclei have a high electric charge (like a strong magnet), they grab onto each other sooner, even if it's still quite hot.
• If they have a low charge, they need to get much colder before they lock together.
• The authors calculated a specific "crystallization temperature" for the outer layers of the star.

3. The Race: Cooling vs. Freezing

The paper tracks a race between two things happening at the star's "surface" (called the neutrinosphere):

1. The Cooling Curve: The temperature of the star dropping over time.
2. The Freezing Line: The specific temperature required for the nuclei to turn solid at that specific density.

The Crust Formation Time is the exact moment when the star's cooling curve dips below the freezing line. That is the moment the first solid patch appears.

The Results: How Long Does It Take?

Using their "recipe" (which includes the star's mass, size, and the type of atoms inside), the authors found that for a typical newborn neutron star:

• The first solid crust usually appears between 100 and 500 seconds after the star is born.
• Heavier stars or smaller stars tend to take longer to form a crust because their "leak" (cooling) is slower.
• Stars with heavier, more charged atoms inside form a crust faster because those atoms stick together more easily.

Why This Matters (According to the Paper)

The authors explain that once this solid crust forms, the star changes character. It goes from being a fluid that can't hold stress to a solid shell that can store "elastic energy" (like a stretched rubber band). This solid shell might also change how the star's magnetic field behaves later on.

Important Note on Limitations:
The authors are careful to say this is a rough estimate, like a weather forecast. They used simplified math (ignoring complex turbulence inside the star) to get a clear, easy-to-use formula. They admit that in reality, the star's interior becomes semi-transparent to neutrinos after about 100 seconds, which makes the math more complicated. However, their formula provides a solid "benchmark" for scientists to understand when this solid shell likely starts to form.

In short: This paper provides a simple stopwatch for the universe, estimating that a newborn neutron star takes about 2 to 8 minutes to grow its first solid skin.

Technical Summary

Technical Summary: Crust Formation Time in Protoneutron Stars

Problem Statement
Neutron stars (NSs) are born as hot, lepton-rich protoneutron stars (PNSs) that cool via neutrino emission. A critical phase in their evolution is the transition from a fluid state to a solid state, marked by the crystallization of heavy ions in the outer layers into a Coulomb lattice. While fully self-consistent numerical simulations with neutrino radiation transport exist, they are computationally expensive and not always suitable for scanning broad parameter dependencies. Consequently, there is a need for an analytic or semi-analytic estimate that explicitly clarifies how macroscopic parameters (such as PNS mass, radius, and composition) control the onset time of crust formation ($t_{\text{crust}}$).

Methodology
The authors develop a simple analytic estimate for the crust-formation timescale during the late, post-convective cooling phase of a PNS. The approach combines two primary components:

1. Neutrino Cooling and Entropy Evolution: Utilizing a diffusion-based analytic solution for neutrino luminosity (Suwa et al. 2021), the authors model the time-dependent thermal trajectory of the PNS. They assume an approximately isentropic interior structure, a condition established by convective mixing in the early cooling phase. This allows for the derivation of time-dependent density ($\rho_\nu$) and temperature ($T_\nu$) at the neutrinosphere.
2. Crystallization Criterion: The onset of the solid phase is determined by the Coulomb crystallization condition for heavy nuclei. This is expressed through the Coulomb coupling parameter ($\Gamma$), where crystallization occurs when the Coulomb interaction energy exceeds thermal agitation. The authors equate the neutrinosphere temperature $T_\nu(t)$ with the critical crystallization temperature $T_c(\rho_\nu)$ evaluated at the neutrinosphere density.

The model treats the PNS as having an isentropic interior and an isothermal outer atmosphere. The authors derive closed-form expressions for the entropy at crystallization ($s_{\text{crust}}$) and the corresponding time ($t_{\text{crust}}$) by matching the cooling trajectory to the crystallization threshold.

Key Results
The study yields explicit scaling relations for the crust-formation time, $t_{\text{crust}}$, which depends on:

• PNS Parameters: Mass ($M_{\text{PNS}}$) and Radius ($R_{\text{PNS}}$).
• Cooling Normalization: An effective diffusion parameter ($g\beta$) and total emitted neutrino energy ($E_{\text{tot}}$).
• Microphysics: Ionic charge ($Z$), heavy-nuclei mass fraction ($x_a$), and the Coulomb coupling parameter at melting ($\Gamma$).

For canonical microphysics parameters ($\Gamma = 175$, $x_a = 0.3$, $Z = 40$) and representative PNS configurations, the authors find that the first solid phase typically appears at $t_{\text{crust}} \sim 100\text{--}500$ seconds.

The derived scaling relations indicate that:

• Larger PNS masses and smaller radii generally delay crystallization due to their impact on the diffusion timescale.
• Higher ionic charges ($Z$) and larger heavy-nuclei mass fractions ($x_a$) accelerate crystallization by raising the melting temperature.

Significance and Claims
The paper claims to provide a "useful late-time analytic benchmark" for the onset of crust formation. The primary value of this work lies in its ability to:

1. Clarify Dependencies: Explicitly show how macroscopic PNS parameters and composition variables control the timing of the fluid-to-solid transition.
2. Enable Rapid Exploration: Allow for the quick scanning of plausible PNS configurations without the computational cost of full radiation-hydrodynamics simulations.
3. Define a Transition Time: Offer a characteristic timescale for the emergence of a solid crust, which may serve as a useful reference for understanding the subsequent evolution of magnetic fields and rotation in newborn neutron stars, as the outer layers transition from fluid-like to solid and highly conducting.

Limitations and Caveats
The authors explicitly note several limitations to maintain modesty regarding their claims:

• Extrapolation: The cooling formulae were calibrated against detailed models up to $\sim$several $\times 10$ seconds; their application to $t \sim 100\text{--}500$ seconds is an extrapolation.
• Diffusion Approximation: By $t \sim 100$ s, the PNS interior becomes increasingly semitransparent to neutrinos, meaning a purely diffusive description is no longer strictly accurate everywhere.
• Microphysical Simplifications: The composition parameters ($Z$, $x_a$, $Z/A$) are treated as effective constants, whereas they likely vary with density and temperature in realistic crusts. The melting prescription relies on a fiducial one-component plasma (OCP) estimate rather than complex finite-temperature crust microphysics.

Therefore, the results are intended to clarify parameter dependencies rather than provide model-specific predictions for individual PNSs. A quantitatively predictive treatment would require multidimensional neutrino-radiation hydrodynamics coupled with improved microphysics.