Wide Ranging Equation of State with Tartarus: a Hybrid Green's Function/Orbital based Average Atom Code

Abstract

Average atom models are widely used to make equation of state tables and for calculating other properties of materials over a wide range of conditions, from zero temperature isolated atom to fully ionized free electron gases. The numerical challenge of making these density functional theory based models work for any temperature, density or nuclear species is formidable. Here we present in detail a hybrid Green's function/orbital based approach that has proved to be stable and accurate for wide ranging conditions. Algorithmic strategies are discussed. In particular the decomposition of the electron density into numerically advantageous parts is presented and a robust and rapid self consistent field method based on a quasi-Newton algorithm is given. Example application to the equation of state of lutetium (Z=71) is explored in detail, including the effect of relativity, finite temperature exchange and correlation, and a comparison to a less approximate method. The hybrid scheme is found to be numerically stable and accurate for lutetium over at least 6 orders of magnitude in density and 5 orders of magnitude in temperature.

References in corpus (9)

• Generation of Projector Augmented-Wave atomic data: a 71 elements validated table in the XML format
• Local Spin-density Approximation Exchange-correlation Free-energy Functional
• Ab initio Exchange-Correlation Free Energy of the Uniform Electron Gas at Warm Dense Matter Conditions
• Nonempirical Semi-local Free-Energy Density Functional for Matter Under Extreme Conditions
• Volume-Collapse Transitions in the Rare Earth Metals
• Gradient corrections to the exchange-correlation free energy
• Generalized Fermi-Dirac Functions and Derivatives: Properties and Evaluation
• dftatom: A robust and general Schrödinger and Dirac solver for atomic structure calculations
• High Temperature Electronic Structure with KKR

Cited by in corpus (14)

• Multiple Scattering Theory for Dense Plasmas
• Reconciling ionization energies and band gaps of warm dense matter derived with ab initio simulations and average atom models
• First-principles derivation and properties of density-functional average-atom models
• Dense plasma opacity via the multiple-scattering method
• Time-Dependent Density Functional Theory Applied to Average Atom Opacity
• Excited States in Warm and Hot Dense Matter
• Predicting Excitation Energies in Warm Dense Matter
• The number of populated electronic configurations in a hot dense plasma
• Real-Space Green's functions for Warm Dense Matter
• Optical and Transport Properties of Plasma Mixtures from Ab Initio Molecular Dynamics
• Charge State Distributions in Dense Plasmas
• Effect of Ionic Disorder on the Principal Shock Hugoniot
• Average Atom Model with Siegert States
• Group Conductivity and Nonadiabatic Born Effective Charges of Disordered Metals, Warm Dense Matter, and Hot Dense Plasma