We acknowledge the support of Natural Science Foundation (NSF) of China (Grant No. 11802092), NSF of Hunan Province (Grant Nos. 2019JJ50221, 2019JJ40127, 2020JJ5260, and 2020JJ4375), the Funding of the Hunan Education Department Project (No. 20A248), the Double first-class construction project of Hunan Agricultural University (No. SYL2019063), and the computation platform of the National Super Computer Center in Changsha (NSCC). We also acknowledge the helpful discussion from X. Y. Liu at The Peac Institute of Multiscale Sciences (PIMS).