A quasi-linear diffusion model including both pitch angle and energy diffusion, adiabatic compression and convective motion with the solar wind flow has been used to investigate the cometary ion pickup process along the Sun-comet line at comet Halley. The total pickup ion densities and magnetic turbulence spectrum levels observed by Giotto were used to constrain the quasi-linear model. Comparisons of the model results were made with energetic ion distributions observed by instruments onboard the Giotto spacecraft. The observed power spectrum index of magnetic turbulence &ggr; is about 2--2.5. However, our simulation shows that when &ggr; was 2, the calculated proton distributions were much more isotropic than the observed ones; hence we have chosen &ggr;=2.5 in our study. Furthermore, we assumed that only about 5% of the total low frequency wave power propagates away from the comet. The numerical solutions of the quasi-linear diffusion equations show that the isotropization of the pickup ion distribution, particularly at the pickup velocity, is not complete even close to the bow shock (but upstream), which agrees with the observations. The calculated and observed proton distributions are somewhat more isotropic than the water group ion distributions in the region far upstream of the shock, but the opposite is true near the bow shock. We find that given the observed turbulence level, quasi-linear theory gives pickup ion energy distributions that agree with the observed ones quite well and easily produces energetic ions with energies up to hundreds of keV. ¿ American Geophysical Union 1993