Because of the availability of fast Fourier transforms when using PWs, their use as basis sets is efficient for GGA DFT calculations. Although the exploitation of symmetry in CRYSTAL17 using GTO basis sets would be helpful for larger unit cells with more molecules involved, we decided on PBE within VASP as a fast and accurate lower level of energy re-evaluation.