The MD-based α DFPT approach yields the best results compared to experiment and will be used as default unless otherwise stated. The DFPT implementation in CP2K is computationally very efficient, with the computational effort not significantly exceeding the time required to perform three additional sets of electric field calculations for each snapshot, as necessary in the induced-field method.