Read the predicted density coefficients and density fitting coefficients with FHI-aims and run one diagonalization of the KS Hamiltonian (\(\rho \rightarrow H_{KS} \rightarrow \text{everything}\))
Parse the output files of AIMS for derived properties, e.g. total/XC/electrostatic energy, forces, etc. Note that some quantities which depend only on the density (e.g. dipoles, electrostatic energy, etc.) do not require the diagonalization above.