The memory effect of test particles interacting with pp-wave Gaussian pulses is investigated for { higher-order polynomial pp–wave modes beyond the standard quadrupolar $+$ and $\times$ profiles}. Massive geodesic equations are solved numerically for several values of the { polynomial index $m$}, allowing the analysis of velocity and energy memory effects. In order to eliminate possible coordinate artifacts, the study is formulated in terms of the relative motion between two test particles. The results show that the relative kinetic energy variation exhibits a quartic dependence on the wave amplitude in the regime of low initial velocities. The coefficient of this scaling is found to depend on the { transverse polynomial structure of the wave, reflecting the spatial gradients associated with the different modes}. It is further shown that the energy memory effect is determined by the integrated tidal field, linking the permanent change in the relative kinetic energy to the history of the spacetime curvature carried by the wave.