ABSTRACT Diffusion coefficients are essential microphysics inputs for modelling white dwarf evolution, as they impact phase separation at crystallization and sedimentary heat sources. Present schemes for computing diffusion coefficients are accurate at weak coupling (Γ ≪ 1), but they have errors up to a factor of 2 in the strongly coupled liquid regime (1 ≲ Γ ≲ 200). With modern molecular dynamics codes, it is possible to accurately determine diffusion coefficients in select systems with per cent-level precision. In this work, we develop a theoretically motivated law for diffusion coefficients that works across the wide range of parameters typical for white dwarf interiors. We perform molecular dynamics simulations of pure systems and two mixtures that respectively model a typical-mass C/O white dwarf and a higher-mass O/Ne white dwarf, and resolve diffusion coefficients for several trace neutron-rich nuclides. We fit the model to the pure systems and propose a physically motivated generalization for mixtures. We show that this model is accurate to roughly 15 per cent when compared to molecular dynamics for many individual elements under conditions typical of white dwarfs, and is straightforward to implement in stellar evolution codes.
Read full abstract