We study the nonlinear interaction of three parallel Alfvén wave packets (AWPs) in an initially uniform plasma using 2.5-dimensional particle-in-cell (PIC) numerical simulations. We aim to help to explain the observation of suprathermal electrons by the collision of multiple Alfvén waves in regions where these waves are trapped like the IAR (Ionospheric Alfvén Resonator), Earth radiation belts or coronal magnetic loops. In the context of the acceleration by the parallel Alfvén waves interactions (APAWI) process that has been described by Mottez (Ann. Geophys., vol. 30, issue 1, 2012, pp. 81–95; J. Plasma Phys., vol. 81, issue 1, 2015, p. 325810104), the interaction of two parallel Alfvén waves (AWs) generates longitudinal density modulations and parallel electric fields at the APAWI crossing region that can accelerate particles effectively in the direction of the background magnetic field. Our simulations show that when a third parallel AWP of different initial position arrives at the APAWI crossing region, it gives rise to a strong parallel electron beam ($V \sim 5\text {--}7 V_{Te}$) at longitudinal cavity density gradients. We suggest that velocity drift from an outgoing AW generates interface waves in the transverse direction, which allows propagating waves to develop parallel electric fields by the phase mixing process when $k_{\perp }^{-1}$ of the wavy density gradient (oblique gradient) is in the range of the electron inertial length $c/\omega _{p0}$.