The location of the H atoms in Ti, Zr, and Hf is crucial to the formation of the hydrides in these metals as it influences the crystal lattice transformation and the hydrogen diffusion involved in the hydride formation process. Although Ti, Zr, and Hf are all of hexagonal close-packed structure with similar lattice parameters, the solute H atom occupies the octahedral interstice in Ti but the tetragonal interstice in Zr and Hf, of which the origin is still mysterious. In the present work, the origin of the distinct site occupation behavior of H atom in Ti and Zr/Hf is investigated through first principles calculations. The calculated solution energies confirm that H prefers the octahedral interstice in Ti but the tetrahedral interstice in Zr and Hf. We ascribe the distinct site occupations of H in Ti and Zr/Hf to the varying Coulomb repulsion between the H (as a screened proton in the metals) and the matrix atoms against the interstitial size. The competition between the H-induced electron accumulation effect and the matrix atom debonding effect might matter as well. We propose that, as a general rule, a H atom prefers the site with a trade-off between a large space and a high electron density in metals.