The chemisorption of atomic hydrogen on metal surface and into subsurface is of great importance to understand the fundamental diffusion mechanism in heterogenous catalysis and hydrogen-induced embrittlement. Using spin-polarized density functional theory, we show that hydrogen prefers the quasi four-fold hollow site near typical Hollow site on Fe(110) and that four-fold hollow site on Fe(100). A very weak surface reconstruction effect induced by hydrogen coverage is reported. Two three-dimensional (3D) potential energy surfaces (PESs) are constructed for modelling hydrogen diffusion on Fe surface and into subsurface by interpolating ab initio energy points (∼1200 for each surface). We appraise the accuracy of PES and plot some contours of potential energies at different adsorption heights including the important subsurface regions. Furthermore, possible minimum energy pathways for hydrogen diffusion on Fe surface and into subsurface are searched out based on these 3D PESs using a mesh method. These pathways are in good agreement with those obtained from the nudged elastic band method. Some trapping regions into subsurface for hydrogen chemisorption are shown and the diffusion coefficient is estimated by classical transition state theory.