AbstractMotivated by the observed high concentration of oxygen ions in the magnetotail during enhanced geomagnetic activity, we investigated the oxygen acceleration in magnetotail reconnection by using 2.5‐D implicit particle‐in‐cell simulations. We found that lobe oxygen ions can enter the downstream outflow region, i.e., the outflow region downstream of the dipolarization fronts (DFs) or the reconnection jet fronts. Without entering the reconnection exhaust, they are accelerated by the Hall electric field. They can populate the downstream outflow region before the DFs arrive there. This acceleration is in addition to acceleration in the exhaust by the Hall and reconnection electric fields. Oxygen ions in the preexisting current sheet are reflected by the propagating DF creating a reflected beam with a hook shape in phase space. This feature can be applied to deduce a history of the DF speed. However, it is difficult to observe for protons because their typical thermal velocity in the plasma sheet is comparable those of the DF and the reflection speed. The oxygen ions from the lobes and the preexisting current sheet form multiple beams in the distribution function in front of the DF. By comparing oxygen concentrations of 50%, 5%, and 0% with the same current sheet thickness, we found that the DF thickness is proportional to the oxygen concentration in the preexisting current sheet. All the simulation results can be used to compare with the observations from the Magnetospheric Multiscale mission.