Schéma itératif rapide pour le calcul des matrices de transfert : application à Mo (100)

La matrice de transfert d'un solide décrite par l'empilement de couches principales est obtenue par une procédure itérative qui prend en compte 2" couches après n itérations, contrairement aux schémas usuels où chaque itération n'inclut qu'une couche de plus. La fonction de Green et la densité d'états à la surface du cristal semi-infini correspondant sont ensuite données par des formules bien connues en termes de la matrice de transfert. Cette méthode, particulièrement pratique près des singularités, est appliquée au calcul des densités d'états spectrales ainsi que totales pour la face (100) du molybdène. L'algorithme de Slater-Koster pour le calcul des paramètres de liaison forte est utilisé avec une base de neuf orbitales par atome (4d, 5s, 5p). Les états de surface et les résonances sont d'abord identifiés puis analysés en composants orbitaux pour trouver leur symétrie dominante. Leur évolution le long des principales lignes de symétrie de la zone de Brillouin bidimensionnelle est donnée explicitement. Le pic des états de surface juste en dessous du niveau de Fermi (bosse de Swanson) n'est pas obtenu. Cela est attribué à la difficulté de placer une condition de bord appropriée à la surface avec le schéma de paramétrisation de liaison forte.