No Arabic abstract
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.
The flow of a matter, accreting onto a magnetized neutron star, is accompanied by an electric current. The closing of the electric current occurs in the crust of a neutron stars in the polar region across the magnetic field. But the conductivity of the crust along the magnetic field greatly exceeds the conductivity across the field, so the current penetrates deep into the crust down up to the super conducting core. The magnetic field, generated by the accretion current, increases greatly with the depth of penetration due to the Hall conductivity of the crust is also much larger than the transverse conductivity. As a result, the current begins to flow mainly in the toroidal direction, creating a strong longitudinal magnetic field, far exceeding an initial dipole field. This field exists only in the narrow polar tube of $r$ width, narrowing with the depth, i.e. with increasing of the crust density $rho$, $rpropto rho^{-1/4}$. Accordingly, the magnetic field $B$ in the tube increases with the depth, $Bpropto rho^{1/2}$, and reaches the value of about $10^{17}$ Gauss in the core. It destroys super conducting vortices in the core of a star in the narrow region of the size of the order of ten centimeters. Because of generated density gradient of vortices they constantly flow into this dead zone and the number of vortices decreases, the magnetic field of a star decreases as well. The attenuation of the magnetic field is exponential, $B=B_0(1+t/tau)^{-1}$. The characteristic time of decreasing of the magnetic field $tau$ is equal to $tausimeq 10^3$ years. Thus, the magnetic field of accreted neutron stars decreases to values of $10^8 - 10^9$ Gauss during $10^7-10^6$ years.
Transiently accreting neutron stars in quiescence (Lx<10^34 erg/s) have been observed to vary in intensity by factors of few, over timescales of days to years. If the quiescent luminosity is powered by a hot NS core, the core cooling timescale is much longer than the recurrence time, and cannot explain the observed, more rapid variability. However, the non-equilibrium reactions which occur in the crust during outbursts deposit energy in iso-density shells, from which the thermal diffusion timescale to the photosphere is days to years. The predicted magnitude of variability is too low to explain the observed variability unless - as is widely believed - the neutrons beyond the neutron-drip density are superfluid. Even then, variability due to this mechanism in models with standard core neutrino cooling processes is less than 50 per cent - still too low to explain the reported variability. However, models with rapid core neutrino cooling can produce variability by a factor as great as 20, on timescales of days to years following an outburst. Thus, the factors of few intensity variability observed from transiently accreting neutron stars can be accounted for by this mechanism only if rapid core cooling processes are active.
The strong magnetic field of neutron stars is intimately coupled to the observed temperature and spectral properties, as well as to the observed timing properties (distribution of spin periods and period derivatives). Thus, a proper theoretical and numerical study of the magnetic field evolution equations, supplemented with detailed calculations of microphysical properties (heat and electrical conductivity, neutrino emission rates) is crucial to understand how the strength and topology of the magnetic field vary as a function of age, which in turn is the key to decipher the physical processes behind the varied neutron star phenomenology. In this review, we go through the basic theory describing the magneto-thermal evolution models of neutron stars, focusing on numerical techniques, and providing a battery of benchmark tests to be used as a reference for present and future code developments. We summarize well-known results from axisymmetric cases, give a new look at the latest 3D advances, and present an overview of the expectations for the field in the coming years.
In quiescent low-mass X-ray binaries (qLMXBs) containing neutron stars, the origin of the thermal X-ray component may be either release of heat from the core of the neutron star, or continuing low-level accretion. In general, heat from the core should be stable on timescales $<10^4$ years, while continuing accretion may produce variations on a range of timescales. While some quiescent neutron stars (e.g. Cen X-4, Aql X-1) have shown variations in their thermal components on a range of timescales, several others, particularly those in globular clusters with no detectable nonthermal hard X-rays (fit with a powerlaw), have shown no measurable variations. Here, we constrain the spectral variations of 12 low mass X-ray binaries in 3 globular clusters over $sim10$ years. We find no evidence of variations in 10 cases, with limits on temperature variations below 11% for the 7 qLMXBs without powerlaw components, and limits on variations below 20% for 3 other qLMXBs that do show non-thermal emission. However, in 2 qLMXBs showing powerlaw components in their spectra (NGC 6440 CX 1 & Terzan 5 CX 12) we find marginal evidence for a 10% decline in temperature, suggesting the presence of continuing low-level accretion. This work adds to the evidence that the thermal X-ray component in quiescent neutron stars without powerlaw components can be explained by heat deposited in the core during outbursts. Finally, we also investigate the correlation between hydrogen column density (N$_H$) and optical extinction (A$_V$) using our sample and current models of interstellar X-ray absorption, finding $N_H ({rm cm}^{-2}) = (2.81pm0.13)times10^{21} A_V$.
The crust of accreting neutron stars plays a central role in many different observational phenomena. In these stars, heavy elements produced by H-He burning in the rapid proton capture (rp-) process continually freeze to form new crust. In this paper, we explore the expected composition of the solid phase. We first demonstrate using molecular dynamics that two distinct types of chemical separation occur, depending on the composition of the rp-process ashes. We then calculate phase diagrams for three-component mixtures and use them to determine the allowed crust compositions. We show that, for the large range of atomic numbers produced in the rp-process ($Zsim 10$--$50$), the solid that forms has only a small number of available compositions. We conclude that accreting neutron star crusts should be polycrystalline, with domains of distinct composition. Our results motivate further work on the size of the compositional domains, and have implications for crust physics and accreting neutron star phenomenology.