In order to evaluate Electron Energy Distribution Function (EEDF) more correctly for radio frequency inductively coupled plasma (RF-ICP) in hydrogen negative ion sources, the Electromagnetic Particle-In-Cell (EM-PIC) simulation code has been improved by taking into account electron-electron Coulomb collision. Binary collision model is employed to model Coulomb collision process and we have successfully modeled it. The preliminary calculation including Coulomb collision has been done and it is shown that Coulomb collision doesn't have significant effects under the condition: electron density n
and high gas pressure pH
= 3 Pa, while it is necessary to include Coulomb collision under high electron density and low gas pressure conditions.