No Arabic abstract
Recently, crust cooling times have been measured for neutron stars after extended outbursts. These observations are very sensitive to the thermal conductivity $kappa$ of the crust and strongly suggest that $kappa$ is large. We perform molecular dynamics simulations of the structure of the crust of an accreting neutron star using a complex composition that includes many impurities. The composition comes from simulations of rapid proton capture nucleosynthesys followed by electron captures. We find that the thermal conductivity is reduced by impurity scattering. In addition, we find phase separation. Some impurities with low atomic number $Z$ are concentrated in a subregion of the simulation volume. For our composition, the solid crust must separate into regions of different compositions. This could lead to an asymmetric star with a quadrupole deformation that radiates gravitational waves. Observations of crust cooling can constrain impurity concentrations.
We calculate the thermal conductivity of electrons for the strongly correlated multi-component ion plasma expected in the outer layers of neutron stars crust employing a Path Integral Monte Carlo (PIMC) approach. This allows us to isolate the low energy response of the ions and use it to calculate the electron scattering rate and the electron thermal conductivity. We find that the scattering rate is enhanced by a factor 2-4 compared to earlier calculations based on the simpler electron-impurity scattering formalism. This findings directly impacts the interpretation of thermal relaxation observed in transiently accreting neutron stars and has implications for the composition and nuclear reactions in the crust that occur during accretion.
Fusion reactions in the crust of an accreting neutron star are an important source of heat, and the depth at which these reactions occur is important for determining the temperature profile of the star. Fusion reactions depend strongly on the nuclear charge $Z$. Nuclei with $Zle 6$ can fuse at low densities in a liquid ocean. However, nuclei with Z=8 or 10 may not burn until higher densities where the crust is solid and electron capture has made the nuclei neutron rich. We calculate the $S$ factor for fusion reactions of neutron rich nuclei including $^{24}$O + $^{24}$O and $^{28}$Ne + $^{28}$Ne. We use a simple barrier penetration model. The $S$ factor could be further enhanced by dynamical effects involving the neutron rich skin. This possible enhancement in $S$ should be studied in the laboratory with neutron rich radioactive beams. We model the structure of the crust with molecular dynamics simulations. We find that the crust of accreting neutron stars may contain micro-crystals or regions of phase separation. Nevertheless, the screening factors that we determine for the enhancement of the rate of thermonuclear reactions are insensitive to these features. Finally, we calculate the rate of thermonuclear $^{24}$O + $^{24}$O fusion and find that $^{24}$O should burn at densities near $10^{11}$ g/cm$^3$. The energy released from this and similar reactions may be important for the temperature profile of the star.
We present mass excesses (ME) of neutron-rich isotopes of Ar through Fe, obtained via TOF-$Brho$ mass spectrometry at the National Superconducting Cyclotron Laboratory. Our new results have significantly reduced systematic uncertainties relative to a prior analysis, enabling the first determination of ME for $^{58,59}{rm Ti}$, $^{62}{rm V}$, $^{65}{rm Cr}$, $^{67,68}{rm Mn}$, and $^{69,70}{rm Fe}$. Our results show the $N=34$ subshell weaken at Sc and vanish at Ti, along with the absence of an $N=40$ subshell at Mn. This leads to a cooler accreted neutron star crust, highlighting the connection between the structure of nuclei and neutron stars.
We study long-term thermal evolution of neutron stars in soft X-ray transients (SXTs), taking the deep crustal heating into account consistently with the changes of the composition of the crust. We collect observational estimates of average accretion rates and thermal luminosities of such neutron stars and compare the theory with observations. We perform simulations of thermal evolution of accreting neutron stars, considering the gradual replacement of the original nonaccreted crust by the reprocessed accreted matter, the neutrino and photon energy losses, and the deep crustal heating due to nuclear reactions in the accreted crust. We test and compare results for different modern theoretical models. We update a compilation of the observational estimates of the thermal luminosities in quiescence and average accretion rates in the SXTs and compare the observational estimates with the theoretical results. Long-term thermal evolution of transiently accreting neutron stars is nonmonotonic. The quasi-equilibrium temperature in quiescence reaches a minimum and then increases toward the final steady state. The quasi-equilibrium thermal luminosity of a neutron star in an SXT can be substantially lower at the minimum than in the final state. This enlarges the range of possibilities for theoretical interpretation of observations of such neutron stars. The updates of the theory and observations leave unchanged the previous conclusions that the direct Urca process operates in relatively cold neutron stars and that an accreted heat-blanketing envelope is likely present in relatively hot neutron stars in the SXTs in quiescence. The results of the comparison of theory with observations favor suppression of the triplet pairing type of nucleon superfluidity in the neutron-star matter.
In this book chapter we review plasma crystals in the laboratory, in the interior of white dwarf stars, and in the crust of neutron stars. We describe a molecular dynamics formalism and show results for many neutron star crust properties including phase separation upon freezing, diffusion, breaking strain, shear viscosity and dynamics response of nuclear pasta. We end with a summary and discuss open questions and challenges for the future.