size from the PC is too small. Furthermore, the possibility of gamma-ray emission from the PC, which would exhibit super-exponential cutoff due to magnetic pair attenuation, was eliminated by the spectral observation such as MAGIC and $\textit{Fermi}$ (eg, [@b9; @b1]). Also, there are not only the study of light curves and averaged spectrum but also that of phase-resolved spectra (eg, [@b10; @b11]).

In PC, SG and OG models, acceleration is caused by an electric field parallel to the magnetic field line, and the accelerating electric field is calculated by the Poisson equation for the non-corotational potential: $$- \nabla^2 \Phi = 4 \pi (\rho - \rho_{\mathrm{GJ}}),$$ where $\boldsymbol{E} = -(\boldsymbol{\Omega} \times \boldsymbol{r})
\times \boldsymbol{B}/c - \nabla \Phi$; $\rho$ and $\rho_{\mathrm{GJ}} = - \boldsymbol{\Omega} \cdot \boldsymbol{B} /
(2 \pi c) [1 - (\boldsymbol{\Omega} \times \boldsymbol{r} / c)^2]$ are the charge density and Goldreich-Julian charge density [@b12], respectively; $\boldsymbol{\Omega}$ is the angular momentum; and $\boldsymbol{r}$ is the position vector. Thus, the field-aligned electric field $E_{||} = -\nabla \Phi \cdot \hat{\boldsymbol{B}}$, with $\hat{\boldsymbol{B}}$ unit vector along the magnetic field line, is formed by the difference between local charge density and Goldreich-Julian density.

The boundary conditions of these models, which are required to solve the 