This work was supported by National Natural Science Funds for Distinguished Young Scholars (Grant 21425728 ), National Key Research and Development Program of China (Grant 2016YFA0203002 ), National Science Foundation of China (Grant 51472100 ), Self-Determined Research Funds of CCNU from the Colleges’ Basic Research and Operation of MOE (Grant CCNU14Z01001 ), and Graduate Education Innovation Funding Project Grant from CCNU ( 2016CXZZ54 ). We also thank the National Supercomputer Center in Jinan for providing high performance computation.