Snapshots of the non-linear evolution of the RTI at t = [318.7 s, 437.9 s, 534.5 s, 598.1 s, 612.9 s] for (x, z)∈[−L₀, L₀]×[−2.5 L₀, 2.5 L₀]. Top panels illustrate the evolution of the in-plane projection of the magnetic field lines overlaid on a grayscale colormap of the mass density of the charges. The bottom panels similarly show the spatial distribution and dynamical evolution of the mass density of the neutrals. The rectangular box marked in the t = 534.5 s panels corresponds to the zoom-in views of a reconnecting current sheet shown in Figs. 3 and B.1 below.