Density functional theory (DFT) calculations were performed with Gaussian 16 using the default UltraFine grid. The geometry optimizations were started from the geometry of the solid-state structures if available using the TPSSh functional and with the Ahlrichs type basis set def2-TZVP as implemented in Gaussian 16. As solvent model, the Polarizable Continuum Model (PCM) was used as implemented in Gaussian 16. As empirical dispersion correction, the D3 dispersion with Becke–Johnson damping was used as implemented in Gaussian 16.