No Arabic abstract
X-ray observations of transiently accreting neutron stars during quiescence provide information about the structure of neutron star crusts and the properties of dense matter. Interpretation of the observational data requires an understanding of the nuclear reactions that heat and cool the crust during accretion, and define its nonequilibrium composition. We identify here in detail the typical nuclear reaction sequences down to a depth in the inner crust where the mass density is 2E12 g/cm^3 using a full nuclear reaction network for a range of initial compositions. The reaction sequences differ substantially from previous work. We find a robust reduction of crust impurity at the transition to the inner crust regardless of initial composition, though shell effects can delay the formation of a pure crust somewhat to densities beyond 2E12 g/cm^3. This naturally explains the small inner crust impurity inferred from observations of a broad range of systems. The exception are initial compositions with A >= 102 nuclei, where the inner crust remains impure with an impurity parameter of Qimp~20 due to the N = 82 shell closure. In agreement with previous work we find that nuclear heating is relatively robust and independent of initial composition, while cooling via nuclear Urca cycles in the outer crust depends strongly on initial composition. This work forms a basis for future studies of the sensitivity of crust models to nuclear physics and provides profiles of composition for realistic crust models.
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.
With the first detections of binary neutron star mergers by gravitational-wave detectors, it proves timely to consider how the internal structure of neutron stars affects the way in which they can be asymmetrically deformed. Such deformations may leave measurable imprints on gravitational-wave signals and can be sourced through tidal interactions or the formation of mountains. We detail the formalism that describes fully-relativistic neutron star models with elastic crusts undergoing static perturbations. This formalism primes the problem for studies into a variety of mechanisms that can deform a neutron star. We present results for a barotropic equation of state and a realistic model for the elastic crust, which enables us to compute relevant quantities such as the tidal deformability parameter. We find that the inclusion of an elastic crust provides a very small correction to the tidal deformability. The results allow us to demonstrate when and where the crust starts to fail during a binary inspiral and we find that the majority of the crust will remain intact up until merger.
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.
In the solid crusts of neutron stars, the advection of the magnetic field by the current-carrying electrons, an effect known as Hall drift, should play a very important role as the ions remain essentially fixed (as long as the solid does not break). Although Hall drift preserves the magnetic field energy, it has been argued that it may drive a turbulent cascade to scales at which Ohmic dissipation becomes effective, allowing a much faster decay in objects with very strong fields. On the other hand, it has been found that there are Hall equilibria, i.e., field configurations that are unaffected by Hall drift. Here, we address the crucial question of the stability of these equilibria through axially symmetric (2D) numerical simulations of Hall drift and Ohmic diffusion, with the simplifying assumption of uniform electron density and conductivity. We demonstrate the 2D-stability of a purely poloidal equilibrium, for which Ohmic dissipation makes the field evolve towards an attractor state through adjacent stable configurations, around which damped oscillations occur. For this field, the decay scales with the Ohmic timescale. We also study the case of an unstable equilibrium consisting of both poloidal and toroidal field components that are confined within the crust. This field evolves into a stable configuration, which undergoes damped oscillations superimposed on a slow evolution towards an attractor, just as the purely poloidal one.
The fastest-spinning neutron stars in low-mass X-ray binaries, despite having undergone millions of years of accretion, have been observed to spin well below the Keplerian break-up frequency. We simulate the spin evolution of synthetic populations of accreting neutron stars in order to assess whether gravitational waves can explain this behaviour and provide the distribution of spins that is observed. We model both persistent and transient accretion and consider two gravitational-wave-production mechanisms that could be present in these systems: thermal mountains and unstable $r$-modes. We consider the case of no gravitational-wave emission and observe that this does not match well with observation. We find evidence for gravitational waves being able to provide the observed spin distribution; the most promising mechanisms being a permanent quadrupole, thermal mountains and unstable $r$-modes. However, based on the resultant distributions alone it is difficult to distinguish between the competing mechanisms.