Charles Edward Starrett

Pseudoatom molecular dynamics

An approach to simulating warm and hot dense matter that combines density-functional-theory-based calculations of the electronic structure to classical molecular dynamics simulations with pair interaction potentials is presented. The method, which we call pseudoatom molecular dynamics, can be applied to single-component or multicomponent plasmas. It gives equation of state and self-diffusion coefficients with an accuracy comparable to orbital-free molecular dynamics simulations but is computationally much more efficient.

Ionic Transport Coefficients of Dense Plasmas without Molecular Dynamics.

We present a theoretical model that allows a fast and accurate evaluation of ionic transport properties of realistic plasmas spanning from warm and dense to hot and dilute conditions, including mixtures. This is achieved by combining a recent kinetic theory based on effective interaction potentials with a model for the equilibrium radial density distribution based on an average atom model and the integral equations theory of fluids. The model should find broad use in applications where nonideal plasma conditions are traversed, including inertial confinement fusion, compact astrophysical objects, solar and extrasolar planets, and numerous present-day high energy density laboratory experiments.

Integral equation model for warm and hot dense mixtures

In a previous work [C. E. Starrett and D. Saumon, Phys. Rev. E 87, 013104 (2013)] a model for the calculation of electronic and ionic structures of warm and hot dense matter was described and validated. In that model the electronic structure of one atom in a plasma is determined using a density-functional-theory-based average-atom (AA) model and the ionic structure is determined by coupling the AA model to integral equations governing the fluid structure. That model was for plasmas with one nuclear species only. Here we extend it to treat plasmas with many nuclear species, i.e., mixtures, and apply it to a carbon-hydrogen mixture relevant to inertial confinement fusion experiments. Comparison of the predicted electronic and ionic structures with orbital-free and Kohn-Sham molecular dynamics simulations reveals excellent agreement wherever chemical bonding is not significant.

Potential of mean force for electrical conductivity of dense plasmas

Abstract The electrical conductivity in dense plasmas can be calculated with the relaxation-time approximation provided that the interaction potential between the scattering electron and the ion is known. To date there has been considerable uncertainty as to the best way to define this interaction potential so that it correctly includes the effects of ionic structure, screening by electrons and partial ionization. Current approximations lead to significantly different results with varying levels of agreement when compared to bench-mark calculations and experiments. We present a new way to define this potential, drawing on ideas from classical fluid theory to define a potential of mean force. This new potential results in significantly improved agreement with experiments and bench-mark calculations, and includes all the aforementioned physics self-consistently.