These Bouka
These Bouka
These Bouka
THESE
Pour l’obtention du Diplôme de Doctorat
Le 08 Mars 2019
JURY
I
Remerciements
Ce travail de thèse a été réalisé à la Faculté des Sciences et Techniques de l’Université
Marien NGOUABI au sein du Groupe de Simulations numériques en Magnétisme et Ca-
talyse (GSMC) avec la collaboration de l’Institut de Chimie et Procédés pour l’Energie,
l’Environnement et la Santé (ICPEES), Strasbourg, France.
II
Second DZABANA HONGUELET, Assistants pour leur soutien.
- du corps professoral du parcours de Physique et tous les enseignants de la Faculté de
Sciences et Techniques pour les enseignements et les conseils qui ont contribué efficace-
ment à mon avancement.
III
Présentation de la structure d’accueil du stage
Le Groupe de Simulations Numériques en Magnétisme et Catalyse (GSMC) est une
unité de recherche associée aux Formations Doctorales en Sciences Physiques et Chi-
miques (FDSPC), et Sciences des Matériaux et Énergies (FDSME) de l’Université Marien
NGOUABI de Brazzaville. Le GSMC regroupe des enseignants parmi lesquels deux Pro-
fesseurs Titulaires (B. M’PASSI MABIALA et P. S. MOUSSOUNDA), deux Maîtres de
Conférences (M. NDOLLO et B. R. MALONDA BOUNGOU), trois Maîtres Assistants
(B. MALOUMBI, L. MOUKETO et D. H. DOUMA), une Assistante (M. D. NKOUA
NGAVOUKA), quatre C.E.T.P (C. C. BOUNGOU, A. L. LOMANGA OKANA, H. J.
J. MAGNOUNGOU et G. M. MOUSSOUNDA MBA) et des étudiants en thèses.
Le GSMC est rattaché à la Faculté des Sciences et Techniques de l’Université Marien
NGOUABI de Brazzaville. Le positionnement scientifique du GSMC s’articule autour
d’une recherche axée sur la physico-chimie des nanostructures. Les études à vocation
théorique sont privilégiées mais les aspects expérimentaux de la recherche ne sont pas
oubliés et sont généralement abordés en liaison avec les applications.
Les recherches sont organisées en trois grands axes :
• Axe 2 : Catalyse : Caractérisation des sites actifs par des propriétés physiques
et chimiques des surfaces des métaux de transition.
Les simulations théoriques s’appuient sur des méthodes ab-initio basées sur la fonction-
nelle de la densité. Le GSMC dispose actuellement 4 approches ab-initio représentées
par les logiciels TB-LMTO, DACAPO, PWscf et TB-DFT.
Les différentes thématiques de recherche sont éffectuées en collaboration étroite avec
d’autres laboratoires de par le monde notamment : l’Institut de Physique et Chimie
des Matériaux de Strasbourg (IPCMS), le Laboratoire des Matériaux, Surfaces et Pro-
cédés pour la Catalyse de Strasbourg (LMSPC), le Laboratoire MDM de Milan, le
Centre International de Physique Théorique (ICTP), le Centre de Physique Atomique
Moléculaire et Optique Quantique de Douala (CEPAMOQ), le Centre de Physique-
IV
Mathématique(CPM) de Mexico, l’Institut pour la Maîtrise de l’Energie (IME) de Ma-
dagascar.
V
TABLE DE MATIERES
Introduction générale 1
VI
2.1.1 Les pseudopotentiels . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.1.2 La base d’ondes planes . . . . . . . . . . . . . . . . . . . . . . . 19
2.1.3 Génération des points-k . . . . . . . . . . . . . . . . . . . . . . . 21
2.1.4 Distribution de Fermi . . . . . . . . . . . . . . . . . . . . . . . . 21
2.1.5 Energie d’échange-corrélation du Gradient Généralisé, GGA . . . 22
2.1.6 Energie de coupure . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.1.7 Calcul périodique, super-cellule ou modèle du slab . . . . . . . . . 22
2.2 Schéma de fonctionnement du code DACAPO . . . . . . . . . . . . . . . 23
2.3 Calcul de la densité d’états locale (LDOS) . . . . . . . . . . . . . . . . . 24
2.4 Calcul du travail de sortie ou Work-function “en anglais” . . . . . . . . . 26
2.5 Calcul et résultats préliminaires . . . . . . . . . . . . . . . . . . . . . . . 26
2.5.1 Paramètres structuraux . . . . . . . . . . . . . . . . . . . . . . . 26
2.5.1.1 Paramètre de maille de Au . . . . . . . . . . . . . . . . . 27
2.5.1.2 Paramètres structuraux des adsorbats . . . . . . . . . . 28
2.5.2 Energies des slabs, du bulk, de surface, de cohésion et du travail
de sortie (work-function) de l’Or . . . . . . . . . . . . . . . . . . . 28
2.5.3 Structure électronique . . . . . . . . . . . . . . . . . . . . . . . . 30
2.5.4 Résultats post calculs préliminaires . . . . . . . . . . . . . . . . . 30
2.6 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
VII
3.1.5.1.a Densité d’états projetée (PDOS) de CH3 libre
sans spin polarisé . . . . . . . . . . . . . . . . . . 39
3.1.5.1.b Densité d’états projetée (PDOS) de CH3 libre
avec spin polarisé . . . . . . . . . . . . . . . . . . 41
3.1.5.2 PDOS des bandes s, p et d de Au libre et adsorbé à 0.25,
0.50 et 0.75 MC sans et avec spin polarisé . . . . . . . . 42
3.1.5.2.a Densité d’états projetée (PDOS) des bandes s,
p et d de Au libre et adsorbé : calculs sans la
polarisation de spin . . . . . . . . . . . . . . . . . 42
3.1.5.2.b Densité d’états projetée (PDOS) des bandes s,
p et d de Au libre et adsorbé : calculs avec la
polarisation de spin . . . . . . . . . . . . . . . . . 44
3.2 Adsorption de l’hydroxyle (OH) : calculs avec et sans la polarisation de
spin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
3.2.1 Contexte bibliographique . . . . . . . . . . . . . . . . . . . . . . . 45
3.2.2 Géométrie de l’hydroxyle (OH) . . . . . . . . . . . . . . . . . . . 46
3.2.3 Sites de préférence . . . . . . . . . . . . . . . . . . . . . . . . . . 47
3.2.4 Energies d’adsorption (avec et sans spin polarisé), travail de sortie
(work function), moment dipolaire et les paramètres géométriques 47
3.2.4.1 Energies d’adsorption avec et sans spin polarisé . . . . . 48
3.2.4.2 Travail de sortie (work function) et moment dipolaire . 51
3.2.4.3 Paramètres géométriques . . . . . . . . . . . . . . . . . 53
3.2.5 Propriétés électroniques de OH et de la surface Au(001) . . . . . 54
3.2.5.1 Densité d’états projetée (PDOS) de O dans OH libre et
adsorbé à 0.25 et 0.50 MC sans et avec spin polarisé . . 54
3.2.5.2 Densité d’états projetée (PDOS) des bandes s, p et d de
Au libre et adsorbé à 0.25 et 0.50 MC sans et avec spin
polarisé . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
3.3 Thermochimie de la décomposition de CH3 et OH . . . . . . . . . . . . . 58
3.4 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
VIII
4.1.2 Géométrie et paramètres structuraux de CH2 OH en phase gazeuse 63
4.1.3 Site d’adsorption de CH2 OH sur la surface Au(001) . . . . . . . . 64
4.1.4 Energies et paramètres structuraux de CH2 OH en phase adsorbée 65
4.1.4.1 Energies d’adsorption . . . . . . . . . . . . . . . . . . . 66
4.1.4.2 Paramètres géométriques . . . . . . . . . . . . . . . . . . 66
4.1.5 Structure électronique . . . . . . . . . . . . . . . . . . . . . . . . 67
4.1.5.1 Site apical . . . . . . . . . . . . . . . . . . . . . . . . . . 67
4.1.5.2 Site ponté . . . . . . . . . . . . . . . . . . . . . . . . . . 70
4.2 Adsorption du formaldéhyde (H2 CO) . . . . . . . . . . . . . . . . . . . . 73
4.2.1 Contexte expérimental et théorique . . . . . . . . . . . . . . . . . 73
4.2.2 Géométrie et paramètres structuraux de H2 CO en phase gazeuse . 74
4.2.3 Site d’adsorption de H2 CO sur Au(001) . . . . . . . . . . . . . . . 75
4.2.4 Energies d’adsorption et paramètres structuraux du formaldéhyde
en phase adsorbée . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
4.2.5 Structure électronique . . . . . . . . . . . . . . . . . . . . . . . . 77
4.2.5.1 Site apical . . . . . . . . . . . . . . . . . . . . . . . . . . 78
4.2.5.2 Site ponté . . . . . . . . . . . . . . . . . . . . . . . . . . 79
4.3 Adsorption du formyle (HCO) sur la surface Au(001) . . . . . . . . . . . 81
4.3.1 Contexte expérimental et théorique . . . . . . . . . . . . . . . . . 81
4.3.2 Géométrie et paramètres structuraux de HCO en phase gazeuse . 82
4.3.3 Site d’adsorption de HCO sur la surface Au(001) . . . . . . . . . 83
4.3.4 Energie d’adsorption et paramètres géométriques du formyle (HCO) 83
4.3.5 Structure électronique . . . . . . . . . . . . . . . . . . . . . . . . 84
4.3.5.1 PDOS du carbone dans HCO . . . . . . . . . . . . . . . 85
4.3.5.2 Structure électronique des bandes s, p et d de Au(001) . 86
4.3.6 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
IX
5.1.5 Structure électronique . . . . . . . . . . . . . . . . . . . . . . . . 92
5.1.5.1 Site apical . . . . . . . . . . . . . . . . . . . . . . . . . . 92
5.1.5.2 Site ponté . . . . . . . . . . . . . . . . . . . . . . . . . . 94
5.1.5.3 Site de symétrie 4 . . . . . . . . . . . . . . . . . . . . . 96
5.2 Adsorption du Carbone (C) sur la surface Au(001) . . . . . . . . . . . . . 99
5.2.1 Contexte expérimental et théorique . . . . . . . . . . . . . . . . . 99
5.2.2 Sites d’adsorption, Energies d’adsorption et paramètres structuraux 99
5.2.3 Structure électronique . . . . . . . . . . . . . . . . . . . . . . . . 100
5.2.3.1 Site apical . . . . . . . . . . . . . . . . . . . . . . . . . . 101
5.2.3.2 Site ponté . . . . . . . . . . . . . . . . . . . . . . . . . . 103
5.2.3.3 Site de symétrie 4 . . . . . . . . . . . . . . . . . . . . . 105
5.3 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
Bibliographie 112
Publication 132
X
Liste des abréviations
XI
LEED : Low-energy electron diffraction (diffraction des électrons à basse énergie)
LMTO : Linear muffin-tin orbital (Orbitales muffin linéarisées)
LSDA : local spin density approximation (Approximation de la Densité Locale de
Spin)
MSD : Metastable De-excitation Spectroscopy (Spectroscopie Métastable de Désex-
citation)
129
Xe NMR : 129 Xenon nuclear magnetic resonance (résonance magnétique du 129 Xenon
nucléaire)
PBE : Perdew-Burke-Ernzerhof
PDOS : projected density of states (densité d’états projetée)
PW91 : Perdew-Wang 91
RAIRS : Reflexion-adsorption infrared spectroscopy (Spectroscopie d’adsorption-
réflexion en infrarouge)
SCAC : including single crystal adsorption calorimetry (adsorption calorimétrique
dans un cristal seul)
SFG : Sum-frequency generation spectroscopy (Spectroscopie par calcul des généra-
tions de fréquences)
SIMS : Secondary ion mass spectrometry (Spectrométrie de masse d’ions secon-
daires)
sp : spin polarisé
ssp : sans spin polarisé
STM : Scanning tunneling microscopy (microscopie à effet tunnel)
TDMS : Thermal desorption mass spectroscopy (spectroscopie de masse par désorp-
tion thermique)
TDS : Thermal desorption spectroscopy (spectroscopie par désorption thermique)
TPD : Temperature programmed desorption (spectroscopie sur la désorption en
température programmée)
TPRS : Temperature programmed reaction spectroscopy (spectroscopie de la réac-
tion à température programmée)
UPS : Ultra violet photoemission spectroscopy (spectroscopie photo émission par
des rayons UV)
XPS : X-ray photoemission spectroscopy (spectroscopie de photoémission par des
rayons X
XII
Liste des tableaux
XIII
4.3 Longueurs de liaison entre C et H (dC−H ), entre C et O (dC−O ), entre O et H (dO−H ),
entre O et Au (dO−Au ), entre C et Au (dC−Au ) et, les angles (C-O-H) et (H-C-H)
obtenus avec un slab de 4 plans après adsorption de CH2 OH sur Au(001). Les résultats
entre crochets sont calculés avec un slab de 5 plans. Les valeurs entre parenthèses
. . . . . . . . . . . . . .
représentent le nombre de longueurs de liaison possibles. 67
4.4 Comparaison entre paramètres géométriques de H2 CO libre. . . . . . . . . . . . . 74
4.5 Energies d’adsorption (Eads ), longueurs de liaison entre C et H (dC−H ), entre C et O
(dC−O ), entre O et Au (dO−Au ) et, les angles (C-O-H) et (H-C-H) obtenus avec un
slab de 4 plans après adsorption de H2 CO sur Au(001). Les résultats entre crochets
sont calculés avec un slab de 5 plans. Les nombres entre parenthèses représentent le
nombre de longueurs de liaison ou angles possibles. . . . . . . . . . . . . . . . . . 76
4.6 Comparaison entre nos paramètres structuraux calculés et plusieurs autres études théo-
riques et expérimentale pour HCO en phase gazeuse. . . . . . . . . . . . . . . . . 83
4.7 Energies d’adsorption (Eads ), longueurs de liaison entre C et H (dC−H ), entre C et O
(dC−O ), entre C et Au (dC−Au ) et l’angle (H-C-O) obtenus après adsorption de HCO
sur Au(001). Les résultats entre crochets sont calculés avec un slab de 5 plans. . . . . 83
XIV
Table des figures
XV
3.10 PDOS de O avec spin polarisé dans OH libre et OH/Au(001) sur le site ponté à 0.25
MC et 0.50 MC. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
3.11 PDOS de la surface de Au libre et adsorbée à 0.25 MC et 0.50 MC sans spin polarisé. 57
3.12 PDOS de la surface de Au libre et adsorbée à 0.25 MC et 0.50 MC avec spin polarisé. 58
3.13 Thermochimie de la décomposition des espèces moléculaires CH3 et OH. La référence
zéro représente l’énergie de l’espèce en phase gazeuse plus l’énergie de la surface libre
de Au(001) sur différentes étapes. Les phases gazeuses sont indiquées par (g) et l’espèce
adsorbée par (a). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
4.1 Géométrie de l’hydroxyméthyle (CH2 OH) en phase gazeuse. La sphère en rouge repré-
sente l’atome d’oxygène, en noir l’atome de carbone et en blanc l’atome d’hydrogène. 63
4.2 Géométrie de l’hydroxyméthyle (CH2 OH) sur la surface de Au(001) avant (a) et après
(b) adsorption sur le site de symétrie 4. La sphère en rouge représente l’atome d’oxy-
gène, en noir l’atome de carbone et en blanc l’atome d’hydrogène. . . . . . . . . . . 64
4.3 Géométrie de l’hydroxyméthyle (CH2 OH) sur la surface de Au(001) avant (a) et après
(b) adsorption sur le site de apical. La sphère en rouge représente l’atome d’oxygène,
en noir l’atome de carbone et en blanc l’atome d’hydrogène. . . . . . . . . . . . . 65
4.4 Géométrie de l’hydroxyméthyle (CH2 OH) sur la surface de Au(001) avant (a) et après
(b) adsorption sur le site de ponté. La sphère en rouge représente l’atome d’oxygène,
en noir l’atome de carbone et en blanc l’atome d’hydrogène. . . . . . . . . . . . . 65
4.5 PDOS sur les orbitales O de CH2 OH libre (en pointillé) par rapport à celles dans
CH2 OH/Au(001) adsorbées (en trait continu) pour le site apical. Le niveau de Fermi
est à 0 eV. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
4.6 PDOS sur les orbitales H de CH2 OH libre (en pointillé) par rapport à celles dans
CH2 OH/Au(001) adsorbées (en trait continu) pour le site apical. Le niveau de Fermi
est à 0 eV. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
4.7 PDOS des bandes de Au libre (en pointillé) par rapport à celles dans CH2 OH/Au(001)
adsorbées (en trait continu) sur le site apical. Le niveau de Fermi est à 0 eV. . . . . 70
4.8 PDOS sur les orbitales O de CH2 OH libre (en pointillé) par rapport à celles dans
CH2 OH/Au(001) adsorbées (en trait continu) pour le site ponté. Le niveau de Fermi
est à 0 eV. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
4.9 PDOS sur l’orbitale H de CH2 OH libre (en pointillé) par rapport à celles dans CH2 OH/Au(001)
adsorbées (en trait continu) pour le site ponté. Le niveau de Fermi est à 0 eV. . . . . 72
4.10 PDOS des bandes de Au libre (en pointillé) par rapport à celles dans CH2 OH/Au(001)
adsorbées (en trait continu) sur le site ponté. Le niveau de Fermi est à 0 eV. . . . . . 73
XVI
4.11 Géométrie du formaldéhyde en phase gazeuse. La sphère en rouge représente l’oxygène,
en noir le carbone et en blanc les hydrogènes. . . . . . . . . . . . . . . . . . . . 75
4.12 PDOS sur l’orbitale O de H2 CO libre (en pointillé) par rapport à celles dans H2 CO/Au(001)
adsorbées (en trait continu) pour le site apical. Le niveau de Fermi est à 0 eV. . . . . 78
4.13 PDOS des bandes de Au libre (en pointillé) par rapport à celles dans H2 CO/Au(001)
adsorbées (en trait continu) sur le site apical. Le niveau de Fermi est à 0 eV. . . . . 79
4.14 PDOS sur l’orbitale O de H2 CO libre (en pointillé) par rapport à celles dans H2 CO/Au(001)
adsorbées (en trait continu) pour le site ponté. Le niveau de Fermi est à 0 eV. . . . . 80
4.15 PDOS des bandes de Au libre (en pointillé) par rapport à celles dans H2 CO/Au(001)
adsorbées (en trait continu) sur le site ponté. Le niveau de Fermi est à 0 eV. . . . . . 81
4.16 Géométrie du formyle (HCO) en phase gazeuse. La sphère en rouge représente l’atome
d’oxygène, en noir l’atome de carbone et en blanc l’atome d’hydrogène. . . . . . . . 82
4.17 PDOS sur l’orbitale de HCO libre via C (en pointillé) par rapport à celles dans
HCO/Au(001) adsorbées (en trait continu) pour le site apical. Le niveau de Fermi
est à 0 eV. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
4.18 PDOS des bandes de Au libre (en pointillé) par rapport à celles dans HCO/Au(001)
adsorbées (en trait continu) sur le site apical. Le niveau de Fermi est à 0 eV. . . . . 87
XVII
5.9 PDOS des bandes de Au libre (en pointillés) par rapport à celles dans C/Au(001)
adsorbées (en trait continu) sur le site apical. Le niveau de Fermi est à 0 eV. . . . . 103
5.10 PDOS sur les orbitales C du carbone libre (en pointillés) par rapport à celles dans
C/Au(001) adsorbées (en trait continu) pour le site ponté. Le niveau de Fermi est à 0
eV. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104
5.11 PDOS des bandes de Au libre (en pointillés) par rapport à celles dans C/Au(001)
adsorbées (en trait continu) sur le site ponté. Le niveau de Fermi est à 0 eV. . . . . . 105
5.12 PDOS sur les orbitales du carbone (C) libre (en pointillés) par rapport à celles dans
C/Au(001) adsorbées (en trait continu) pour le site de symétrie 4 . Le niveau de Fermi
est à 0 eV. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 106
5.13 PDOS des bandes de Au libre (en pointillés) par rapport à celles dans C/Au(001)
adsorbées (en trait continu) sur le site de symétrie 4. Le niveau de Fermi est à 0 eV. . 107
XVIII
Introduction générale
Les découvertes dans le domaine très vaste des sciences expérimentale et théorique
en général et dans le domaine très particulier de la physique et la chimie combinatoires
ont permis de mettre en exergue un nombre incommensurable de réactions physico-
chimiques qui peuvent être bénéfiques ou alarmantes pour l’humanité. Parmi les si-
tuations alarmantes, on peut citer l’exploitation et la combustion des énergies fossiles
(pétrole, charbon, huile de schiste, etc) qui génèrent des polluants tels que les oxydes de
soufre, d’azote et de carbone. Ceux-ci sont potentiellement des sources à effet de serre,
qui une fois rejetés dans l’atmosphère, modifient l’équilibre environnemental mondial.
L’enjeu majeur est de trouver d’autres sources alternatives aux énergies fossiles. On
peut y arriver en valorisant les molécules hydrocarbonées, hydro-oxygénées ou hydro-
sulfurées issues de la biomasse. Parmi celles-ci, nous retiendrons les fragments molé-
culaires tels que : CH3 , OH, CH2 OH, H2 CO et HCO. Concrètement, ceux-ci peuvent
être de bons candidats pour l’extraction de l’hydrogène qui peut être utilisé comme une
nouvelle source d’énergie dite propre ; car sa combustion ne produit que la vapeur d’eau.
Quoiqu’il en soit, pour atteindre cet objectif, il faut au préalable l’analyse des pro-
priétés physico-chimiques des fragments moléculaires lorsqu’ils sont en contact avec l’en-
vironnement.
De manière générale, l’étude des adsorbats sur les surfaces des métaux est d’une
grande importance en science des surfaces grâce aux applications industrielles en catalyse
[1-3]. Mais qu’est-ce que la catalyse ? Le terme « catalyse » désigne un processus qui vise à
modifier la cinétique chimique en abaissant la barrière énergétique à franchir. L’objectif
de la catalyse est d’augmenter la vitesse à laquelle se produit une réaction chimique
ou, parfois, de privilégier une réaction plutôt qu’une autre. Généralement la réaction
catalytique se fait en présence d’un catalyseur, c’est-à-dire une substance liquide, solide
ou gazeuse. Dans notre étude, on se limitera à la catalyse hétérogène pour laquelle les
phases du catalyseur (ici le solide) et des réactifs sont différentes.
L’enjeu de cette thèse est d’utiliser l’Or comme substrat. Mais pourquoi l’Or ? Long-
temps considéré comme métal inerte en catalyse, l’intérêt pour les propriétés catalytiques
1
des nanoparticules d’or a émergé à la fin des années 1980, lorsque Haruta et ses collabo-
rateurs ont reporté que ces particules sont des catalyseurs efficaces pour l’oxydation du
monoxyde de carbone (CO) aux basses températures [4, 5]. Dans la même décennie, une
autre découverte clé a été faite par Hutchings, quand il a rapporté que les nanoparticules
d’or sont de très bons catalyseurs sélectifs pour l’hydrochloration de l’acétylène [6]. Tou-
tefois, la commercialisation des catalyseurs en Or n’a commencé qu’une décennie plus
tard, vers 1992, lorsqu’ils ont été prouvés pour décomposer l’ammoniac et la triméthyla-
mine à la température ambiante, et ont été utilisés principalement pour décomposer les
odeurs [7]. Une caractéristique particulière qui a rendu ces catalyseurs si efficaces pour
cette application est que l’humidité améliore leur activité. Depuis cette commercialisa-
tion, les nanoparticules d’or ont également montré des activités très élevées de sélectivité
pour de nombreuses autres réactions, en particulier pour les réactions d’oxydation sélec-
tive et d’hydrogénation. Certaines des réactions dans lesquelles des catalyseurs d’or ont
été prouvés utiles incluent : l’époxydation du propylène [8, 9], l’oxydation du glycérol
[10-12], l’oxydation de l’éthylène glycol [13], la synthèse directe du peroxyde d’hydrogène
et de l’oxygène [14-16], l’hydrogénation de l’acétylène à l’éthylène [17] et réduction de
NO par H2 [18].
Les développements considérables des recherches sur les catalyseurs en Or sont reliés
à des intérêts surtout appliqués. Parallèlement au progrès des techniques expérimentales,
à la fois de production et d’analyse catalytiques, les méthodes théoriques pour étudier
les interactions adsorbats-catalyseur en Or ont connu un essor considérable ces dernières
années ; notamment il est vrai grâce au développement des outils informatiques de plus
en plus puissants. Les méthodes ab initio de la chimie quantique ou de la physique de
la matière condensée de type fonctionnelle de la densité sont aujourd’hui bien adaptées
pour étudier les interactions adsorbats-surface d’or.
Ce mémoire est consacré à « l’adsorption des espèces atomique (C), moléculaire (CO)
et des fragments moléculaires (CH3 , OH, CH2 OH, H2 CO et HCO) sur la surface Au(001)
». L’étude des fragments moléculaires est motivée par le fait qu’on peut extraire de l’hy-
drogène utilisable comme carburant idéal pour le futur, comme nous l’avons déjà évoqué
ci-haut. Quant aux espèces atomique et moléculaire, leur étude permet d’élucider les dif-
férentes étapes dans la réaction de dissociation du CO2 , qui est un gaz potentiellement
à effet de serre.
Comme, nous l’avons déjà dit, le métal utilisé sera l’or (Au) qui cristallise dans un
réseau cubique à faces centrées. Les études expérimentales et théoriques sur l’adsorption
des espèces atomiques et moléculaires ont été essentiellement menées sur la face (111)
2
de Au du fait qu’elle est plus stable et plus dense que les deux autres faces Au(110) et
Au(001) [19]. Pendant qu’il a été observé que la surface Au(001) se reconstruit selon deux
structures décrites, soit par un arrangement d’atomes à courte portée représenté par une
structure (5 x N) [20, 21], soit par un arrangement d’atomes à longue distance décrit par
une superstructure c(28 x 28) [22], notre étude ne tient compte que des réactions sur la
surface non-reconstruite Au(001).
Ce travail est donc une contribution théorique à l’étude de la réactivité de la sur-
face Au(001) lorsque sont adsorbés des espèces atomique (C), moléculaire (CO) et des
fragments moléculaires (CH3 , OH, CH2 OH, H2 CO, et HCO). Les simulations sont me-
nées dans le cadre de la théorie de la fonctionnelle de la densité, en ayant recours au
code DACAPO [Dacapo (version 2.7.3) Danish Code Ab-initio Pseudopotential] qui est
un programme de calcul élaboré à CAMP (Center for Atomic-scale Materials Physics)
au Département de Physique de l’Université Technique de Danemark. Ce programme,
fondé sur les principes de la mécanique quantique et de la physique du solide, calcule des
données essentielles du système à étudier et liées à son « énergétique » (énergie totale,
énergie de liaison, affinité, énergie d’ionisation, énergie d’atomisation...), à son moment
magnétique, à la densité électronique, aux fréquences de vibration, au travail de sortie
et aux charges des atomes qui forment le système. De plus, il permet l’optimisation de
la structure atomique par minimisation des contraintes entre les atomes. Dacapo utilise
une base d’ondes planes pour représenter les orbitales atomiques. L’action des électrons
profonds, ainsi que celle du noyau des atomes, qui sont considérés comme gelés dans
le programme, sont régies par des pseudo-potentiels « ultrasoft » de Vanderbilt. Les
équations mono-électroniques de Kohn-Sham sont résolues de façon auto-cohérente. Le
système, à étudier, est représenté sous forme de film mince (slab) dans le calcul et re-
produit, selon la technique de la supercellule, pour compenser la cassure de la symétrie
suivant la direction perpendiculaire à la surface du slab [23, 24].
Dans ce manuscrit, nous présentons dans le premier chapitre le cadre théorique des
méthodes ab initio que nous avons utilisées. Nous abordons d’une part, les approxi-
mations de Born-Oppenheimer, de Hartree-Fock-Slater, et d’autre part, la théorie de
la fonctionnelle de densité (DFT). Enfin, nous terminerons par les approximations des
fonctionnelles d’échange-corrélation et des pseudopotentiels.
Dans le deuxième chapitre, nous présentons le code Dacapo et son fonctionnement.
Nous présentons ensuite les résultats préliminaires sur le paramètre de maille de Au
en volume, l’énergie de surface, l’énergie de cohésion et la représentation de la surface
Au(1x1) (001) avec les sites d’adsorption : apical (top), ponté (bridge) et de symétrie 4
3
(hollow).
Le chapitre 3 est consacré à l’adsorption du méthyle (CH3 ) et de l’hydroxyle (OH) sur
la surface Au(001). Nous présentons d’abord les résultats sur CH3 /Au(001) avec l’analyse
des énergies d’adsorption avec et sans polarisation de spin, des paramètres géométriques,
les structures électroniques (travail de sortie, moment dipolaire et densité des états
projetés-PDOS). Dans une deuxième partie, nous analysons les résultats sur l’adsorption
de OH/Au(001), et nous terminons ce chapitre par l’analyse de la thermochimie pour la
décomposition de CH3 et OH.
Le chapitre 4 développe les études sur les adsorptions sur l’hydroxyméthyle (CH2 OH),
du formaldéhyde (H2 CO) et du formyle (HCO). Les discussions des résultats reprennent
la même démarche élaborée au chapitre précédent.
Dans le dernier chapitre, nous exposons les résultats sur l’adsorption du carbone (C)
et du monoxyde de carbone (CO). Ici aussi, nous étudions la réactivité de la surface
de Au(001) à travers l’analyse des énergies d’adsorption et les densités d’état locales
(PDOS).
Enfin, une conclusion générale résumera les principaux résultats de cette thèse et
ensuite, nous terminerons par les perspectives.
4
Chapitre 1
Le système que nous allons étudier est constitué du substrat (Au) et des adsorbats
[c.à.d. des espèces atomique (C), moléculaire (CO) et des fragments moléculaires (CH3 ,
OH, CH2 OH, H2 CO et HCO)] constitués de plusieurs particules. L’étude des propriétés
physiques des atomes, des molécules, des nanostructures à la surface nécessite une grande
connaissance des états électroniques mis en jeu. En théorie de la matière condensée, la
description quantique de tel système se fait par le biais de la résolution de l’équation de
Schrödinger.
Dans ce chapitre, nous supposons que le système est formé de L électrons et K
noyaux en interaction. Les équations seront écrites en notation du système d’unités
atomiques pour lesquelles : me = e = ~ = 4πεo = 1 ( me = masse électronique, e =
h
charge électronique et } = 2π = constante de réduite de Planck).
→
− →
−
Hψs (→
−
r i , R α ) = Es ψs (→
−
r i, R α) (1.1)
où H est l’opérateur hamiltonien du système décrivant l’ensemble des interactions qui s’y
produisent et constitué de K noyaux de masse Mα , de numéro atomique Zα , de position
5
→
−
R α , et de L électrons de position →
−
r i ; Es est l’énergie totale du système (K noyaux +
L électrons). Pour un système donné, l’hamiltonien total du système s’écrit :
H =T +V (1.2)
T = Te + TN (1.3)
V = Vee + VN e + VN N (1.4)
l’hamiltonien devient :
H = Te + TN + Vee + VN e + VN N (1.5)
avec :
1X
Te = − ∆i (1.6)
2 i
1 X ∆α
TN = − (1.7)
2 α Mα
X 1
Vee = (1.8)
i<j
| ri − →
→
− −
rj |
X Zα
VN e = − −
→ → (1.9)
i,α |Rα − −
ri |
6
X Zα Zβ
VN N = −
→ − → (1.10)
α<β |Rα − Rβ |
où les termes ci-dessus représentent respectivement : énergie cinétique des électrons (1.6) ;
énergie cinétique des noyaux (1.7) ; énergie d’interaction coulombienne électron-électron
(1.8) ; énergie d’interaction électrostatique noyau-électrons (1.9) et énergie d’interaction
noyau-noyau (1.10).
L’hamiltonien H total du système de L électrons et K noyaux [24] s’écrit donc :
1X X 1 X Zα 1 X ∆α X Zα Zβ
H=− ∆i + →
− →
− − −
→ − + → (1.11)
2 i i<j
| ri − rj | →
−
i,α |Rα − ri |
2 α Mα α<β |−
→ −
Rα − Rβ |
L’ équation de Schrödinger (1.1) avec un tel hamiltonien (1.11), ne peut être résolue
qu’en utilisant certaines approximations que nous exposons dans ce qui suit.
1.1.2 Approximation de Born Oppenheimer
Cette approximation, dite adiabatique basée sur la simplification de la résolution de
l’équation de Schrödinger (1.1), consiste à séparer la partie électronique de la partie
−
→
nucléaire dans la fonction d’onde ψs (→
−
ri , Rα ). En effet, l’énorme différence entre la masse
électronique me et la masse nucléaire Mα (Mα = 1836 me ) nous conduit à supposer que
le mouvement des noyaux est très lent par rapport à celui des électrons.
La fonction d’onde du système peut alors s’écrire comme le produit de deux fonctions :
→
− −
→
ψs (→
−
ri , R α ) = φ− →
−
→ ( r i ).χ(Rα )
Rα
(1.12)
→
− −→
où φ− → ( r i ) et χ(Rα ) sont des fonctions d’ondes décrivant respectivement le mouvement
Rα
des électrons et des noyaux.
La résolution de l’équation de Schrödinger du système revient donc à la recherche
de la fonction d’onde découplée (1.12) vérifiant l’équation (1.1), alors que les fonctions
→
− −→
d’onde φ− → ( r i ) et χ(Rα ) vérifient respectivement les équations de Schrödinger
Rα
!
1X X Zα X 1
− ∆i − + φRα (r) = Eel (R)φRα (r) (1.13)
2 i α,i
|Rα − ri | i<j |ri − rj |
7
et
!
1 X ∆α X Zα Zβ
− + + Eel χ(R) = Es χ(R) (1.14)
2 α Mα α<β |Rα − Rβ |
→
−
avec : r = →
−r i , R = R α , Eel = énergie électronique du système. On peut en outre
souligner que la contribution du terme d’énergie cinétique des noyaux à l’énergie to-
tale du système est négligeable (noyaux supposés fixes dans l’approximation de Born-
Oppenheimer). Ainsi, l’énergie totale du système pour une configuration R des noyaux
s’écrit :
X Zα Zβ
Es ' Eel + (1.15)
α<β
|Rα − Rβ |
où le terme
X Zα Zβ
= cte (noyauxf ixes). (1.16)
α<β
|Rα − Rβ |
On peut conclure alors que le calcul de l’énergie totale et d’autres propriétés physico-
chimiques du système dans le cadre de l’approximation de Born-Oppenheimer se réduit
à la résolution de l’équation de Schrödinger suivante :
" #
1X X Zα X 1
− ∆i − + φRα (r) = Eel (R)φRα (r) (1.17)
2 i α,i
|Rα − ri | i<j
|r i − rj |
on remarque que le dernier terme du membre de gauche entre crochets de (1.17) est un
terme bi-électronique alors que les deux premiers sont mono-électroniques ; ce qui pose
une difficulté majeure pour le traitement du système électronique. Ainsi on fait recours
à d’autres approximations supplémentaires.
1.1.3 Approximation de Hartree et Hartree-Fock
1.1.3.1 Approximation de Hartree
L’équation (1.17) est très complexe à résoudre du fait d’un grand nombre d’électrons
mis en jeu. L’approximation de Hartree consiste donc à réduire le problème de L corps à
celui d’une particule. Dans cette approximation, la fonction d’onde totale φ(→
−
r1 , →
−
r2 , ..., →
−
rL )
est écrite sous forme d’un produit de fonctions mono-électroniques :
8
φ(→
−
r 1, →
−
r 2 , ..., →
−
r L ) = ϕ1 (→
−
r 1 )ϕ2 (→
−
r 2 )...ϕL (→
−
r L) (1.18)
L’équation de Schrödinger à une particule ou (équation de Hartree) est définie donc par :
hi ϕi (→
−
ri ) = εi ϕi (→
−
ri ) (1.19)
avec
1
hi = − ∆i + Vext (→
−
r ) + Vi (→
−
r) (1.20)
2
où hi est l’hamiltonien à un électron ; Vext décrit à la fois le potentiel des interactions élec-
trostatiques noyau-noyau et celles des autres électrons-noyaux c.à.d les (L - 1) électrons-
noyaux dans le système ; Vi (→−r ) le potentiel de Hartree pour le ieme électron définie par
[24] :
→
−
− ρi ( r0 )
→
Z
Vi (→
−
r)= 3 0
d r →
− (1.21)
|→
−
r − r0 |
avec
→
− X →
−
ρi ( r0 ) = |ϕj ( r0 )|2 (1.22)
i6=j
→
−0 2
→
−
1 →
− →
−
XZ
3 0 |ϕj ( r )| →
− →
−
− ∆i + Vext ( r ) ϕi ( r ) + d r −0 ϕi ( r ) = εi ϕi ( r )
→ (1.23)
2 →
−
|r − r |
i6=j
9
Cette formulation de Hartree est très difficile à mettre en œuvre, surtout si le nombre
d’électrons L est grand. On admet de plus que la densité mono-électronique ρi (→ −r ) est
→
−
égale à la densité électronique ρ( r ) :
ρ(→
−
r ) = ρi (→
− |ϕj (→
−
X
r)= r )|2 (1.24)
i6=j
Mais la tâche est encore plus compliquée du fait que le potentiel subit par chaque
électron est le même, on introduit donc une interaction de chaque électron avec lui même
ce qui est incorrect pour les systèmes localisés tels que les atomes.
1.1.3.2 Approximation de Hartree-Fock
En 1930, Fock [26] a montré que la fonction d’onde de Hartree (1.18) viole le prin-
cipe d’exclusion de Pauli parce qu’elle ne respecte pas l’antisymétrisation de la fonction
d’onde électronique par permutation de deux électrons quelconques. La fonction d’onde
monoélectronique doit être changée par un déterminant de Slater [27] :
ϕ1 (→ −
r1 ) ϕ1 (→ −
r2 ) · · · ϕ1 (→ −
rL )
→
− →
−
ϕ2 ( r1 ) ϕ2 ( r2 ) · · · ϕ2 (rL )→
−
1
φ(→
−
r1 , →
−
r2 , ..., →
−
rL ) = √ .. .. .. .. (1.25)
L!
. . . .
ϕL ( r1 ) ϕL ( r2 ) · · · ϕL (→
→
− →
− −
rL )
→
−0 2 !
1 →
− |ϕ ( )|
Z
r
− ∆i + Vext (→
− →
− →
−
X j
r )+ d3 r0 −0 − Vx ϕi ( r ) = εi ϕi ( r )
→ (1.26)
2 →
−
|r − r |
i6=j
→
−
XZ − ϕ∗i ( r0 )ϕj (→
→ −r) →
Vx ϕi (→
−
r)=− 3 0
d r →
− ϕi (−
r) (1.27)
→
−
|r − r | 0
i6=j
10
est appliquée à un métal puisque la méthode de Hartree-Fock prédit une densité d’états
nulle à la surface de Fermi d’un métal. Ce désaccord s’explique par un phénomène
physique important absent dans la théorie de Hartree-Fock : la corrélation électronique.
1.1.4 Approximation de Hartree-Fock-Slater
Contrairement à l’idée de Hartree-Fock, Slater [27] fait une approximation du terme
du potentiel d’échange en supposant qu’il possède un caractère local et postule la forme
de ce potentiel :
3ρ(→
−
1
r) 3
Vx (→
−
r ) = −6α (1.28)
4π
11
- Théorème 1 :
Z
E[ρ] = d3 →
−
r Vext (→
−
r )ρ(→
−
r ) + F [ρ] (1.29)
avec F [ρ] = T [ρ] + Eee [ρ] , fonction de Hohenberg-Kohn ; T [ρ] : fonctionnelle énergie
cinétique ; Eee [ρ] : fonctionnelle énergie d’interaction électron-électron.
- Théorème 2 :
Ce théorème est basé sur la recherche de l’énergie minimun de E[ρ]. Cette énergie
est minimun lorsque la densité électronique quelconque ρi (→
−
r ) correspond à la densité
→
−
électronique de l’état fondamental ρEF ( r ). On a :
Z
d3 →
−
r Vext (→
−
r )ρEF (→
−
r ) + F [ρEF ] = EEF (1.30)
→
−
− ρ(→
→ −r )ρ( r0 )
Z Z
E[ρ(→
−
r )] = Ts [ρ(→
−
r )] + d r Vext (→
3→
− −
r )ρ(→
−
r )+ 3→
− 3 0
d rd r − +
→ (1.31)
|→
−
r − r0 |
Exc [ρ(→
−
r )]
avec
- Ts [ρ(→
−
r )], fonctionnelle de l’énergie cinétique d’un gaz d’électrons non interagissants ;
R 3→
- d−r Vext (→
−
r )ρ(→
−r ), fonctionnelle de l’énergie d’interaction avec le potentiel exté-
12
rieur ;
- Exc [ρ(→−
r )], fonctionnelle de l’énergie d’échange-corrélation qui contient non seule-
ment les effets d’échange et corrélation, mais aussi les contributions cinétiques et cou-
lombiennes liées à l’interaction entre les électrons.
→
−
δE[ρ(→
−r )] δTs [ρ(→
−
r )] − ρ( r0 )
→ δExc [ρ(→
−
r )]
Z
= →
− 3 0
+ Vext ( r ) + d r + =µ (1.32)
→
−
δρ( r ) →
−
δρ( r ) →
− →
− →
−
δρ( r )
| r − r0 |
δE[ρ(→
−r )] δTs [ρ(→
−
r )]
→
− = →
− + Vef f (→
−
r)=µ (1.33)
δρ( r ) δρ( r )
→
−
− ρ( r0 )
→ δExc [ρ(→
−
r )]
Z
Vef f (→
−
r ) = Vext (→
−
r )+ 3 0
d r →
− + →
− (1.34)
|→
−
r − r0 | δρ( r )
et
δExc [ρ(→
−
r )]
Vxc (→
−
r)= →
− (1.36)
δρ( r )
13
L
ρ(→
− |ϕi (→
−
X
r)= r )|2 (1.37)
i=1
→
− 1 →
− → − KS →
−
HKS ϕi ( r ) = − 2 ∆i + Vef f ( r ) ϕi ( r ) = εi ϕi ( r )
(1.38)
ρ(→
−
r ) = i fi |ϕi (→
−
r )|2
P
avec
- εi KS : valeur propre mono électronique de l’hamiltonien, HKS , de Kohn et Sham ;
- ϕi (→
−r ) : fonction d’onde mono électronique ;
- fi : nombre d’occupation de l’orbitale i
avec : N
P occ
i=1 fi = L et Nocc nombre d’états occupés.
Pour les systèmes polarisés en spin,
L ↑ L↓
ρ(→
−
r ) = ρ↑ (→
−
r ) + ρ↓ (→
− |ϕi ↑ (→
− |ϕi ↓ (→
−
X X
r)= r )|2 + r )|2 (1.39)
i=1 i=1
1
où fi = 1 car chacun des spins contribue pour 2
;
ρ↑ (→
−
r ) : densité d’états de spin up (↑) ;
ρ↓ (→
−
r ) : densité d’états de spin down (↓) ;
L = L↑ + L↓
→
− ρ↑ (→
−
r ) − ρ↓ (→
−
r)
η( r ) = →
− (1.40)
ρ( r )
14
du système devient :
→
−
L
− ρ(→
→ −r )ρ( r0 )
Z Z
E[ρ(→
− d r Vext (→
3→
− −
r )ρ(→
− 3→
− 3 0 →
−
X
r )] = εi − r )− d rd r −0 + Exc [ρ( r )] (1.41)
→
→
−
|r − r |
i=1
Dans cette expression de l’énergie totale de Kohn-Sham, les termes de l’énergie ci-
nétique et leur potentiel sont déterminables sauf celui d’échange-corrélation qui n’est
pas connu exactement et, est source de problème. Nous devons donc recourir à d’autres
approximations.
Z
LDA
Exc [ρ(→
−
r )] = d3 →
−
r εLDA →
−
xc [ρ( r )] (1.42)
où εLDA →
−
xc [ρ( r )] est l’énergie d’échange-corrélation par particule d’un système d’électrons
homogène de densité ρ(→ −r ).
On peut l’exprimer en deux termes :
LDA
Exc [ρ(→
−
r )] = ExLDA [ρ(→
−
r )] + EcLDA [ρ(→
−
r )] (1.43)
avec :
1 − 1
- ExLDA [ρ(→
−
r )] = − 43 π3 3 ρ(→
r ) 3 , est le terme d’énergie d’échange déduite du modèle de
Thomas-Fermi-Dirac [33] ;
- EcLDA [ρ(→
−
r )] exprime l’énergie de corrélation plus complexe à déterminer ; cependant
15
elle peut être évaluée de différentes façons, par exemple à l’aide des calculs Monte-Carlo
quantiques [25]. Mais, de nos jours, on peut la connaître par la théorie des perturbations.
Les approximations sur la fonctionnelle d’échange-corrélation sont améliorées si on
tient compte du spin électronique. D’où l’appellation de l’approximation de la densité
locale de spin (L.S.D.A, “local spin density approximation”). Il faut aussi noter qu’il
existe dans la littérature plusieurs paramétrisations de la fonctionnelle d’échange corré-
lation dans ce cadre de la L.S.D.A parmi lesquelles, on peut citer : Perdew and Zunger
[34], Vosko et al. [35]. Elle s’écrit :
Z
LSDA
Exc [ρ(→
−
r )] = d3 →
−
r εLSDA
xc [ρ ↑ (→
−
r ) + ρ ↓ (→
−
r )] (1.44)
avec
ρ(→
−
r ) = ρ ↑ (→
−
r ) + ρ ↓ (→
−
r)
où ρ ↑ (→
−
r ) et ρ ↓ (→
−
r ) représentent respectivement les densités d’électrons associées aux
états de spin up(↑) et down(↓)
1.2.2 Fonctionnelle semi-locale ou Approximation du Gradient Généralisé
(GGA)
La théorie sur la LSDA a donné des bons résultats et a eu un succès intéressant dans
la description des propriétés de l’état fondamental d’un grand nombre de matériaux et
convient même à l’étude des systèmes isotropes. Mais elle a aussi ses limites lorsque la
densité ρ(→
−r ) n’est pas uniforme. Pour ce, on fait appel à l’approximation du Gradient
Généralisé (GGA, “Generalised Gradient Approximation”) donnée par :
Z
GGA
Exc [ρ(→
−
r )] = d3 →
− GGA
r fxc [ρ(→
−
r ), 5ρ(→
−
r )] (1.45)
GGA
où fxc est une fonction dépendante de la densité locale et du gradient de la densité.
1.2.3 Les pseudo-potentiels
Les pseudo-potentiels représentent une partie non négligeable des méthodes de calcul
en physique de la matière et la plupart des propriétés physiques de la matière découle des
propriétés d’ électrons de valence ; par exemple dans l’étude des métaux de transition. De
plus les pseudo-potentiels apportent un modèle simple des interactions entre l’ion formé
par le noyau et les électrons de coeur et les électrons de valence qui eux, participent
16
activement aux réactions chimiques [36]. Cette approximation permet de réduire un
grand nombre d’équations à résoudre.
1.2.4 Construction des pseudo-potentiels
Dans la construction des pseudo-potentiels, on en note deux méthodes : la méthode
de la conservation de la norme et la méthode de la non conservation de la norme.
- La méthode de la conservation de la norme : dans cette méthode le pseudo-potentiel est
basé sur les pseudo-fonctions d’ondes lisses et consiste en la détermination des énergies
et des états propres de l’atome isolé dans son état fondamental. Il doit vérifier un certain
nombre de propriétés par sa qualité, sa transférabilité c’est à dire pour des environne-
ments, chimiques différents on doit pouvoir utiliser le même potentiel et son efficacité.
- La méthode de la non conservation de la norme : contrairement à la première mé-
thode, les pseudo-fonctions d’ondes ne sont plus lisses à cause de la non conservation de
la norme.
1.3 Conclusion
Dans ce chapitre, nous avons passé en revue les points importants qui fixent le cadre
théorique des calculs ab initio de notre étude. Nous avons rappelé l’équation de Schrö-
dinger qui est l’équation de base du physico-chimiste et les méthodes qui résultent sa
résolution. Dans cette résolution, nous avons abordé plusieurs approximations jusqu’au
développement de la théorie de la fonctionnelle de la densité (DFT). Cette théorie est
basée sur deux théorèmes de Hohenberg et Kohn dépendant de l’énergie de la fonction
d’onde et de la densité de l’état fondamental. Ensuite, nous avons utilisé les équations
de Kohn-Sham développant cette théorie en tenant compte de l’approche de Thomas-
Fermi. Nous avons aussi abordé la description de l’énergie d’échange et corrélation avec
quelques approximations comme l’Approximation du Gradient Généralisé (G.G.A) et
enfin sur les pseudo-potentiels qui permettent de décrire les interactions ions-électrons.
17
Chapitre 2
Le Code Dacapo ayant été utilisé dans plusieurs travaux précédents [23, 24], nous
allons ici résumer son architecture.
18
électrons de valence. Cet ion est considéré par DACAPO comme gelé et représenté par
un pseudo-potentiel. Parmi les pseudo-potentiels, on distingue entre autres : les pseudo-
potentiels de Trouillier et Martin [37], Hamann [38] et Vanderbilt [39]. DACAPO est
un code qui utilise cependant le pseudo-potentiel ultra-doux (ultrasoft) de Vanderbilt
[39]. Ce pseudo-potentiel utilisé par DACAPO fait une description nette des paramètres
énergétiques et géométriques du système étudié.
2.1.2 La base d’ondes planes
Kohn et Sham proposèrent en 1965 une équation qui ramène la détermination de
la fonctionnelle d’énergie d’un système d’équations autocohérentes à un système mo-
noélectronique (cf. Chapitre 1). Ce système monoélectronique est représenté par un
hamiltonien hi qui satisfait à l’équation aux valeurs propres :
hi ϕi = εi ϕi (2.1)
Les équations de Kohn et Sham bien que exactes, leur résolution pose problème.
Pour surmonter cette difficulté, il est préférable que l’étude se fasse dans les systèmes
périodiques qui respectent le théorème de Bloch [40]. Ce théorème de Bloch permet de
représenter la fonction d’onde sous la forme :
−
→− →
ϕnb,k (→
−
r ) = ei k . r Φnb,k (→
−
r) (2.2)
avec :
ϕnb : fonction d’onde du système périodique,
→
−
k : vecteur de l’espace réciproque du cristal,
Φnb : fonction périodique de même périodicité que le système étudié, associée à une
bande nb .
L’application de ce théorème de Bloch suppose que le calcul se fasse dans une cellule
qui, elle même est une partie du réseau de l’espace réel étudié. Le réseau réciproque est
périodique et sa cellule de base est la première zone de Brillouin. Chaque fonction peut
s’écrire d’une façon explicite sous forme d’une série de Fourier. Le théorème décompose
alors la fonction Φnb en ondes planes de la forme :
19
−
→− →
Φnb,k (→
− Cnb,k (→
−
X
r)= g )ei g . r (2.3)
−
→
g
avec :
nb = 1, ..., L
→
−g : vecteur de translation du réseau réciproque tel que →
−
g .→
−
r = 2mπ, avec m un entier
→
−
k : vecteur de la zone de Brillouin
→
−r : position des noyaux
En portant (2.3) dans (2.2), on a l’expression :
→ −
− → − →
ϕnb,k (→
− Cnb,k (→
−
X
r)= g )ei( g + k ). r (2.4)
−
→
g
→
− − 2
|k +→
g | < Ecut (2.5)
Des approximations ont été apportées pour surmonter cette difficulté telles que :
- tout l’espace doit être représenté de façon uniforme à l’aide des ondes planes,
- les ondes planes doivent être transformées en fonction d’ondes planes à l’aide de la
transformée de Fourier rapide (FFT). Ce qui leur confère une grande efficacité d’utilisa-
tion,
- en se basant sur le nombre d’ondes planes, la précision et la convergence des calculs
sur les propriétés physiques sont contrôlées,
- les ondes planes sont indépendantes des positions atomiques ; ce qui permet de calculer
facilement les forces F agissant sur les atomes par le théorème de Hellmann-Feynman
[41] où les forces sont données par l’expression :
20
∂H
F =− ϕ ϕ (2.6)
∂R
où
- H : hamiltonien du système ;
- R : position atomique ;
- ϕ : fonction d’onde.
2.1.3 Génération des points-k
Les fonctions d’ondes trouvées dans l’espace réciproque (cf. éq. 2.4) sont définies
pour un nombre infini de points-k (k-points de la zone de Brillouin) [25]. Chaque point-
k contribue au potentiel électronique dont on déduit l’énergie totale. Afin de déterminer
numériquement les intégrales portant sur les états occupés de chacun des points-k, il
convient de restreindre le calcul à une grille de points-k de taille finie. Pour générer cette
grille de points-k, la méthode de Monkhorst et Pack [42, 43] est utilisée.
Nous résumons ici le principe de détermination des points-k tel que l’auteur de la
référence [23] a rapporté. Pour déterminer la densité des points-k nécessaire à la conver-
gence de l’énergie totale, cet auteur a utilisé une maille (1x1) de 4 plans et 5 plans
de Au(001) et relaxé le plan de surface en faisant varier le nombre de points-k. L’aug-
mentation de la densité des points-k augmente le temps de calcul dans la recherche de
l’autocohérence mais une bonne précision sur la valeur de l’énergie totale demande une
densité de points-k relativement importante. Les tests de convergence sur le nombre de
points-k ont montré que la variation d’énergie n’est que de l’ordre de 0.009 eV. N’dollo
[23] a trouvé une convergence raisonnable de l’énergie pour un nombre de points-k de
(12x12x1) pour une maille (1x1). Cette convergence de l’énergie a été ensuite vérifiée
pour une maille (2x2) avec un nombre de points-k de (6x6x1). Ces densités de points-k
relatives aux mailles (1x1) et (2x2) ont été conservées dans tous les calculs présentés
dans cette thèse.
2.1.4 Distribution de Fermi
D’après la prévision de la théorie de la mécanique quantique, les électrons se placent
par bande. Pour pouvoir trouver le nombre des électrons dans les états occupés, le
code DACAPO utilise la méthode de Fermi-Dirac. Au voisinage du niveau de Fermi,
le code utilise la distribution de Fermi pour accélérer la convergence du calcul. Dans
nos calculs, on considère kB T = 0.1 eV (kB , constante de Boltzmann). Les résultats
21
sont ensuite extrapolés vers 0 K (Température correspondant à l’état fondamental du
système étudié).
2.1.5 Energie d’échange-corrélation du Gradient Généralisé, GGA
DACAPO utilise la méthode de la théorie de la fonctionnelle de la densité dans la-
quelle la fonctionnelle énergie d’échange-corrélation est évaluée. Cette énergie d’échange-
corrélation n’a pas jusqu’à maintenant une formulation exacte. Dans nos simulations,
nous avons utilisé une fonctionnelle GGA proposée par Perdew et Wang (PW91) [44].
Cette approximation se justifie car dans plusieurs cas les résultats obtenus sont proches
de la réalité physique.
2.1.6 Energie de coupure
La base des ondes planes utilisée pour la représentation des fonctions d’onde est à
priori infinie. Pour les calculs numériques, il est indispensable de limiter cette base. Pour
DACAPO, cela est réalisé en fixant l’énergie maximale des ondes planes. Les tests de
convergence sur l’énergie de coupure (Ecut ) montrent que entre 400 eV et 600 eV, les
variations d’énergies sont de l’ordre de 0.03 eV [23]. Dans leurs travaux, Baraldi et al.
[45] montrent qu’une valeur de 340 eV est suffisante pour des pseudopotentiels “ultra-
doux“ dans le cas des métaux de transitions. Dans le cas de notre travail, elle a été fixée
à 400 eV pour faciliter la rapidité et la convergence des calculs.
2.1.7 Calcul périodique, super-cellule ou modèle du slab
Le code DACAPO utilise le modèle des fonctions d’ondes dont le système étudié
respecte le théorème de Bloch donc, périodique. Il étudie le système réel en utilisant la
technique de la super cellule. Le système étudié par le code DACAPO forme une couche
mince d’atomes superposés les uns sur les autres ; périodique suivant les directions (OX)
et (OY). C’est la méthode du slab. Dans ce modèle, on introduit dans la direction (OZ)
un vide entre deux slabs avec une distance correcte de telle sorte que les slabs (couches)
se reproduisent périodiquement.
Dans le cas de notre solide (Au), ce modèle est représenté sur la figure 2.1. Dans le
cas de notre travail, le vide a été fixé à 12 Å de telle sorte que les interactions entre slabs
soient négligeables.
22
Figure 2.1: Schéma du modèle de slab
23
Figure 2.2: Visualisation de l’algorithme de convergence du code DACAPO
24
X
ρ(ε) = hϕn |ϕn i δ(ε − εn ) (2.7)
n
X X
ρ(ε) = ρi (ε); ρi (ε) = hϕn |ii hi|ϕn i δ(ε − εn ) (2.8)
i n
Z
d→
−
r ρ(→
−
r , ε); ρ(→
− hϕn |→
−
r i h→
−
X
ρ(ε) = r , ε) = r |ϕn i δ(ε − εn ) (2.9)
n
Notons que :
Z
d→
−
r |→
−
r i h→
−
X
|ii hi| = 1 et r|=1 (2.10)
i
Par définition, on appelle densité d’états locale (LDOS, local density of states), la
quantité ρi (ε) :
X
ρi (ε) = hϕn |ii hi|ϕn i δ(ε − εn ) (2.11)
n
La densité d’états projetée (PDOS, projected density of states) est définie par :
ρ(→
− hϕn |→
−
r i h→
−
X
r , ε) = r |ϕn i δ(ε − εn ) (2.12)
n
Enfin, notons qu’en intégrant la PDOS multipliée par la distribution de Fermi nF (ε),
on retrouve la densité électronique :
25
Z
dε nF (ε)ρ(→
−
r , ε) = n(→
−
r) (2.13)
où φsyst. est le travail de sortie du système (adsorbat-substrat) et φsurf ace celui corres-
pondant de la surface nue du substrat.
26
2.5.1.1 Paramètre de maille de Au
La connaissance du paramètre de maille de notre substrat est donc primordiale. Nous
avons fait un test de calcul sur le paramètre de maille a de Au en volume pour une densité
de point-k = (14, 14, 14) et en le faisant varier autour de la valeur expérimentale (4.08
Å) [28]. Nous avons traçé la variation de l’énergie totale en fonction de a comme le
montre la figure 2.3.
27
Tableau 2.1: Illustration du paramètre de maille de Au en volume.
1
Esurf = [Eslab − NAu Ebulk ] (2.15)
2A
28
Tableau 2.2: Tableau comparatif sur l’énergie par unité de surface.
L’analyse sur l’énergie de surface montre que pour la maille (1x1) et (2x2), cette
énergie de surface est presque la même. Ceci révèle que la variation de l’énergie de
surface est négligeable entre la maille (1x1) et (2x2). Notre valeur de ∼ 0.88 J/m2
peut être comparée aux valeurs expérimentale et théorique. Elle est presque la moitié
de la valeur expérimentale qui est de 1.540 J/m2 [53]. Théoriquement, Zólyomi et al.
[54] trouvent une valeur de 0.864 J/m2 en utilisant la DFT-PBE-GGA ; cette valeur est
légèrement inférieure à la notre. Wang and Wang [55] trouvent dans leur calcul sur la
DFT-LDA une valeur de 1.359 J/m2 . Vitos et al. [56] par contre trouvent 1.627 J/m2 avec
la fonctionnelle FCD-LMTO (potential linear muffin-tin orbital). Les valeurs trouvées
par Wang and Wang et Vitos et al. sont nettement plus grandes à la valeur de notre
maille. Cette grande différence peut se justifier sur les méthodes utilisées.
L’énergie de cohésion est déterminée par la relation suivante :
Ebulk
Ecoh = Eat − (2.16)
N
où Eat est l’énergie d’un atome de Au isolé, calculée dans une boîte cubique d’arête
15 Å et Ebulk l’énergie totale du cristal (massif) élémentaire contenant N atomes. Avec
cette définition, une valeur positive de l’énergie de cohésion indique une interaction
interatomique.
Nous avons calculé une énergie de cohésion de 4.65 eV. La valeur expérimentale est
de 3.81 eV [28] ; elle est 0.84 eV inférieure à la nôtre. Gómez et al. [49], et Janthon et
al. [46] calculent des valeurs respectives de 3.07 eV et 2.99 eV, lesquelles ne sont pas
en accord quantitatif avec notre détermination. Le désaccord peut s’expliquer par les
différences entre les méthodes de calculs (corrections relativistes dans [49] vs calculs «
purs » dans notre étude).
Nous avons aussi calculé le travail de sortie (work-function) de la surface nue de l’Or
pour une maille (1x1) et (2x2) pour un slab de 4 plans. Nous trouvons respectivement
29
5.098 eV pour une maille (1x1) et 5.080 eV pour une maille (2x2). A 5 plans, nous
trouvons aussi respectivement 5.092 eV et 5.074 eV pour la maille (1x1) et (2x2). Notre
valeur de 5.080 eV se rapproche de celle de l’expérience [57] qui est de 5.220 eV. Notre
valeur est en très bon accord avec celle trouvée (5.080 eV) par N’dollo [23] qui utilise
une méthode similaire. Elle est aussi cependant en bon accord avec celle trouvée par
d’autres auteurs [48] qui est de 5.140 eV.
2.5.3 Structure électronique
Dans la suite de notre travail, nous allons parler de la structure électronique et donc
de la densité des états projétée (PDOS) qui traduira les interactions de la surface de
l’Or (substract) sur sa face (001) avec les différents fragments (adsorbats).
2.5.4 Résultats post calculs préliminaires
Ces résultats porteront sur l’adsorption des différents fragments (CH3 , OH, CH2 OH,
H2 CO, HCO, CO et C) sur la face (001) de Au. La stabilité relative des sites d’adsorp-
tion sur chacun des fragments sera déterminée ainsi que les énergies d’adsorption, les
structures géométriques, électroniques et le travail de sortie du système. Comme nous
l’avons déjà évoqué au sous-paragraphe 2.5.1, l’énergie en phase gazeuse de chacun des
fragments a été déterminée à partir des calculs réalisés en les plaçant seul dans une boîte
cubique d’arête de 15 Å ou 20 Å au point gamma. L’adsorption de chaque adsorbat sera
faite sur trois sites de la surface Au(001) comme le montre la figure 2.4.
Figure 2.4: Schéma de la surface libre de Au(001) pour une maille (1x1). En bleu, les atomes de Au
en plan de surface et en jaune, l’atome du 2nd plan de surface.
30
Nous avons les sites suivants :
- le site apical (top, en anglais), le radical est placé au dessus d’un atome de Au ;
- le site ponté (bridge, en anglais), le radical est placé entre deux atomes plus proches
voisins de Au ;
- le site de symétrie quatre (hollow, en anglais), le radical est placé entre quatre
atomes plus proches voisins de Au.
2.6 Conclusion
Dans cette partie, nous avons décrit le fonctionnement du code DACAPO basé sur
la théorie de la fonctionnelle de la densité. Ce code permet d’extraire toute l’énergie liée
au phénomène d’adsorption. Par ce code, nous avons obtenu les propriétés du métal à
savoir : les énergies de surface, de cohésion, le paramètre de maille de Au. Nous avons
enfin présenté les différents sites d’adsorption.
31
Chapitre 3
32
ces études est la réactivité des surfaces permettant de situer l’importance de l’activité
catalytique.
3.1.2 Géométrie du méthyle
Dans notre étude, le radical méthyle est un fragment qui est issu des molécules
hydrocarbonées de la biomasse c’est-à-dire de la rupture d’une liaison C-OH dans la
molécule du méthanol (CH3 OH). Le méthyle est composé d’un atome de carbone lié à
trois atomes d’hydrogène. La géométrie du méthyle (CH3 ) libre dans le plan en phase
gazeuse est donnée sur la figure 3.1.
L’optimisation du calcul sans spin polarisé (ssp) et avec spin polarisé (sp), nous a
permis de remarquer que le CH3 libre est magnétique avec un moment magnétique total
de ∼ 1.00 µB . Nous ne disposons pas de données théorique et expérimentale pour des
éventuelles comparaisons. Cette magnétisation n’affecte pas les paramètres géométriques.
Les longueurs de liaisons C-H et les angles ne sont réduits que de 0.001 Å et 0.03o
respectivement par rapport au calcul sans spin polarisé. Les valeurs de longueurs de
liaisons et angles trouvées sont représentées dans le tableau 3.1.
33
Tableau 3.1: Comparaison entre paramètres structuraux calculés et expérimentaux de CH3 libre. Les
nombres entre parenthèses représentent les nombres de longueurs de liaison et angles possibles.
En se référant au tableau 3.1, nous remarquons que nos valeurs sont en bon accord
avec les résultats expérimentaux [78]. Elles le sont aussi avec les résultats théoriques
trouvés par Moussounda et al. [76], Michaelides et al. [77] et Lai et al. [69].
3.1.3 Site préférentiel d’adsorption
Dans nos calculs, nous avons utilisé les taux de couverture de 0.25, 0.50, 0.75 et
1.00 MC. Les taux de couverture de 0.25, 0.50 et 0.75 MC sont obtenus en adsorbant,
respectivement, un, deux ou trois radicaux de CH3 sur la surface (2x2) de Au(001).
Signalons que pour le taux de couverture 1.00 MC, nous avons utilisé une maille (1x1)
dans laquelle nous avons fait adsorber un radical CH3 .
Dans tous les calculs, l’énergie de coupure a été fixée à 400 eV. De plus nous avons
utilisé un slab formé de 4 et 5 plans de Au. Nous avons ensuite après adsorption, déter-
miné les paramètres géométriques, les énergies d’adsorption du système CH3 /Au(001)
avec et sans spin polarisé, la work function et le moment dipolaire.
Notons qu’après adsorption, le moment magnétique total de CH3 est approximati-
vement réduit à zéro (∼ +25.10−5 µB ) c’est-à-dire que la magnétisation de CH3 est
considérablement réduite. Nous remarquons que le CH3 magnétique et la surface Au ne
sont pas couplés magnétiquement parce que, les interactions entre CH3 et la surface de
Au annihilent complètement le moment magnétique de spin.
Enfin, le radical méthyle (CH3 ) a été adsorbé sur un seul site de la surface Au(001).
Dans la phase d’optimisation du système CH3 /Au(001), nous avons constaté que les
sites ponté (bridge) et de symétrie 4 (hollow) sont instables. Le radical méthyle (CH3 )
dans les deux cas migre vers le site apical. Ceci n’est pas nouveau car Moussounda et
al. [76], par un calcul DFT-GGA, avaient trouvé aussi que les sites ponté et de symétrie
4 étaient instables lors de l’adsorption de CH3 sur le Pt(100).
34
3.1.4 Energies d’adsorption (avec et sans spin polarisé), travail de sortie
(work function), moment dipolaire et paramètres géométriques
Pour le site apical, les énergies (avec et sans spin polarisé), le travail de sortie (work
function), le moment dipolaire et les paramètres géométriques calculés en fonction du
nombre de plans du slab et des taux de couverture après adsorption sont donnés dans
le tableau 3.2.
Tableau 3.2: Energies d’adsorption (Eads (sp) et Eads (ssp)), travail de sortie (work-function) (∆Φ),
moment dipolaire (µ), longueurs de liaisons entre C et H (dC−H ), entre C et Au (dC−Au ), et angles
(H-C-H) obtenus après adsorption de CH3 sur le site apical de Au(001) avec un slab de 4 plans. Les
résultats entre crochets sont calculés avec un slab de 5 plans.
1
Eads = ECH3 /slab − Eslab − N ECH3 (3.1)
N
où ECH3 /slab est l’énergie totale du slab Au(001) après adsorption de N méthyle sur une
seule face, Eslab l’énergie totale de Au(001) nue et ECH3 l’énergie de CH3 libre (avec
ou sans spin polarisé). Notons qu’une grande valeur négative de Eads correspond à une
grande stabilité du système.
Les énergies d’adsorption obtenues avec et sans spin sur le site apical pour un slab à 4
et 5 plans sont données dans le tableau 3.2. Nous remarquons que le spin polarisé a affecté
considérablement les énergies d’adsorption (par rapport à celles sans spin) et celles-ci
sont réduites de ∼ 0.614 eV. Nos résultats sur les énergies d’adsorption avec l’effet de
35
polarisation confirment ceux trouvés par Wang et al. [79] sur le système CH3 O/Ni(111).
Le tableau 3.2, ci-dessus, montre que les énergies d’adsorption changent avec la varia-
tion du taux de couverture. Ce faisant, les énergies d’adsorption décroissent en général
quand le taux de couverture augmente. Notons que la stabilité de CH3 est plus favorable
à 5 plans de Au(001) sur 0.25 MC et à 4 plans sur 0.50 MC. Cependant, les énergies d’ad-
sorption augmentent lorsque le nombre de plans augmente c’est-à-dire, lorsqu’il passe
de 4 à 5 plans pour les taux de couvertures 0.25, 0.50 et 0.75 MC. Par contre pour
le taux de couverture de 1.00 M.C, l’énergie d’adsorption décroit lorsque le nombre de
plans augmente. Nous pouvons faire remarquer que la variation de l’énergie d’adsorption
avec le nombre de plans est dans le rapport (Eads )5plans /(Eads )4plans = 1.27, 1.14, 1.12 et
0.96 pour 0.25 MC, 0.50 MC, 0.75 MC et 1.00 MC respectivement pour le calcul sans
spin ; les valeurs correspondantes pour les calculs avec spin sont de 1.40, 1.21, 1.20 et
0.85 respectivement. Nous pouvons dire que, à l’exception du taux de couverture de 1.00
MC, lorsque le nombre d’atomes du slab augmente, l’interaction entre l’adsorbat (CH3 )
et le slab augmente.
L’évolution des énergies d’adsorption avec spin (a) et sans spin (b) en fonction du
taux de couverture pour les slabs de 4 (en rouge) et 5 (en noir) plans est représentée sur
la figure 3.2.
Figure 3.2: Courbe représentative des énergies d’adsorption en fonction du taux de couverture pour
4 (en rouge) et 5 (en noir) plans.
Comme on peut le voir sur la figure 3.2, pour une couverture donnée et quelque soit
le type de calcul examiné ici, l’énergie d’adsorption à 5 plans est supérieure à celle de
36
4 plans. La seule exception existe pour une couverture de 1.00 MC pour laquelle nous
constatons une quasi dégénérescence en énergie d’adsorption. Le haut taux de couverture
se caractérise ainsi par des répulsions stériques entre adsorbats.
Poursuivons notre discussion par une étude comparative de nos énergies d’adsorption
avec d’autres travaux antérieurs. Signalons tout de suite, qu’il n’existe pas dans la lit-
térature des travaux expérimentaux et théoriques sur l’adsorption de CH3 sur Au(001).
Cependant, plusieurs études ont été faites sur d’autres métaux comme le Ni(100) [69],
Pt(100) [76] et Au(111) [19]. En titre d’illustration, notre valeur d’énergie d’adsorption
de -1.903 eV à 4 plans et à 0.25 MC peut être comparée avec celles obtenues dans les
études antérieures de l’adsorption du CH3 sur les métaux de transition. Moussounda [24]
étudie l’adsorption de CH3 sur un slab de 4 plans de Pt(100) pour un taux de couverture
de 0.25 MC et il trouve une valeur de -2.60 eV qui est de 36.6 % plus grande que la
nôtre sur le même site (apical). Lai et al. [69] et wang et al. [107] trouvent respective-
ment une énergie de -2.05 eV et -1.62 eV avec un slab de 5 plans de l’adsorption de CH3
sur Ni(100) et sur Ni(111). Ces valeurs sont respectivement 15 % et 33 % plus faible
que notre présent résultat de -2.423 eV à 5 plans. Santiago-Rodriguez sur Au(111) [face
plus dense que Au(001)] [19] trouve une valeur de -1.19 eV. Les désaccords constatés ici
peuvent s’expliquer par les natures et les faces cristallographiques différentes entre les
surfaces étudiées [Au(001) vs Ni(100), Ni(111), Au(111)].
3.1.4.2 Travail de sortie (work function) et moment dipolaire
Le travail de sortie (work-function) est défini comme l’énergie minimale qu’il faut
fournir pour extraire ou transférer un électron d’une surface à une autre. Pour com-
prendre ce phénomène dans cette réaction catalytique, nous avons réalisé une étude sur
la variation de la work-function ∆Φ mais aussi du moment dipolaire µ. Cette variation
est ainsi marquée par la différence de la work-function entre le système (substrat +
adsorbat) et la surface nue de Au (cf. Eq. 2.14 du chapitre 2). La détermination de ∆Φ
conduit ainsi au moment dipolaire défini par la relation de Helmholtz [80] :
A∆Φ
µ= (3.2)
12πθ
37
function de Au diminue lorsque le méthyle est adsorbé sur la surface. Ce qui montre que
l’électron quitte plus facilement la surface adsorbée. En conséquence, une valeur négative
de ∆Φ est trouvée, c’est-à-dire une réduction de la work-function. En observant le tableau
3.2, nous remarquons que la variation de la work-function (∆Φ) n’est pas proportionnelle
au taux de couverture θ à l’exception de θ = 0.50 M.C pour lequel ∆Φ ∼ -0.50 eV. En
outre, nous notons une décroissance de la variation de la work-function lorsque le taux de
couverture augmente lors de l’adsorption de CH3 , c’est-à-dire pour θ = 0.25, 0.75 et 1.00
M.C ; cependant, il y a une exception à θ = 0.50 M.C et à 5 plans où ∆Φ est constant
pour les deux slabs de 4 et 5 plans. Ce genre de situation, c’est-à-dire la diminution de
la work-function est possible pour les surfaces des métaux mais pas toujours évident.
C’est pourquoi la question sur la relation de la work-function et du transfert de charge
vers le substrat apparait. Mais, cette question a été analysée par Pettersson et Bagus
sur l’étude du système Cl/Cu(001) [81], Bagus et Illas sur O/Cu(001) [82], Migani et al.
[83] sur Cl/MT(001) (avec MT = Cu, Ag, Au, Rh, Pd, Pt), Pick sur CH3 O/Cu [84], Lin
et Groβ sur H2 O/Au [85]. Sur cette question, ils concluent que selon la nature atomique
des surfaces, la présence des charges sur ces métaux donne une réponse en terme de
polarisation. En conséquence, il est plus correct de parler de polarisation des surfaces
que de large transfert de charge. Nous notons que la valeur inattendue de ∆Φ ∼ -0.50 eV
(pour θ = 0.50 M.C) en utilisant un slab de 4 et 5 plans de Au proviendrait de la même
polarisation provoquée par l’adsorption. Enfin, nous avons remarqué que l’augmentation
du nombre de plans influence faiblement ∆Φ.
Pour terminer ce sous-paragraphe, notons que le signe du moment dipolaire µ est
négatif, ce qui signifie que le vecteur moment dipolaire pointe vers le vide constituant
la super cellule. En plus, le moment dipolaire µ augmente avec l’augmentation du taux
de couverture. Ceci est une conséquence direct de l’adsorption de CH3 qui provoque une
redistribution de la densité électronique.
3.1.4.3 Paramètres géométriques
Après la discussion sur les énergies d’adsorption du système CH3 /Au(001) avec et
sans spin polarisé pour les mailles (2 x 2) et (1 x 1), nous avons déterminé les para-
mètres géométriques. Les résultats obtenus révèlent que le spin polarisé n’affecte pas les
paramètres géométriques (longueurs de liaison et angles). Ces résultats sont rassemblés
dans le tableau 3.2.
Nous remarquons aussi que l’augmentation du nombre de plans dans le slab et du taux
de couverture n’influence pas les paramètres géométriques. Une exception est cependant
38
observée à θ = 1.00 M.C avec une légère modification de l’angle H-C-H qui diminue de
110.4o à 108.4o lorsque le nombre de plans augmente.
A θ = 0.25 M.C, nos longueurs de liaison dC−Au = 2.123 [2.121] Å sont en bon accord
avec celles de Moussounda et al. [76] et de Santiago-Rodriguez et al. [19] qui trouvent des
valeurs respectives de 2.08 Å et 2.13 Å dans les adsorptions respectives de CH3 /Pt(100)
et CH3 /Au(111). Lai et al. [69] dans leur étude du système CH3 /Ni(100) trouvent une
valeur de 1.95 Å.
Par ailleurs, nous constatons que les longueurs de liaisons dC−H , ne sont pas trop
affectées par l’adsorption par comparaison avec la phase libre de CH3 . Cependant les
angles H-C-H de CH3 sont réduits de 120o (pour le radical CH3 libre) à 111.6o (pour
le système CH3 /Au(001) à 0.25 MC et, entre 111.1o et 108.4o pour les autres taux de
couverture. Ainsi, nous pouvons conclure que l’interaction entre CH3 et la surface de
Au(001) fait que les hydrogènes s’écartent de la surface. Ceci montre que CH3 passe
d’un état d’hybridation sp2 à sp3 lorsqu’il est adsorbé.
3.1.5 Propriétés électroniques de CH3 et de la surface Au(001)
Pour mieux élucider la nature de l’interaction entre le radical méthyle (CH3 ) et la
surface Au(001), nous allons analyser l’évolution de la structure électronique à travers
la densité d’états projetée (PDOS). Cette étude de la PDOS est faite sans et avec spin
polarisé.
3.1.5.1 PDOS de CH3 libre et adsorbé à 0.25, 0.50 et 0.75 MC sans et avec
spin polarisé
3.1.5.1.a Densité d’états projetée (PDOS) de CH3 libre sans spin polarisé
Sur la figure 3.3, nous donnons la PDOS des orbitales moléculaires de valence de CH3
libre et adsorbé sur Au(001).La PDOS est déterminée avec les taux de couvertures 0.25
M.C, 0.50 M.C et 0.75 M.C. Sur cette figure, la PDOS de CH3 libre présente trois pics
situés à -7.53 eV, -1.27 eV et +4.39 eV équivaut respectivement aux orbitales (2a1 ), (1e)
et (3a1 ), selon la notation de Schönflies. Le premier pic (2a1 ) est formé par la combinaison
de l’orbitale 2s du carbone et des orbitales 1s des trois atomes hydrogènes. Le deuxième
pic (1e) est formé d’une paire d’orbitales de liaison C-H ou chacune des orbitales est
constituée d’une orbitale p du carbone et d’une orbitale s de l’hydrogène. Le dernier pic
(3a1 ) est la combinaison des orbitales 2p du carbone. Nous pouvons voir que les pics
de 2a1 et 1e sont des états moléculaires occupés alors que l’état 3a1 est inoccupé. La
PDOS de CH3 libre est dominée par l’orbitale 1e. L’orbitale moléculaire 3a1 est l’orbitale
la plus fortement occupée. Nous pouvons voir que sous adsorption, l’augmentation du
39
taux de couverture a une influence notable sur les PDOS. En effet, lorsque le taux de
couverture augmente, les intensités des pics diminuent à l’exception du pic situé à -
7.53 eV pour 0.25 M.C. Sous adsorption, avec tous les taux de couverture examinés
ici, l’orbitale moléculaire 3a1 disparait car elle est transformée en continuum sur
l’intervalle [+2.73 eV, +6.00 eV] par rapport à celle sur CH3 libre. Ceci s’explique par
un fort mélange entre l’orbitale 3a1 de CH3 et les bandes s, p et d de la surface Au
comme nous allons le voir au sous-paragraphe 3.1.5.2.b. Toutes ces observations sont
qualitativement en accord avec celles obtenues par Michaelides et Hu [77] sur le système
CH3 /Ni(111) et Moussounda [24] sur CH3 /Pt(100).
Figure 3.3: PDOS sur les orbitales du groupe CH3 libre (en pointillé) par rapport à celles dans
CH3 /Au(001) adsorbées (en trait continu ) pour 0.25, 0.50 et 0.75 M.C sur le site apical. Le niveau de
Fermi est à 0 eV.
40
3.1.5.1.b Densité d’états projetée (PDOS) de CH3 libre avec spin polarisé
La PDOS de CH3 avec spin polarisé est représentée sur la figure 3.4. En observant
la phase gazeuse, nous remarquons que les états de spin up et de spin down présentent
trois pics. Ces pics sont localisés à -8.78 eV, -2.11 eV et +2.53 eV pour le spin up et,
-7.77 eV, -1.91 eV et +5.77 eV pour le spin down. Après adsorption, les états de spins
up et down décroissent quand le taux de couverture augmente, mais nous notons une
exception à 0.25 M.C pour les états localisés à -8.78 eV et -7.77 eV. Nous remarquons
également que sous adsorption, les états de spin up se décalent vers les hautes énergies
mais la situation n’est pas la même sur les états de spin down par rapport aux pics
situés à -7.77 eV et -1.91 eV. Lorsque le taux de couverture augmente (états de spin
down), les pics se décalent vers les basses énergies par rapport à celui situé à -7.77 eV
(phase gazeuse) et vers les hautes énergies par rapport à -1.91 eV (phase gazeuse). Cette
inversion de situation sur les deux pics (-7.77 eV et -1.91 eV) peut être expliquée par
la magnétisation qui disparait lors de l’adsorption de CH3 sur la surface de Au(001).
En conséquence, sous adsorption les pics des états de spins up et down sont localisés
à ∼ -8.17 eV et -1.51 eV. Notons un autre point important aussi sur cette PDOS : les
deux états inoccupés positionnés à +2.53 eV (pour les états de spin up) et +5.77 eV
(pour les états de spin down) au-delà du niveau de Fermi sont très affectés par la surface
c’est-à-dire sont transformés en continuum. Ceci peut être expliqué par une hybridation
entre les états de CH3 et les bandes de Au.
Cette étude est pionnière car nous n’avons pas trouvé dans la littérature de travaux
expérimentaux ou théoriques pour d’éventuelles comparaisons.
41
Figure 3.4: PDOS sur les orbitales du groupe CH3 libre (en pointillé) par rapport à celles dans
CH3 /Au(001) pour 0.25, 0.50 et 0.75 M.C sur le site apical. Le niveau de Fermi est à 0 eV.
Sur la figure 3.5, nous comparons les PDOS des bandes s, p et d de la surface de Au
libre par rapport à celles adsorbées pour les taux de couvertures de 0.25, 0.50 et 0.75
M.C. Sous adsorption du méthyle sur Au, nous remarquons l’apparition d’un nouveau
pic de résonance totalement absent sur la surface libre de Au, localisé à -7.53 eV au
niveau de la PDOS des bandes s, p et d de Au. Ce pic de résonance est exactement
situé sur la même position que l’orbitale 2a1 sur la PDOS de CH3 (Figure 3.3) ; elle
décroit avec l’augmentation du taux de couverture. Ceci résulte d’un fort mélange entre
les bandes s, p et d de Au et l’orbitales 2a1 de CH3 .
La figure 3.5a représente la PDOS de la bande s de Au. En observons la PDOS de
cette figure, nous remarquons que le pic le plus dominant est situé à -1.07 eV et ne
42
change presque pas avec l’augmentation du taux de couverture à l’exception de celui sur
0.25 MC qui augmente. Enfin, sur les deux intervalles [-0.26 eV, +3.37 eV] et [+4.38 eV,
+5.39 eV], les états des pics augmentent et décroissent respectivement.
Sur la figure 3.5b représentant la PDOS de la bande p, nous remarquons dans l’inter-
valle d’énergie [-2.48 eV, +6.48 eV], la présence de plusieurs pics de résonance avec une
augmentation significative quand le taux de couverture diminue par rapport à la PDOS
de la surface libre. Toutes ces transformations sont interprétées par un fort mélange
entre les orbitales moléculaires 1e et 3a1 de CH3 et la bande p de Au.
Enfin, la figure 3.5c représente la PDOS de la bande d de Au. Nous remarquons
qu’après adsorption, la PDOS décroit avec l’augmentation du taux de couverture par
rapport à la surface libre de Au. Nous remarquons aussi que le centre de gravité de la
bande d se décale de 0.30 eV vers les basses énergies par rapport à surface libre. Nous
notons deux domaines importants d’énergie. Le premier intervalle [-1.26 eV, +1.15 eV]
nous révèle la présence de trois nouveaux états localisés à -1.07 eV, -0.26 eV et +0.75
eV respectivement. Le second intervalle [+3.17 eV, +6.20 eV] montre une décroissance
des états inoccupés avec le taux de couverture. Ces changements sont dus à l’interaction
entre les orbitales moléculaires 1e et 3a1 de CH3 et la bande d de Au.
43
Figure 3.5: PDOS des bandes s, p et d de la surface nue de Au (en pointillés) par rapport à celles
dans CH3 /Au(001) sur 0.25, 0.50 et 0.75 MC pour le site apical. Le niveau de Fermi est 0 à eV.
44
d décroissent lorsque le taux de couverture augmente.En effet dans le domaine d’énergie
[-3.12 eV, +5.97 eV], la PDOS de la bande p de Au augmente ; ceci est une conséquence
de l’hybridation avec les orbitales de CH3 . Concernant les bandes s et d de Au, nous
remarquons qu’elles sont faiblement affectées par l’adsorption.
Figure 3.6: PDOS des bandes s, p et d de la surface nue de Au (en pointillés) par rapport à celles
dans CH3 /Au(001) sur 0.25, 0.50 et 0.75 MC pour le site apical. Le niveau de Fermi est à 0 eV.
45
tion), la Spectroscopie des pertes d’énergie à électrons (EELS, electron energy Loss
spectroscopy), la Distribution Angulaire d’ion de Désorption Stimulée par électron (ES-
DIAD, Electron-Stimulation Desorption Ion Angular Distribution), par spectroscopie
photo-électronique des rayons UV (UPS) et par Spectroscopie Métastable de Désexcita-
tion (MSD, Metastable De-excitation Spectroscopy).
Sur le plan théorique, plusieurs études ont été réalisées sur l’adsorption de OH sur les
surfaces de Pt(111) [91], Ni(100) [92], Ni(111) [93], Cu(111) [94], Ag(111) [95], Ag(110)
et Ag(100) [96]. Dans ses études, Anderson [91] montre que les interactions entre les
surfaces métalliques et OH sont dues à un mélange de leurs états, c’est-à-dire, mélange
entre les orbitales de OH et les bandes du métal.
3.2.2 Géométrie de l’hydroxyle (OH)
L’hydroxyle est un fragment qui peut être issu : soit de la décomposition de la mo-
lécule du méthanol (issue de la biomasse) en CH3 + OH, soit de la décomposition de
H2 O en H + OH. Il est donc composé d’un atome d’oxygène lié à un atome d’hydro-
gène. Sa géométrie en phase gazeuse est l’une des étapes dans l’étude de l’adsorption de
OH/Au(001). Elle est représentée sur la figure 3.7.
Figure 3.7: Visualisation du fragment OH libre. En rouge l’atome d’oxygène et en blanc l’atome
d’hydrogène.
Nous avons vérifié que les propriétés structurales de OH sont bien reproduites. Pour
cela, nous avons calculé l’énergie totale de OH en phase gazeuse en la plaçant dans
une boîte cubique de 15 Å. Le tableau 3.3 compare les longueurs de liaison, dOH , des
différentes études expérimentales et théoriques.
46
Tableau 3.3: Comparaison des longueurs de liaison (dO−H ) théorique et expérimentale de OH.
Nous remarquons que notre valeur est quantitativement en accord raisonnable avec
le résultat expérimental [78]. Elle l’est aussi avec les résultats théoriques trouvés par
d’autres chercheurs [79, 97, 98] comme le montre le tableau 3.3. Signalons que le spin
polarisé n’a pas affecté la valeur de la longueur de liaison lors du calcul.
3.2.3 Sites de préférence
Comme nous l’avons déjà évoqué plus haut sur l’adsorption de CH3 /Au(001), les
taux de couverture de 0.25, 0.50, 0.75 et 1.00 M.C sont aussi utilisés dans l’adsorption
de OH/Au(001). L’énergie de coupure de la base d’ondes planes a été aussi fixée à 400
eV. Un slab formé de 4 et 5 plans de Au a été utilisé.
Nous avons après adsorption, déterminé les paramètres géométriques et les énergies
d’adsorption du système OH/Au(001) sans et avec spin polarisé. Le radical hydroxyle
s’adsorbe sur les trois sites de la surface Au(001) à savoir : apical (top), ponté (bridge)
et de symétrie 4 (hollow). Cependant, nous avons remarqué que le site ponté est le plus
stable (site préféré) pour des taux de couverture 0.25, 0.50 et 0.75 MC.
Notons que sur la phase gazeuse avec spin polarisé, le moment magnétique de OH est
de (∼ 0.999 µB ). Après adsorption, il est approximativement réduit à zéro (de l’ordre
de ∼ 10−5 µB ) c’est-à-dire la magnétisation de OH est réduite considérablement par la
surface de Au. Nous remarquons aussi que le couplage entre OH magnétique et Au n’est
pas possible car, les interactions entre OH et la surface de Au annihilent complètement
l’action du spin. Signalons que nous n’avons pas eu des résultats expérimentaux et
théoriques dans la littérature pour d’éventuelles comparaison.
3.2.4 Energies d’adsorption (avec et sans spin polarisé), travail de sortie
(work function), moment dipolaire et les paramètres géométriques
Nous avons déterminé les énergies (avec et sans spin polarisé), le travail de sortie
(work function), le moment dipolaire et les paramètres géométriques pour des slabs de
4 et 5 plans et pour les taux de couverture de 0.25, 0.50, 0.75 et 1.00 MC.
47
3.2.4.1 Energies d’adsorption avec et sans spin polarisé
Les énergies d’adsorption (avec et sans spin polarisé) du système OH/Au(001) ont
été calculées en utilisant la formule (3.1) et les résultats sont donnés dans le tableau 3.4.
Nous notons une grande différence d’énergie d’adsorption entre les deux types de
calcul c’est à dire avec et sans la polarisation de spin. Les énergies d’adsorption avec
spin polarisé sont inférieures à celles sans spin polarisé ; et elles varient entre 0.390 et
0.791 eV.
Contrairement à l’adsorption de CH3 qui n’est stable que sur le seul site apical
(comme nous l’avons dit précédemment), l’adsorption de OH est stable sur les sites
apical, ponté et de symétrie 4. Sur ces trois sites d’adsorption, nous remarquons que les
énergies d’adsorption augmentent lorsque le nombre de plans de slab augmente (excepté
à 1.00 MC sur les sites apical et ponté) (cf Tableau 3.4).
48
Tableau 3.4: Energies d’adsorption (Eads (sp) et Eads (ssp)) obtenues pour le système OH/Au(001)
avec un slab de 4 et 5 plans. Les résultats avec 5 plans sont donnés entre crochets. Notons que Eads (sp)
et (Eads (ssp)) sont des énergies d’adsorption calculées avec et sans spin polarisé par calcul de la DFT.
Pour les deux types de calcul, c’est-à-dire, pour les calculs sans spin et avec spin,
l’analyse sur les énergies d’adsorption montre que le site ponté est le plus stable pour les
taux de couverture 0.25, 0.50 et 0.75 MC avec des valeurs respectives de 2.681 [3.170],
(2.289 [2.777]), 2.870 [3.134] eV et (2.480) [2.742]), 2.270 [2.450] et 1.877 [2.057]) eV
comme nous le montre le tableau 3.4.
La différence d’énergies entre le site préféré (ponté) et le plus métastable énergétique-
ment (apical) est constante c’est-à-dire indépendant du type de calcul de polarisation ;
les différences d’énergie d’adsorption sont ainsi de 0.470 [0.476] eV, 0.641 [0.630] eV et
0.950 [0.947] eV, respectivement pour 0.25, 0.50 et 0.75 MC. Nous les avons aussi cal-
culées entre les sites ponté et de symétrie 4 (le moins métastable) et nous trouvons des
49
valeurs de ∼ 0.248 [0.263] eV, 0.524 [0.516] eV et 0.163 [0.098] eV, respectivement.
Nous remarquons que, l’augmentation du nombre de plans du slab influence faible-
ment les différences d’énergie d’adsorption entre les sites préféré et métastables, ex-
ception trouvée pour le taux de couverture de 0.75 MC pour lequel nous avons 0.163
[0.098].
Enfin, pour θ = 1.00 MC, l’hydroxyle est plus stable sur le site de symétrie 4. Les
différences d’énergie d’adsorption entre les sites apical et ponté sont respectivement
de 0.673 [1.220] (0.668 [0.822]) eV et 0.259 [0.729] (0.250 [0.367]) (les nombres entre
parenthèses représentent les valeurs des calculs avec spin polarisé).
Dans cette dernière partie sur la discussion de nos énergies d’adsorption, nous allons
les comparer, à titre illustratif, avec celles d’autres auteurs qui ont travaillé sur les métaux
de transitions et nobles. En utilisant un modèle périodique de slab, Pessoa et al. [99] ont
obtenu les énergies d’adsorption de 1.97, 3.07 et de 3.11 eV pour les sites apical, ponté
et de symétrie 4 respectivement. Ces résultats sont quantitativement supérieurs aux
nôtres. Cependant, en utilisant un modèle de clusters, Pessoa et al. [99] trouvent pour
les sites ponté et de symétrie 4 des valeurs des énergies d’adsorption de 2.53 et 2.77 eV,
respectivement. Par conséquent, nos énergies d’adsorption de OH sur les sites ponté et
de symétries 4 sont proches de celles trouvées par Pessoa et al. [99]. Santiago-Rodriguez
et al. [19] trouvent une valeur de 1.70 eV sur le site ponté du système OH/Au(111).
Cette valeur est faible par rapport à la nôtre ; ceci justifie la faible interaction entre OH
et Au(111) [Au(111) est plus dense que Au(001)]. Zhou et al. [100] à 1/6 MC obtiennent
pour les sites ponté et de symétrie 4, les énergies d’adsorption de 3.093 et 3.367 eV
respectivement pour le système OH/Ni(100). Des désaccords existent entre ces résultats
et les nôtres. Pour le système OH/Ni(111) et à 1/6 MC, Wang et al. [79] dans leurs
travaux trouvent des énergies d’adsorption de 3.071 et 3.185 eV pour les sites apical
et ponté, respectivement. Leurs résultats théoriques sont élevés par rapport aux nôtres.
Cette différence peut s’expliquer par le fait que, Wang et al. [79] dans leur méthode
de calcul, ont utilisé la fonctionnelle GGA-PBE (non-self-consistent) pour un slab de 3
plans de Ni(111) et une énergie de coupure de 340 eV ; par contre, nous avons utilisé la
GGA-PW91 (self-consistent) pour un slab de 4 et 5 plans de Au(001), avec une énergie
de coupure de 400 eV. Par contre, notre valeur de 2.210 eV à 4 plans et pour θ = 0.25
MC sur le site apical est en bon accord avec celle de 2.13 eV trouvée par Ford et al. [98]
pour le même site à 3 plans de Pt(111). Enfin, Hu et Nakatsuji [96] dans leur calcul sur le
sytème OH/Ag(100) en utilisant le modèle de cluster, trouvent des énergies d’adsorption
de 4.211, 4.432 et 4.705 eV respectivement pour les sites apical, ponté et de symétrie 4.
50
Ces valeurs sont presque deux fois plus élevées que les nôtres.
3.2.4.2 Travail de sortie (work function) et moment dipolaire
La variation de la work-function ∆Φ entre le système OH/Au(001) et la surface libre
de Au(001) ainsi que le moment dipolaire µ sont listés dans le tableau 3.5.
Tableau 3.5: Variation de la work-function (∆Φ) ; moment dipolaire (µ) obtenus pour le système
OH/Au(001) avec un slab de 4 et 5 plans. Les résultats avec 5 plans sont donnés entre crochets.
θOH Sites ∆Φ µ
(MC) (eV) (Debye)
apical +0.617 +0.572
[+0.430] [+0.399]
0.25 ponté +0.050 +0.046
[-0.158] [-0.147]
symétrie 4 -0.662 -0.614
[-0.921] [-0.854]
apical +1.041 +0.483
[+1.041] [+0.483]
0.50 ponté +0.365 +0.169
[+0.335] [+0.155]
symétrie 4 -1.032 -0.479
[-1.035] [-0.480]
apical -0.412 -0.127
[-0.418] [-0.129]
0.75 ponté +0.153 +0.047
[+0.114] [+0.035]
symétrie 4 -1.150 -0.356
[-1.198] [-0.370]
apical -0.904 -0.210
[-0.917] [-0.213]
1.00 ponté -1.229 -0.285
[-1.309] [-0.304]
symétrie 4 -1.251 -0.290
[-1.271] [-0.295]
51
pour θ = 0.25, 0.50, 0.75 et 1.00 MC pour lesquels nous obtenons ∆Φ = -0.412 [-0.418]
(-0.904 [-0.917]) eV et -0.662 [-0.921] ; -1.032 [-1.035] ; -1.150 [-1.198] ; -1.251 [-1.271] eV,
respectivement. Les résultats similaires ont été obtenus par Moussounda et al. [101] et
Zeng et al. [102] lors des études des systèmes NO/Pt(100) [101] et NO/Pt(111) [102] où
ils trouvent des signes opposés de ∆Φ pour les sites ponté et apical [101], et sur les sites
de symétrie 4 et apical [102]. Probablement, ces résultats proviennent des différentes re-
distributions de charge provoquée par l’adsorption. En effet, sur la figure 3.8, nous avons
visualisé la redistribution de la densité d’électrons localisée principalement au dessus du
plan de surface et pour OH adsorbé à 0.25 MC. Pour OH adsorbé sur le site de symétrie
4 (Figure 3.8b), nous remarquons une forte accumulation de charges autour OH et une
faible accumulation autour des quatre de Au proches atomes voisins. Pour le site apical
(Figure 3.8a), il y a aussi une accumulation de charges de densité d’électrons autour de
OH et une faible densité de charge au niveau des premiers proches atomes de Au, alors
qu’il y’a une accumulation et une faible densité de charge autour des seconds atomes
proches voisins de Au. Ainsi, cette différence significative de distribution de charges entre
les sites de symétrie 4 (hollow) et apical (top) à 0.25 MC semble être responsable du
signe opposé de la variation de la work-function ∆Φ.
Figure 3.8: Vue de profil des isosurfaces de OH adsorbée sur les sites apical (a) et de symétrie 4 (b)
sur la surface Au(001) à 0.25 MC. L’accumulation des charges est en rouge et la déplétion en bleu.
Continuant notre analyse des résultats rapportés dans le tableau 3.5, nous remar-
quons que le moment dipolaire µ présente plusieurs comportements. Premièrement, sur
le site de symétrie 4 nous observons l’augmentation du moment dipolaire avec l’augmen-
tation du taux de couverture. Deuxièmement, la valeur positive de µ est obtenue sur
l’apical (ponté) à 0.25 et 0.50 MC avec des moments dipolaires de 0.572 D et 0.483 D,
respectivement. Les valeurs correspondantes pour le site ponté sont 0.046 D et 0.169 D,
respectivement. Troisièmement, pour 0.25 à 0.75 MC, µ décroit significativement lors de
52
l’adsorption de OH sur le site apical. La dépolarisation du système OH/Au peut expli-
quer cette observation. En outre, pour une couverture de 1.00 MC, le moment dipolaire
µ a approximativement la même valeur. Enfin, pour θ = 0.25 et 0.50 MC, les valeurs de
µ sont dans l’ordre suivant : apical > ponté > symétrie 4.
3.2.4.3 Paramètres géométriques
Nous rappelons que les simulations du système OH/Au(001) ont été faites avec et
sans la polarisation du spin. Les paramètres géométriques optimisés sont rassemblés dans
le tableau 3.6.
Tableau 3.6: Longueur de liaison de OH (dO−H ), longueur de liaison entre O et Au (dO−Au ), angles
(H-O-Au) obtenus pour le système OH/Au(001) avec un slab de 4 et 5 plans. Les résultats avec 5 plans
sont donnés entre crochets.
Nous avons remarqué que la polarisation du spin n’influence pas les paramètres géo-
53
métriques. De même, l’augmentation du nombre de plans du slab et le taux de couverture
ont une influence très limitée sur les paramètres géométriques. En observant le tableau
3.6, nous remarquons que les longueurs de liaison O-H sont situées entre 0.984 Å et
1.015 Å ; par exemple, on assiste soit à une compression, soit à un étirement de dO−H
par rapport à la phase gazeuse (0.993 Å). Santiago-Rodriguez et al. [19] et Ford et al.
[98] ont étudié l’adsorption de OH sur Au(111) et Pt(111), respectivement. Ils s’ac-
cordent pour une longueur de liaison dO−H de 0.990 Å, qui est située dans le domaine
de variation de nos résultats. En utilisant la méthode des clusters combinée avec les
approximations de Hartree-Fock (HF) et de Moller-Plesset 2 (MP2), Hu et Nakatsuji
[96] étudient l’adsorption de OH/Ag(100) et calculent les longueurs de liaison de 0.94 Å,
0.94 Å et 0.95 Å sur les sites apical, ponté et de symétrie 4, respectivement. Ces valeurs
peuvent être favorablement comparées aux nôtres. La distance dO−Au augmente avec la
coordination, c’est-à-dire, nous trouvons des valeurs respectives de ∼ 2.084 Å, 2.219 Å
et 2.503 Å sur les sites apical, ponté et de symétrie 4. Cette tendance à l’augmentation
de la distance dO−Au avec la coordination a été aussi confirmée par Hu et Nakatsuji [96]
sur OH/Ag(100) qui trouvent des valeurs 2.137 Å, 2.310 Å et 2.513 Å sur les sites apical,
ponté et de symétrie 4. Zhou et al. [100] étudient l’adsorption de OH/Ni(100) sur le site
de symétrie 4 (hollow) et trouvent une distance dO−N i de 2.090 qui est cependant 19.76
% plus petite par rapport à notre valeur.
3.2.5 Propriétés électroniques de OH et de la surface Au(001)
Dans cette sous-section, nous analysons l’évolution de la structure électronique à
travers la densité d’états projetée (PDOS). Ceci nous permet de comprendre la nature
de l’interaction entre OH et la surface Au(001). Nous limiterons cette analyse sur le site
préférentiel d’adsorption (le site ponté) avec deux taux de couverture de 0.25 et 0.50
MC pour lesquels les PDOS sont les plus perturbées. Nous analyserons les PDOS dans
les calculs avec et sans la polarisation de spin.
3.2.5.1 Densité d’états projetée (PDOS) de O dans OH libre et adsorbé à
0.25 et 0.50 MC sans et avec spin polarisé
Sur la figure 3.9 nous présentons la PDOS totale de l’oxygène correspondant à
OH/Au(001) et de O dans OH en phase gazeuse sans spin polarisé. Le niveau de Fermi
se trouve à 0 eV. En phase gazeuse, nous voyons deux pics (états inoccupés) localisés
respectivement à +1.91 et +5.14 eV. Sous adsorption, nous remarquons que l’augmen-
tation du taux de couverture n’influence pas la PDOS de O de OH à l’exception du pic
dominant situé à -1.53 eV. L’effet le plus remarquable est le changement de la PDOS
54
sous adsorption par rapport à celle de la phase gazeuse pour les états occupés. En effet,
les deux pics dominants au dessus du niveau de Fermi se trouvent dispersés de façon
continue entre 0 eV et +7 eV. Ces changements s’expliquent par un mélange entre les
orbitales s et p de l’oxygène avec les bandes s, p et d de la surface Au.
Figure 3.9: PDOS de O dans OH libre et OH/Au(001) sur le site ponté à 0.25 MC et 0.50 MC.
Par contre, sur la figure 3.10 nous analysons, la PDOS de l’oxygène dans OH libre
et adsorbé à 0.25 et 0.50 MC dans un calcul de polarisation du spin. Nous remarquons
que sur la phase gazeuse, les états du spin up ont deux pics localisés à +0.89 et +3.52
eV. Par contre, les états de spin down ont deux pics situés à +1.50 et +5.34 eV. Sous
adsorption, nous notons que la PDOS de O est très affectée. Ainsi, les états sont trans-
formés en continuum et s’étendent sur l’intervalle [-2.55 eV, +7.56 eV]. En conséquence,
les interactions entre les orbitales de O et les bandes de Au annihilent le spin et font
apparaitre à 0.25 MC, un nouveau pic situé à -1.53 eV sur les deux états du spin (up et
down).
55
Figure 3.10: PDOS de O avec spin polarisé dans OH libre et OH/Au(001) sur le site ponté à 0.25 MC
et 0.50 MC.
56
Figure 3.11: PDOS de la surface de Au libre et adsorbée à 0.25 MC et 0.50 MC sans spin polarisé.
Nous montrons sur la figure 3.12, la PDOS des bandes s, p et d de la surface Au avec
spin polarisé. Nous remarquons un faible changement de la PDOS sur les deux taux de
couvertures (0.25 et 0.50 MC) par rapport à celle de la surface libre. Les états occupés
ainsi qu’inoccupés sont similaires pour les spins up et down. Nous voyons aussi que la
PDOS de la bande p de Au est plus perturbée que celle des bandes s et d. Ces effets
peuvent être expliqués par un mélange entre les orbitales de l’hydroxyle et les bandes de
la surface de Au.
57
Figure 3.12: PDOS de la surface de Au libre et adsorbée à 0.25 MC et 0.50 MC avec spin polarisé.
58
•OH(g) −→ OH(a) −→ O(a) + H(a) (3.6)
59
Figure 3.13: Thermochimie de la décomposition des espèces moléculaires CH3 et OH. La référence
zéro représente l’énergie de l’espèce en phase gazeuse plus l’énergie de la surface libre de Au(001) sur
différentes étapes. Les phases gazeuses sont indiquées par (g) et l’espèce adsorbée par (a).
3.4 Conclusion
Dans ce chapitre, nous avons fait l’étude par la DFT-GGA (PW91) de l’adsorption
de CH3 et OH sur la surface Au(001) pour les taux de couverture de θ = 0.25, 0.50, 0.75
et 1.00 MC. Les simulations ont été menées avec des mailles (2x2) et (1x1) construites
pour deux slabs de 4 et 5 plans.
Nous avons fait l’étude avec et sans la polarisation de spin. Nous avons constaté que
la polarisation de spin réduit considérablement les énergies d’adsorption. Cependant, les
paramètres structuraux des adsorbats ne sont pas affectés par le calcul magnétique.
60
Dans la phase d’optimisation de CH3 /Au(001), nous avons trouvé que les sites ponté
(bridge) et de symétrie 4 (hollow) étaient instables. Dans les deux cas, le radical CH3
convergeait vers le site apical (top). Les énergies d’adsorption dépendent du taux de
couverture et du nombre de plans de slab de Au(001). Les paramètres géométriques ne
dépendent pas du nombre de plans. La variation de la work function ∆Φ diminue avec
l’augmentation du taux de couverture par contre le moment dipolaire µ augmente. Ceci
est dù à la polarisation de la surface Au(001). Les changements sur les PDOS ont révélé
un fort mélange entre les orbitales de CH3 et les bandes s, p et d de Au(001).
Nous avons constaté que le fragment OH est stable sur les trois sites d’adsorption.
L’augmentation du nombre de plans de slab impacte les énergies d’adsorption. A l’excep-
tion du taux de couverture de 1.00 MC, la préférence du site d’adsorption est dans cet
ordre : ponté > symétrie 4 > apical. La différence d’énergie entre le site le plus stable
(ponté) et les métastables (symétrie 4 et apical) a été évaluée. Nous avons constaté des
signes contraires de la variation de la work-function ∆Φ ; la redistribution de charges
est à l’origine de ces faits. Enfin l’analyse sur la PDOS de l’oxygène de OH a montré un
mélange entre les orbitales de OH et les bandes de Au.
Finalement, en étudiant la thermochimie de l’adsorption et de la décomposition di-
recte de CH3 et OH sur la surface de Au(001), la désorption a été trouvée favorable pour
les deux adsorbats.
Notons, pour terminer, que tous ces résultats ont fait l’objet d’une publication
scientifique paru dans Surface Review and Letters (Bouka-Pivoteau et al., Surface Re-
view and Letters, 1850198 (16 pages) c World Scientific Publishing Company DOI :
10.1142/S0218625X18501986) [103]
61
Chapitre 4
Dans cette partie, nous avons travaillé sur des slabs de 4 et 5 plans de la surface de
Au et pour le calcul des énergies d’adsorption, nous avons utilisé la formule suivante :
où Esyst est l’énergie totale du système formée du slab Au(001) et du fragment après
adsorption ; Eslab , l’énergie totale de Au(001) nue et Ef ragment , l’énergie du fragment
libre.
Notons que, les fragments étudiés dans ce chapitre sont : l’hydroxyméthyle (CH2 OH),
le formaldéhyde (H2 CO) et le formyle (HCO). Ces études ont un intérêt majeur parce
que ces adsorbats sont des fragments hydrocarbonés dont il est possible d’en extraire de
l’hydrogène qui peut être utilisé comme carburant propre pour les véhicules.
62
fonctionnelle de la densité basée sur le modèle du slab.
Expérimentalement plusieurs techniques ont été aussi réalisées sur l’étude de ce frag-
ment (CH2 OH) parmi lesquelles : la spectroscopie sur la désorption en température
programmée (TPD, Temperature programmed desorption), la spectroscopie à haute ré-
solution de perte d’énergies d’électrons (HREELS, High resolution electron energy loss
spectroscopy) [79], la spectroscopie de photo émission par des rayons UV (UPS, Ultra
violet photoemission spectroscopy), la spectroscopie de photoémission par des rayons
X (XPS, X- ray photoemission spectroscopy), la spectroscopie des rayons Auger (AES,
Auger electron spectroscopy) [104]. Ces études ont été menées pour la recherche de la
stabilité des sites d’adsorption, des énergies d’adsorption ainsi que les structures géomé-
triques et électroniques.
4.1.2 Géométrie et paramètres structuraux de CH2 OH en phase gazeuse
Rappelons que, le radical hydroxyméthyle peut être issus de la rupture d’une liaison
(C-H) du groupement CH3 dans la molécule du méthanol. En phase gazeuse, nous avons
déterminé la géométrie (Figure 4.1) et des paramètres structuraux théoriques (Tableau
4.1) de l’hydroxyméthyle.
Figure 4.1: Géométrie de l’hydroxyméthyle (CH2 OH) en phase gazeuse. La sphère en rouge représente
l’atome d’oxygène, en noir l’atome de carbone et en blanc l’atome d’hydrogène.
63
Ce tableau 4.1 montre que le résultat obtenu pour la longueur de liaison, dC−O , est
en bon accord raisonnable avec les études de Wang et al.[79]
4.1.3 Site d’adsorption de CH2 OH sur la surface Au(001)
Nous avons réalisé ce travail pour un taux de couverture de 0.25 MC obtenu après
adsorption de CH2 OH sur la maille (2x2) de la surface Au(001) ; le slab utilisé étant
composé de 4 et 5 plans. Dans les phases d’optimisation, seuls les sites d’adsorption
apical (top) et ponté (bridge) se sont avérés stables. Le site de symétrie 4 (hollow) est
instable, le radical CH2 OH dans ce cas migre vers le site de ponté et l’hydrogène lié à
l’atome d’oxygène se détache du groupement (Figure 4.2b). Notons aussi que pour les
deux sites stables, la configuration initiale (c’est-à-dire avant optimisation) représentant
l’adsorption de CH2 OH/Au(001) s’est faite par l’oxygène (O). Nous avons remarqué que
les configurations finales (c.-à-d. après optimisation) se caractérisent par une adsorption
qui ne se fait plus par l’oxygène mais plutôt par l’hydrogène comme on peut le voir sur les
figures 4.3b et 4.4b. Ceci n’est pas nouveau, car Wang et al. [79] dans leurs travaux sur
Ni(111) ont signifié que l’adsorption de CH2 OH était favorable par le carbone. Greeley
et al. [105] et, Kua et Goddard [106] ont aussi indiqué que l’adsorption de CH2 OH était
plus favorable par le carbone que par l’oxygène.
Figure 4.2: Géométrie de l’hydroxyméthyle (CH2 OH) sur la surface de Au(001) avant (a) et après (b)
adsorption sur le site de symétrie 4. La sphère en rouge représente l’atome d’oxygène, en noir l’atome
de carbone et en blanc l’atome d’hydrogène.
64
Figure 4.3: Géométrie de l’hydroxyméthyle (CH2 OH) sur la surface de Au(001) avant (a) et après (b)
adsorption sur le site de apical. La sphère en rouge représente l’atome d’oxygène, en noir l’atome de
carbone et en blanc l’atome d’hydrogène.
Figure 4.4: Géométrie de l’hydroxyméthyle (CH2 OH) sur la surface de Au(001) avant (a) et après (b)
adsorption sur le site de ponté. La sphère en rouge représente l’atome d’oxygène, en noir l’atome de
carbone et en blanc l’atome d’hydrogène.
65
4.1.4.1 Energies d’adsorption
Pour calculer les énergies d’adsorption, nous avons utilisé la formule 4.1 ; les résultats
sont présentés dans le tableau 4.2.
Tableau 4.2: Energies d’adsorption (Eads ) obtenues avec un slab de 4 plans après adsorption de
CH2 OH sur Au(001). Les résultats entre crochets sont calculés avec un slab de 5 plans.
apical -0.972
[-1.612]
0.25 ponté -1.899
[-2.484]
66
(top) ; ils s’accordent pour les longueurs de liaison dC−H (1.098 Å) et dC−O (1.476 Å). Ces
valeurs sont quantitativement dans le bon ordre de grandeur aux nôtres. Zhou et al. [100]
trouvent une longueur de liaison dC−O de 1.471 Å dans leur étude de CH2 OH/Ni(100)
sur le site apical (top). Cette valeur est légèrement supérieure à la nôtre.
Tableau 4.3: Longueurs de liaison entre C et H (dC−H ), entre C et O (dC−O ), entre O et H (dO−H ),
entre O et Au (dO−Au ), entre C et Au (dC−Au ) et, les angles (C-O-H) et (H-C-H) obtenus avec un slab
de 4 plans après adsorption de CH2 OH sur Au(001). Les résultats entre crochets sont calculés avec un
slab de 5 plans. Les valeurs entre parenthèses représentent le nombre de longueurs de liaison possibles.
67
Figure 4.5: PDOS sur les orbitales O de CH2 OH libre (en pointillé) par rapport à celles dans
CH2 OH/Au(001) adsorbées (en trait continu) pour le site apical. Le niveau de Fermi est à 0 eV.
La figure 4.6 représente la PDOS de H (du groupement OH) de CH2 OH avant et après
adsorption. Nous avons une situation presque similaire à la figure 4.5. Après adsorption,
les PDOS se décalent vers les basses énergies par rapport à celles correspondantes à la
phase gazeuse. Les faits les plus marquants sont : (i) l’apparition d’un pic de résonance
localisé à ∼ -2.70 eV ; (ii) la transformation des pics d’intensité situés au-delà du niveau
de Fermi (notamment à +3.56 eV et +6.79 eV) en un continuum qui s’étend dans un
domaine d’énergie [0 eV, +10 eV]. Toutes ces transformations des PDOS observées sous
adsorption s’expliquent par l’interaction forte entre le radical hydroxyméthyle et les
atomes de surface de Au.
68
Figure 4.6: PDOS sur les orbitales H de CH2 OH libre (en pointillé) par rapport à celles dans
CH2 OH/Au(001) adsorbées (en trait continu) pour le site apical. Le niveau de Fermi est à 0 eV.
Sur la figure 4.7, nous présentons les PDOS des bandes s, p et d des atomes Au de
surface dans CH2 OH/Au(001) et la surface nue. Globalement, les états occupés de ces
bandes ne sont perturbés par l’adsorption. Notons une augmentation des états inoccupés
sur les bandes s et p de Au, et un déplacement du centre de gravité de la bande d vers
les basses énergies. Ces changements s’expliquent par un mélange entre les orbitales de
CH2 OH et les bandes s, p et d de Au.
69
Figure 4.7: PDOS des bandes de Au libre (en pointillé) par rapport à celles dans CH2 OH/Au(001)
adsorbées (en trait continu) sur le site apical. Le niveau de Fermi est à 0 eV.
70
Figure 4.8: PDOS sur les orbitales O de CH2 OH libre (en pointillé) par rapport à celles dans
CH2 OH/Au(001) adsorbées (en trait continu) pour le site ponté. Le niveau de Fermi est à 0 eV.
Sur la figure 4.9, nous avons représenté la PDOS de H dans CH2 OH/Au(001) et
CH2 OH libre. Nous notons un effet contraire à la PDOS du même objet examiné à la
figure 4.6. Contrairement à la figure 4.6, les PDOS de la figure 4.9 se décalent légèrement
vers les hautes énergies. Après adsorption, nous remarquons que la PDOS de H (figure
4.9) est marquée par une diminution en intensité du pic situé à -8.96 eV et une légère
augmentation de ceux situés à -5.32 eV et -0.68 eV par rapport à la phase libre. Enfin,
notons que les états inoccupés ont complètement disparu et sont alors dispersés de façon
quasi-continue entre +2.5 eV et +10.0 eV. Le même argument concernant l’explication
de tous ces changements des PDOS reste valable à savoir l’interaction entre CH2 OH et
la surface Au(001).
71
Figure 4.9: PDOS sur l’orbitale H de CH2 OH libre (en pointillé) par rapport à celles dans
CH2 OH/Au(001) adsorbées (en trait continu) pour le site ponté. Le niveau de Fermi est à 0 eV.
Sur la figure 4.10, nous présentons l’évolution des PDOS des bandes s, p et d de Au
avant et après adsorption. L’évolution de ces PDOS est similaire à celle observée sur la
figure 4.7. On note surtout une augmentation des états inoccupés des bande s et p des
atomes Au de surface. La bande d est peu perturbée, mais nous notons un léger décalage
de son centre de gravité vers les basses énergies moins prononcé que celui observé sur la
figure 4.7.
72
Figure 4.10: PDOS des bandes de Au libre (en pointillé) par rapport à celles dans CH2 OH/Au(001)
adsorbées (en trait continu) sur le site ponté. Le niveau de Fermi est à 0 eV.
73
niques en sciences des surfaces. La spectroscopie des électrons à basse énergie (EELS,
energy loss spectroscopy) utilisée dans la décomposition du méthanol sur les métaux de
transition conduit à la formation des espèces intermédiaires H2 CO et HCO [108, 109].
L’adsorption du H2 CO a été aussi réalisée par spectroscopie de photo émission des rayons
UV sur la surface Ru(110) [110] et par spectroscopie à haute résolution (HREELS, high-
resolution electron energy loss spectroscopy) sur l’oxyde du nickel (NiO). La spectro-
scopie infrarouge par transformation de Fourrier (FTIR, Fourrier transform infrared
spectroscopy) a été aussi réalisée dans l’étude de l’adsorption de H2 CO/TiO2 [111]. Liu
et al. [112] ont rapporté que sur Pt, pour des faibles taux de couverture, le méthanol se
décomposait en méthoxyle (CH3 O) à des températures de 130 à 170 K, lequel aussi se
décompose en formaldéhyde (H2CO) à des températures entre 170 et 200 K.
Sur le plan théorique, l’adsorption du formaldéhyde (H2 CO) a aussi suscité un grand
intérêt. Certaines études ont été réalisée sur la surface Cu(111) par Gomes et al. [113] qui
utilisent la théorie de la fonctionnelle de densité (DFT) basée sur le modèle de clusters.
Des calculs utilisant la théorie de la fonctionnelle de densité dans l’approche des slabs
ont été réalisés sur les surfaces Ni(100) [100], Ni(111) [79], Pt(111) [114], Au(111) [104]
et sur le graphène [115].
4.2.2 Géométrie et paramètres structuraux de H2 CO en phase gazeuse
Le radical formaldéhyde (H2 CO) peut être issu de la réaction de décomposition de
l’hydroxyméthyle (CH2 OH) suivant le schéma ci-dessous :
74
Figure 4.11: Géométrie du formaldéhyde en phase gazeuse. La sphère en rouge représente l’oxygène,
en noir le carbone et en blanc les hydrogènes.
Nous remarquons que nos résultats sont en accord raisonnable avec les valeurs expé-
rimentales [78] et théoriques [104, 116, 117].
4.2.3 Site d’adsorption de H2 CO sur Au(001)
Dans l’étude de l’adsorption du formaldéhyde, nous avons utilisé un taux de couver-
ture de 0.25 MC réalisé lorsqu’un fragment de celui-ci est déposé sur la surface Au(001).
L’adsorption s’est faite dans la configuration η 1 (O), c.-à-d., via l’atome d’oxygène (O)
avec la liaison dC−O complètement perpendiculaire au plan de la surface Au(001). Dans
la phase d’optimisation, nous avons constaté que le site de symétrie 4 était instable ;
le radical H2 CO migre complètement vers le site ponté (bridge). Ainsi donc, l’adsorp-
tion de H2 CO/Au(001) a été examinée sur les sites apical (top) et ponté (bridge). Enfin
l’adsorption a été faite sur une seule face de deux slabs de 4 et 5 plans.
4.2.4 Energies d’adsorption et paramètres structuraux du formaldéhyde en
phase adsorbée
Les énergies et les paramètres structuraux obtenus après adsorption de H2 CO/Au(001)
sont rassemblés dans le tableau 4.5. Les énergies d’adsorption ont été calculées en utili-
sant la formule 4.1.
75
Tableau 4.5: Energies d’adsorption (Eads ), longueurs de liaison entre C et H (dC−H ), entre C et O
(dC−O ), entre O et Au (dO−Au ) et, les angles (C-O-H) et (H-C-H) obtenus avec un slab de 4 plans après
adsorption de H2 CO sur Au(001). Les résultats entre crochets sont calculés avec un slab de 5 plans. Les
nombres entre parenthèses représentent le nombre de longueurs de liaison ou angles possibles.
• Energies d’adsorption
Nous commençons à exploiter les résultats du tableau 4.5 par l’analyse de l’énergie
d’adsorption. Le premier constat est que l’augmentation du nombre de plans a une
influence significative sur l’énergie d’adsorption du système ; en effet, le rapport E5plans
ads
/ E4plans
ads ≈ 68. Le deuxième constat est que l’adsorption de H 2 CO/Au(001) a une très
faible valeur positive de Eads (+ 9.6 meV) dans le cas d’un slab à 4 plans pour le site
d’adsorption apical (top). Donc, le site apical combiné avec un slab de 4 plans génère une
réaction (désorption) endothermique du système H2 CO/Au(001) par rapport au radical
H2 CO gazeux et la surface Au(001) nue. Un autre constat est que le site préférentiel
d’adsorption de H2 CO est le site ponté pour un slab de 5 plans de Au.
Nous comparons nos énergies d’adsorption avec celles déterminées sur d’autres mé-
taux. Chi et Zhao [115] explorent l’adsorption de H2 CO sur le graphène et calculent
de faibles énergies d’adsorption de -32 meV. Remediakis et al. [121] utilisant l’approche
DFT-slab avec le code Dacapo, calculent une énergie d’adsorption de -0.14 eV lors de
l’adsorption de H2 CO/Ni(111). En utilisant la méthode DFT dans l’approche des clus-
ters sur Cu(111), Gomes et al. [113] calculent des énergies d’adsorption contenues dans
le domaine [-0.228 eV, -0.259 eV]. Sur Au(111), Chen et al. [104] rapportent une énergie
d’adsorption de -0.245 eV. Ces derniers résultats sont en moyenne 3 fois plus faibles
aux nôtres. Cahyanto et al. [117] trouvent pour une même conformation, c.-à-d., une
adsorption par O sur Pt(111) et sur le site apical une valeur de -0.65 eV. Ce résultat est
quantitativement en accord aux nôtres. Les valeurs élevées calculées par Zhou et al. [100]
sur Ni(100) et Ni(111) (sur le site apical) sont de -1.54 eV et -1.025 eV, respectivement.
Ce désaccord profond entre nos résultats et ceux de Zhou et al. [100] peut être expliqué
par : (i) la nature différente entre les deux substrats (Au vs Ni), et (ii) les conformations
différentes dans les adsorptions (nous avons une adsorption par O contre une adsorption
76
par C dans [100]).
• Paramètres structuraux
Les paramètres structuraux obtenus après optimisation du système H2 CO/Au(001)
sont résumés dans le tableau 4.5. Nous remarquons que contrairement à l’énergie d’ad-
sorption, l’augmentation du nombre de plans n’influence pas les paramètres structuraux.
Pour les deux sites d’adsorption examinés ici, les paramètres structuraux ont les mêmes
valeurs. Cependant, comparée à la phase gazeuse, l’adsorption ne modifie pas les lon-
gueurs de liaison C-H et C-O ; par contre elle a tendance à rapprocher les deux atomes
d’hydrogène vers la surface, ce qui tend alors à ouvrir l’angle H-C-H [qui passe 116.2o
(phase gazeuse) à 121.8o (phase adsorbée)] et à fermer l’angle H-C-O [qui passe 121.9o
(phase gazeuse) à 116.4o (phase adsorbée)]. Sur le site ponté, le plus stable, la longueur
de liaison, dO−Au , optimisée est d’environ 3.510 Å. Cependant, elle est d’environ 3.110
Å dans le cas du site apical, le métastable. Comme nous l’avons dit dans nos conclusions
antérieures, l’augmentation de la distance avec la coordination au site d’adsorption est
conforme au comportement général constaté dans les études expérimentales et théoriques
de l’adsorption atomique.
Chen et al. [104] ont étudié l’adsorption de H2 CO/Au(111) sur le site apical (top) ;
ils rapportent une distance, dO−Au , de 3.868 Å, qui est 0.758 Å plus élevée par rapport
à notre résultat. En outre, ces auteurs [104] trouvent deux hydrogènes inéquivalents
avec des longueurs de liaison, dC−H , de 1.108 et 1.162 Å. Dans nos calculs, les deux
hydrogènes sont équivalents avec une longueur de liaison, dC−H , de ∼ 1.115 Å. Enfin,
ces auteurs [104] prédisent un angle H-C-H de 121.58o , qui est plus ouvert, comparé avec
notre résultat (116.4o ).
Sur Cu(111), Gomes et al. [113] rapportent les paramètres structuraux dC−O , dC−H
et H-C-H et H-C-O de 1.121 Å, 1.109 Å, 119.8o et 121.6o , respectivement. Les valeurs
des longueurs de liaison dC−O et dC−H sont en accord avec les nôtres. Cependant, les
différences existent entre les valeurs des angles H-C-H et H-C-O sur Au(001) (121.8o ,
116.4o ) et celles sur Cu(111) (119.8o , 121.6o ).
4.2.5 Structure électronique
Pour élucider la nature de l’interaction entre H2 CO et la surface Au, nous allons
analyser les changements des PDOS de l’oxygène dans H2 CO et des atomes de surface
de Au avant et après adsorption de H2 CO dans les deux sites d’adsorption apical (top)
et ponté (bridge) examinés ici.
77
4.2.5.1 Site apical
Sur la figure 4.12, nous avons représenté la PDOS de O dans H2 CO libre et adsorbée.
En phase gazeuse, la PDOS de O dans H2 CO est caractérisée par des pics localisés à
-7.41, -3.77, -2.58, -1.70, +2.99 et +5.52 eV. Sous adsorption, nous remarquons que les
PDOS se décalent de ∼ 0.6 eV vers les basses énergies. Les pics situés à -7.41, -3.77 et
+2.99 eV augmentent légèrement en intensité ; par contre, ceux localisés à -2.58, -1.70
et +5.52 eV diminuent considérablement en intensité. Tous ces changements des PDOS
peuvent s’expliquer par l’interaction entre les orbitales de O dans H2 CO et les bandes
de Au.
Figure 4.12: PDOS sur l’orbitale O de H2 CO libre (en pointillé) par rapport à celles dans
H2 CO/Au(001) adsorbées (en trait continu) pour le site apical. Le niveau de Fermi est à 0 eV.
La figure 4.13 représente les PDOS des bandes s, p et d des atomes de surface,
après adsorption. Les PDOS correspondantes à la surface nue ont été incorporées pour
d’éventuelles comparaisons. En général, on observe qu’après adsorption, les PDOS des
bandes s, p et d de Au ne sont pas affectées. Cependant, on note une légére augmentation
de l’intensité du pic de la bande s située à +5.16 eV comme le montre la figure 4.13(a)
mais aussi, celui de la bande p situé presque au même endroit (Figure 4.13(b)). Ceci
montre qu’il y’a une très faible interaction entre les bande s, p et d avec les orbitales de
78
O de H2 CO.
Figure 4.13: PDOS des bandes de Au libre (en pointillé) par rapport à celles dans H2 CO/Au(001)
adsorbées (en trait continu) sur le site apical. Le niveau de Fermi est à 0 eV.
79
Figure 4.14: PDOS sur l’orbitale O de H2 CO libre (en pointillé) par rapport à celles dans
H2 CO/Au(001) adsorbées (en trait continu) pour le site ponté. Le niveau de Fermi est à 0 eV.
Dans la figure 4.15, nous comparons les PDOS des bandes s, p et d de la surface nue
de Au avec celles de la surface adsorbée. Sous adsorption, l’évolution des PDOS de ces
bandes est similaire à celle observée dans le cas du site apical. Ainsi, nous remarquons
que les bandes s et p sont faiblement affectées ; par contre, les bandes d ne sont pas
affectées par cette interaction. Ceci montre que l’adsorption de H2 CO ne perturbe pas
d’une façon significative les bandes de Au.
80
Figure 4.15: PDOS des bandes de Au libre (en pointillé) par rapport à celles dans H2 CO/Au(001)
adsorbées (en trait continu) sur le site ponté. Le niveau de Fermi est à 0 eV.
81
programmée en température (TPD) sur la phase gazeuse de l’oxydation du méthanol sur
Ni(100) [118] et par spectroscopie infrarouge (IR) sur Ni(100) [119]. Aux basses tempéra-
tures, Michell et al. [120] ont utilisé la technique de spectroscopie à haute résolution des
électrons à perte d’énergie (HREELS) pour étudier l’adsorption de HCO sur la surface
Ru(001).
Sur le plan théorique, l’adsorption du formyle a aussi fait l’objet de plusieurs études.
Santiago-Rodrìguez et al. [19] utilisent la théorie de la fonctionnelle de la densité (DFT-
GGA) pour étudier l’adsorption du formyle sur la surface Au(111). D’autres études ont
été menées par d’autres chercheurs. En utilisant l’approche DFT-slab (GGA-PW91), Xu
et al. [122] étudient l’adsorption de HCO/Fe(110) et rapportent une énergie d’adsorption
de -2.75 eV. Morgan et al. [123] et Remediakis et al. [121] ont étudié l’adsorption de HCO
sur Ru(0001) et Ni(111), respectivement.
4.3.2 Géométrie et paramètres structuraux de HCO en phase gazeuse
Le radical formyle (HCO) est aussi l’un des fragments hydrocarbonés généré par la
décomposition du méthanol qui conduit à la formation du radical formaldéhyde (H2 CO)
qui se décompose lui aussi selon la réaction :
Figure 4.16: Géométrie du formyle (HCO) en phase gazeuse. La sphère en rouge représente l’atome
d’oxygène, en noir l’atome de carbone et en blanc l’atome d’hydrogène.
82
Tableau 4.6: Comparaison entre nos paramètres structuraux calculés et plusieurs autres études théo-
riques et expérimentale pour HCO en phase gazeuse.
Notons que ces paramètres structuraux de HCO libre sont en accord raisonnable avec
les résultats théoriques antérieurs [117, 121, 124] ; mais ils sont légèrement supérieurs
(respectivement inférieurs) de 2.06 % (respectivement 3.14 %) aux longueurs de liaison
(respectivement aux angles) expérimentales [78].
4.3.3 Site d’adsorption de HCO sur la surface Au(001)
L’adsorption du radical HCO a été étudiée sur les trois sites d’adsorption apical
(top), ponté (bridge) et de symétrie 4 (hollow) de la surface (2x2) de Au(001). Nous
avons travaillé avec un taux de couverture de 0.25 obtenu en adsorbant HCO sur la
cellule (2x2) de 4 et 5 plans de Au, pour une densité de points-k = (6 x 6 x 1).
Dans les phases d’optimisation du radical formyle (HCO), nous avons trouvé que
les sites ponté et de symétrie 4 sont instables. Le fragment HCO sur ces sites se dirige
toujours vers le site apical (top).
4.3.4 Energie d’adsorption et paramètres géométriques du formyle (HCO)
Nous avons déterminé l’énergie d’adsorption et les paramètres géométriques du ra-
dical formyle (HCO), pour le site apical de haute symétrie de la surface Au(001). Ces
grandeurs sont rassemblées dans le tableau 4.7.
Tableau 4.7: Energies d’adsorption (Eads ), longueurs de liaison entre C et H (dC−H ), entre C et O
(dC−O ), entre C et Au (dC−Au ) et l’angle (H-C-O) obtenus après adsorption de HCO sur Au(001). Les
résultats entre crochets sont calculés avec un slab de 5 plans.
83
• Energie d’adsorption
Dans le tableau 4.7, on note clairement que l’augmentation du nombre de plans a
une influence considérable sur l’énergie d’adsorption. L’énergie d’adsorption augmente
de ∼ 0.60 eV quand on passe de 4 à 5 plans.
A notre connaissance, l’adsorption de HCO/Au(001) n’a jamais été examinée ni expé-
rimentalement ni théoriquement. Des études théoriques similaires sur d’autres métaux
de transition ou noble existent dans la littérature. Zhou et al. [100] utilisent le code
STATE (Simulation Tool for Atom Technology) pour étudier l’adsorption de HCO sur
Ni(111) et Ni(100) dans l’approximation DFT-slab (GGA-PBE). Ils rapportent des éner-
gies d’adsorption respectives de -2.41 eV et -3.15 eV. Sur Pt(111), Cahyanto et al. [117]
calculent une énergie d’adsorption de -2.89 eV. Notons que les résultats de Zhou et al.
[100] ont été obtenus sur une maille p(2x3) de 3 plans de Ni(111) et Ni(100). Par contre
Cahyanto et al. [117] utilisent une maille (4x4) de 3 plans de Pt(111). Comparées à notre
résultat (-1.748 eV) de 4 plans, les valeurs de ces auteurs [100, 117] sont en désaccord
avec les nôtres. Remediakis et al. [121] utilisant le code Dacapo dans l’approximation
GGA-RPBE, calculent, avec un slab de 3 plans de Ni(111), une énergie d’adsorption de
-1.79 eV à comparer avec notre résultat (-1.748 eV) obtenu avec un slab de 4 plans de
Au.
• Paramètres structuraux
Les paramètres structuraux obtenus après optimisation du système HCO/Au(001)
sont aussi résumés dans le tableau 4.7. Nous remarquons que l’augmentation du nombre
de plans du slab à une influence très limitée sur les paramètres structuraux. Sous ad-
sorption, nous remarquons qu’en général, les paramètres géométriques se rétrécissent
par rapport à la phase gazeuse, à l’exception de la longueur de liaison dC−O pour la-
quelle nous avons un effet contraire. Ainsi donc, la longueur dC−H est réduite de 0.022
Å, tandis que dC−O augmente de 0.014 Å par rapport au formyle libre. Ceci suppose
que la rupture de la liaison dC−O est plus favorable dans la réaction de décomposition
de HCO. Cahyanto et al. [117] ont étudié l’adsorption de HCO/Pt(111). Ils rapportent
la distance dC−O de 1.20 Å, la distance dC−H de 1.14 Å et la distance dC−P t de 2.02 Å.
Ces valeurs sont en accord raisonnable avec nos calculs.
4.3.5 Structure électronique
Comme nous l’avons déjà souligné ci-haut, l’adsorption de HCO/Au(001) s’est faite
sur le site apical (top). Nous nous contentons ici d’étudier la structure électronique par
l’analyse de la PDOS de HCO en phases gazeuse et adsorbée. Rappelons que l’adsorption
84
de HCO s’est faite via l’atome de carbone.
4.3.5.1 PDOS du carbone dans HCO
Sur la figure 4.17, nous présentons la PDOS du carbone (C) dans HCO/Au(001).
Celle de la phase gazeuse a été incorporée pour une comparaison. Après adsorption,
nous remarquons qu’il y a augmentation des intensités des pics situés à -6.46 eV et -3.23
eV, et une diminution de ceux situés à -2.22 eV et -1.42 eV par rapport aux pics sur la
phase libre. Au dessus du niveau de Fermi, les pics (avant adsorption) situés à +4.23
eV ; +6.05 eV et +8.48 eV sont totalement transformés en continuum dans un
domaine d’énergie [0 eV, +10 eV]. Ceci s’explique par un fort mélange entre les orbitales
du carbone de HCO et les bandes de Au.
Figure 4.17: PDOS sur l’orbitale de HCO libre via C (en pointillé) par rapport à celles dans
HCO/Au(001) adsorbées (en trait continu) pour le site apical. Le niveau de Fermi est à 0 eV.
85
4.3.5.2 Structure électronique des bandes s, p et d de Au(001)
Sur la figure 4.18, nous avons représenté les PDOS des bandes s, p et d des atomes
de Au de surface avant et après adsorption.
Après adsorption, nous remarquons que la bande s est fortement modifiée par rapport
à celle de la surface nue (Figure 4.18a). On note principalement, l’apparition des pics
de résonance situés à -6.46 eV ; -3.23 eV et -2.42 eV. Ceux-ci sont aux mêmes endroits
dans la PDOS de C dans HCO (Figure 4.17). Ensuite, nous notons une diminution en
intensité du pic situé à -0.81 eV après adsorption par rapport à celui sur la phase libre.
Enfin, au-delà du niveau de Fermi, c.-à-d., dans l’intervalle [+0.19 eV ; +3.83 eV], nous
notons une augmentation des états inoccupés. Sur la figure 4.18b représentant la PDOS
de la bande p, nous remarquons, comme sur la figure 4.18a, l’apparition des nouveaux
états représentés par des pics de résonance situés à -6.46 eV à -3.23 eV (totalement
absents dans la PDOS de la surface nue). Enfin, une augmentation des états inoccupés
dans le domaine d’énergie [+2.5 eV, +7.5 eV]. La figure 4.18c (montrant la PDOS de d)
présente une évolution de la PDOS similaire aux figures 4.18a et 4.18b. En effet, nous
notons l’apparition de deux pics de résonance, plus intenses, localisés à -6.46 eV et -3.23
eV. Ces deux résonances sont présentes à la fois dans les PDOS de la figure 4.18. Un
dernier fait important à signaler est le déplacement du centre de gravité de la bande vers
les basses énergies. Tous ces changements s’expliquent par un mélange entre les orbitales
de HCO et les bandes s, p et d des atomes Au de surface.
86
Figure 4.18: PDOS des bandes de Au libre (en pointillé) par rapport à celles dans HCO/Au(001)
adsorbées (en trait continu) sur le site apical. Le niveau de Fermi est à 0 eV.
4.3.6 Conclusion
Dans ce chapitre, nous avons apporté notre contribution à l’étude de l’adsorption de
l’hydroxyméthyle (CH2 OH), du formaldéhyde (H2 CO) et du formyle (HCO). L’adsorp-
tion de ces trois fragments s’est faite sur les trois sites (apical, ponté et de symétrie 4)
de haute symétrie de la surface Au(001), pour un taux de couverture de 0.25 MC obtenu
en déposant un fragment sur la maille unitaire (2x2) de la surface Au(001).
L’adsorption de l’hydroxyméthyle (CH2 OH) a commencé par nous révéler des sur-
prises. Dans la phase d’optimisation, seuls les sites d’adsorption apical (top) et ponté
(bridge) se sont avérés stables. Pour ces deux sites d’adsorption, la configuration initiale
87
via l’oxygène se termine par une adsorption via l’hydrogène. Nous avons obtenu des éner-
gies qui augmentent lorsque le nombre de plans dans le slab augmente de 4 à 5 plans.
Nous avons remarqué que la préférence au site d’adsorption est dans cet ordre : ponté
> apical. Les paramètres structuraux sont peu modifiés avec l’augmentation du nombre
de plans. L’étude des PDOS a révélé que celles des bandes s, p et d des atomes de Au
de surface sont modifiées. Par contre celles de O et H dans CH2 OH se décalent vers les
basses énergies par rapport à leurs positions dans la phase gazeuse. Ces changements des
PDOS s’expliquent par un mélange des électrons des orbitales de O ou H dans CH2 OH
avec les bandes s, p et d de Au de surface.
L’adsorption du formaldéhyde (H2 CO) nous a montré que le site de symétrie 4 est
instable. Dans la phase d’optimisation, le fragment sur ce site se dirige toujours vers le
site ponté. Les énergies d’adsorption calculées avec un slab de 4 plans sont très faibles
(∼ 97 meV). Notons qu’avec un slab de 4 plans, une désorption en énergie d’adsorption
est observée. Nous avons remarqué que l’augmentation du nombre de plans influence
considérablement les énergies d’adsorption car on passe de ∼ -97 meV (4 plans) à ∼
-0.658 eV (5 plans). Par contre, les paramètres géométriques ne sont pas affectés par
l’augmentation du nombre de plans. Sur l’analyse des PDOS, nous avons observé que
les évolutions des PDOS de O dans H2 CO et des bande s, p et d des atomes de Au de
surface sont similaires à ceux du système CH2 OH/Au(001).
L’adsorption du formyle (HCO) a montré que seul le site apical est stable. Ici aussi,
nous avons constaté que les énergies d’adsorption augmentent lorsque le nombre de plans
dans le slab augmente. Mais cette augmentation du nombre de plans a une influence très
limitée sur les paramètres géométriques. L’analyse des PDOS de C dans HCO et des
bandes s, p et d fait apparaître des nouveaux états représentés par des pics de résonance.
Ces nouveaux états s’expliquent par l’interaction entre HCO et la surface Au(001).
Enfin, les énergies d’adsorption calculées en utilisant la fonctionnelle GGA PW91
décroissent dans l’ordre suivant : H2 CO > CH2 OH > HCO.
88
Chapitre 5
89
sur Pt(111) a été étudiée par des techniques SCAC [132], LEED [133, 134], TDS [135],
EELS [135], RAIRS [136], SFG [137, 138] et calorimétriques [139]. De même, l’adsorption
de CO sur Pd(111) a été étudiée selon les techniques TDS [140, 141], SFG [141, 142],
LEED [143], AES [143], RAIRS [144] et STM [146-148] et micro-calorimétriques [142].
L’adsorption et l’énergétique de CO/Rh(111) a fait aussi l’objet de plusieurs études
expérimentales. La technique LEED [149-151], a été utilisée pour déterminer les sites
d’adsorption de CO aux basses et hautes couvertures. De même, l’adsorption de CO sur
Ir(111) [152-156], Fe(110) [157-168], Au(111) [169-171] et Ru(0001) [172-182] a été étu-
diée en utilisant les techniques FT-IRAS, TPD, LEED, HAS, TDS, XPS, EELS, FTIR,
DRIFTS, STM, 129 Xe NMR, FT-ICR, SFG, TDMS et UPS.
De nombreuses études ab initio et DFT ont été réalisées sur l’adsorption de CO
sur Pt(111) [98, 139][183-185], Pd(111) [185], Rh(111) [186, 187], Ir(111) [72], Cu [188],
Fe(110) [122], Ni(111) [79, 100, 189], Au(111) [19, 48], Re(0001) [190] et Ru(0001) [191].
Tous ces travaux ont permis d’étudier la réactivité de CO sur ces surfaces conduisant
à l’évaluation des sites d’adsorption les plus stables et les énergies d’adsorption corres-
pondantes, les fréquences de vibration et les effets des surfaces.
5.1.2 Géométrie et paramètres structuraux de CO en phase gazeuse
La géométrie de CO en phase gazeuse a été obtenue en le plaçant verticalement dans
une boîte cubique d’arête 15 Å, avec un atome de carbone positionné en position basse
par rapport à l’oxygène. La géométrie optimisée de CO en phase gazeuse est représentée
sur la figure 5.1.
Figure 5.1: Géométrie du monoxyde de carbone en phase gazeuse. En rouge l’atome d’oxygène et en
noir l’atome de carbone.
Nous avons calculé une longueur de liaison, dC−O , de 1.160 Å. Cette valeur est en
bon accord avec les résultats théorique (1.140 Å) [189] et expérimental (1.128 Å) [78].
90
5.1.3 Sites d’adsorption de CO sur la surface Au(001)
Dans nos simulations, nous avons étudié l’adsorption de CO sur les trois sites de
haute symétrie de la surface Au(001) : apical (top), ponté (bridge) et de symétrie 4
(hollow). Signalons que cette étude est faite avec un taux de couverture de 0.25 obtenu
en adsorbant une molécule de CO sur la maille (2x2) de la surface Au(001) représentée
par deux slabs de 4 et 5 plans.
5.1.4 Energies d’adsorption et paramètres géométriques du monoxyde de
carbone après adsorption sur Au(001)
L’énergie d’adsorption de la molécule de CO a été déterminée en utilisant la formule
suivante :
où les trois termes du membre de droite de l’équation (5.1) représentent les énergies
totales de la surface avec la molécule adsorbée, de la surface nue, et du CO en phase
gazeuse, respectivement. Avec cette définition, une valeur négative de Eads correspond à
une situation stable de CO sur la surface.
Les molécules CO ont été adsorbées dans les trois sites d’adsorption précitées afin de
déterminer le site préférentiel d’adsorption. La molécule CO a été adsorbée perpendicu-
lairement à la surface avec l’atome de carbone lié à la surface, comme rapportée dans
un certain nombre d’expériences [193-195]. Les énergies d’adsorption et les paramètres
structuraux de CO/Au(001) sont résumés dans le tableau 5.1
Tableau 5.1: Energies d’adsorption (Eads ), longueurs de liaison entre C et O (dC−O ) et entre C et
Au (dC−Au ) obtenues avec un slab de 4 et 5 plans après adsorption de CO sur Au(001). Les résultats
entre crochets sont calculés avec un slab de 5 plans.
Le tableau 5.1 montre que les énergies d’adsorption augmentent avec l’augmentation
91
du nombre de plans sur tous les sites. Nous trouvons que le site ponté (bridge) est le site
préférentiel d’adsorption de CO/Au(001). Ceci est en accord avec les travaux de Meh-
mood et al. [48] qui ont étudié l’adsorption de CO/Au(100) avec une couverture de 0.50
MC. L’adsorption de CO a été examinée théoriquement sur d’autres métaux de transi-
tion et noble. Santiago-Rodrìguez et al. [19] ont étudié l’adsorption de CO/Au(111) avec
4 plans de Au et une couverture de 0.25 MC. Ils rapportent sur le site bridge a une éner-
gie d’adsorption de -0.36 eV qui est 2 fois plus petite que la nôtre. La valeur de l’énergie
d’adsorption expérimentale de CO/Au(111) obtenue en utilisant la technique TPD est
de -0.43 eV [195] . Yim et ses collaborateurs [196] ont publié les énergies d’adsorption
expérimentales qui varient dans le domaine de [-0.46 eV, -0.56 eV], à comparer à notre
valeur pour le site ponté (bridge). Des études similaires ont été menées sur d’autres
surfaces, par exemple Pt(111) [98], Rh(111) [187], Pd(111) [185] sur lesquels CO a des
énergies d’adsorption respectives de -1.89 eV, -1.89 eV et -1.77 eV sur le site ponté.
Quant aux paramètres géométriques, nous constatons que l’augmentation du nombre
de plans n’a aucun effet sur ceux-ci. Cependant, la longueur de liaison, dC−O , passe de
1.160 Å (phase gazeuse) à 1.190 Å (phase adsorbée). Cette valeur extrême (1.190 Å) est
en accord avec le résultat de Santiago-Rodriguez et al. [19]. Sur le site ponté, Mehmood
et al. [48] étudient l’adsorption de CO/Au(100) et trouvent des distances, dC−O , et dC−Au
de 1.16 Å et 2.03 Å. Ces valeurs sont en bon accord avec les nôtres (1.18 Å et 2.13 Å).
5.1.5 Structure électronique
Comme dans les chapitres précédents, nous analysons les structures électroniques de
la surface Au(001) et de la molécule CO à travers l’examen des densités d’état projetées
(PDOS) des orbitales moléculaires de CO et des bandes s, p et d des atomes de Au de
surface. Pour ce taux de couverture de 0.25, nous allons examiner les changements des
PDOS des orbitales moléculaires de CO et des atomes de surface de Au avant et après
adsorption de CO dans les trois sites.
5.1.5.1 Site apical
Ici, nous analysons les PDOS des orbitales moléculaires de CO et des bandes s, p et
d de Au.
Les PDOS des orbitales moléculaires de CO/Au(001) et de CO libre sont représentées
sur la figure 5.2. Les orbitales moléculaires de CO libre sont nommées 4σ, 1π, 5σ et
2π, et elles sont positionnées par rapport à leur localisation en énergie sur la figure
5.2. Ces états moléculaires sont localisés à -2.71 eV, -0.25 eV, +2.36 eV, +9.03 eV,
respectivement. A partir de la figure 5.2, nous pouvons voir que l’état moléculaire 2π
92
est le plus affecté par l’adsorption. En effet, cet état 2π, vide dans la phase gazeuse
est totalement transformé en « continuum » après adsorption. Nous remarquons aussi
que les états moléculaires occupés (4σ et 1π) et inoccupés (5σ) dans la phase gazeuse,
se décalent respectivement vers les basses énergies de 1.0 eV, 1.01 eV et 1.39 eV, par
rapport à la phase adsorbée. Notons enfin que les intensités des états moléculaires 1π,
5σ diminuent sensiblement tandis que celle de l’orbitale moléculaire 4σ ne change pas.
Tous ces changements s’expliquent par l’interaction entre ces orbitales moléculaires et
la bande p de la surface.
Figure 5.2: PDOS sur les orbitales de CO libre (en pointillés) par rapport à celles dans CO/Au(001)
adsorbées (en trait continu) pour le site apical. Le niveau de Fermi est à 0 eV.
Les PDOS des bandes s, p et d des atomes de Au de la surface nue et adsorbée sont
représentées sur la figure 5.3. Il apparaît sur la figure 5.3 que les états électroniques
des bandes s et d ne sont pas modifiés par l’interaction avec la molécule CO. Finale-
ment, seule la bande p de Au est légèrement perturbée avec une augmentation des états
inoccupés à +1.04 eV et au dessus de +5.68 eV. Nous pensons que ces faibles change-
ments des PDOS des bandes de la surface Au sont corrélés à la faible adsorption (énergie
d’adsorption de 42 meV) calculée sur le site ponté avec un slab de 4 plans.
93
Figure 5.3: PDOS des bandes de Au libre (en pointillés) par rapport à celles dans CO/Au(001)
adsorbées (en trait continu) sur le site apical. Le niveau de Fermi est à 0 eV.
94
Figure 5.4: PDOS sur les orbitales de CO libre (en pointillés) par rapport à celles dans CO/Au(001)
adsorbées (en trait continu) pour le site ponté. Le niveau de Fermi est à 0 eV.
Sur la figure 5.5, nous avons représenté les PDOS de la bande s, p et d de Au avant
et après adsorption. Sous adsorption, nous remarquons tout de suite un fort changement
de la PDOS de la bande s de Au par rapport à celle de la phase libre (Figure 5.5a).
En effet, nous notons l’apparition des pics de résonance situés à -4.50 eV et -1.49 eV.
Ces nouveaux états, totalement absents dans la PDOS de la surface nue de Au, sont aux
mêmes endroits que les états moléculaires 4σ et 1π de CO sur la figure 5.4. De plus, nous
notons une diminution en intensité du pic localisé à +0.31 eV par rapport à celui sur la
phase libre ; celui-ci se décale de 0.16 eV vers les hautes énergies. Enfin, dans l’intervalle
[+6.19 eV ; +7.50 eV], nous notons aussi une augmentation des états inoccupés.
La figure 5.5b montre l’évolution de la PDOS de la bande p de Au avant et après
adsorption de CO. Tout d’abord, nous voyons l’apparition des mêmes pics de résonance
situés à -4.50 eV et -1.49 eV, exactement aux mêmes endroits que sur la PDOS de
CO/Au(001) (cf. figure 5.4). Ces nouveaux états sur la surface s’expliquent par les mé-
langes entre les orbitales moléculaires 4σ et 1π de CO et la bande p de Au. Enfin, de 0
eV à +7.50 eV, on note une augmentation des états inoccupés par rapport à ceux de la
surface nue.
Sur la figure 5.5c montrant les PDOS de la bande d, nous remarquons qu’après ad-
sorption, la bande d de Au est significativement modifiée. Nous notons aussi l’apparition
95
de ces deux pics de résonance localisés à -4.50 eV et -1.49 eV. Enfin, nous pouvons aussi
noter que le centre de gravité de la PDOS de la bande d de Au se déplace vers les basses
énergies. L’origine de ces changements de la PDOS de la bande d de Au sous adsorption
s’explique par le mélange entre les états moléculaires de CO et la bande d de Au.
Figure 5.5: PDOS des bandes de Au libre (en pointillés) par rapport à celles dans CO/Au(001)
adsorbées (en trait continu) sur le site ponté. Le niveau de Fermi est à 0 eV.
96
dans un domaine d’énergie [0 eV, +10 eV]. Notons une forte diminution de l’intensité
des orbitales moléculaires 4σ et 1π. Nous remarquons aussi que les états moléculaires 4σ
et 1π se décalent vers les basses énergies de 1.36 eV et 0.90 eV, respectivement. Tous ces
changements de la PDOS de CO/Au(001) traduisent une forte interaction entre états
moléculaires de CO et les bandes des atomes du slab d’or.
Figure 5.6: PDOS sur les orbitales de CO libre (en pointillés) par rapport à celles dans CO/Au(001)
adsorbées (en trait continu) sur le site de symétrie 4. Le niveau de Fermi est à 0 eV.
Nous avons par contre comparé sur la figure 5.7, la PDOS des bandes de Au dans
le système CO/Au avec celle de la surface nue. La figure 5.7a montre que la bande
s de Au dans CO/Au change énormément par rapport à celle sur la surface nue. Nous
constatons tout d’abord, l’apparition de deux pics de résonance situés à -4.08 eV et -1.26
eV, totalement absents sur la surface nue. Ces nouveaux états sont localisés aux mêmes
positions d’énergie que les états moléculaires 4σ et 1π de CO/Au(001) sur la figure 5.6.
Nous notons aussi une diminution de l’intensité du pic situé au niveau de Fermi. Enfin,
nous remarquons une augmentation des états inoccupés au-delà de +5.80 eV. Toutes ces
modifications de la PDOS de la bande s de Au s’expliquent par un fort mélange entre
cette bande et les états moléculaires de CO.
Sur la figure 5.7b représentant la PDOS de p de Au, nous remarquons, comme sur
la figure 5.7a, la présence de deux nouveaux états localisés à -4.08 eV et -1.26 eV,
exactement aux mêmes endroits que les états moléculaires 4σ et 1π de CO/Au(001) sur
97
la figure 5.6. Au-delà de +6.36 eV, nous notons une augmentation substantielle des états
inoccupés. Le même argument expliquant ces changements des PDOS reste valable.
La PDOS de la bande d de Au est représentée sur la figure 5.7c. Cette bande d est
la moins perturbée par l’adsorption de CO. Nous notons toutefois l’apparition de deux
pics situées à -4.08 eV et -1.26 eV, mais également un léger déplacement de son centre
de gravité vers les basses énergies. Tous ces changement traduisent l’interaction entre
les orbitales moléculaires de CO avec les bandes de la surface de Au.
Figure 5.7: PDOS des bandes de Au libre (en pointillés) par rapport à celles dans CO/Au(001)
adsorbées (en trait continu) sur le site de symétrie 4. Le niveau de Fermi est à 0 eV.
98
5.2 Adsorption du Carbone (C) sur la surface Au(001)
L’étude du carbone atomique (C) est nécessaire si on souhaite élucider les mécanismes
de décomposition du CO en C + O. En réalité, cette décomposition du CO est l’une des
réactions possibles dans la longue chaîne des réactions de décomposition du méthanol
[79].
5.2.1 Contexte expérimental et théorique
Le carbone atomique (C) est certainement l’un des systèmes les plus simples étudiés
en catalyse. L’adsorption de C sur les métaux de transition et noble a fait l’objet de
plusieurs études théoriques et expérimentales.
Sur le plan expérimental, les expériences de caractérisation du carbone adsorbé sur
les métaux sont très rares. Les quelques expériences trouvées dans la littérature ont été
menées en utilisant les techniques LEED [198-203], UPS [157], AES [199, 200, 203] et
STM [200, 203].
Sur le plan théorique, plusieurs investigations ab initio et DFT ont aussi été faites sur
l’adsorption de C sur les métaux de transition et nobles. Krekelberg et al. [72], utilisant
le code Dacapo, calculent, avec un slab de 4 plans, les énergies d’adsorption de C/Ir(111)
de -5.30 eV, -6.21 eV, -6.54 eV et -6.71 eV sur les sites top, bridge, hollow fcc et hollow
hcp, respectivement. De même, Ford et al. [98] rapportent pour C/Pt(111) des énergies
d’adsorption de -4.48 eV, -6.75 eV et -6.65 eV sur les sites top, hollow fcc et hollow hcp,
respectivement. D’autres valeurs ont été obtenues par Herron et al. [185] pour les sites
hollow fcc et hcp sur Pd (111) ; ils rapportent des énergies d’adsorption de -6.27 eV et
-6.30 eV, respectivement. Pour la surface Rh(111), Mavrikakis et ses collaborateurs [187]
ont publié des énergies d’adsorption de -5.10 eV, -6.51 eV, -6.78 eV et -7.11 eV pour les
quatre sites de haute symétrie top, bridge, hollow fcc et hollow hcp, respectivement.
D’autres études ont été menées sur Au(111) [19], Fe(110) [122][205-207], Re(0001) [190],
Ru(0001) [191] et Ni(111) [208-212].
Etant donné l’importance significative du catalyseur à base d’or dans les applica-
tions industrielles, une connaissance fondamentale et comportementale de l’adsorption
atomique du carbone sur la surface Au(001) est nécessaire ; ceci permet de renforcer le
champ de recherche fondamentale en sciences des surfaces.
5.2.2 Sites d’adsorption, Energies d’adsorption et paramètres structuraux
L’étude que nous avons menée sur le carbone atomique a porté sur les sites apical
(top), ponté (bridge) et de symétrie 4 (hollow). Sur chaque site, les énergies d’adsorption
99
ont été déterminées ainsi que les distance dC−Au . Cette étude a été faite sur un taux de
couverture de 0.25 sur un slab de 4 et 5 plans de Au. L’énergie d’adsorption a été calculée
seulement sur l’une des faces de deux surfaces exposées de chaque slab, en utilisant la
formule 5.1.
Le tableau 5.2 présente les résultats de notre étude.
Tableau 5.2: Energies d’adsorption (Eads ), longueurs de liaison entre C et Au (dC−Au ) obtenus avec
un slab de 4 et 5 plans après adsorption de C sur Au(001). Les résultats entre crochets sont calculés
avec un slab de 5 plans.
Dans le tableau 5.2, on note que l’augmentation du nombre de plans dans le slab a
une influence considérable sur l’énergie d’adsorption. Quand on passe de 4 à 5 plans,
l’énergie d’adsorption augmente de 0.548 eV, 0.602 eV et 0.595 eV pour les sites apical,
ponté et de symétrie 4, respectivement. Le site préférentiel d’adsorption de C est le
site de symétrie 4 avec une énergie d’adsorption de -5.773 [-6.368] eV (le nombre entre
crochets représente la valeur pour un slab de 5 plans). Nous avons aussi trouvé deux
autres sites métastables : apical et ponté ; cependant les énergies d’adsorption sur ces
sites sont faibles (-3.479 [-4.027] eV et -4.820 [-5.422] eV, respectivement). Krekelberg et
al. [72], Ford et al. [98] et Mavrikakis et al. [187] ont étudié l’adsorption de C sur Ir(111),
Pt(111) et Rh(111) sur le site top ; ils rapportent des énergies d’adsorption respectives
de -5.30 eV, -4.48 eV et -5.10 eV. La plus faible valeur (-2.12 eV) a été obtenue par
Santiango-Rodriguez et al. [19] sur le site top de Au(111).
Pour un site donné, on note que l’augmentation du nombre plans du slab n’influence
pas la distance, dC−Au ; cependant elle augmente avec la coordination. Sur le site préfé-
rentiel d’adsorption (site de symétrie 4), l’atome de carbone est à une distance de 2.164
Å au-dessus de la surface Au.
5.2.3 Structure électronique
Nous analysons des densités d’état projetées (PDOS) de la surface Au(001) et du
carbone avant et après adsorption. Nous examinons successivement les sites apical, ponté
100
et de symétrie 4 traités dans cette sous-section.
5.2.3.1 Site apical
Nous avons tracé sur la figure 5.8, la PDOS du carbone libre et celle du carbone dans
le système C/Au. La PDOS du carbone libre est caractérisée par deux pics d’intensité
localisés à -2.92 eV et +5.69 eV. Le pic à -2.92 eV est un état occupé à caractère s ;
tandis que celui à +5.63 eV est un état vide à caractère p, c.-à-d., une combinaison de
(px+py+pz).
Après adsorption, le pic d’intensité de l’état occupé situé à -2.92 eV se déplace de
2.57 eV vers les basses énergies par rapport à celui sur le carbone libre. D’autre part, on
note une très forte diminution en intensité du pic situé à +5.69 eV (avant adsorption) et
qui, se transforme presque en « continuum » dans un domaine d’énergie [-1.52 eV, +6.74
eV]. Ces changements s’expliquent par l’interaction entre les orbitales de C et celles des
atomes de Au, comme nous allons le voir dans ce qui suit.
Figure 5.8: PDOS sur les orbitales du carbone (C) libre (en pointillés) par rapport à celles dans
C/Au(001) adsorbées (en trait continu) pour le site apical. Le niveau de Fermi est à 0 eV.
101
Sur la figure 5.9, nous avons représenté les PDOS des bandes s, p et d de la surface
de Au libre et adsorbé. La PDOS de la bande s de la surface nue est cantonnée dans
un domaine d’énergie [-3.38 eV, +6.82 eV] avec un pic d’intensité localisé à -0.18 eV
(cf. Figure 5.9a) et une densité d’états de 0.11 électrons au niveau de Fermi. Après
adsorption, nous observons sur la figure 5.9a, une diminution drastique de la PDOS de
la bande s au niveau de Fermi (0.021 électrons) et l’apparition d’un pic de résonance
(nouvel état), totalement absent sur la surface nue, et localisé à -10.79 eV. Ce pic de
résonance est positionné à une distance en énergie environ 2 fois plus grande que la
position du pic sur le carbone adsorbé (Figure 5.8). Il est fort possible que ce nouvel
état provienne du mélange des orbitales du carbone et de la bande s de Au. Nous notons
enfin un déplacement de la bande s vers les basses énergies par rapport à la phase
non-adsorbée.
Sur la bande de p (figure 5.9b), nous remarquons la même situation que sur la figure
5.9a, avec l’apparition du pic de résonance situé à -10.79 eV. Nous pouvons aussi signaler
une légère augmentation de la densité des états au niveau de Fermi qui passe de 0.012
électrons (phase libre) à 0.015 électrons (phase adsorbée). Notons enfin que cette bande
p de Au après adsorption se déplace vers les basses énergies.
La bande d (figure 5.9c) quant à elle, fait apparaître aussi un nouvel état plus intense
à -10.79 eV, une largeur de la bande de ∼ 4 eV et une densité importante des états (2.53
électrons) au niveau de Fermi. Notons enfin que son centre de gravité se déplace de 0.46
eV vers les basses énergies par rapport à la phase adsorbée.
Tous ces changements observés sur les PDOS des bandes s, p et d de Au, notamment
l’apparition d’un nouvel état localisé à ∼ -10.79 eV et le déplacement des PDOS vers les
basses énergies sont clairement la conséquence de la stabilisation de ces bandes ; ceci est
due à la chimisorption qui provoque une passivation chimique de la surface Au.
102
Figure 5.9: PDOS des bandes de Au libre (en pointillés) par rapport à celles dans C/Au(001) adsorbées
(en trait continu) sur le site apical. Le niveau de Fermi est à 0 eV.
103
Figure 5.10: PDOS sur les orbitales C du carbone libre (en pointillés) par rapport à celles dans
C/Au(001) adsorbées (en trait continu) pour le site ponté. Le niveau de Fermi est à 0 eV.
Les variations de la PDOS des bandes s, p et d des atomes de surface Au(001) sont
présentées dans la figure 5.11. Sur toutes les courbes, on remarque l’apparition d’un pic
de résonance localisé à -5.88 eV. Ce nouvel état totalement absent sur la surface nue
se justifie par les mélanges des orbitales de C avec les bandes s, p et d des atomes de
surface Au(001). En outre, on note que la densité d’états des bandes p (Figure 5.11b)
au dessus du niveau de Fermi, dans le cas du carbone adsorbé, est supérieure à celle de
la surface seule : cela implique une perte de charge globale des bandes p pour les Au de
la surface avec l’adsorbat. Enfin, l’autre modification importante des bandes d (Figure
5.11c) avec l’adsorption du carbone est liée au déplacement du centre de gravité vers les
basses énergies.
Toutes ces modifications sont une conséquence de l’interaction due à un fort mélange
entre les orbitales du carbone et des bandes s, p et d de la surface de Au.
104
Figure 5.11: PDOS des bandes de Au libre (en pointillés) par rapport à celles dans C/Au(001)
adsorbées (en trait continu) sur le site ponté. Le niveau de Fermi est à 0 eV.
105
situé à -2.92 eV.
Figure 5.12: PDOS sur les orbitales du carbone (C) libre (en pointillés) par rapport à celles dans
C/Au(001) adsorbées (en trait continu) pour le site de symétrie 4 . Le niveau de Fermi est à 0 eV.
La figure 5.13 représente les PDOS des bandes s, p et d des atomes de Au de surface.
Comme nous l’avons vu dans les cas précédents, on note l’apparition d’un pic de réso-
nance situé à -6.49 eV. Ce pic à -6.49 eV correspond à un pic de PDOS sur le carbone
adsorbé (Figure 5.12), donc on peut l’attribuer à l’interaction carbone-surface Au. En
outre, on note au niveau de Fermi, une diminution de la densité d’états de la bande s.
Enfin, le centre de gravité des bandes d se décale légèrement vers les basses énergies.
Pour ce dernier cas, les mouvements des bandes d vers les basses énergies au cours de
l’adsorption du carbone sur les sites apical et ponté dominent celui de l’adsorption de C
sur le site de symétrie 4.
106
Figure 5.13: PDOS des bandes de Au libre (en pointillés) par rapport à celles dans C/Au(001)
adsorbées (en trait continu) sur le site de symétrie 4. Le niveau de Fermi est à 0 eV.
5.3 Conclusion
Dans ce chapitre, nous avons fait l’étude de l’adsorption du monoxyde de carbone
(CO) et du carbone (C) sur les sites apical, ponté et de symétrie 4 de la surface (001)
de Au pour un taux de couverture 0.25 selon des slabs sur 4 et 5 plans.
Selon la définition que nous avons choisie pour le calcul des énergies d’adsorption,
nous avons trouvé que celles-ci sont négatives sur tous les sites d’adsorption ; ce qui
signifie qu’il y’a bien adsorption des deux adsorbats (CO et C) sur la surface Au(001).
Nous avons remarqué que lorsque nous passons de 4 à 5 plans, l’énergie d’adsorption
augmente avec le nombre de plans.
L’adsorption de CO a révélé que le site préférentiel d’adsorption est le site ponté. Les
107
énergies d’adsorption et les paramètres structuraux ont été confrontés avec les résultats
expérimentaux et théoriques d’autres métaux de transition et noble. L’analyse des PDOS
a révélé que les états moléculaires 5σ et 2π de CO sont les plus perturbés par l’adsorption.
Quant à l’adsorption du carbone (C), nous avons trouvé qu’elle est plus favorable sur
le site de symétrie 4. L’analyse des PDOS des bandes s, p et d des atomes de Au a révélé
l’apparition des nouveaux états et le déplacement des PDOS vers les basses énergies ;
ceci est une stabilisation de ces bandes qui résulte de la chimisorption provoquant ainsi
une passivation chimique de la surface Au.
108
Conclusion générale et perspectives
Au cours de ce travail de thèse, nous avons étudié les propriétés énergétiques, élec-
troniques et structurales des espèces atomique (C), moléculaire (CO) et des fragments
moléculaires (CH3 , OH, CH2 OH, H2 CO et HCO) sur la surface Au(001). Nous avons
mené une étude théorique de toutes ces espèces avec les objectifs suivants. D’une part,
l’étude des fragments moléculaires hydrocarbonés se justifie parce que ceux-ci sont, en
perspectives, de bons candidats à l’extraction de l’hydrogène qui peut être utilisé comme
une nouvelle source d’énergie propre. D’autre part, il fallait étudier les espèces atomique
(C) et moléculaire (CO) dans le but de comprendre, en perspectives, la réaction de
désoxygénation du CO qui est un gaz potentiellement à effet de serre.
Nous avons utilisé comme substrat la surface de Au. En effet, les travaux théoriques
et expérimentaux récents portant sur l’adsorption de molécules légères sur la surface de
Au ont montré les possibilités catalytiques de ce métal noble. C’est ainsi que notre choix
s’est porté sur la face (001) de Au qui n’a fait l’objet d’aucune étude.
Nous avons commencé par les études préliminaires en déterminant le paramètre de
maille théorique de Au, en construisant les slabs de 4 et 5 plans de Au et en déterminant
les paramètres géométriques de ces adsorbats en phase gazeuse. A ce stade, nos résultats
structuraux sont en bon accord avec ceux expérimentaux et théoriques trouvés dans la
littérature.
Nous avons fait ensuite une étude systématique de l’adsorption de ces adsorbats sur
la surface Au(001) en considérant plusieurs configurations géométriques sur les trois sites
de haute symétrie (apical, ponté et de symétrie 4) de la surface Au(001).
Dans le chapitre 3, nous avons étudié l’adsorption de CH3 et OH sur Au(001) avec
et sans spin polarisé, pour des couvertures de 0.25, 0.50, 0.75 et 1.00 MC. Nous avons
constaté que les calculs avec la polarisation du spin réduisent considérablement les éner-
gies d’adsorption.
Dans le cas de CH3 /Au(001), nous avons constaté que CH3 est instable dans les sites
d’adsorption ponté et de symétrie 4. Concernant le système OH/Au(001), nous avons
constaté que le fragment OH s’adsorbe sur les trois sites (apical, ponté et de symétrie 4).
109
La préférence de stabilité énergétique suit cet ordre (excepté pour le taux de couverture
de 1.00) : ponté > symétrie 4 > apical. Les énergies d’adsorption dépendent du taux de
couverture et de l’augmentation du nombre de plans. Sur CH3 la variation de la work-
function ∆Φ diminue lorsque le taux de couverture augmente ; par contre, le moment
dipolaire µ augmente avec l’augmentation du taux de couverture. Sur OH, Nous avons
constaté des signes contraires de la variation de la work-function ∆Φ et du moment
dipolaire µ. Ceci a été expliqué en termes de redistribution de charges. L’analyse des
PDOS a révélé un mélange entre les orbitales de CH3 mais aussi de OH et les bandes s,
p et d de Au.
Enfin l’étude de la thermochimie par décomposition directe de CH3 et OH sur la
surface de Au(001) a démontré que la dissociation de ces deux adsorbats est défavorable
sur la surface Au(001).
Dans le chapitre 4, nous avons étudié l’adsorption de CH2 OH, H2 CO et HCO sur
Au(001) pour une couverture de 0.25 MC.
Pour l’hydroxyméthyle (CH2 OH), nous avons constaté que l’adsorption initiale via
l’oxygène n’est pas favorable car elle se termine par l’hydrogène. Sur l’adsorption du
formaldéhyde (H2 CO), nous avons noté une désorption en énergie avec un slab de 4 plans
et une chimisorption avec un slab de 5 plans dans le site apical. Nous avons noté aussi
que le site de symétrie 4 est instable pour les deux cas. L’adsorption du formyle (HCO) a
révélé que seul le site d’adsorption apical est stable. L’analyse des PDOS de l’hydrogène
et de l’oxygène (dans CH2 OH) a révélé que ces objets ont subi des modifications, celle de
C dans HCO révèle l’apparition des nouveaux états représentés par des pics de résonance
traduisant l’interaction entre ces adsorbats et les bandes s, p et d du substrat. Celle de
H2 CO via O avec la surface de Au ne perturbe pas les bandes s, p et d des atomes de
Au de surface.
Dans le dernier chapitre subdivisé en deux parties, nous avons étudié l’adsorption du
monoxyde de carbone (CO), d’une part, et du carbone atomique (C), d’autre part.
L’adsorption de CO a révélé que le site d’adsorption le plus favorable est le site ponté
et celui de l’adsorption de C, le site de symétrie 4. L’analyse des PDOS a révélé que les
états moléculaires 5σ et 2π de CO sont les plus perturbés par l’adsorption. Celle des
PDOS des orbitales du carbone et des bandes s, p et d des atomes de Au de surface,
met en évidence l’apparition de nouveaux états.
L’intérêt majeur de ces travaux a été de démontrer l’efficacité et la réactivité de la
surface Au(001) pour son utilisation comme catalyseur efficient. Ainsi, en perspectives,
il nous paraît important de poursuivre ce travail premièrement selon deux axes : (i)
110
l’étude des fragments moléculaires (CH3 , OH, CH2 OH, H2 CO et HCO) en présence d’un
adatome comme l’oxygène atomique, avec la polarisation de spin ; ce qui nous permettrait
de tester l’importance du magnétisme dans ces réactions catalytiques ; (ii) déterminer
avec et sans adatome les barrières d’activation dans les réactions de déshydrogénation
de ces fragments ; ce qui permettrait de conclure dans quel contexte le catalyseur à base
de Au est actif en déshydrogénation de ces fragments. Enfin, le poursuivre sur les taux
de couvertures 0.5, 0.75 et 1.00 pour les espèces CH2 OH, H2 CO, HCO, CO et C sans et
avec spin de polarisation puis, l’étendre sur les autres faces (111) et (110) de Au.
111
Bibliographie
[4] Haruta M., Kobayashi T., Sano H. and Yamada N. (1987), Novel gold catalysts
for the oxidation of carbon monoxide at a temperature far below 0o C, ”Chemistry
Letters“, no 2, vol. 16, p. 405-408.
[5] Haruta M., Yamada N., Kobayashi T. and Lijima S. (1989), Gold catalysts pre-
pared by coprecipitation for low-temperature oxidation of hydrogen and of carbon
monoxide, ”Journal of Catalysis“, no 2, vol. 115, p. 301-309.
[8] Sinha A. K., Seelan S., Tsubota S. and Haruta M. (2004), A tree-dimensional me-
soporous titanosilicate support for gold nanoparticles : vapor-phase epoxidation of
propene with high conversion, “Angewandte Chemie international edition”, no 12,
vol. 43, p. 1546-1548.
112
[9] Hayashi T., Tanaka K. and Haruta M. (1998), Selective vapor-phase epoxidation of
propylene over Au/TiO2 catalysts in the presence of oxygen and hydrogen, ”Journal
of Catalysis”, no 2, vol. 178, p. 566-575.
[10] Carrettin S., McMorn P., Johnston P., Griffin K., Kiely C. J., Attard G. A. and Hut-
chings G. J. (2004), Oxydation of glycerol using supported gold catalysts, “Topics
in Catalysis“, no 1-4, vol. 27, p. 131-136.
[11] Porta F. and Prati L. (2004), Selective oxidation of glycerol to sodium glycerate with
gold-on-carbon catalyst : an insight into reaction selectivity, ”Journal of Catalysis”,
no 2, vol. 224, p. 397-403.
[12] Demirel-Gülen S., Lucas M. and Claus P. (2005), Liquid phase oxidation of glycerol
over carbon supported gold catalysts, “Catalysis Today”, no , vol. 102-103, p. 166-
172.
[13] Porta F., Prati L., Rossi M., Coluccia S. and Martra G. (2000), Metal sols as a
useful tool for heterogeneous gold catalyst preparation : reinvestigation of a liquid
phase oxidation, “Catalysis Today”, no 1-4, vol. 61, p. 165-172.
[14] Li G., Edwards J., Carley A. F. and Hutchings G. J. (2006), Direct synthesis of
hydrogen peroxide from H2 and O2 using zeolite-supported Au catalysts, “Catalysis
Today”, no 4, vol. 114, p. 369-371.
[15] Landon P., Collier P. J., Papworth A. J., Kiely C. J. and Hutchings G. J. (2002),
Direct formation of hydrogen peroxide from H2 /O2 using a gold catalyst, “Chemical
Communications“, no 18, vol. 0, p. 2058-2059.
[16] Edwards J., Solsona B., Landon P., Carley A., Herzing A., Kiely C. and Hut-
chings G. (2005), Direct synthesis of hydrogen peroxide from H2 and O2 using
TiO2 -supported Au-Pd catalysts, ”Journal of Catalysis”, no 1, vol. 236, p. 69-79.
[17] Choudhary T. V., Sivadinarayana C., Datye A. K., Kumar D. and Goodman D. W.
(2003), Acetylene hydrogenation on Au-based catalysts, “Catalysis Letters”, no 1-3,
vol. 86, p. 1-8.
[18] Salama T. M., Ohnishi R., Shido T. and Ichikawa M. (1996), highly selective cataly-
tic reduction of NO by H2 over Auo and Au(I) impregnated in nay zeolite catalysts,
”Journal of Catalysis”, no 2 , vol. 162, p. 169-178.
113
[19] Santiago-Rodriguez Y., Herron J. A., Curet-Arana M. C. and Mavrikakis M. (2014),
Atomic and molecular adsorption on Au(111), “Surface Science”, no , vol. 627 , p.
57-69.
[20] Havu P., Blum V., Havu V., Rinke P. and Scheffler M. (2010), Large-scale surface
reconstruction energetics of Pt(100) and Au(100) by all-electron density functional
theory, “Physical Review B”, no 16, vol. 82, p. 161418 (1-4).
[21] Pieczyrak B., Trembulowicz A., Antczak G. and Jurczyszyn L. (2017), Nature of
monovacancies on quasi-hexagonal structure of reconstructed Au(100) surface, “Ap-
plied Surface Science“, no 1, vol. 407, p. 345-352.
[22] Hammer R., Sander A., Förster S., Kiel M., Meinel K. and Widdra W. (2014),
surface reconstruction of Au(001) : high-resolution real-space and reciprocal-space
inspection, “Physical Review B”, no 3, vol. 90, p. 035446 (1-13).
[28] Kittel C. (2005), Introduction to solid state physics, 8th edition : Hoboken, J. Wiley,
p. 680.
114
[30] Fermi E. (1927), Un Metodo Statistico per la determinazione di alcune proprietà
dell atomo “Rendiconti Accademia dei Lincei “, no 6, vol. 6, p. 602-607.
[31] Hohenberg P. and Kohn W. (1964), Inhomogeneous Electron Gas, ”Physical Review
B“, no 3B, vol. 136, p. B864-B871.
[32] Kohn W. and Sham L. J. (1965), Self-Consistent Equations Including Exchange and
Correlation Effects, ”Physical Review A“, no 4A, vol. 140, p. A1133-A1138.
[35] Vosko S. H., Wilk L. and Nursair M. (1980), Accurate spin-dependent electron
liquid correlation energies for local spin density calculations : a critical analysis,
”Canadian Journal of Physical”, no 8, vol. 58, p. 1200-1211.
[37] Trouiller N. and Martins J. L. (1991), Efficient pseudopotential for plane-wave cal-
culation, ”Physical Review B“, no 3, vol. 43, p. 1993-2006.
[42] Monkhorst H. J. and Pack J. D. (1976), Special points for Brillouin-zone integra-
tions, ”Physical Review B“, no 12, vol. 13, p. 5188-5192.
115
[43] Monkhorst H. J. and Pack J. D. (1978), Special points for Brillouin-zone integra-
tions, ”Physical Review B“, no 10, vol. 18, p. 5897-5899.
[44] Perdew J. P. and Wang Y. (1992), Accurate and simple analytic representation of
the electron-gas correlation energy, ”Physical Review B“, no 23, vol. 45, p. 13244-
13249.
[45] Baraldi A., Lizzit S., Comelli G., Kiskinova M., Rosei R., Honkala K. and Nørkov
J. K. (2004), spectroscopic Link between Adsorption Site Occupation and Local
Surface Chemical Reactivity, ”Physical Review and Letter“, no 4, vol. 93, p. 046101
(1-4).
[46] Janthon P., Kozlov S. M., Vines F., Limtrakul J. and Illas F. (2013), Establishing
the Accuracy of Broadly Used Density Functionals in Describing Bulk Properties
of Transition Metals, ”Journal of Chemical Theory and Computation“, no 3 , vol.
9, p. 1631-1640.
[47] Xu B., Haubrich J., Baker T. A., Kaxiras E. and Friend C. M. (2011), Theoreti-
cal Study of O-Assisted Selective Coupling of Methanol on Au(111), ”Journal of
Physical Chemistry C, no 5, vol. 115, p. 3703-3708.
[48] Mehmood F., Kara A., Rahman T. S. and Henry C. R. (2009), Comparative study
of CO adsorption on flat, stepped, and kinked Au surfaces using density functional
theory, “Physical Review B”, no 7, vol.79, p. 075422 (1-6).
[49] Gómez E. V., Amaya-Roncancio S., Avalle L. B., Linares D. H. and Gimenez M.
C. (2017), DFT Study of Adsorption and Diffusion of Atomic Hydrogen on Metal
Surfaces, ”Applied Surface Science“, no , vol. 420, p. 1-8.
[50] Shi H. and Stampfl C. (2007), First principles investigation of the structure and
stability of oxygen adsorption and surface oxide formation at Au(111), ”Physical
Review B.“, no 7, vol. 76, p. 075327 (1-14).
[51] Ladmann M., Rauls E and Schmidt W. G. (2009), First principles calculations of
clean Au(110) surfaces and chimisorptions of atomic oxygen, ”Physical Review B.“,
no 4, vol. 79, p. 045412 (1-9).
[52] Daigle A. D. and Belbruno J. J. (2011), Density functional theory study of the ad-
sorption of oxygen atoms on gold (111), (100) and (211) surfaces, ”Surface Science”,
no 13-14, vol. 605, p. 1313-1319.
116
[53] Tyson W. R. and Miller W. A. (1977), Surface free energies of solid metals : Esti-
mation from liquid surface tension measuremente, ”Surface Science“, no 1, vol. 62,
p. 267-276 .
[54] Zólyomi V., Vitos L., Kwon S. K. and Kollár J. (2009), Surface relaxation and stress
for 5d transition metals, “Journal of Physics Condensed Matter”, no 9, vol. 21, p.
095007 (1-5).
[55] Wang J. and Wang S.-Q. (2014), Surface energy and work function of fcc and bcc
crystals : Density functional study, “Surface Science”, no , vol. 630, p. 216-224.
[56] Vitos L., Ruban A. V., Skriver H. L. and Kollár J. (1998), The surface energy of
metals, “Surface Science“, no 1-2, vol. 411, p. 186-202.
[57] Hansson G. V. and Flodström S. A. (1978), Photoemission study of the bulk and
surface electronic structure of single crystal of gold, “Physical Review B condensed
matter”, no 4 , vol. 18, p. 1572-1585.
[58] Haroun M. F., Moussounda P. S., Légaré P. and Parlebas J. C. (2010), Adsorp-
tion and co-adsorption of CH3 and H on flat and defective ,nickel (111) surface,
“European Physical Journal B”, no 3, vol. 7, p. 353-358.
[59] Pang X. Y., Wang C., Yu-Huazhou, Zhao J-M. and Wang G.-C. (2010), DFT Study
of the Structure Sensitivity for the Adsorption of Methyl, Methoxy, and Formate on
Ni(111), Ni(100), and Ni(110) Surfaces, ”Journal of molecular structure : THEO-
CHEM“, no 1-3, vol. 948, p. 1-10.
[60] Steinbach F., Kiss J., Krall R. (1985), Identification and stability of CH3 , CH2 , and
CH species on Co and Ni surfaces, investigation, “Surface Science”, no 2-3, vol. 157,
p. 401-412.
[62] Ceyer S. T., Beckerle J. D., Lee M. B. , Tang S. L., Yang Q. Y., Hines M. A. (1987),
Effet of translational and vibrational energy on adsorption : The dynamics of mole-
cular and dissociative chemisorption, ”Journal of Vaccum Science and Technology
A : vaccum, Surfaces, and Films“, no 4, vol. 5, p. 501-507.
117
[63] Watson D. T. P., Titmuss S. and King D. A. (2002), Surface products of the disso-
ciative adsorption of methane on Pt(110)(2x2), ”Surface Science“, no , vol. 505 p.
49-57.
[65] Gavin R. M., Rentt J. and Muetterties E. L. (1986), Procedings of the National
Academy of Sciences of the United States of America, 78, 3981.
[66] Minot C., Van Hove M. A. and Somorjai G. A. (1983), A covalent model for the
bonding of adsorbed hydrocarbon fragments on the (111) face of platinum, ” Surface
Science“, no 3, vol. 127, p. 441-460.
[67] Petersen M. A., Jenkins S. J., King D. A (2004), Theory of Methane Dehydroge-
nation on Pt(110)(1x2). Part I : Chemisorption of CHX (X= 0-3), “ Journal of
Physical Chemistry B”, no 19, vol. 108, p. 5909-5919.
[68] Siegbahn P. E. M., Panas I., (1990), A theoretical study of CHX chemisorption on
the Ni(100) and Ni(111) surfaces, “Surface Science”, no 1-3, vol. 240, p. 37-49.
[69] Lai W., Xie D., Zhang D. H. (2005), First-Principles Study of Adsorption of Me-
thyl, Coadsorption of Methyl and hydrogen, and Methane dissociation on Ni(100),
“Surface Science”, no 1-3, vol. 594, p. 83-92.
[70] Ciobîcà I. M., Frechard F., Van Santen R. A., Kleyn A. W., Hafner J. (1999), A
theoretical study of CHX chemisorption on the Ru(0001) surface, ”Chemical Physics
and letters”, no 3-4, vol. 311, p. 185-192.
[71] Klinke D. J. II, Wilke S., Broadbelt L. J. (2002), A theoretical study of carbone
chemisorption on Ni(111) and Co(0001) surface “Journal of Catalysis”, no 2, vol.
178, p. 540-554.
[72] Krekelberg W. P., Greeley J. and Mavrikakis M. (2004), Atomic and Molecular
Adsorption on Ir(111), “Journal of Physical Chemistry B”, no 3, vol. 108, p. 987-
994.
118
[74] Liu H., Yan R., Zhang R., Wang B. and Xie K. (2011), A DFT theoretical study
of CH4 dissociation on gold-alloyed Ni(111) surface, ”Journal of Natural Gas Che-
mistry“, no 6, vol. 20, p. 611-617.
[75] Kratzer P., Hammer B. and Nørskov J. K. (1996), A theoretical study of CH4
dissociation on pure and gold-alloyed Ni(111) surfaces, ”the Journal of Chemical
Physics“, no 13, vol. 105, p. 5595-5604.
[76] Moussounda P. S., Haroun M. F., Rakotovelo G. and Légaré P. (2007), A Theoretical
Study of CH4 Dissociation on Pt(100), “Surface Science“, no 18, vol. 601, p. 3697-
3701.
[78] Lide D. R. (2009), ” Handbook of chemistry and Physics, 89th edition, CRC Press.
[79] Wang G.-C., Zhou Y.-H., Morikawa Y., Nakamura J., Cai Z-S and Zhao X.-Z.
(2005), Kinetic Mechanism of Methanol Decomposition on Ni(111) surface : A Theo-
retical Study, “Journal of Physical Chemistry”, no 25, vol.109, p. 12431-12442.
[80] Soon A., Todorova M., Delley B., and Stampfl C. (2006), Oxygen adsorption and
stability of surface oxides on Cu(111) : A first-principles investigation, “Physical
Review B”, no 16, vol. 73, p. 165424 (1-12).
[82] Bagus P. S. and Illas F. (1990), Theoretical analysis of the bonding of oxygen to
Cu(100), ”Physical Review B”, no 17, vol. 42, p. 10852-10857.
[83] Migani A., Sousa C. and Illas F. (2005), chemisorption of atomic chlorine on me-
tal surface and the interpretation of the induced work function change, “Surface
Science”, no 2-3, vol. 574, p. 297-305.
119
[85] Lin X. and Groβ A. (2012), First-principles study of the water structure on flat
and stepped gold surfaces, “Surface Science”, no 11-12, vol. 606, p. 886-891.
[87] Ljunström S., Hall J., Kasemo B., Rosèn A. and Wahnström T. (1987), A compa-
rative study of OH radical desorption in the H2 + O2 reaction on Pt, Pd, Rh, Ir
and Ni, ”Journal of Catalysis“, no 2, vol. 107, p. 548-556.
[89] Fujimoto G. T., Selwyn G. S., Kelser J. T., Lin M. C. (1983), Temperature effect
on the removal of hydroxyl radicals by a polycrystalline platinum surface, ”Journal
of Physical Chemistry“, no 11, vol. 87, p. 1906-1910.
[90] Bange K., Grider D. E., Madey T. E. and Sass J. K. (1984), The surface chemistry
of H2 O on clean and oxygen-covered Cu(110), ”Surface Science“, no 1, vol. 137, p.
38-64.
[91] Aderson A. B. (1981), Reactions and structures of water on clean and oxygen co-
vered Pt(111) and Fe, ”Surface Science“, no 1, vol. 105 p. 159-176.
[93] Yang H. and Whitten J. L. (1989), The adsorption of water and hydroxyl on Ni(111),
”Surface Science“, no 1-2, vol. 223 p. 131-150.
[95] Carter E. A. and Goddard III W. A. (1989), Chemisorption of oxygen, chlorine, hy-
drogen, hydroxide, and ethylene on silver cluster : A model for the olefin epoxidation
reaction, ”Surface Science“, no 1-2, vol. 209, p. 243-289.
120
[96] Hu Z.-M. and Nakatsuji H. (1999), Adsorption and disproportionation reaction of
OH on Ag surfaces : dipped adcluster model study, “Surface Science“, no 2-3, vol.
425, p. 296-312.
[97] Todorova M., Reuter K. and Scheffler M. (2004), Oxygen overlayer on Pd(111)
studies by density functional theory, “Journal of Physical and Chemistry B”, no 38,
vol. 108, p. 14477-14483.
[98] Ford D. C., Xu Y. and Mavrikakis M. (2005), Atomic and molecular adsorption on
Pt(111), “Surface Science“, no 3, vol. 587, p. 159-174.
[99] Pessoa A. M., Fajín J. L. C., Gomes J. R. B. and Cordeiro M. N. D. S. (2012), Ionic
and radical adsorption on the Au(hkl) surfaces : A DFT study, “ Surface Science”,
no 1-2, vol. 606, p. 69-77.
[100] Zhou Y.-H., Hong Lv P. and Wang G.-C. (2006), DFT Studies of methanol decom-
position on Ni(100) Surface : Compared with Ni(111) Surface, “Journal of Molecular
Catalysis A : Chemical”, no 1-2, vol. 258, p. 203-2015.
[101] Moussounda P. S., N’dollo M., Dintzer T., Mayéla Nkouka J. C. and Garin F.
(2013), Adsorption of NO on Pt(100) from first principles, “Physica Scripta”, no 5,
vol. 87, p. 055601 (1-7).
[102] Zeng Z.-H, Da Silva J. L. F., Deng H.-Q. and Li W.-X., (2009), Density functional
theory study of the energetics, electronic structure, and core-level shifts of NO
adsorption on the Pt(111) surface, “Physical Review B”, no 20, vol. 79, p. 205413
(1-13).
[103] Bouka-Pivoteau G. B., N’dollo M., Malonda-Boungou B. R., Maloumbi B., Mous-
sounda P. S. and Dintzer T. (2018), DFT Study of Methyl (CH3 ) and Hydroxyl
(OH) adsorption on a gold (001) surface, “Surface Review and Letters”, c (World
Scientific Publishing Company DOI : 10.1142/S0218625X18501986), p. 1850198 (1-
16)
[104] Chen W. K., Liu S. H., Cao M. J., Yan Q. G. and Lu C. H. (2006), Adsorption and
Dissociation of Methanol on Au(111) Surface : A First-Principles Periodic Density
Functional Study, “Journal of Molecular Structure (THEOCHEM), no 1-3, vol. 770,
p. 87-91.
121
[105] Greeley J. and Mavrikakis M. (2004), Competitive paths for methanol decompo-
sition on Pt(111), “Journal of the American Chemical Society”, no 12, vol. 126, p.
3910-3918.
[106] Kua J. and Goddard III W. A. (1999), Oxidation of Methanol on 2nd and 3rd Row
Group VIII Transition Metals (Pt, Ir, Os, Pd, Rh, and Ru) : Application to Direct
Methanol Fuel Cells, “Journal of the American Chemical Society”, no 47, vol. 121,
p. 10928-10941.
[107] Wang Z., Kan E.-J. and Yang J.-L. (2012), First-Principles Study of Adsorption
and Dissociation of Methanol on the Pt(100) Surface, “Chinese Journal of Chemical
Physics”, no 2, vol. 25, p. 199-203.
[109] Mai C.-F., Shue C.-H., Yang Y.-C., Yang L.-Y. O., Yau S.-L. and Itaya K. (2005),
Adsorption of formaldehyde on Pt(111) and Pt(100) Electrodes : cyclic voltammetry
and scanning tunneling microscopy, ”Langmuir“, no 11, vol. 21, p. 4964-4970.
[111] Raskó J., Kecskés T. and Kiss J. 2004, Adsorption and reaction of formaldehyde on
TiO2- supported Rh catalysts studied by FTIR and mass spectrometry, “ Journal
of Catalysis“, no 1, vol. 226, p. 183-191.
[112] Liu Z., Sawada T., Takagi N., Watanabe K. and Matsumoto Y. (2003), reaction
intermediates in the oxidation of methanol on Pt(111)-(2x2) O surface, ”Journal of
Chemical Physics“, no 9, vol. 119, p. 4879-4886.
122
studies of Pt- Zn/X-zeolite for isobutane dehydrogenation, ”Journal of Catalysis and
Letters, no 1-2, vol. 74, p. 17-25.
[115] Chi M. and Zhao Y. P. (2009), Adsorption of formaldehyde molecule on the in-
trinsic and Al-doped graphene :A first principle study, “Computational Materials
Science“, no 4, vol. 46, p. 1085-1090.
[116] Wang L. and Chen W. K. (2010), Density functional theory for adsorption of
HCHO on the FeO(100) surface, “Journal of Natural Gas Chemestry”, no 1, vol.
19, p. 21-24.
[117] Cahyanto W. T., Shukri G., Agusta M. K. and Kasai H. (2013), Adsorption of
formaldehyde and formyl intermediates on Pt, PtRu-, and PtRuMo-alloy surfaces :
A density functional study, ”Applied Surface Science”, no , vol. 266, p. 405-409.
[118] Yates J. T., Jr., Goodman D. W. , Madey T. E. (1977), Proceeding of the 7th
International Vaccum Congress and 3rd International Conference On Solid Surfaces,
Vienna, p. 1133.
[119] Baudais F. L. , Borschke A. J., Fedyk J. D., Dignam M. J. (1980), The decompo-
sition of methanol on Ni(100), “Surface Science”, no 1, vol. 100, p. 210-224.
[120] Mitchell W. J., Xie J., Wang Y., Weinberg W. H. (1993), Carbone monoxide
hydrogenation on the Ru(001) surface at low temperature using gas-phase atomic
hydrogen 1, “Journal of Electron Spectroscopy and Related Phenomena”, no , vol.
64, p. 427-433.
[122] Xu L., Kirvassilis D., Bai Y. and Mavrikakis M. (2017), Atomic and Molecular
Adsorption on Fe(110), ” Surface Science“, vol. 667, p. 54-65.
[123] Morgan Jr G. A., Sorescu D. C., Zubkov T., Yates Jr J. T. (2004), The forma-
tion and stability of adsorbed formyl as a possible intermediate in Fisher-Tropsch
chemistry on Ruthenium, “Journal of Physical Chemistry B”, no 11, vol. 108, p.
3614-3624.
123
[124] Gomes J. R. B. and Gomes J. A. N. F. (2000), Adsorption of the formyl species on
transition metal surfaces, “Journal of Electroanalytical Chemistry no 1-2, vol. 483,
p. 180-187.
[125] Takenaka S., Shimizu T. and Otsuka K. (2004), Complete removal of carbon mo-
noxyde in hydrogen-rich gas stream through methanation over supported metal
catalysts ” International Journal of Hydrogen Energy“, no 10, vol. 29, p. 1065-1073.
[127] Solymosi F., Erdöhelyi A. and Kocsis M. (1981), Methanation of CO2 on supported
Ru catalysts, ”Journal of the Chemical Society, Faraday Transactions I“, no 5, vol.
77, p. 1003-1012.
[129] Ekerdt J. G. and Bell A. T. (1980), Evidence for intermediates involved in Fischer-
Tropsch synthesis over Ru, ”Journal of Catalysis“, no 1, vol. 62, p. 19-25.
[130] Alayoglu S., Nilekar A. U., Mavrikakis M. and Eichhorn B. (2008), Ru-Pt core-shell
nanoparticules for preferential oxydation of carbon monoxide in hydrogen, ”Nature
Materials“, vol. 7, p. 333-338.
[131] Park E. D., Lee D., Lee H. C. (2009), Recent progress in selective CO removal in
a H2 -rich stream, ”Catalysis Today“, no 4, vol. 139, p. 280-290.
[132] Yeo Y. Y., Vattuone L. and King D. A. (1997), calorimetric heats for CO and
oxygen adsorption and for the catalytic CO oxidation reaction on Pt(111) , ”Journal
of Chemical Physics“, no 1, vol. 106, p. 392-401.
[133] Ertl G., Neumann M. and Streit K. M. (1977), Chemisorption of CO on the Pt(111)
surface, ”Surface Science“, no 2, vol. 64, p. 393-410.
[134] Wheeler M. C., Reeves C. T., Seets D. C. and Mullins C. B. (1998), Experimental
study of CO oxidation by an atomic oxygen beam on Pt(111), Ir(111), and Ru(001),
”Journal of Chemical Physics“, no 7, vol. 108, p. 3057-3063.
124
[135] Steininger H., Lehwald S. and Ibach H. (1982), On the adsorption of CO on
Pt(111), ”Surface Science“, no 2-3, vol. 123, p. 264-282.
[136] Persson B. N. J., Tüshaus M. and Bradshaw A. M. (1990), On the nature of dense
CO adlayers, ”Journal of Chemical Physics“, no 8, vol. 92, p. 5034-5046.
[137] Kung K. Y., Chen P., Wei F., Shen Y. R. and Somorjai G. A. (2000), Sum-
frequency generation spectroscopic study of CO adsorption and dissociation on
Pt(111) at high pressure and temperature ”Surface Science“, no 2, vol. 463 , p.
L627-L633.
[139] Podkolzin S. G., Shen J. Y., de Pablo J. J. and Dumesic J. A. (2000), Equilibrated
adsorption of CO on silica-supported Pt catalysts, ”Journal of Physical Chemistry
B“, no 17, vol. 104, p. 4169-4180.
[140] Guo X. C. and Yates J. T. (1989), Dependence of effective desorption kinetic para-
meters on surface coverage and adsorption temperature : CO on Pd(111), ”Journal
of Chemical Physics“, no 11, vol. 90, p. 6761-6766.
[141] Morkel M., Rupprechter G. and Freund H. J. (2003), Ultrahigh vacuum and high-
pressure coadsorption of CO and H2 on Pd(111) : A combined SFG, TDS, and
LEED study, ”Journal of Chemical Physics“, no 20, vol. 119, p. 10853-10866.
[143] Morkel M., Unterhalt H., Salmeron M., Rupprechter G., Freund H.-J. (2003), SFG
spectroscopy from 10−8 to 1000 mbar : less-ordered CO structures and coadsorption
on Pd(111), ”Surface Science“, vol. 532-535, p. 103-107.
125
[144] Bradshaw A. M. and Hoffmann F. M. (1978), The chemisorption of carbon mo-
noxide on palladium single crystal surfaces : IR spectroscopic evidence for localised
site adsorption ”Surface Science“, no 3, vol. 72, p. 513-535.
[145] Sautet P., Rose M. K., Dunphy J. C., Behler S. and Salmeron M. (2000), Adsorp-
tion and energetics of isolated CO molecules on Pd(111), ”Surface Science“, no 1-3,
vol. 453, p. 25-31.
[146] Rose M. K., Mitsui T., Dunphy J., Borg A., Ogletree D. F., Salmeron M. and
Sautet P. (2002), Ordered structures of CO on Pd(111) studied by STM, ”Surface
Science“, no 1-2, vol. 512, p. 48-60.
[147] Mitsui T., Rose M. K., Fomin E., Ogletree D. F. and Salmeron M. (2005), Diffusion
and pair interactions of CO molecules on Pd(111), ”Physical Review Letters“, no 3,
vol. 94, p. 036101 (1-4).
[148] Dubois L. H., and Somorjai G. A. (1980), The chemisorption of CO and CO2
on Rh(111) studied by high resolution electron energy loss spectroscopy, ”Surface
Science, no 2-3, vol. 91, p. 514-532.
[149] Beutler A., Lundgren E., Nyholm R., Andersen J. N., Setlik B. and Heskett D.
(1997), On the adsorption sites for CO on the Rh(111) single crystal surface, “Sur-
face Science”, no 2-3, vol. 371, p. 381-389.
[150] Beutler A., Lundgren E., Nyholm R., Andersen J. N., Setlik B. and Heskett D.
(1998), Coverage and temperature-dependent site occupancy of carbon monoxyde
on Rh(111) studied by high-resolution corelevel photoemission, “Surface Science”,
no 1-3, vol. 396, p. 117-136.
[151] Lerotholi J. T., Held G., and King D. A. (2007), phase mixing and phase separation
accompanying the catalytic oxydation of CO on Ir(100), “Surface Science”, no 5,
vol. 601, p. 1285-1295.
[152] Küppers J., Plagge A. J., (1976), Interaction of CO and O2 with Ir(111) surface,
“Journal of Vaccum Science and Technology“, no 1, vol. 13, p. 259-263.
[153] Lauterbach J., Boyle R. W., Schick M., Mitchell W. J., Meng B. and Weinberg W.
H. (1996), The adsorption of CO on Ir(111) investigated with, FT-IRAS, ”Surface
Science“, no 1-3, vol. 350, p. 32-44.
126
[154] Comrie C. M. and Weinberg W. H. J., (1976), The chemisorption of carbon mo-
noxyde on the Iridium III surface, ”the Journal Chemical Physics”, no 1, vol. 364,
p. 250-259.
[155] Zhdan P. A., Boreskov G. K., Baronin A. I., Egelhoff W. F. and Weinberg W.
H. (1976), Au XPS and UPS investigation of the chemisorption of CO on Ir(111),
“Chemical of Physics Letters”, no 3, vol. 44, p. 528-532.
[157] Brodén G., Gafner G. and Bonzel H. P. (1977), A UPS and LEED/Auger study
of adsorbates on Fe(110), “Applied Physics“, no 3, vol. 13, P. 333-342.
[158] Gafner G. and Feder R. (1976), A low-energy electron diffraction study of Fe(110) :
I. Experi- mental, ”Surface Science“, no 1, vol. 57, p. 37-44.
[161] Feder R. and Gafner G. (1976), A low-energy electron diffraction study of Fe(110) :
II. theory and comparison with the experiment, “Surface Science”, no 1, vol. 57, p.
45-52.
[166] Wedler G. and Ruhmann H. (1982), Laser induced thermal desorption of carbon
monoxide from Fe(110) surfaces, “Surface Science”, no 3, vol. 121, p. 464-486.
[163] Gonzalez L., Miranda R. and Ferrer S. (1982), A thermal desorption study of
the adsorption of CO on Fe(110) ; enhancement of dissociation by surface defects,
“Surface Science”, no 1, vol. 119, p. 61-70.
[164] Maruyama T., Sakisaka Y., Kato H., Aiura Y. and Yanashima H. (1994), CO on
Fe(110) : an angle-resolved photoemission study, “Surface Science”, no 3, vol. 304,
p. 281-290.
127
[165] Ertl G., Küppers J., Nitschke F. and Weiss M. (1977), Photoelectron spectroscopic
study of surface complexes on Fe(110), “Chemical Physical and Letters”, no 2, vol.
52, p. 309-313.
[166] Wedler G. and Ruhmann R. (1983), Study of the adsorption and decomposition
of CO on Fe(110) with an improved vibrating capacitor, “Applied Surface Science”,
no 2, vol. 14, p. 137-148.
[169] Lemire C., Meyer R., Shaikhutdinov S. K. and Freund H.-J., (2004), CO adsorp-
tion on oxide supported gold : from small clusters to monolayer islands and three-
dimensional nanoparticles, “Surface Science”, no 1-3, vol. 552, p. 27-34.
[170] Piccolo L., Loffreda D., Cadete Santos Aires F. J, Deranlot C., Jugnet Y., Sautet P.
and Bertolini J. C. (2004), The adsorption of CO on Au(111) at elevated pressures
studied by STM, RAIRS and DFT calculations, “Surface Science”, no (Part 2), vol.
566-568, p. 995-1000.
[171] Pfnür H., Feulner P. and Menzel D. (1983), The influence of adsorbate interac-
tions of kinetics and equilibrium for CO on Ru(001). II. Desorption kinetics and
equilibrium, “The Journal of Chemical Physics”, no 9, vol. 79, p. 4613-4623.
[172] Hoffmann F. M., Weisel M. D. and Peden C. H. F. (1991), In-situ FT-IRAS study
of the CO oxidation reaction over Ru(001) : II. coadsorption of carbone monoxide
and oxygen, “Surface Science”, no 1-3, vol. 253, p. 59-71.
[173] Thomas G. E. and Weinberg W. H. (1979), The vibrational spectrum and adsorp-
tion site of CO on the Ru(001) surface, “The Journal of Chemical Physics”, no 3,
vol. 70, p. 1437-1439.
[174] Pfnür H., Menzel D., Hoffmann F. M., Ortega A. and Bradshaw A. M. (1980),
High resolution vibrational spectroscopy of CO on Ru(001) : The importance of
lateral interactions, “Surface Science”, no 2-3, vol. 93, p. 431-452.
128
[175] Michalk G., Moritz W., Pfnür H. and Menzel D. (1983), A LEED determination
√ √
of the structures of Ru(001) and of CORu(001)- 3X 3 R30o , “Surface Science”,
no 1, vol. 129, p. 92-106.
[177] Bonn M., Hess C., Funk S., Miners J. H., Persson B. N. J., Wolf M., Ertl G.
(2000), femtosecong surface vibrational spectroscopy of CO adsorbed on Ru(001)
during desorption, “Physical Review Letters”, no 20, vol. 84, p. 4653-4656.
[178] Fuggle J. C., Madey T. E., Steinkilberg M. and Menzel D. (1975), Photoelectron
spectroscopic studies of adsorption of CO and oxygen on Ru(001), “Surface Science”,
no 3, vol. 52, p. 521-541.
[179] Fuggle J. C., Umbach E., Feulner P. and Menzel D. (1977), Electron induced
dissociation of CO on Ru(001) : A study by thermal desorption and electron spec-
troscopies, “Surface Science”, no 1, vol. 64, p. 69-84.
[180] Ciobica I. M., Kleyn A. W. and Van Santen R. A. (2003), Adsorption and coad-
sorption of CO and H on Ruthenium surface, “ The Journal of Physical Chemistry
B”, no 1, vol. 107, p. 164-172.
[181] Gonzalez S. and Illas F. (2005), CO adsorption on monometallic Pd, Rh, Cu and
bimetallic PdCu and RhCu monolayers supported on Ru(0001), “ Surface Science”,
no 1-3, vol. 598, p. 144-155.
[182] Eichler A. and Hafner J. (1999), Reaction channels for the catalytic oxidation of
CO on Pt(111), “ Physical Review B”, no 8, vol. 59, p. 5960-5967.
[184] Kandoi S., Gokhale A. A., Grabow L. C., Dumesic J. A. and Mavrikakis M. (2004),
Why Au and Cu are more selective than Pt for preferential oxidation of CO at low
temperature, ”Catalysis Letters”, no 1-2, vol. 93, p. 93-100.
[185] Herron J. A., Tonelli S. and Mavrikakis M. (2012), Atomic and molecular adsorp-
tion on Pd(111), “Surface Science”, no 21-22, vol. 606, p. 1670-1679.
129
[186] Zhang C. J., Hu P. and Lee M. H. (1999), A density functional theory study on
the interaction between chemisorbed CO and S on Rh(111), “Surface Science”, no
3, vol. 432, p. 305-315.
[187] Mavrikakis M., Rempel J., Greeley J., Hansen L. B. and NørskovJ. K. (2002),
Atomic and molecular adsorption on Rh(111), “Journal of Chemical Physics” no
14, vol. 117, p. 6737-6744.
[188] Mehmood F., Kara A., rahman T. S. and Bohnen P. K. (2006), Energetics Of
CO on stepped an kinked Cu surfaces : A comparative theoretical study (2006),
“Physical Review B”, no 15, vol. 74, p. 155439 (1-9).
[189] Shah V., Li T., Baumert K. L., Cheng H. and Sholl D. S. (2003), A comparative
study of CO chemisorption on flat and stepped Ni surfaces using density functional
theory, “Surface Science“, no 3, vol. 537, p. 217-227.
[191] Herron J. A., Tonelli S. and Mavrikakis M. (2013), Atomic and molecular adsorp-
tion on Ru(0001), “Surface Science”, vol. 614, p. 64-74.
[193] Kresse G. and Furthmüller J. (1996), Efficient iterative schemes for ab initio total-
energy calculations using a plane-wave basis set, “Physical Review B”, no 16, vol.
54, p. 11169-11186.
[194] Kresse G. and Furthmüller J. (1996), Efficient of ab-initio total energy calcula-
tions for metals and semiconductors using a plane-wave basis set, “Computational
Materials Science”, no 1, vol. 6, p. 15-50.
[196] Yim W.-L., Nowitzki T., Necke M., Schnars H., Nickut P., Biener J., Biener M. M.,
Zielasek V., Al-Shamery K., Klüner T. and Bäumer M. (2007), Universal phenomena
130
of CO adsorption on gold surfaces with low-coordinated sites, “Journal of Physical
Chemistry C”, no 1, vol. 111, p. 445-451.
[197] Kelemen S. R., Kaldor A. (1981), The interaction of surface sulfur with carbon on
Fe(110), “Journal of Chemical Physics”, no 3, vol. 75, p. 1530-1537.
[198] Buckley D. H. (1972), Adhesion of metals to a clean iron surface studied with
LEED and auger emission spectroscopy, ”Wear”, no 1, vol. 20, p. 89-103.
[200] Klink C., Stensgaard I., Besenbacher F. and Laegsgaard E. (1995), An STM study
of carbon induced structures on Ni(111)-Evidence for a carbidic-phase clock recons-
truction, ”Surface Science“, no 1-3, vol. 342, p. 250-260.
[201] Gardin D. E., Batteas J. D., Vanhove M. A., Somorjai G. A. (1993), Carbon,
nitrogen, and sulfur on Ni(111)-Formation of complex structures and consequences
for molecular decomposition, ”Surface Science“, no 1, vol. 296, p. 25-35.
[202] Mccarrol J., Edmonds T., Pitkethl R. (1969), Interpretation of a complex low
energy electron diffraction pattern : Carbonaceous and sulphur-containing struc-
tures on Ni(111), ”Nature”, no 5212, vol. 223, p. 1260-1262.
[204] Arabczyk W. and Rausche E. (1996), Model cluster studies of the restructuring
of the Fe(110) surface induced by the adsorption of carbon, nitrogen or oxygen,
”Polonica Journal of Chemical”, no 7, vol. 70, p. 910-974.
[205] Jiang D. E., Carter E. A. (2005), Carbon atom adsorption on and diffusion into
Fe(110) and Fe(100) from first principles, “Physical Review B Condensed Matter
and Materials Physics” no 4, vol. 71, p. 045402 (1-6).
131
[207] Watwe R. M., Bengaard H. S., Rostrup-Nielsen J. R., Dumesic J. A. and Nørskov J.
K. (2000), Theoretical studies of stability and reactivity of CHx species on Ni(111),
”Journal of Catalysis”, no 1, vol. 189, p. 16-30.
[208] Zhu Y. A., Dai Y. C., Chen D. and Yuan W. K. (2007), First-principles study of C
chemisorption and diffusion on the surface and in the subsurfaces of Ni(111) during
the growth of carbon nanofibers, “Surface Science”, no 5, vol. 601, p. 1319-1325.
[210] Wang S. G., Liao X. Y., Cao D. B., Li Y. W., Wang J. G. and Jiao H. J. (2007),
Formation of carbon species on Ni(111) : Structure and stability, ”Journal of Phy-
sical Chemistry C“, no 29, vol. 111, p. 10894-10903.
[211] Bai Y., Kirvassilis D., Xu L. and Mavrikakis M. (2018), Atomic and Molecular
Adsorption on Ni(111), ”Surface Science“, no , vol. 679 , p. 240-253.
132
Publication
133
Surface Review and Letters, 1850198 (16 pages)
°c World Scienti¯c Publishing Company
DOI: 10.1142/S0218625X18501986
§
mbblaise2002@yahoo.fr
¶paulsand_moussounda@yahoo.com
by WSPC on 05/30/18. For personal use only.
T. DINTZER
Institut de Chimie et Procedes pour l'Energie,
l'Environnement et la Sante (ICPEES),
Ecole europeenne de Chimie, Polymeres et Materiaux,
25 rue Becquerel, 67087 Strasbourg Cedex 2, France
thierry.dintzer@unistra.fr
Received 4 March 2017
Revised 17 March 2018
Accepted 13 April 2018
Published 18 May 2018
We used the DACAPO code with the GGA-PW91 approximation to study the adsorption of methyl
(CH3) and hydroxyl (OH) for four- and ¯ve-layer gold (Au) (001) slabs. We have determined for
each species the best binding site, adsorption energy, the change in the work function, surface
energy, surface dipole moment, geometrical parameters and projected density of states (PDOS).
We performed spin-unpolarized and spin-polarized DFT calculations for free and adsorbed CH3 and
OH species. The most important point is that the spin polarization diminishes the adsorption
energies but does not change the geometrical parameters. For the CH3 species, only the top site was
found to be stable for di®erent coverages. We found that during the optimization phase, the hollow
and bridge sites were found to be unstable. In both cases the CH3 species moves toward the top site.
We observe that the adsorption energy decreases when increasing the coverage. However, the OH
species was stable in all investigated sites (top, bridge and hollow). We notice that the adsorption
energy is dependent on the number of slab layers and the bridge is the best site in adsorption energy.
The analysis of the calculated O PDOS of OH radical shows a mixing between the O orbitals and the
Au bands.
†Corresponding author.
1850198-1
G. B. Bouka-Pivoteau et al.
hydrochlorination of acetylene.11 Thus, over the last authors have investigated the Au(001) surface and
few decades, investigations into the catalysis of me- ascertained that it reconstructs into ð5 1Þ,68,69
by WSPC on 05/30/18. For personal use only.
tallic gold have been intensi¯ed. ð5 20Þ70–72 and cð26 68Þ73 surface cells. Grønlund
The CH3 adsorption on the surface of the transi- and Højlund Nielson74 have demonstrated using
tion metals (TMs) has been intensively studied both HREELS that the Au(001) surface reconstructs into
experimentally and theoretically. Experimentally, two ð5 20Þ orthogonal domains, while Reider and
several techniques such as secondary ion mass spec- coworkers75 using He di®raction techniques con-
troscopy (SIMS) and high-resolution electron energy ¯rmed that it reconstructs into ð5 1Þ surface cells.
loss spectroscopy (HREELS) on Ni(111) surface,12–16 Using ultra-high vacuum transmission electron mi-
temperature-programmed reaction spectroscopy croscopy (UHVTEM) Yamazaki and Takayanagi76
(TPRS)17 and re°exion–adsorption infrared spec- found a ð28 5Þ reconstructed cell for the Au(001)
troscopy (RAIRS) on Pt(111) surfaces18 were used. surface. More recently, Havu et al.77 and Pieczyrak
Theoretically, some studies concerning the adsorption et al.78 reported that a short-range atoms arrange-
of methyl species on transition metal surfaces and ment of the quasihexagonal \hex" surface recon-
clusters have been reported. A DFT method was used struction of the Au(001) surface was described by a
by Liu et al.,19 Kratzer et al.20 and Siegbahn and ð5 NÞ unit cell, where N 20–30, i.e. > 100 atoms
Panas21 to study the adsorption of CH3 species on are to be considered in each layer. Hammer and
Au/Ni(111), Ni(111) and Ni(100) surfaces, respec- coworkers79 reported that the long-range arrange-
tively. Many additional theoretical studies of methyl ment of the Au(001) surface can be described by the
adsorption on di®erent metal surfaces have been cð28 48Þ superstructure. It has been observed that
performed.22–35 the Au(001) surface catalytically reconstructs as
The adsorption of hydroxyl species was studied Pt(001).77,80,81 Nevertheless, surface adsorbates are
experimentally on Ni,36–39 Pd,39–41 Ru,42,43 Pt,44–49 known to lift reconstructions of the (001) surface.
Rh,39,50 Cu51–53 and Ag.53–59 Theoretically, OH was These take account of CO,82,83 NO,84 ethane,85
investigated on Pt(111),60 Fe(100),60 Ni(100),61,62 Ni hydrocarbon83 sulfur and thiolates.86–88 Hence, when
(111),62–64 Cu(111),65 Ag.66,67 the Au(001) surface is in interaction with its envi-
The Au(001) surface is of intermediate atomic ronment, it should correspond to a stable con¯gura-
density between Au(111) and Au(110) surfaces, i.e. tion. So it makes sense to investigate the adsorbate
the Au(001) surface is more open than the dense (CH3 or OH) adsorption on the unreconstructed Au
surfaces of Au(111) and Au(110), which are thermo- (001) surface on account of the di±culties associated
dynamically most stable. Yet, to the best of our with the description of two types of structures on
1850198-2
DFT Study of CH3 and OH Adsorption on an Au (001) Surface
the reconstructed Au(001) surface within a periodic corresponding to an electronic temperature of 0.1 eV.
supercell approach.89 All energies were extrapolated to zero temperature.
In our precedent works,8,9 we have reported the The adsorbates (CH3 and OH) and the topmost of the
DFT calculations of the adsorption of CH3OH mole- metal-slab were relaxed. We used the Broyden–
cule, CH3 O radical and atomic hydrogen (H), and the Fletcher–Goldfard–Shanno (BFGS) method to opti-
decomposition of CH3OH with CH3O and H on the mize all geometries. In our works, the calculations
(001) surface of gold. First, the goal of this paper is were stopped when the force was less than 0.05 eV/ A.
to study the adsorption of CH3 and OH species on an In our previous works, the optimized lattice-bulk of
Au(001) surface to determine the surface energy, Au was of 4.18 A.8,9 Our result is 10% higher
adsorption energy and electronic structure [work than the experimental value (4.08 A).93 However, our
function and projected density of states (PDOS)]. result is in agreement with the values of 4.17 A and
Second, our study is motivated by the following most
4.19 A calculated by Xu et al. 114
and Mehmood
important point: due to the signi¯cant relevance of et al.,94 respectively. The ð2 2Þ unit cell of Au(001),
gold catalysts in industrial applications, fundamental occupied by one, two and three adsorbates, was used
understanding of the adsorption behavior of CH3 and to obtain the coverages of 0.25, 0.50 and 0.75 ML,
OH over the gold surface can be a high topic for many respectively. The ð1 1Þ unit cell, containing one
investigations in surface science. In addition, in par- adsorbate (CH3 or OH), was used for a coverage of
Surf. Rev. Lett. Downloaded from www.worldscientific.com
cal oxidation of glucose and carbon oxidation. dence of the CH3 or OH adsorption energy on the
Consequently, due to the importance of OH and CH3 number of layers. Our study showed that if the
species, with respect to the catalytic activity of gold- number of Au layers is more than ¯ve the adsorption
based electrodes, it seems crucial to understand energy is stabilized. We have reported the adsorption
their adsorption properties. Third, the objective of energies for CH3 and OH using four and ¯ve Au layers
this study is to give a database of the benchmarks to show that the di®erences are signi¯cant. But, the
that can be accessible to many researchers of gold di®erences using 5–7 Au layers were very weak
catalysis. ( 0:03–0.05 eV).
The adsorption energy ðEads ) of adsorbate was
calculated using the equation
2. Theoretical Methods 1
Eads ¼ ½E Eslab NEadsorbate ; ð1Þ
90
In this work we used DACAPO package to simulate N system
the interaction between two adsorbates (CH3 and
where N is the total number of adsorbates (CH3 or
OH) and the surface of Au(001). The basis of plane
OH) in the surface unit cell at the considered cover-
waves was used to describe the Kohn–Sham equations
age, Esystem is the total energy of the Au(001) slab
using the generalized gradient approximation (GGA)
with the adsorbate, Eadsorbate the energy of the ad-
of Perdew and Wang 91 (GGA-PW91)91 and the
sorbate in gas phase and Eslab the total energy of the
ultra-soft pseudo-potentials of Vanderbilt.92 The
free Au(001) surface. By this de¯nition, exothermic
kinetic energy (cuto®) was limited to 400 eV. We used
adsorption is indicated by a positive Eads value; i.e. a
a ð1 1Þ unit cell of four or ¯ve slabs of Au(001) and a
more positive Eads value is associated with enhanced
convergence in energy was obtained with the k-points
binding strength, i.e. is associated to the most stable
of ð12 12 1Þ. Afterwards, this convergence in en-
or favorable site.
ergy was tried for a ð2 2Þ unit cell using the k-points
The surface energy ðEsurf Þ was evaluated by using
of ð6 6 1Þ. The method of Monkhorst and Pack
the equation
was used to generate the grid of the k-points in the
Brillouin zone. The region of vacuum separating two 1
Esurf ¼ ½E NAu Ebulk ; ð2Þ
successive slabs was equivalent to six-layer slab in the 2A slab
z-direction. The numeric stability of the ground where Eslab , Ebulk , NAu and A represents the energy of
state was obtained using the Fermi distribution a slab calculation, energy of a single atom in bulk, the
1850198-3
G. B. Bouka-Pivoteau et al.
number of atoms in the supercell and the area per 3. Results and Discussion
ð1 1Þ or ð2 2Þ surface unit cell, respectively.
We ¯rst consider the clean metal surface. The surface
We have also examined our results by analyzing
energies Esurf calculated with the ð1 1Þ and ð2 2Þ
the work-function change calculated as
surface unit cells are 0.882 J/m2 and 0.881 J/m2,
¼ ðadsorbate Þ clean ; ð3Þ respectively. It reveals that the surface energy is not
very sensitive to the choice of surface unit cells, and
where ðadsorbate Þ and clean are the work functions of
the e®ects on the surface energies are also negligibly
the Au(001) slab with the adsorbate and free Au slab.
small. Our surface energy ( 0:88 J/m2) only
The Helmholtz equation95 was used to calculate
accounts for half of the surface energy determined
the surface dipole moment
from measurements (1.54 J/m2).99,100 At this point, it
A is very di±cult to decide whether the di®erences arise
¼ : ð4Þ
12 from de¯ciencies in the experimental method. We
To include the charge density di®erences ðrÞ, we compare also our results to the Esurf values calculated
performed all calculations as implemented in the for the Au(001) surface. Zólyomi et al.101 reported the
plane-wave self-consistent ¯eld (PWSCF) package of Esurf of 0.864 J/m2 through DFT-PBE-GGA, which is
the QUANTUM-ESPRESSO distribution.96 The somewhat less than ours. Wang and Wang102 and
ðrÞ was calculated by the usual relation
Surf. Rev. Lett. Downloaded from www.worldscientific.com
Table 1. Theoretical and experimental bond lengths, dCH for CH3 and
dOH for OH, and angles H–C–H for CH3.
dCH (
A) 1.09 1.08 1.085, b 1.09 c
H–C–H angle ( ) 120 120 120, b 120 c
dOH (
A) 0.993 0.971 0.972 b
1850198-4
DFT Study of CH3 and OH Adsorption on an Au (001) Surface
Furthermore, these properties are also calculated by species moves in both cases to the top site. This
using a ¯ve-layer Au slab and are provided in the text phenomenon has already been observed by Mous-
and the tables enclosed by square brackets next to the sounda et al.83 for CH3/Pt(100).
four-layer Au slab values. Before studying the CH3 adsorption, we have
calculated a total magnetic moment of 1:00 B , for
the free CH3. No experimental and theoretical data
3.1. CH3 adsorption on Au(001) are available for eventual comparison. Upon adsorp-
3.1.1. Adsorption energies and structural tion, the total magnetic moment of CH3 is approxi-
parameters mately zero ð þ25:10 5 B Þ. Consequently, the CH3
magnetization is reduced by the Au(001) surface
Figure 1 shows the adsorption site of CH3 on Au(001)
dramatically. We can conclude that the CH3 species is
surface. The top site is stable for di®erent coverages.
not coupled magnetically to the Au atoms because
We found that during the optimization phase, the the interactions between CH3 species and surface
hollow and bridge sites were unstable. The CH3
metal annihilate the net spin on each CH3 species.
The adsorption energies obtained by the spin-
unpolarized and spin-polarized DFT calculations are
H hollow B bridge and T top sites summarized in Table 2. This table reveals three most
Surf. Rev. Lett. Downloaded from www.worldscientific.com
Table 2. Adsorption energies [Eads ðspÞ and Eads ðunpÞ], work-function change
(), dipole moment (), bond length between C and H (dCH Þ, bond length
between C and Au (dCAu ), and angles (H–C–H) for CH3/Au(001) at the top site
with a four-layer slab. The calculated results with a ¯ve-layer slab are given in
square brackets. Note that the Eads ðspÞ [Eads ðunpÞ] refers to the adsorption energies
calculated using the spin-polarized (spin-unpolarized) DFT calculations.
1850198-5
G. B. Bouka-Pivoteau et al.
[1.574], 0.840 [1.011] and 0.237 [0.202] eV, respec- adsorption energies between 1.86 eV and 2.08 eV for
tively. We can conclude that the adsorption energy adsorption on top which are in agreement with our
decreases when the coverage increases. We think that value for four-layer slab. We note that Santiago-
the increased concentration of the adsorbates on the Rodriguez et al.105 also calculated an adsorption en-
Au(001) surface results in enhanced interaction be- ergy of 1.19 eV in their DFT calculations of CH3
tween the adsorbates and hence weakens their bond adsorbed only at the top site of Au(111). In summary,
with the substrate as shown in Table 2. In Sec. 3.1.2, our adsorption energies are in general higher than
we will discuss the surface electronic modi¯cation due those of CH3/TM [TM ¼ Ni(111) and Au(111)23,105].
to the dependence of the CH3 adsorption energies on Note that this is due to (i) the di®erence of local
the CH3 coverage. Moreover, as shown in Table 2, the electronic structures between these transition and
stability of CH3 is more favorable for 0.25 ML with a noble metal surfaces (Au and Ni) and (ii) the more
¯ve-layer Au(001) slab and 0.50 ML with a four-layer open (001) surface with lower coordination atoms
Au(001) slab. Third, in general, we notice that there which is characterized by a higher adsorption energy.
is a dependence of the adsorption energy on the The structural parameters of CH3/Au(001) are
number of layers. The variations of the adsorption summarized in Table 2. Figure 2 illustrates the cal-
energy with the number of layers in the ratio are Eads culation of these structural parameters. We notice
(¯ve-layers)/Eads (four-layersÞ ¼ 1:27, 1.14, 1.12 and that the increase in number of layers and the spin-
Surf. Rev. Lett. Downloaded from www.worldscientific.com
0.96 for 0.25, 0.50, 0.75 and 1.00 ML, respectively, for polarized DFT calculation have no in°uence on the
the spin-unpolarized DFT calculations; the corre- geometrical parameters. For example, the coverage
by WSPC on 05/30/18. For personal use only.
sponding values for the spin-polarized DFT calcula- does not a®ect the bond lengths of C–Au. For
tions are 1.40, 1.21, 1.20 and 0.85, respectively. So, we 0.25 ML coverage our results are 2.123 [2.121] A, in
see that, with the exception for the coverage of good agreement with the data of Moussounda et al.
1.00 ML, the adsorption energy increases when the (2.08
A)83 and Santiago-Rodriguez et al. (2.13 A)105
number of layers increases. In e®ect, as the number of for CH3 adsorbed on the top sites of Pt(100) and Au
atoms in the slab increases the interaction of the CH3 (111) surfaces, respectively, except that Lai et al.97
species with the slab increases. Many interpretations found a bond length between C and Au of 1.95 A in
of the e®ects of the number of layers are also possible. their adsorption of CH3 on Ni(100). The C–H dis-
First, one can consider that the mechanism such as tances vary between 1.073 A and 1.096 A which are
dipole interactions become dominant. Second, the
electronic structures of the clean surface and
the methyl species can be interpreted by the analysis
of the corresponding electronic structure of the
adsorbed CH3 species. No experimental and theoret-
ical studies of the CH3 adsorption on Au(001) have
been reported. However, many workers have exam-
ined the CH3 adsorption on several transition metals
such as Ni,23,97 Pt(100)98 and Au(111).105 Then, we
may compare brie°y our results with those reported
by these authors which have been performed without
spin-polarized calculations. That goes without say-
ing, we have the conscience that the adsorption en-
ergies could not be similar. Lai et al.97 and
Moussounda98 reported adsorption energies of
2.05 eV and 2.60 eV for CH3/Ni(100) and CH3/Pt
(100) at the top site for a coverage of 0.25 ML, re-
spectively. The di®erence between our results and
previous DFT studies97,98 is due to the di®erences
in computational conditions. Schule et al.,23 using the Fig. 2. (Color online) Side view for the optimized struc-
cluster models [for CH3/Ni(111)], published the ture of CH3/Au(001) at the top site for 0.25 ML coverage.
1850198-6
DFT Study of CH3 and OH Adsorption on an Au (001) Surface
not a®ected when the number of Au layers increases. 0:50 eV (for ¼ 0:50 ML) using four and ¯ve
We notice that the CH3 angles (H–C–H) are reduced Au layers might come from the same polarization
from 120:0 (for the free CH3) to 111:6 [for CH3/Au induced by adsorption. Second, the increase of the
(001)] and between 110:0 and 111:1 for other cov- number of layers in°uences weakly the work-function
erages. We conclude that the interaction between the change, excluding for the coverage of 0.25 ML.
methyl species and the gold surface causes a repulsion Consequently, when the coverage is higher than
of hydrogen atoms by the gold surface and that 0.25 ML, the electrostatic-potential change (’ ¼
consequently the methyl has changed from the sp 2 ’adsorbate=slab ’slab ) is the same for the two types of
hybridization state to the sp 3 hybridization in the slabs examined in this study, e.g. 31.5 mV (0.50 ML),
CH3 adsorption. 63.5 mV (0.75 ML) and 113.5 mV (1.00 ML). Third,
we have observed that the sign of is negative, i.e.
points upward to the vacuum. Fourth, the dipole
3.1.2. Electronic structures moment increases with the coverage. This is a direct
3.1.2.1. Work function and dipole moment consequence of the CH3-induced redistribution of the
electronic density. In addition in all cases experi-
Now, we discuss the electronic structure of the
mental and theoretical data are not available.
CH3–Au(001) system. Table 2 lists the work function
change and dipole moment . First, we note that
Surf. Rev. Lett. Downloaded from www.worldscientific.com
constant for both 4 and 5-Au layers. For metal sur- 0.25 ML CH3/Au minority states
0.50 ML CH3/Au majority states
faces, such behavior is possible but not obvious. In 4 0.50 ML CH3/Au minority states
0.75 ML CH3/Au majority states
this circumstance, the question about the relation of 0.75 ML CH3/Au minority states
PDOS (states/eV)
1850198-7
G. B. Bouka-Pivoteau et al.
located at 8:78, 2:11 and þ2:53 eV, while the 0.05 0.25 ML CH3/Au majority states
0.25 ML CH3/Au minority states
minority-spin states have also three energy peaks 0.50 ML CH3/Au majority states
0.5 ML CH3/Au minority states
0
positioned at 7:77, 1:91 and þ5:77 eV. Upon ad-
0.75 ML CH3/Au majority states
0.75 ML CH3/Au minority states
E - EF (eV)
positioned at 8:17 eV and 1:51 eV. The most im-
portant change of the PDOS is that both the unoc- Fig. 4. (Color online) Spin-PDOSs for the clean surface
cupied antibonding states located at þ2:53 eV (for and adsorbed surface at 0.25, 0.50 and 0.75 ML coverages,
the majority-spin states) and þ5:77 eV (for the mi- with majority (minority)-(minority)-s, p and d metal states.
nority-spin states) are mostly a®ected by the surface,
i.e. are transformed in continuum below the Fermi Second, we analyze now the PDOS without the
level. These results should be explained by a strong spin-polarized DFT calculation. Figure 5 shows the
hybridization between the CH3 states and the s, p and CH3 PDOSs for 0.25, 0.50 and 0.75 ML coverages.
d metal states. The CH3 PDOS for the isolated CH3 is included for
Now we examine the PDOS of the Au(001) surface. comparison. The Fermi level is positioned to 0 eV and
It is known that the Au metal is not magnetic. the adsorbed states are calibrated to their isolated
Consequently, the PDOS majority-spin states and states. For the isolated states of CH3, we observe
minority-spin states have the same evolution. Figure 4 three peaks located at 7:53, 1:27 and þ4:39 eV,
illustrates the evolution of the PDOS of majority
(minority)-spin s, p and d states. The PDOS of clean unpolarized spins EF
15
surface is included for comparison. For 0.25, 0.50 and free CH3 states
0.25 ML CH3/Au states
0.75 ML coverages, we can notice a discrepancy be- 0.50 ML CH3/Au states
1e
0.75 ML CH3/Au states
tween the clean- and adsorbed-surface Au s, p and d
states. For the adsorbed surface, we see a new state 10
PDOS (states/eV)
is hybridized with the CH3 orbitals. In contrast, the s Fig. 5. (Color online) Projected densities of states for the
and d metal states are weakly a®ected by the CH3 free CH3 and CH3/Au(001) at 0.25, 0.50 and 0.75 ML
adsorption. coverages, using the spin-unpolarized DFT calculation.
1850198-8
DFT Study of CH3 and OH Adsorption on an Au (001) Surface
and termed (2a1 ), (1e) and (3a1 ), respectively, using peak at 7:53 eV in the PDOSs of the s, p and d metal
Sch€on°ies notation. We can see that (2a1 ) and bands. This is a peak of resonance having exactly the
(1e) are occupied molecular states while (3a1 ) is the same position as (2a1 ) in the PDOS for CH3/Au(001)
unoccupied molecular state in the free CH3 species. (see Fig. 5) and it decreases with the increase of the
We can notice that the coverage in°uences the coverage. Therefore, there is a strong hybridization
PDOS. When the coverage increases the intensities of between s, p and d metal bands and (2a1 ) orbital.
the peaks decrease with the exception of the peak at Looking at the PDOS of Au-s band, we can see that
7:53 eV for 0.25 ML coverage. The most important the most intense peak is located at 1:07 eV. Its in-
point is that the (3a1 ) orbital disappears compared tensity weakly increases when the coverage increases.
to the isolated CH3; this orbital mixes strongly Other change in the PDOS of Au-s band is that in
with the s, p and d metal bands as emphasized by both intervals [0:26 eV, þ3:37 eV] and [þ4:38 eV,
Lai et al.97 in CH3/Ni(100) study. In reality, this þ5:39 eV], the number of states per energy increases
(3a1 ) orbital is transformed into continuum of states and decreases, respectively. In the case of PDOS of
in the interval [þ2:73 eV, þ6:00 eV]. All these Au-p band, there is an increase of the intensity of
observations are qualitatively consistent with those PDOS in the interval [2:48 eV, þ6:48 eV] accom-
obtained by Michaelides and Hu25 for CH3/Ni(111) panied by several peaks of resonance, which are
and Moussounda98 for CH3/Pt(100). interpreted by the hybridization between the CH3 (1e
Surf. Rev. Lett. Downloaded from www.worldscientific.com
Figure 6 shows the PDOSs of the s, p and d metal and 3a1 ) orbitals and Au-p band. To ¯nish, we can see
bands at the coverages of 0.25, 0.5 and 0.75 ML (the that the PDOSs of the Au-d band in the clean surface,
by WSPC on 05/30/18. For personal use only.
corresponding PDOSs in the clean surface are incor- the occupied states and the unoccupied states are not
porated for comparison). As shown in Fig. 6, on ad- equal. Upon adsorption, we notice that the PDOS
sorption of the methyl radical, there appears a new decreases when the coverage increases. Another in-
teresting e®ect is that the center of gravity of PDOS
EF
has moved by 0.30 eV towards the negative energies,
unpolarized spins
0.2 with respect to its localization at Au clean surface.
Au (s) PDOS (states/eV)
0.04
(b)
1850198-9
G. B. Bouka-Pivoteau et al.
Table 3. Adsorption energies [Eads ðspÞ and Eads ðunpÞ], work-function change (), dipole
moment (), bond length of OH (dOH ), bond length between O and Au (dOAu ) and angles
(H–O–Au) obtained for OH/Au(001) with a four-layer slab. The calculated results with a ¯ve-
layer slab are given in square brackets. Note that the Eads ðspÞ ½Eads ðunpÞ values refer to the
adsorption energies calculated using the spin-polarized (spin-unpolarized) DFT calculations.
between these two types of calculations is large: the the two types of calculations; the di®erences are 0.470
adsorption energies of the spin-polarized calculations [0.476], 0.641 [0.630] and 0.950 [0.947] eV, respec-
are about 0.392–0.791 eV lower than the unpolarized tively, for 0.25, 0.50 and 0.75 ML coverages. The
calculations. corresponding values between the bridge and inter-
The adsorption of OH on Au(001) is stable on mediate metastable (hollow) sites are 0.248 [0.263],
hollow, bridge and top sites contrary to the CH3 ad- 0.525 [0.516] and 0.163 [0.098] eV, respectively. We
sorption which is stable only in the top site (see can deduce that the number of layers in the slab
Sec. 3.1.1). For these types of calculations, our in°uences weakly the energy di®erence calculated
Table 3 reveals that the adsorption energies increase from the preferred site to the metastable site, with the
when the number of layers in the slab increases (ex- exception of 0.75 ML coverage for the latter case. To
cept for 1.00 ML coverage for the top and bridge ¯nish, the hydroxyl is most stable in the hollow site
sites), i.e. Eads (¯ve-layers) > Eads (four-layers). The for ¼ 1:00 ML with energy di®erences with respect
analysis of the spin-unpolarized (spin-polarized) cal- to the top and bridge sites of 0.673 [1.220] (0.668
culations shows that the bridge is the most stable for [0.822]) eV and 0.259 [0.729] (0.250 [0.367]) eV,
0.25, 0.50 and 0.75 ML coverages with binding ener- respectively (numbers in parentheses refer to the
gies of 2.681 [3.170] (2.289 [2.777]), 2.870 [3.134] calculated values for spin polarization).
(2.480 [2.742]) and 2.270 [2.450] (1.877 [2.057]) eV. Using our spin-unpolarized DFT calculation, we
The energy di®erence between the preferred (bridge) can compare our adsorption energies with similar
and lowest metastable (top) sites is independent of studies recently performed for other transition and
1850198-10
DFT Study of CH3 and OH Adsorption on an Au (001) Surface
noble metals, such as for Au(100),112 Au(111),105 Ni hollow sites are of 2.084, 2.219 and 2.503 A, respec-
(100),7 Ni(111),4 Pt(111) Ref. 34 and Ag(100).67 tively. On the bridge site, our results are similar to
Using a periodic-slab approach, Pessoa et al.112 those obtained by Hu and Nakatswi67 on the Ag(100)
obtained adsorption energies of 1.97, 3.07 and 3.11 eV surface. Calculated H–O–Au angle has mean values of
for the top, bridge and hollow sites, respectively. 105:9 , 104:9 and 122:2 for OH for the top, bridge
These values are quantitatively higher than our and hollow sites, respectively.
results. However, using a cluster model, Pessoa
et al.112 obtained adsorption energies of 2.53 eV and
2.77 eV for the bridge and hollow sites, respectively. 3.2.2. Electronic structures
Thus, our adsorption energies for the OH species on
3.2.2.1. Work function and dipole moment
the bridge and hollow sites are close to the theoretical
values reported by Pessoa et al.112 Santiago-Rodriguez In Table 3, the work-function change and dipole
et al.105 obtained a value of 1.70 eV corresponding to moment are listed. A positive is expected due to
the adsorption of OH/Au(111) at the bridge site. This the di®erence in electronegativities between hydrogen
smaller adsorption energy on Au(111) surface indi- (2.2), oxygen (3.44) and gold (2.54). In fact, this is
cates a weaker interaction between OH and Au(111). found for OH in the bridge (top) site at 0.25, 0.50 and
For OH/Ni(100) at 1/6 ML and for the bridge and 0.75 ML, for example, ¼ 0:050 (0.617 [0.430]),
Surf. Rev. Lett. Downloaded from www.worldscientific.com
hollow sites, Zhou et al.7 obtained adsorption energies 0.365 [0.335] (1.041 [1.041]) and 0.153 [0.114] eV, re-
of 3.093 eV and 3.367 eV, respectively, which are in spectively. But, this is not the case for OH in the
by WSPC on 05/30/18. For personal use only.
disagreement with our results. Wang and coworkers4 top site at 0.75 ML (1.0 ML) and the hollow site
in their works on the adsorption on OH/Ni(111) for 0.25, 0.50, 0.75 and 1.00 ML for which we
reported (at ¼ 1=6 ML) adsorption energies of obtained ¼ 0:412½0:418ð0:904½0:917Þ eV
3.071 eV and 3.185 eV for top and bridge sites, re- and 0:662½0:921, 1:032½1:035, 1:150½1:198
spectively. These results are higher than our theo- and 1:251½1:271 eV, respectively. Similar results
retical values. These discrepancies could be due to the have been obtained for NO/Pt(100)112 and NO/Pt
fact that, in their calculations, Wang et al.4 used the (111),113 where these authors112,113 obtained the op-
GGA-PBE (non-self-consistent) functional, three- posite sign of for the bridge and top sites as well
layer Ni(111) slab and a cuto® energy of 340 eV for as the hollow and top sites, respectively. It looks like
wave functions while in our calculations, we use the these ¯ndings come from the di®erent charge redis-
GGA-PW91 (self-consistent) functional, four-layer tributions induced by adsorption. In e®ect, in
(or ¯ve-layer) Au(001) slab and a cuto® of 400 eV. Figs. 7(a) and 7(b), we visualized the electron-density
For ¼ 0:25 ML, our OH adsorption energy at the redistributions localized principally in the topmost
top site (2.210 eV) for four-layer slab is consistent surface layer and adsorbed OH in the case of 0.25 ML.
reasonably with the result of Ford et al.34 at top site For OH in the hollow site [Fig. 7(b)], we noted a
(2.13 eV) for three-layer slab of Pt(111). To ¯nish, strong charge accumulation around OH and a de-
Hu,67 in their cluster model, reported adsorption en- pletion in the region of the four nearest Au neighbors.
ergies of 4.705, 4.432 and 4.211 eV for OH/Ag(100) in For the top site [Fig. 7(a)], there is also a charge
the hollow, bridge and top sites, respectively. These accumulation of electron density around the hydroxyl
values are about two times than our corresponding and a depletion of charge density for the ¯rst nearest
values. Au, whereas there is an accumulation and depletion
Table 3 presents also the O–H bond lengths, O–Au around the second nearest Au neighbor. Thus, these
distance and H–O–Au angle. The O–H bond lengths signi¯cant di®erences of the charge redistribution
are situated between 0.984 A and 0.993 A, and they between the hollow and top sites at 0.25 ML seem to
are in agreement with the values obtained (0.990 A) be responsible for the opposite sign of .
by Santiago-Rodriguez et al.105 and Ford et al.34 for From Table 3, we can be seen that the dipole
OH adsorbed on Au(111) and Pt(111), respectively. moment presents several behaviors. First, for the
For a given site, the O–Au distance is not dependent hollow site a negative increasing of the dipole moment
on the number of slab layers and the OH coverage. is observed with increasing coverage. Second, a posi-
The calculated O–Au distances in the top, bridge and tive is obtained for OH in the top (bridge) at
1850198-11
G. B. Bouka-Pivoteau et al.
EF polarized spins
free OH majority states
free OH minority states
4 0.25 ML OH/Au majority states
0.25 ML OH/Au minority states
0.50 ML OH/Au majority states
0.50 ML OH/Au minoriry states
PDOS (states/eV)
0
-2
(a) (b) -4
1850198-12
DFT Study of CH3 and OH Adsorption on an Au (001) Surface
clean Au states
0.25 ML OH/Au states
0.04 0.15 0.50 ML OH/Au states
0.02
0 0.1
-0.02
-0.04 0.05
-0.06
0
(b)
Au (p) PDOS (states/eV)
(b)
-0.02
0.25 ML OH/Au majority states 0.005
Au (d) PDOS (states/eV)
1
-2
-3 0 3 6 9
E - EF (eV) 0
by WSPC on 05/30/18. For personal use only.
-3 0 3 6 9
E - EF (eV)
Fig. 9. (Color online) Spin-PDOSs for the clean surface
and the adsorbed surface for 0.25 ML and 0.50 ML cov- Fig. 11. (Color online) Spin-PDOSs for the clean surface
erages, using the spin-polarized DFT calculation. and the adsorbed surface for 0.25 ML and 0.50 ML cov-
erages, using the spin-unpolarized DFT calculation.
OH. The hybridization between the s and p orbitals of
O and the s, p and d metal bands revealed that the
at 1:73 eV appears on the O PDOS, which is in
initial unoccupied orbitals are dispersed on a domain
the same energy position as the s, p and d metal bands
of energy ( 7 eV). Notice that one new state located
(see Fig. 11). This shows a strong mixing between O
orbitals and the s, p and d bands of gold. In addition,
EF the Au-p band is particularly most a®ected by
unpolarized spins
10 the adsorption, because the high PDOS value is ob-
free OH states
0.25 ML OH/Au states served between þ0:08 eV and þ4:33 eV; this signi¯es
0.50 ML OH/Au states
that the Au-p band hybridizes strongly with the O
8
orbitals.
PDOS (states/eV)
6 4. Conclusion
DFT calculations are used to investigate the ad-
4 sorption of two adsorbates (CH3 and OH) on Au(001)
surface at = 0.25, 0.50, 0.75 and 1.00 ML. For these
calculations, a ð2 2Þ or ð1 1Þ unit cell is used to
2 construct a four-layers (¯ve-layer) Au(001) slab. We
performed spin-unpolarized and spin-polarized DFT
calculations for free and adsorbed CH3 and OH spe-
0
-3 0 3 6 9 cies. The spin-polarized calculations signi¯cantly re-
E - EF (eV)
duce the adsorption energies.
Fig. 10. (Color online) Spin-PDOSs of O in OH/Au(001) For CH3 species, we found that during the opti-
on a bridge site for 0.25 ML and 0.50 ML coverages, using mization phase, the hollow and bridge sites were
the spin-unpolarized DFT calculation. found to be unstable. In both case the CH3 species
1850198-13
G. B. Bouka-Pivoteau et al.
moves towards the top site. The adsorption energies 6. H. A. C. M. Hendrickx, A. Hoek and B. E. Nieu-
calculated are dependent on the coverage and n-layer wenhuys, Surf. Sci. 135 (1983) 81.
7. Y.-H. Zhou, P.-H. Lv and G.-C. Wang, J. Mol. Catal.
Au(001) slab ðn ¼ 4; 5Þ. The work-function change
A, Chem. 258 (2006) 203.
, increases noticeably as the coverage is reduced. 8. M. N'dollo, P. S. Moussounda, T. Dintzer, B.
The changes of the PDOS revealed a strong mixing of M'Passi-Mabiala and F. Garin, Surf. Interface Anal.
CH3 orbital with s, p and d metal bands. 45 (2013) 1410.
For OH species, we notice that the number of 9. M. N'dollo, P. S. Moussounda, T. Dintzer and
layers in the slab in°uences the adsorption energy. F. Garin, J. Mod. Phys. 4 (2013) 409.
10. M. Haruta, N. Yamada, T. Kobayashi and S. Lijima,
The OH species can bind on either the hollow, top or
J. Catal. 115 (1989) 301.
bridge site of Au(001) with the preferred adsorption 11. G. J. Hutchings and M. Haruta, Appl. Catal. A 291
site in the order (except for the 1.00 ML coverage) (2005) 2.
bridge > hollow > top. The energy di®erence be- 12. M. P. Kaminsky, N. Winograd, G. L. Geo®roy and
tween the preferred (bridge) and lowest metastable M. A. Vannice, J. Am. Chem. Soc. 108 (1986) 1315.
13. M. B. Lee, Q. Y. Yang, S. L. Tang and S. T. Ceyer,
(top and hollow) sites was estimated. The site de-
J. Chem. Phys. 85 (1986) 1693.
pendence of the di®erent charge redistributions 14. S. T. Ceyer, J. D. Beckerle, M. B. Lee, S. L. Tang,
explains the opposite sign of and induced by Q. Y. Yang and M. A. Hines, J. Vac. Sci. Technol. A
OH adsorption in the bridge, hollow and top sites. 5 (1987) 501.
Surf. Rev. Lett. Downloaded from www.worldscientific.com
The analysis of the oxygen (in the OH species) pro- 15. Q. Y. Yang and S. T. Ceyer, J. Vac. Sci. Technol. A 6
jected density of states shows a strong mixing of the (1988) 851.
by WSPC on 05/30/18. For personal use only.
1850198-14
DFT Study of CH3 and OH Adsorption on an Au (001) Surface
35. G. Papoian, J. K. Nørskov and R. Ho®mann, J. Am. 65. K. Hermann, M. Witko, L. G. M. Pettersson and
Chem. Soc. 122 (2000) 4129. P. Siegbahn, J. Chem. Phys. 99 (1993) 610.
36. J. T. Keiser, M. A. Ho®bauer and M. C. Lin, J. Phys. 66. E. A. Carter and W. A. Goddard III, Surf. Sci. 209
Chem. 89 (1985) 2635. (1989) 243.
37. F. P. Netzer and T. E. Madey, Phys. Rev. Lett. 47 67. Z.-M. Hu and H. Nakatsuji, Surf. Sci. 425 (1999) 296.
(1981) 928. 68. D. G. Fedak and N. A. Gjostein, Acta Metall. 15
38. C. Benndorf, C. N€ obl and T. E. Madey, Surf. Sci. 138 (1967) 827.
(1984) 292. 69. D. G. Fedak and N. A. Gjostein, Phys. Rev. Lett. 16
39. C. Nyberg and P. Uvdal, J. Chem. Phys. 84 (1986) (1966) 171.
4631. 70. D. G. Fedak and N. A. Gjostein, Surf. Sci. 8 (1967)
40. L. C. Anderson, C. E. Mooney and J. H. Lunsford, 77.
Chem. Phys. Lett. 196 (1992) 445. 71. P. N. J. Dennis and P. J. Dobson, Surf. Sci. 33 (1972)
41. P. A. Thiel, F. M. Ho®mann and W. H. Weinberg, 187.
J. Chem. Phys. 75 (1981) 5556. 72. G. E. Rhead, J. Phys. F 3 (1973) L53.
42. K. Kretschmar, J. K. Sass, A. M. Bradshaw and 73. M. A. Van Hove, R. J. Koestner, P. C. Stair, J. P.
S. Holloway, Surf. Sci. Lett. 115 (1982) 183. Biberian, L. L. Kesmodel, I. Bartos and G. A.
43. G. T. Fujimoto, G. S. Selwyn, J. T. Keiser and M. C. Somorjai, Surf. Sci. 103 (1981) 189.
Lin, J. Phys. Chem. 87 (1983) 1906. 74. F. Grønlund and P. E. Hojlund Nielsen, J. Appl.
44. L. V. Novakoski and D. S. Y. Hsu, J. Chem. Phys. 92 Phys. 43 (1972) 391.
(1990) 1999. 75. K. H. Reider, T. Engel and R. H. Swendsen, Surf. Sci.
Surf. Rev. Lett. Downloaded from www.worldscientific.com
1850198-15
G. B. Bouka-Pivoteau et al.
93. D. R. Lide (ed.), Handbook of Chemistry and Physics 104. G. V. Hansson and S. A. Flodstr€om, Phys. Rev. B 18
(CRC Press, Boca Raton, 2005). (1978) 1572.
94. F. Mehmood, A. Kara, T. S. Rahman and C. R. 105. Y. Santiago-Rodriguez, J. A. Herron, M. C.
Henry, Phys. Rev. B 79 (2009) 075422. Curet-Arana and M. Mavrikakis, Surf. Sci. 627
95. M. Todorova, K. Reuter and M. Sche®ler, J. Phys. (2014) 57.
Chem. B 108 (2004) 14477. 106. L. G. M. Pettersson and P. S. Bagus, Phys. Rev. Lett.
96. Quantum Espresso Foundation, Quantum Espresso 56 (1986) 500.
(2017), http://www.quantum-espresso.org. 107. P. S. Bagus and F. Illas, Phys. Rev. B 42 (1990)
97. W. Lai, D. Xie and D. H. Zhang, Surf. Sci. 594 (2005) 10852.
83. 108. A. Migani, C. Sousa and F. Illas, Surf. Sci. 574 (2005)
98. P. S. Moussounda, Adsorption et activation du 297.
methane et du methanol sur la surface (100) du pla- 109. Š. Pick, J. Phys., Condens. Matter 22 (2010) 395002.
tine: Une etude par la fonctionnelle de la densite 110. Z.-H. Zeng, J. L. F. Da Silva, H.-Q. Deng and W.-X.
Ph.D. thesis, Universite Louis Pasteur, Strasbourg, Li, Phys. Rev. B 79 (2009) 205413.
France (2006). 111. S. Schnur and A. Groβ, New J. Phys. 11 (2009)
99. W. R. Tyson and W. A. Miller, Surf. Sci. 62 (1977) 125003.
267. 112. A. M. Pessoa, J. L. C. Fajín, J. R. B. Gomes and
100. H. B. Michaelson, J. Appl. Phys. 48 (1977) 4729. M. N. D. S. Cordeiro, Surf. Sci. 606 (2012) 69.
101. V. Zólyomi, L. Vitos, S. K. Kwon and J. Kollar, 113. P. S. Moussounda, M. N'dollo, T. Dintzer, J. C.
J. Phys., Condens. Matter 21 (2009) 095007. Mayela Nkouka and F. Garin, Phys. Scr. 87 (2013)
Surf. Rev. Lett. Downloaded from www.worldscientific.com
102. J. Wang and S.-Q. Wang, Surf. Sci. 630 (2014) 216. 055601.
103. L. Vitos, A. V. Ruban, H. L. Skriver and J. Kollar, 114. B. Xu, J. Haubrich, T. A. Baker, E. Kaxiras and
by WSPC on 05/30/18. For personal use only.
Surf. Sci. 411 (1998) 186. M. Friend, J. Phys. Chem. C 115 (2011) 3703.
1850198-16
Résumé
La magnitude des énergies d’adsorption (en valeur absolue) pour toutes les espèces analysées peut être
ordonnée de la manière suivante : H2 CO < CO < HCO < CH2 OH < CH3 < OH < C. Le carbone
atomique a été adsorbé préférentiellement sur le site de symétrie 4. Une préférence pour le site apical
pour les espèces CH3 et HCO a été observée. Par contraste, les espèces CO, OH, H2 CO et H2 COH
s’adsorbent sur le site apical à l’exception de OH (symétrie 4 avec une couverture de 1.00 MC). Les
énergies des espèces adsorbées (CH3 et OH) ont été utilisées pour calculer les énergies thermochimiques
des potentiels de surface pour la décomposition de ces fragments moléculaires (CH3 et OH), donnant
ainsi l’aperçue de la thermochimie de ces réactions sur la surface Au(001). Cette thermochimie démontre
que la décomposition complète de CH3 ou OH est thermodynamiquement défavorable. Finalement, les
PDOS des adsorbats et substrat ont été analysées.
Mots clés : calculs par la théorie de la fonctionnelle de la densité, méthyle, hydroxyméthyle, formal-
déhyde, formyle, hydroxyle, monoxyde de carbone, carbone, l’Or, adsorption, catalyse.
Abstract
We used the Dacapo code with the GGA-PW91 approximation to study the adsorption of one ato-
mic (C) and molecular (CO) species, and molecular fragments (CH3 , CH2 OH, H2 CO, HCO, OH) on the
Au(001) surface. We have determined for each species the best binding site, adsorption energy, geometri-
cal parameters and projected density of states (PDOS). The magnitude of the adsorption energy (in ab-
solute value) for all the analyzed species can be ordered as follows : H2 CO < CO < HCO < CH2 OH < CH3
< OH < C. The atomic carbon preferred to absorb at the hollow site. A preference for the top site for
the CH3 and HCO species was observed. In contrast, CO, OH, H2 CO and H2 COH species favor top
sites with the exception of OH (hollow with a coverage of 1.0 monolayer (ML)). The energies of the
adsorbed (CH3 and OH) species were used to calculate the thermochemical potential energy surfaces
for decomposition of these molecular fragments (CH3 and OH), giving insight into the thermochemis-
try of these reactions on the Au(001) surface. This thermochemistry demonstrated that the complete
decomposition of CH3 or OH is thermodynamically unfavorable. Finally, the PDOS of adsorbates and
substrate were analyzed.
Keys words : Density functional theory calculations, methyl, hydroxymethyl, formaldehyde, formyl,
hydroxyl, carbon monoxyde, carbon, gold, adsorption, catalysis.