The diffusion of elements is a key process in understanding the unusual surface composition of white dwarfs stars and their spectral evolution. The diffusion coefficients of Paquette et al. (1986) have been widely used to model diffusion in white dwarfs. We perform new calculations of the coefficients of inter-diffusion and ionic thermal diffusion with 1) a more advanced model that uses a recent modification of the calculation of the collision integrals that is more suitable for the partially ionized, partially degenerate and moderately coupled plasma, and 2) classical molecular dynamics. The coefficients are evaluated for silicon and calcium in white dwarf envelopes of hydrogen and helium. A comparison of our results with Paquette et al. shows that the latter systematically underestimates the coefficient of inter-diffusion yet provides reliable estimates for the relatively weakly coupled plasmas found in nearly all types of stars as well as in white dwarfs with hydrogen envelopes. In white dwarfs with cool helium envelopes (Teff < 15000K), the difference grows to more than a factor of two. We also explored the effect of the ionization model used to determine the charges of the ions and found that it can be a substantial source of discrepancy between different calculations. Finally, we consider the relative diffusion time scales of Si and Ca in the context of the pollution of white dwarf photospheres by accreted planetesimals and find factor of > 3 differences between calculations based on Paquette et al. and our model.

Models of diffusion in white dwarfs universally rely on the coefficients calculated by Paquette et al. (1986). We present new calculations of diffusion coefficients based on an advanced microscopic theory of dense plasmas and a numerical simulation approach that intrinsically accounts for multiple collisions. Our method is validated against a state-of-the-art method and we present results for the diffusion of carbon ions in a helium plasma.

The thermal conductivity of a dense {\it multi-component} plasma is critical to the modeling of accreting neutron stars. To this end, we perform large-scale molecular dynamics simulations to calculate the static structure factor of the dense multi-component plasma in the neutron star crust from near the photosphere-ocean boundary to the vicinity of the neutron drip point. The structure factors are used to validate a microscopic linear mixing rule that is valid for arbitrarily complex plasmas over a wide range of Coulomb couplings. The microscopic mixing rule in turn implies and validates the linear mixing rule for the equation of state properties and also the linear mixing rule for the electrical and thermal conductivities of dense multi-component plasmas. To make our result as useful as possible, for the specific cases of electrical and thermal conductivities, we provide a simple analytic fit that is valid for arbitrarily complex multi-component plasmas over a wide range of Coulomb couplings. We compute the thermal conductivity for a representative compositional profile of the outer crust of an accreting neutron star in which hundreds of nuclear species can be present. We utilize our results to re-examine the so-called impurity parameter formalism as used to characterize impure plasmas.