We have investigated the dissolution and diffusion behaviors of hydrogen (H) in copper (Cu) based on first-principles calculations in combination with the classical thermodynamics models. A single H atom energetically prefers to occupy the octahedral interstitial site (OIS) instead of the tetrahedral interstitial site (TIS). This can be confirmed by the charge density results. The dissolved equilibrium concentration of H is 8.98×10 −7 at a typical temperature 600K. We demonstrate that the OIS→TIS→OIS path is the optimal diffusion path of H in Cu with diffusion barrier of 0.23eV. Double H atoms tend to be paired up at the two neighboring OIS’s along the 〈110〉 direction with a distance of 2.59Å and a binding energy of 0.07eV. This suggests a weak attractive interaction between H atoms, with the implication that self-trapping of H and thus formation of the H 2 molecules are quite difficult in the bulk Cu.