Density functional theory calculations of copper cluster Cu n with n = 2-10 are analyzed with respect to their molecular orbitals. It is shown that shell type orbitals of s, p, and d character govern the electronic structure growth and have an influence on the geometric structure. They appear as highest occupied molecular orbital, lowest unoccupied molecular orbital, or as orbitals energetically close to these. Their occurrence can explain the magic numbers of shell models. The bonding in such structures can be explained by an analysis of the critical points of the molecular electrostatic potential. Six types of critical points are observed and related to bonds, rings, cages, lone pairs, and their connections.