<p>The Strongly Constrained and Appropriately Normed (SCAN) meta-GGA (generalized gradient approximation) density functional and its regularized derivatives (e.g., SCAN and r²SCAN) have been proposed as a post-standard exchange-correlation functional, and are widely believed to replace conventional GGA functionals&#xa0;(e.g. PBE-GGA) owing to greatly improved electronic structures of strongly correlated systems and overall accuracy of total energies. While these improvements have been widely demonstrated for various systems, we report a significant failure of SCAN functionals related to erroneous stability of multivalent states of copper: SCAN and its derivatives (r²SCAN) critically fail to predict the relative stability of copper in oxidation states Cu^+ 1(<i>d</i>¹⁰) and Cu^+ 2(<i>d</i>⁹), excessively stabilizing Cu^+ 2 over Cu^+ 1, which leads to wrong relative stability of Cu₂O and CuO. This spurious bias also results in unphysical oxygen defect structures of YBa₂Cu₃O_7–δ for small <i>δ</i>. While the PBE-GGA functional can be fixed with a simple Hubbard-<i>U</i> correction (PBE + <i>U</i>) to predict both the spectral and thermochemical properties of copper compounds correctly, this is shown to not be the case for SCAN functionals. Our work advocates careful consideration of SCAN meta-GGA functionals when they are applied for cuprate superconductors, catalysis, and defect studies of copper compounds.</p>