17 Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, A. Bansil, D. Grauer, Y. S. Hor, R. J. Cava et al. , Nature Phys. 5 , 398 (2009), 10.1038/nphys1274. 18 Y. L. Chen, J. G. Analytis, J. H. Chu, Z. K. Liu, S. K. Mo, X. L. Qi, H. J. Zhang, D. H. Lu, X. Dai, Z. Fang et al. , Science 325 , 178 (2009). 19 C.-X. Liu, H. Zhang, B. Yan, X.-L. Qi, T. Frauenheim, X. Dai, Z. Fang, and S.-C. Zhang, arXiv.org:0908.3654 (accepted by Phys. Rev. B, 2010 in press). 20 T. Jungwirth, J. Sinova, J. Maˇsek, J. Kuˇcera, and A. H. MacDonald, Rev. Mod. Phys. 78 , 809 (2006). 21 T. Dietl, H. Ohno, F. Matsukura, J. Cibert, and D. Ferrand, Science 287 , 1019 (2000). 22 T. Dietl, H. Ohno, and F. Matsukura, Phys. Rev. B 63 , 195205 (2001) 23 H. Ohno, Science 281 , 951 (1998). 24 G. Zhang, H. Qin, J. Teng, J. Guo, Q. Guo, X. Dai, Z. Fang, and K. Wu, Applied Physics Letters 95 , 053114 (pages 3) (2009). 25 Y.-Y. Li, G. Wang, X.-G. Zhu, M.-H. Liu, C. Ye, X. Chen, Y.-Y. Wang, K. He, L.-L. Wang, X.-C. Ma et al. , arXiv.org:0912.5054 (2009). 26 Y. Zhang, K. He, C.-Z. Chang, C.-L. Song, L. Wang, X. Chen, J. Jia, Z. Fang, X. Dai, W.-Y. Shan et al. , arXiv.org:0911.3706 (2009). 27 T. Zhang, P. Cheng, X. Chen, J.-F. Jia, X. Ma, K. He, L. Wang, H. Zhang, X. Dai, Z. Fang et al. , Phys. Rev. Lett. 103 , 266803 (2009). 28 V. Kul bachinskii, A. Kaminskii, K. Kindo, Y. Narumi, K. Suga, P. Lostak, and P. Svanda, JETP Letters 73 , 352 (2001), 10.1134/1.1378118. 29 J. S. Dyck, P. H´ ajek, P. Loˇst’´ ak, and C. Uher, Phys. Rev. B 65 , 115212 (2002). 30 Y. J. Chien, Transition Metal-Doped Sb$_{2}$Te$_{3}$ and Bi$_{2}$Te$_{3}$ Diluted Magnetic Semiconductors , Ph.D dissertation, University of Michigan, USA (2007). 31 J. H. Van Vleck, The Theory of Electronic and Magnetic Susceptibilities (Oxford University Press, London, 1932). 32 The parameters in Eq.(2) are determined by fitting to the ab initio results $^{16}$, which gives M$_{0}$ =- 0.28eV, A$_{xy}$ =2.2eV ˚ A, A$_{z}$ =4.1eV ˚ A, B$_{xy}$ =10eV ˚ A $^{2}$, B$_{z}$ =56.6eV ˚ A $^{2}$, C =-0.0068eV, D$_{xy}$ =1.3eV ˚ A $^{2}$, D$_{z}$ =19.6eV ˚ A $^{2}$. 33 The plane wave pseudopotential calculations were done using the BSTATE package in the framework of the Perdew-Burke-Ernzerhof-type (PBE-type) generalized gradient approxima- 13