We investigate the residual zonal flow (Rosenbluth and Hinton 1998 Phys. Rev. Lett. 80 724) in the presence of fusion product -particles in tokamak plasmas using the modern gyrokinetic approach in the electrostatic limit. The residual zonal flow is predicted to be enhanced considerably for ITER plasmas in the radial wave-number regime of , where is the Larmor radius of a thermal ion of background plasma consisting of Deuterium and Tritium. This is a consequence of the fact that larger Larmor radius (banana orbit width, ) of an energetic -particle leads to an enhancement (a reduction) of the classical (neoclassical) polarizability at that wave-number regime compared to the case without -particles. This enhancement is slightly more pronounced for the slowing-down distribution compared to the Maxwellian if keV for MeV. In addition, we find that the Rosenbluth–Hinton formula for the residual zonal flow level that has been derived for the Maxwellian equilibrium ion distribution remains valid in the long wavelength () and high aspect ratio limit for any well-behaved ion distribution function which is isotropic in velocity space.