No Arabic abstract
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.
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.
Above-barrier fusion cross-sections for an isotopic chain of oxygen isotopes with A=16-19 incident on a $^{12}$C target are presented. Experimental data are compared with both static and dynamical microscopic calculations. These calculations are unable to explain the $sim$37% increase in the average above-barrier fusion cross-section observed for $^{19}$O as compared to $beta$-stable oxygen isotopes. This result suggests that for neutron-rich nuclei existing time-dependent Hartree-Fock calculations underpredict the role of dynamics at near-barrier energies. High-quality measurement of above-barrier fusion for an isotopic chain of increasingly neutron-rich nuclei provides an effective means to probe this fusion dynamics.
We compute the binding energy of neutron-rich oxygen isotopes and employ the coupled-cluster method and chiral nucleon-nucleon interactions at next-to-next-to-next-to-leading order with two different cutoffs. We obtain rather well-converged results in model spaces consisting of up to 21 oscillator shells. For interactions with a momentum cutoff of 500 MeV, we find that 28O is stable with respect to 24O, while calculations with a momentum cutoff of 600 MeV result in a slightly unbound 28O. The theoretical error estimates due to the omission of the three-nucleon forces and the truncation of excitations beyond three-particle-three-hole clusters indicate that the stability of 28O cannot be ruled out from ab-initio calculations, and that three-nucleon forces and continuum effects play the dominant role in deciding this question.
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.
The Modular Neutron Array (MoNA) was used in conjunction with a large-gap dipole magnet (Sweeper) to measure neutron-unbound states in oxygen isotopes close to the neutron dripline. While no excited states were observed in 24O, a resonance at 45(2) keV above the neutron separation energy was observed in 23O.