Structure of accreted neutron star crust

Using molecular dynamics simulations, we determine the structure of neutron star crust made of rapid proton capture nucleosynthesis material. We find a regular body-centered cubic lattice, even with the large number of impurities that are present. Low-charge-Z impurities tend to occupy interstitial positions, while high-Z impurities tend to occupy substitutional lattice sites. We find strong attractive correlations between low-Z impurities that could significantly increase the rate of pycnonuclear (density driven) nuclear reactions. The thermal conductivity is significantly reduced by electron impurity scattering. Our results will be used in future work to study the effects of impurities on mechanical properties such as the shear modulus and breaking strain.