Thermal and transport properties of the neutron star inner crust

We review the nuclear and condensed matter physics underlying the thermal and transport properties of the neutron star inner crust. These properties play a key role in interpreting transient phenomena such as thermal relaxation in accreting neutron stars, superbursts, and magnetar flares. We emphasize simplifications that occur at low temperature where the inner crust can be described in terms of electrons and collective excitations. The heat conductivity and heat capacity of the solid and superfluid phase of matter is discussed in detail and we emphasize its role in interpreting observations of neutron stars in soft X-ray transients. We highlight recent theoretical and observational results, and identify future work needed to better understand a host of transient phenomena in neutron stars.

💡 Research Summary

The paper provides a comprehensive review of the microscopic physics governing the thermal and transport properties of the neutron‑star inner crust and explains why these properties are crucial for interpreting a variety of transient astrophysical phenomena. The inner crust, located between densities of roughly 4 × 10¹¹ g cm⁻³ and the nuclear saturation density, consists of a Coulomb lattice of neutron‑rich nuclei, a sea of dripped neutrons, and a highly degenerate electron gas. At temperatures below ∼10⁹ K the electrons behave as a nearly ideal Fermi liquid, the nuclei support quantized lattice vibrations (phonons), and the dripped neutrons become superfluid through BCS‑type pairing. This separation of degrees of freedom enables a series of simplifications: electron‑phonon scattering, electron‑impurity scattering, and electron‑neutron scattering can be treated independently, and the superfluid gap suppresses neutron‑related scattering at low temperature.

The authors first derive the electron thermal conductivity κₑ using the Boltzmann transport equation. κₑ is expressed in terms of the electron relaxation time τₑ, which is the inverse sum of the scattering rates from phonons, impurities, and superfluid neutrons. In the low‑temperature regime Umklapp processes are frozen out, so κₑ scales roughly as T⁻¹ when impurity scattering dominates, but reverts to a steeper T⁻³ dependence when electron‑phonon scattering becomes significant. The presence of lattice defects, impurity clusters, or exotic “pasta” phases (spaghetti, lasagna) introduces additional scattering channels that can dramatically lower κₑ.

Heat capacity is dissected into three additive contributions. The electronic part Cₑ is linear in temperature (Cₑ ∝ T) because the electrons are degenerate. The phonon part follows the Debye law (C_ph ∝ T³) up to the Debye temperature, after which it saturates. The superfluid neutron contribution C_sf is exponentially suppressed (C_sf ∝ exp(−Δ/k_BT)) below the critical temperature T_c ≈ 10⁸ K, where Δ is the pairing gap. Consequently, at temperatures well below T_c the total heat capacity is dominated by electrons and phonons, while near T_c the neutron term can become comparable.

The review then connects these microphysical inputs to observable phenomena. In soft X‑ray transients (SXTs) the post‑outburst cooling curve is governed by the product κ · C of the crust. A high κ (clean lattice, few impurities) yields rapid thermal relaxation on timescales of 10⁴–10⁵ s, whereas a low κ (disordered or pasta‑rich crust) can extend relaxation to 10⁶ s or longer. Superbursts and magnetar flares inject large amounts of heat into the crust; the subsequent temperature decay is sensitive to the size of the neutron pairing gap because a larger gap reduces the neutron heat capacity, leading to faster cooling. The authors demonstrate that current X‑ray observations can be used to constrain key microscopic parameters such as the impurity parameter Q_imp, the Debye temperature, and the superfluid gap Δ.

Finally, the authors outline several avenues for future work. First, realistic band‑structure calculations for electrons in pasta phases are needed to quantify the electron‑phonon coupling in these highly anisotropic environments. Second, the role of strong magnetic fields on spin‑orbit coupling and possible spin‑flip scattering channels must be explored, especially for magnetar flare modeling. Third, multi‑scale simulations that couple microscopic transport coefficients to global neutron‑star thermal evolution codes are required to reduce systematic uncertainties. By addressing these challenges, the community will be able to extract more precise information about the inner‑crust composition and state of matter from the growing body of transient neutron‑star observations.