The research is partially supported by the RGC grant of Hong Kong and by the NNSF of China (Nos. 11371383 , 11271148 , 11171100 and 11371382 ), and the first author is also supported by the Fundamental Research Funds for the Central Universities and Guangdong Province Key Laboratory of Computational Science at the Sun Yat-Sen University.