Automated optimization and uncertainty quantification of convergence parameters in plane wave density functional theory calculations

First principles approaches have revolutionized our ability in using computers to predict, explore, and design materials. A major advantage commonly associated with these approaches is that they are fully parameter-free. However, numerically solving the underlying equations requires to choose a set of convergence parameters. With the advent of high-throughput calculations, it becomes exceedingly important to achieve a truly parameter-free approach. Utilizing uncertainty quantification (UQ) and linear decomposition we derive a numerically highly efficient representation of the statistical and systematic error in the multidimensional space of the convergence parameters for plane wave density functional theory (DFT) calculations. Based on this formalism we implement a fully automated approach that requires as input the target precision rather than convergence parameters. The performance and robustness of the approach are shown by applying it to a large set of elements crystallizing in a cubic fcc lattice.