To partially heal the underestimation of the energy gap typical of PBE when applied to transition metals, such as Cu, we applied to the relaxed structures a Hubbard U correction to the Cu-3d orbital (UCu = 11.50 eV) obtained following the procedure reported in ref.