The spatial asymmetries in the α and β Kohn-Sham MOs require special attention as they are not identical to each other. Furthermore, their relative orders are also not guaranteed to be the same. To address this challenge, an MO reordering scheme based on the maximum orbital overlap between α and β orbitals was implemented, ensuring that the α and β pair orbitals align accordingly.