The potential energy surface of the systems was described with a deep neural network (DNN) potential trained on the density functional SCAN. SCAN calculations were performed using periodic boundary conditions, as implemented in the Quantum-ESPRESSO package. Wavefunctions and electron density were plane-wave expanded with kinetic energy cutoffs of 200 and 800 Ry, respectively. Pseudopotentials of Troulier-Martins 32 type replaced explicit core-valence electron interactions.