Simulating the Novel Phase Separation of a Rapid Proton Capture Ash Composition

Nucleosynthesis in the oceans of accreting neutron stars can produce novel mixtures of nuclides, whose composition is dependent on the exact astrophysical conditions. Many simulations have now been done to determine the nucleosynthesis products in the ocean, but the phase separation at the base of the ocean, which determines the composition of the crust, has not been as well studied. In this work, we simulate the phase separation of a composition, which was predicted to produce a crust enriched in light nuclei, in contrast with past work which predicts that crust is enriched in heavy nuclei. We perform molecular dynamics simulations of the phase separation of this mixture using the methods of Horowitz et. al. (2007). We find good agreement with the predictions of Mckinven et al. (2016) for the phase separation of this mixture. Moreover, this supports their method as a computationally efficient alternative to molecular dynamics for calculating phase separation for a wider regime of astrophysical conditions.