Density Functional Theory ToolKit (DFTTK) to automate first-principles thermodynamics via the quasiharmonic approximation

The Helmholtz energy is a key thermodynamic quantity representing available energy to do work at a constant temperature and volume. Despite a well-established methodology from first-principles calculations, a comprehensive tool and database are still lacking. To address this gap, we developed an open-source Density Functional Theory Tool Kit (DFTTK), which automates first-principles thermodynamics using the quasiharmonic approximation (QHA) for Helmholtz energy predictions. This Python-based package provides a solution for automating the calculation and analysis of various contributions to Helmholtz energy, including the static total energy contributions at 0 K in terms of DFT-based energy-volume curves, vibrational contributions from the Debye-Grüneisen model and phonons, and thermal electronic contributions via the electronic density of states. The QHA is also implemented to calculate the Gibbs energy and associated properties at constant temperature and pressure. The present work demonstrates DFTTK's capabilities through case studies on a simple FCC Al and various collinear magnetic configurations of Invar Fe₃Pt, where DFTTK enumerates all unique configurations and their associated multiplicities. DFTTK is freely available on GitHub, and its modular design allows for the easy addition of new workflows.