All simulations are based on DFT at the GGA level, employing the RPBE 24 exchangecorrelation functional with the D3 dispersion correction scheme. Electronic states are described using the projector-augmented wave (PAW) method and a plane wave basis set, as implemented within the VASP package. Further details on all computational settings can be found in the SI, Section S1.