Direct molecular dynamics simulation of liquid-solid phase equilibria for a three-component plasma.
The neutron-rich isotope ²²Ne may be a significant impurity in carbon and oxygen white dwarfs and could impact how the stars freeze. We perform molecular dynamics simulations to determine the influence of ²²Ne in carbon-oxygen-neon systems on liquid-solid phase equilibria. Both liquid and solid phases are present simultaneously in our simulation volumes. We identify liquid, solid, and interface regions in our simulations using a bond angle metric. In general we find good agreement for the composition of liquid and solid phases between our MD simulations and the semianalytic model of Medin and Cumming. The trace presence of a third component, neon, does not appear to strongly impact the chemical separation found previously for two-component carbon and oxygen systems. This suggests that small amounts of ²²Ne may not qualitatively change how the material in white dwarf stars freezes. However, we do find systematically lower melting temperatures (higher Γ) in our MD simulations compared to the semianalytic model. This difference seems to grow with impurity parameter Q_{imp} and suggests a problem with simple corrections to the linear mixing rule for the free energy of multicomponent solid mixtures that is used in the semianalytic model.
Department of Physics and Nuclear Theory Center, Indiana University, Blooomington, Indiana 47405, USA. jhughto@astro.indiana.edu
J
J
Hughto
C J
CJ
Horowitz
A S
AS
Schneider
Zach
Z
Medin
Andrew
A
Cumming
D K
DK
Berry