Crystal plasticity finite element analysis method considering the accumulation of geometrically necessary (GN) dislocations was applied to monotonic loading of pure magnesium bi-crystal. The deformation mechanisms considering in the present analysis method are basal slip , prismatic slip , 1st order pyramidal slip , 2nd order pyramidal slip and tensile twinning . Tensile twinning is incorporated into crystal plasticity analysis assuming that twinning plane and direction of shear by twinning are equivalent to slip plane and slip direction, respectively. Critical resolved shear stresses (CRSSs) for each slip system in the literatures were used. Analysis model is designed to investigate the influence of grain boundary on the activation of slip systems. That is, one grain consisting of bi-crystal (grain A) had the crystal orientation whose Schmid factor for prismatic slip is 0.5. The crystal orientation of the other grain (grain B) was slightly deviated from that of grain A. The result of the calculation of tensile loading of the bi-crystal showed that both grains are deformed by the multiple slip of basal slip system, which resulted in the formation of GN dislocation bands.