Boundary plasma physics plays an important role in tokamak confinement, but is difficult to simulate in a gyrokinetic code due to the scale-inseparable nonlocal multi-physics in magnetic separatrix and open magnetic field geometry. Neutral particles are also an important part of the boundary plasma physics. In the present paper, noble electrostatic gyrokinetic techniques to simulate the flux-driven, low-beta electrostatic boundary plasma is reported. Gyrokinetic ions and drift-kinetic electrons are utilized without scale-separation between the neoclassical and turbulence dynamics. It is found that the nonlinear intermittent turbulence is a natural gyrokinetic phenomenon in the boundary plasma in the vicinity of the magnetic separatrix surface and in the scrape-off layer.