Shear viscosity for dense plasmas by equilibrium molecular dynamics in asymmetric Yukawa ionic mixtures.

We present molecular dynamics (MD) calculations of shear viscosity for asymmetric mixed plasma for thermodynamic conditions relevant to astrophysical and inertial confinement fusion plasmas. Specifically, we consider mixtures of deuterium and argon at temperatures of 100-500 eV and a number density of 10^{25} ions/cc. The motion of 30,000-120,000 ions is simulated in which the ions interact via the Yukawa (screened Coulomb) potential. The electric field of the electrons is included in this effective interaction; the electrons are not simulated explicitly. Shear viscosity is calculated using the Green-Kubo approach with an integral of the shear stress autocorrelation function, a quantity calculated in the equilibrium MD simulations. We systematically study different mixtures through a series of simulations with increasing fraction of the minority high-Z element (Ar) in the D-Ar plasma mixture. In the more weakly coupled plasmas, at 500 eV and low Ar fractions, results from MD compare very well with Chapman-Enskog kinetic results. In the more strongly coupled plasmas, the kinetic theory does not agree well with the MD results. We develop a simple model that interpolates between classical kinetic theories at weak coupling and the Murillo Yukawa viscosity model at higher coupling. This hybrid kinetics-MD viscosity model agrees well with the MD results over the conditions simulated, ranging from moderately weakly coupled to moderately strongly coupled asymmetric plasma mixtures.

Diffusivity in asymmetric Yukawa ionic mixtures in dense plasmas.

In this paper we present molecular dynamics (MD) calculations of the interdiffusion coefficient for asymmetric mixed plasma for thermodynamic conditions relevant to astrophysical and inertial confinement fusion plasmas. Specifically, we consider mixtures of deuterium and argon at temperatures of 100-500 eV and a number density ∼10(25) ions/cm(3). The motion of 30,000-120,000 ions is simulated in which the ions interact via the Yukawa (screened Coulomb) potential. The electric field of the electrons is included in this effective interaction; the electrons are not simulated explicitly. The species diffusivity is then calculated using the Green-Kubo approach using an integral of the interdiffusion current autocorrelation function, a quantity calculated in the equilibrium MD simulations. Our MD simulation results show that a widely used expression relating the interdiffusion coefficient with the concentration-weighted sum of self-diffusion coefficients overestimates the interdiffusion coefficient. We argue that this effect due to cross-correlation terms in velocities is characteristic of asymmetric mixed plasmas. Comparison of the MD results with predictions of kinetic theories also shows a discrepancy with MD giving effectively a larger Coulomb logarithm.

Self-diffusivity and interdiffusivity of molten aluminum-copper alloys under pressure, derived from molecular dynamics.

We use molecular dynamics (MD) to simulate diffusion in molten aluminum-copper (AlCu) alloys. The self-diffusivities and Maxwell-Stefan diffusivities are calculated for AlCu mixtures using the Green-Kubo formulas at temperatures from 1000 to 4000 K and pressures from 0 to 25 GPa, along with additional points at higher temperatures and pressures. The diffusivities are corrected for finite-size effects. The Maxwell-Stefan diffusivity is compared to the diffusivity calculated from the self-diffusivities using a generalization of the Darken equation. We find that the effects of cross-correlation are small. Using the calculated self-diffusivities, we have assessed whether dilute hard-sphere and dilute Lennard-Jones models apply to the molten mixture. Neither of the two dilute gas diffusivities describes the diffusivity in molten Al and Cu. We report generalized analytic models for the self-diffusivities and interdiffusivity (mutual diffusivity) that fit the MD results well. The MD-derived transport coefficients are in good agreement with the available experimental data. We also report MD calculations of the viscosity and an analytic fit to those results. The ionic thermal conductivity is discussed briefly.