In this paper, a permanent magnet of multicusp ion source is calculated in theory, using the binary collision model to deal with coulomb collision between electrons. Collision between electrons and hydrogen species is treated with the null-collision model. Therefore, the three-dimensional PIC-MCC simulation algorithm is developed. Two optimization design models of multicusp ion source, ET-60U, are simulated by the software. Also, the spatial distribution characteristics and the relevant issues of the volume negative hydrogen ion production rate of the two kinds of ion sources are explored. Finally, The basic idea about the design of the negative hydrogen ion source is proposed. indicating that appropriately regulating the magnetic field distribution of the multicusp ion source can produce uniform ion beam, and that appropriately adjusting the size of the extraction magnetic field and the ion source structure can obtain both the ion beam spatial uniformity and the high production rate.