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.

[1]  C. G. Broyden.  A Class of Methods for Solving Nonlinear Simultaneous Equations , 1965 .

[2]  R. Piron,et al.  Variational-average-atom-in-quantum-plasmas (VAAQP) code and virial theorem: equation-of-state and shock-Hugoniot calculations for warm dense Al, Fe, Cu, and Pb. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[3]  J. Dufty,et al.  Nonempirical Semilocal Free-Energy Density Functional for Matter under Extreme Conditions. , 2016, Physical review letters.

[4]  N. Mermin.  Thermal Properties of the Inhomogeneous Electron Gas , 1965 .

[5]  Johnson,et al.  Modified Broyden's method for accelerating convergence in self-consistent calculations. , 1988, Physical review. B, Condensed matter.

[6]  W. Johnson.  Atomic structure theory , 2007 .

[7]  C. Starrett.  A Green's function quantum average atom model , 2015 .

[8]  N. M. Gill,et al.  Tartarus: A relativistic Green’s function quantum average atom code , 2017 .

[9]  B. Cichocki,et al.  Variational theory of average-atom and superconfigurations in quantum plasmas. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[10]  B. Rozsnyai.  Relativistic Hartree-Fock-Slater Calculations for Arbitrary Temperature and Matter Density , 1972 .

[11]  R. More.  Pressure Ionization, Resonances, and the Continuity of Bound and Free States , 1985 .

[12]  N. Metropolis,et al.  Equations of State of Elements Based on the Generalized Fermi-Thomas Theory , 1949 .

[13]  J. Chihara.  Unified description of metallic and neutral liquids and plasmas , 1991 .

[14]  Sonnad,et al.  Purgatorio—a new implementation of the Inferno algorithm , 2005 .

[15]  W. Kohn,et al.  Self-Consistent Equations Including Exchange and Correlation Effects , 1965 .

[16]  Ishikawa,et al.  Pressure ionization in the spherical ion-cell model of dense plasmas and a pressure formula in the relativistic Pauli approximation. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[17]  Generalized Fermi–Dirac functions and derivatives: properties and evaluation , 2001, astro-ph/0102329.

[18]  P. A. Loboda,et al.  Transport and dielectric properties of dense ionized matter from the average-atom RESEOS model , 2016 .

[19]  Xavier Gonze,et al.  A brief introduction to the ABINIT software package , 2005 .

[20]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[21]  C. Starrett.  High-temperature electronic structure with the Korringa-Kohn-Rostoker Green's function method. , 2018, Physical review. E.

[22]  Allan H. MacDonald,et al.  A relativistic density functional formalism , 1979 .

[23]  Perrot.  Dense simple plasmas as high-temperature liquid simple metals. , 1990, Physical review. A, Atomic, molecular, and optical physics.

[24]  D. Liberman.  Self-consistent field model for condensed matter , 1979 .

[25]  G. Bertsch,et al.  Optical properties of plasmas based on an average-atom model , 2006 .

[26]  D. Saumon,et al.  Electronic and ionic structures of warm and hot dense matter. , 2013, Physical review. E, Statistical, nonlinear, and soft matter physics.

[27]  M. Dharma-wardana,et al.  Density-functional theory of hydrogen plasmas , 1982 .

[28]  Donald G. M. Anderson.  Iterative Procedures for Nonlinear Integral Equations , 1965, JACM.

[29]  G. Zimmerman,et al.  A new quotidian equation of state (QEOS) for hot dense matter , 1988 .

[30]  V. Eyert.  A Comparative Study on Methods for Convergence Acceleration of Iterative Vector Sequences , 1996 .

[31]  F. Malone,et al.  Ab initio Exchange-Correlation Free Energy of the Uniform Electron Gas at Warm Dense Matter Conditions. , 2017, Physical review letters.

[32]  Travis Sjostrom,et al.  Accurate homogeneous electron gas exchange-correlation free energy for local spin-density calculations. , 2014, Physical review letters.

[33]  M. Trzhaskovskaya,et al.  Atomic structure data based on average-atom model for opacity calculations in astrophysical plasmas , 2018 .

[34]  D. Young,et al.  Theory and experiment for ultrahigh pressure shock Hugoniots , 2001 .

[35]  Stephanie B. Hansen,et al.  Equation of state, occupation probabilities and conductivities in the average atom Purgatorio code , 2006 .

[36]  A. Zunger,et al.  Self-interaction correction to density-functional approximations for many-electron systems , 1981 .

[37]  D. Liberman.  Inferno: A better model of atoms in dense plasmas , 1982 .

[38]  J. Daligault,et al.  Gradient corrections to the exchange-correlation free energy , 2014, 1408.1380.

[39]  G. E. Matthews,et al.  A Projector Augmented Wave (PAW) code for electronic structure calculations, Part I: atompaw for generating atom-centered functions , 2001 .

[40]  Peter Weinberger,et al.  Electron Scattering in Solid Matter: A Theoretical and Computational Treatise , 2004 .

[41]  Stefan Goedecker,et al.  ABINIT: First-principles approach to material and nanosystem properties , 2009, Comput. Phys. Commun..

[42]  J. Korringa,et al.  On the calculation of the energy of a Bloch wave in a metal , 1947 .

[43]  Jirí Vackár,et al.  dftatom: A robust and general Schrödinger and Dirac solver for atomic structure calculations , 2012, Comput. Phys. Commun..

[44]  Marc Torrent,et al.  Generation of Projector Augmented-Wave atomic data: A 71 element validated table in the XML format , 2013, Comput. Phys. Commun..

[45]  M. Pénicaud.  An average atom code for warm matter: application to aluminum and uranium , 2009, Journal of physics. Condensed matter : an Institute of Physics journal.

[46]  D. Saumon,et al.  A simple method for determining the ionic structure of warm dense matter , 2014 .