*

γ-CsPbI₃ solar cells have achieved promising efficiencies, yet the quantitative understanding of their defect properties is limited due to severe computational challenges of hybrid functionals. We have discovered an algorithm to improve the convergence speed by a combination of structural relaxation with strongly constrained and appropriately normed (SCAN) Meta-generalized-gradient approximation (Meta-GGA) functional and further ionic and electronic calculations with Heyd-Scuseria-Ernzerhof (HSE) hybrid functional. The static HSE calculation with SCAN results as inputs is qualitatively reliable in defect calculations, different from one-ionic step HSE calculation based on GGA inputs. Contradictory to previous GGA defect results, a suppressed bipolar conductivity by p-type V_Cs and V_Pb, and n-type Cs_I is found. Additionally, stable bipolar defects I_int and Cs_Pb, featured by strong bond orbital coupling or structural deformation, detrimentally serve as carrier-traps. Such a strengthened bond orbital coupling in γ-CsPbI₃ causes more defect charge states than organic perovskites with larger lattice constants.