Colossal negative thermal expansion (CNTE) materials are of significant importance for applications in precision engineering and related fields. However, simultaneously achieving a colossal negative thermal expansion coefficient (CTE) and a broad temperature range remains a central challenge. Using a combination of synchrotron x-ray diffraction, Raman spectroscopy, and density-functional-theory calculations, we find that the (, Mg, Zn) series exhibits significant CNTE during the α-β phase transition. Structural analysis indicates that during the phase transition, the pronounced contraction of the axial bonds drives the lattice volume reduction. Meanwhile, temperature-dependent Raman spectroscopy reveals that the 105 vibrational mode, assigned to antisymmetric Cu–O bond stretching, undergoes softening accompanied by significant broadening, indicating strong electron–lattice coupling in the series. Combined with energy pathway calculations, the presence of the pseudo-Jahn–Teller effect destabilizes the high-symmetry β phase and drives its transformation into the low-symmetry α phase, thereby leading to significant volume contraction upon heating. Further analysis indicates that the shallow and flattened potential energy surface at the Cu site in enables pronounced anti-off-centering displacement of Cu atoms with increasing temperature, which in turn leads to significant NTE before the phase transition, with a volumetric CTE of . Furthermore, substituting V at the P site leads to the formation of , which exhibits a colossal negative CTE () and a broad CNTE temperature range of 100–900 K. This study not only reveals a physical driving mechanism of CNTE in the series, but also provides important theoretical guidance for the design of novel functional materials that combine colossal negative thermal expansion coefficients with a broad temperature range.