Molecular Dynamics Simulations of Non-accreting Neutron Star Crusts

Neutron stars which do not accrete material can still have impurities in their crust, due to nuclear reactions in the crust. These impurities in the crust could affect the mechanical properties of the neutron star crust. In order to investigate the properties of the crust of a non-accreting neutron star we are performing molecular dynamic simulations of the crust. We are using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS), where simulations are run at fixed energy and volume, with the isotopes interacting via a repulsive Yukawa potential. Here we are presenting the preliminary results of the equilibrium structure of the solid neutron star crust using abundances of a non-accreting source.