“…We take advantage of the 2.5D PIC simulation code of GPIC (Xiong, Huang, Yuan, Jiang, Wei, et al, 2023;. The double Harris current sheet reconnection model has been used in previous studies (e.g., Xiong et al, 2022aXiong et al, , 2022bXiong et al, , 2022cHuang et al, 2014Huang et al, , 2024Huang, Zhou, et al, 2015;Zhou et al, 2012Zhou et al, , 2014. For the initial condition of the magnetic field, we set two localized perturbations in the same current sheet, and it forms multiple X points to generate turbulent outflow due to Kelvin-Helmholtz instability (e.g., Huang et al, 2017;Huang, Lu, et al, 2015;Lu et al, 2023).…”