Multicomponent mutual diffusion in the warm, dense matter regime
In this work, we present the formulation, simulations, and results for multicomponent mutual diffusion coefficients in the warm, dense matter regime. While binary mixtures have received considerable attention for mass transport, far fewer studies have addressed ternary and more complex systems. We therefore explicitly examine ternary systems utilizing the Maxwell-Stefan formulation that relates diffusion to gradients in the chemical potential. Onsager coefficients then connect the macroscopic diffusion to microscopic particle motions, evinced in trajectories characterized by positions and velocities, through various autocorrelation functions (ACFs). These trajectories are generated by molecular dynamics (MD) simulations either through the Born-Oppenheimer approximation, which treats the ions classically and the electrons quantum-mechanically by an orbital-free density-functional theory, or through a classical MD approach with Yukawa pair-potentials, whose effective ionizations and electron screening length derive from quantal considerations. We employ the reference-mean form of the ACFs and determine the center-of-mass coefficients through a simple reference-frame-dependent similarity transformation. The Onsager terms in turn determine the mutual diffusion coefficients. We examine a representative sample of ternary mixtures as a function of density and temperature from those with only light elements (D-Li-C, D-Li-Al) to those with highly asymmetric mass components (D-Li-Cu, D-Li-Ag, H-C-Ag). We also follow trends in the diffusion as a function of number concentration and evaluated the efficacy of various approximations such as the Darken approximation.