A large number of numerical experiments has been performed in order to study the interaction of interstellar pickup protons and helium ions with quasiparallel collisionless shocks. The shocks are modeled by a one-dimensional hybrid simulation method which treats the ions as macroparticles and the electrons as a massless fluid. Solar wind alpha particles and pickup protons are included self-consistently. In addition, the particle splitting method is used for the solar wind ions so that the distribution function can be followed over more than 10 orders of magnitude. A large part of the pickup ion distribution is reflected; the reflection efficiency is very high, and can reach in cases where the pickup ion density is low as much as 50%–60%. The reflection efficiency is almost independent of magnetic field–shock normal angle. This indicates that magnetic mirroring is unimportant and does not lead to larger reflection efficiencies. The reflection efficiency of pickup protons rapidly decreases when the pickup ion density exceeds a few percent of the solar wind density. An addition of 25% pickup protons decreases the reflection coefficient for these ions to ∼10%. This represents the fact that a quasiparallel shock cannot be considered as being uncoupled from the upstream region: at high additions of pickup ions the shock structure is changed in such a way as to reflect less pickup ions. The intensity of diffuse ions upstream of a quasiparallel shock does not depend on the temperature of the core distribution. Within the framework of the present model even solar wind distributions with a hard power law tail do not produce higher intensities of diffuse ions. It is argued that this can be understood by the fact that the intrinsic self-consistency between the processes in the upstream region and at the shock transition determines the injection and reflection properties of the core solar wind distribution.
