A general electron model have been developed to study the low frequency electromagnetic modes in magnetized plasmas. Coupled to the gyrokinetic ions, the global gyrokinetic particle simulation of tearing modes have been developed and verified in the gyrokinetic toroidal code (GTC). GTC linear simulations in the fluid limit of the kink-tearing, resistive tearing modes and collisonless tearing mode in the cylindrical geometry agree well with the magnetohydrodynamic eigenvalue and initial value codes. Ion kinetic effects are found to reduce the radial width of the modes structure and the linear growth rate of tearing mode in cylinder geometry. GTC simulations of the resistive tearing modes in the toroidal geometry find that the toroidicity reduces the growth rates[2].