All computations have been carried out based on density functional theory (DFT) as implemented in the CASTEP code. The pseudo-potential is adopted by the norm conserving, and the exchange correlation functional is adopted by the Perdew-Burke-Ernzerhof in local generalized gradient approximation (GGA).