The authors would like to thank Dr Yong Xiao, Dr Miaohui Li for fruitful discussions. The work was supported by the National Natural Science Foundation of China under Grant no. 11775108, and the National Key R&D Program of China under Grant no. 2017YFE0300406. The numerical computations were performed on the ShenMa High Performance Computing Cluster in the Institute of Plasma Physics, Chinese Academy of Sciences.