The segregation energy distributions of hydrogen in four α-Fe symmetric tilt grain boundaries(GBs) are analyzed by using molecular statics(MS), and then the shear responses of four GBs embedded with different number of hydrogen atoms at the room temperature are investigated by using molecular dynamics(MD) methods. To facilitate our quantitative analysis, the hydrogen density ρ is defined as the ratio between the number of hydrogen atoms and the GB area. At different hydrogen densities, the variations of ini...