atoMEC: An open-source average-atom Python code

Abstract

Average-atom models are an important tool in studying matter under extreme conditions, such as those conditions experienced in planetary cores, brown and white dwarfs, and during inertial confinement fusion. In the right context, average-atom models can yield results with similar accuracy to simulations which require orders of magnitude more computing time, and thus can greatly reduce financial and environmental costs. Unfortunately, due to the wide range of possible models and approximations, and the lack of open-source codes, average-atom models can at times appear inaccessible. In this paper, we present our open-source average-atom code, atoMEC. We explain the aims and structure of atoMEC to illuminate the different stages and options in an average-atom calculation, and to facilitate community contributions. We also discuss the use of various open-source Python packages in atoMEC, which have expedited its development.

arXiv comment9 pages, 8 figures. Published in Proceedings of the 21st Python in Science Conference (SciPy 2022)

References in corpus (5)

• Finding Density Functionals with Machine Learning
• Self-averaging stochastic Kohn-Sham density functional theory
• Stochastic Density Functional Theory at Finite Temperatures
• Accelerating Finite-temperature Kohn-Sham Density Functional Theory with Deep Neural Networks
• First-principles derivation and properties of density-functional average-atom models