Particle-in-cell is a popular technique for a global five dimensional numerical simulation of the neoclassical plasma phenomena in a toroidal plasma. In this paper, we briefly review the physical and mathematical aspects of the modern neoclassical particle simulation methodology for a plasma edge simulation and present representative results recently obtained from XGC (X-point included Guiding Center) code. The strength and weakness in the modern neoclassical particle simulation techniques will also be discussed. (c) 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.