* Corresponding authors

^a Department of Physics and Hebei Advanced Thin Film Laboratory, Hebei Normal University, Shijiazhuang 050024, China
E-mail: jwang@hebtu.edu.cn

^b Beijing Computational Science Research Center, Beijing 100193, China

^c Department of Applied Physics and Applied Mathematics, Columbia University, New York, New York 10027, USA

^d National Key Laboratory for Materials Simulation and Design, Beijing 100083, China