$\bar{m}$ and $\bar{q}$, time steps of integration, angular velocity $\bar{\Omega}$, frequency of pair creation, critical electric field for pair creation $E_{\mathrm {cr}}$, pair multiplicity $\alpha$, frequency of injection from the stellar surface, and softening parameter of Coulomb force. It is notable that the total number of particles in the simulation box is determined self-consistently but can be controlled by $\bar{q}$ and the parameters of pair creation. We do not have any grid points for the electromagnetic field. We give typical scales of time and length for the simulation, normal pulsars, and MSPs in Table \[tab:00\] and \[tab:01\], where $\gamma_1$ and $\gamma_2$ are the Lorentz factors in units of $10$ and $10^2$, respectively; $P_{-3}$ and $P_{1}$ are the stellar rotation periods in units of $10^{-3}$ sec and $1$ sec, respectively; $B_{8}$ and $B_{12}$ are the magnetic field strength in units of $10^{8}$ G and $10^{12}$ G, respectively; $R_6$ is the stellar radius in units of $10^6$ cm; and $M_2$ is the multiplicity which is defined as $M \equiv n / n_{\mathrm {GJ}}$ in units of $10^2$. The time step was taken to be much smaller than the Larmor period to reproduce drift motions due to the electric field, the gradient of 