Ab initio Calculations of Multilayer Relaxations of Stepped Cu Surfaces

We present trends in the multilayer relaxations of several vicinals of Cu(100) and Cu(111) of varying terrace widths and geometry. The electronic structure calculations are based on density functional theory in the local density approximation with norm-conserving, non-local pseudopotentials in the mixed basis representation. While relaxations continue for several layers, the major effect concentrates near the step and corner atoms. On all surfaces the step atoms contract inwards, in agreement with experimental findings. Additionally, the corner atoms move outwards and the atoms in the adjacent chain undergo large inward relaxation. Correspondingly, the largest contraction (4%) is in the bond length between the step atom and its bulk nearest neighbor (BNN), while that between the corner atom and BNN is somewhat enlarged. The surface atoms also display changes in registry of upto 1.5%. Our results are in general in good agreement with LEED data including the controversial case of Cu(511). Subtle differences are found with results obtained from semi-empirical potentials.