An all-electron scalar relativistic calculation on Cu n H (n = 1-13) clusters has been performed by using density functional theory with the generalized gradient approximation at the PW91 level. Our results reveal that the hydrogen atom prefers to occupy the two fold coordination site for Cu n H (n = 2, 4-6, 8, 10-13) clusters, the single fold coordination site for Cu n H (n = 1, 3, 7) and the three fold coordination site for Cu 9 H cluster.