Subscriber access provided by NEW YORK UNIV
Article
Theoretical Modeling of Tribochemical Reaction on Pt and Au Contacts: Mechanical Load and Catalysis Yubo Qi, Jing Yang, and Andrew M. Rappe ACS Appl. Mater. Interfaces, Just Accepted Manuscript • DOI: 10.1021/acsami.5b12350 • Publication Date (Web): 24 Feb 2016 Downloaded from http://pubs.acs.org on March 1, 2016
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
ACS Applied Materials & Interfaces is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 21
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Applied Materials & Interfaces
Theoretical Modeling of Tribochemical Reaction on Pt and Au Contacts: Mechanical Load and Catalysis Yubo Qi, Jing Yang, and Andrew M. Rappe∗ 1
The Makineni Theoretical Laboratories, Department of Chemistry, University of Pennsylvania, Philadelphia, PA 19104-6323 USA E-mail:
[email protected] Phone: 215-898-8313. Fax: 215-573-2112
Abstract
2
3
Micro–electro–mechanical system and nano–electro–mechanical system (MEMS and
4
NEMS) transistors are considered promising for size–reducing and power–maximizing
5
electronic devices. However, the tribopolymer which forms due to the mechanical load
6
to the contacts affects the conductivity dramatically. This is one of the challenging
7
problems that prevent the widespread practical use of these otherwise promising de-
8
vices. Here, we use density functional theory (DFT) to investigate the mechanisms of
9
tribopolymer formation, including normal mechanical load, the catalytic effect, as well
10
as the electrochemical effect of the metal contacts. We select benzene as the background
11
gas, because it is one of the most common and severe hydrocarbon contaminants. Two
12
adsorption cases are considered: one is benzene on the reactive metal surface, Pt(111),
13
and the other is benzene on the noble metal, Au(111). We demonstrate that the for-
14
mation of tribopolymer is induced by both the mechanical load and the catalytic effect
15
of the contact. First, benzene molecules are adsorbed on the Pt surfaces. Then, due ∗
To whom correspondence should be addressed
1
ACS Paragon Plus Environment
ACS Applied Materials & Interfaces
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
16
to the closure of the Pt contacts, stress is applied to the adsorbates, making the C–H
17
bonds more fragile. As the stress increases further, H atoms are pressed close to the Pt
18
substrate and begin to bond with Pt atoms. Thus Pt acts as a catalyst, accelerating
19
the dehydrogenation process. When there is voltage applied across the contacts, the
20
catalytic effect is enhanced by electrochemistry. Finally, due to the loss of H atoms,
21
C atoms become more reactive and link together or pile up to form tribopolymer.
22
By understanding these mechanisms, we provide guidance on designing strategies for
23
suppressing tribopolymer formation.
24
KEYWORDS: tribopolymer, Pt and Au contacts, mechanical load, catalysis, reaction path
25
Introduction
26
Metal oxide semiconductor field effect transistors (MOSFETs) are the fundamental compo-
27
nents in electronic logic devices. Moore’s law predicts that the number of transistors in a
28
defined size electronic device doubles every two years. Effort toward shrinking the MOSFET
29
scale is always ongoing, but a vital challenge is that the leakage current becomes unac-
30
ceptably large as the MOSFET size decreases. 1 Besides, the minimum subthreshold swing
31
S = 60 meV/decade sets a lower limit for energy dissipation in the MOSFET operation. 2,3
32
To overcome these problems, microelectromechanical system (MEMS) switches are being
33
developed. 4–8 The source and drain electrodes are separated by an air gap, and they can also
34
be connected by an adjustable conductive bar. Such an electronic device has nearly zero
35
leakage current, and the minimum subthreshold swing also breaks the S = 60 meV/decade
36
limitation. However, a prime concern about this approach is the stability of the electrical
37
contact resistance (ECR): a hydrocarbon contaminant forms on the surfaces of the electrodes
38
and on the conductive bar after many switching cycles. This tribopolymer may increase the
39
ECR beyond the tolerance of logic applications. 9–13
40
The study of the reliability of MEMS devices has been a fast–moving field. 14 In MEMS
41
devices, Pt and Au are two common electrodes in use. The benefits of using Au electrodes 2
ACS Paragon Plus Environment
Page 2 of 21
Page 3 of 21
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Applied Materials & Interfaces
42
are corrosion retardation, low resistivity, and ease of deposition. 15–17 However, Au is a soft
43
metal with low cohesive energy, so surface wear and material transfer are the main drawbacks
44
of Au contacts. Also, some contamination can still form after electrical cycling. 9,18 As for
45
Pt contacts, previous experiments showed that tribopolymer is more prone to form on Pt
46
than on Au. Wabiszewski et al. used atomic force microscopy (AFM) to mimic the Pt/Pt
47
electrical contacts. The ECR after cycling was six orders of magnitude higher than before
48
. 19 Brand et al. observed considerable amounts of tribopolymer formation on Pt contacts,
49
leading to a dramatic increase of ECR. 11,12
50
For Pt contacts, the rate of formation and the amount of tribopolymer are related to how
51
the MEMS switches are operated. 11,12 Three types of switching-cycle mechanisms are consid-
52
ered in their experimental setup: cold switching, hot switching, and mechanical switching.
53
In cold switching, voltage is applied after the two contacts are closed, and removed before
54
the contacts open. In hot switching, voltage is applied across the device the entire time,
55
whether the contacts are closed or open. In mechanical switching, in order to better un-
56
derstand the role of electrochemistry relative to other effects, the contacts undergo cycles of
57
closing and opening without applying any voltage. When a voltage of 5 V was applied across
58
the contacts during the cold–switching cycle, a large amount of tribopolymer was generated.
59
However, interestingly, even in mechanical switching, the tribopolymer is still formed, but
60
less compared to the cold switching protocol. In this work, we applied a computational com-
61
pression experiment to study the effects of mechanical stress and catalysis on tribopolymer
62
formation.
63
Methodology
64
It is believed that the tribopolymer results from the polymerization of hydrocarbon gases
65
from the atmosphere or from packaging. Among them, benzene is the one which causes the
66
most severe contamination. 13 Therefore, benzene is selected as the background gas in our
3
ACS Paragon Plus Environment
ACS Applied Materials & Interfaces
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
67
study. In order to investigate tribopolymer formation on Pt contacts in a benzene atmo-
68
sphere, we conducted compression computational experiments with first–principles density
69
functional theory (DFT) calculations. Au contacts are also studied as a control system to
70
illustrate the tribopolymer formation mechanism. All calculations are performed with the
71
plane–wave DFT package Quantum–espresso 20 using the generalized gradient approxi-
72
mation (GGA) exchange correlation functional. A plane–wave cutoff energy of Ecut ≈ 680 eV
73
and an 8×8×1 Monkhorst–Pack k–point mesh 21 are used for energy and force calculations.
74
All atoms are represented by norm-conserving pseudopotentials generated by the OPIUM 22
75
code with a plane-wave cutoff, ~2 qc2 /2me ≈ 680 eV. The maximum threshold for atomic force
76
convergence is set as 5.0×10−3 eV/Å per atom.
77
A supercell geometry is used with periodically separated metal slabs and two benzene
78
molecules per surface between the slabs to model metal contacts in a benzene atmosphere.
79
A four–layer metal slab is used to simulate top and bottom contacts. FIG. 1 shows the
80
initial setup of the supercell structures. The space between the top and bottom contacts
81
is 9.3 Å. The vertical distance between the two parallel benzene molecules is selected to be
82
4.1 Å, so that the two benzene molecules have little interaction. In the initial structure, the
83
benzene molecules are intact without any deformation. The distance between one benzene
84
molecule and its closest metal contact is 2.6 Å. At this distance, the benzene molecule has an
85
attractive adsorption force toward the contact, and no stress is required to press the molecule
86
to such a position. After the relaxation, the two benzene molecules are strongly chemisorbed
87
on the top and bottom contacts. Once the adsorption system is fully optimized, the value
88
of the distance between each pair of neighboring metal atomic layers, a, is decreased by ∆a.
89
Since there are three intervals in the four–layer slabs, the height of the supercell is reduced
90
by 3∆a in total in each compression step. At the beginning, the value of ∆a is selected
91
as 0.2 Å. Once any significant deformation of the adsorbates is observed, we decrease ∆a
92
to 0.02 Å or less, in order to have a subtle observation of the structural changes correlated
93
with the height of the supercell. Then, we optimize the compressed system, and after the
4
ACS Paragon Plus Environment
Page 4 of 21
Page 5 of 21
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Applied Materials & Interfaces
94
relaxation, we repeat this compression and optimization process. Our manual adjustment of
95
the supercell is analogous to applying stress to the adsorbed benzene. As the compression
96
process goes on, the normal load (stress) increases. The average stress is calculated from the
97
energy of the system and the dimension of the supercell
σ=−
1 ∂E S ∂z
(1)
98
where σ is the stress along the z direction of the system, E is the energy of the system, and
99
S is the area of the surface. Since we have optimized the structure to equilibrium, the stress
100
in the supercell is approximately uniform and the stress on the molecules is the same as the
101
average one.
102
Energy barriers are calculated with the nudged elastic band (NEB) method implemented
103
in the Quantum–espresso package. The optimization scheme is the quasi–Newton Broy-
104
den’s second method. Path optimization stops when the norm of the force orthogonal to the
105
path is less than 0.05 eV/Å.
106
Results
107
In Pt–benzene systems, we choose three different registries as shown in FIG. 1. In registry
108
(1), there is one carbon atom overlapping between the two benzene molecules from the top
109
view. In the second registry, there is no carbon atom overlap. In registry (3), the two
110
benzene molecules have a relatively large overlap.
111
During the relaxation, the compressed electrodes will recover to a new equilibrium struc-
112
ture, and hence will expand and squeeze the center vacuum space. The adsorbed molecules
113
at the center will sense the stress applied by the two closing contacts. We conduct compres-
114
sion computational experiments and observe that for all three different registries, when the
115
height of the supercell is reduced to 11.06 Å, chemical reaction occurs: some hydrogen atoms
116
dissociate from the adsorbed benzene molecules, making them less saturated. The carbon 5
ACS Paragon Plus Environment
ACS Applied Materials & Interfaces
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
(a)
Page 6 of 21
(b) Top View of three different registries
Benzene molecules are placed at different locations
{ {
Compressed Compressed
{
Not Compressed
{
Compressed
a (1) One C atom overlap 2.6
4.1 (2) No C atom overlap
z
y (3) A large overlap
x
y
z
x
Figure 1: Ball–and–stick model of the initial registries of two benzene molecules. The lattice vectors of the supercell are v1 = (8.117, 0.000, 0.000), v2 = (4.058, 7.029, 0.000) and v3 = (0.000, 0.000, 16.461). (a) Side view of the compression supercell setup with 4–layer Pt (111) slab and two benzene molecules. Compression is initially applied to the metal interlayer spacings, leaving the center vacuum uncompressed (Pt = light grey, C = brown, H = light pink). (b) Top view of three different benzene initial registries on upper and lower Pt slabs. The left panel shows the entire supercell. On the right panel, all benzene molecules sit on the hollow site (three C atoms on top of Pt atoms and three in the Pt 3–fold hollow) of the Pt surfaces with different horizontal intermolecular distances: (1) two carbon directly overlapped, (2) the furthest, no overlap at all, (3) significantly overlapped.
6
ACS Paragon Plus Environment
Page 7 of 21
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Applied Materials & Interfaces
117
atoms which lose their bonded hydrogen atoms form a C–C bond making a biphenyl–like
118
structure, as shown in FIG. 2.
Side View
Top View
(1)
(2)
(3)
Figure 2: Side and top views of the biphenyl–like structure formed after the height of the supercell is reduced to 11.06 Å for the three different registries described in Fig. 1. Each view consists of 2×2×1 supercells. For clarity, in the top view, only the adsorbates and the first layer of Pt atoms beneath are shown.
119
In the following part, we take a more detailed look at the formation of the biphenyl–like
120
structure. For the structure without any compression, the benzene molecules adsorb on the
121
surface and distort to become non–planar. The adsorption is at the hollow site, as it is the
122
most favorable position, 23 as shown in FIG. 4 (a). The C–C bonds are around 1.5 Å, which is
123
approximately the length of C–C single bond. These indicate the chemisorption of benzene
124
on Pt(111) surface. The bonding between C and Pt induces the H atoms to move away from
125
the surface, indicating sp3 C character. The C–Pt bond distances are approximately 2.1 Å.
126
We should also note that before the compression (no external stress applied), the distance
127
between one H atom and its closest Pt atom is approximately 2.7 Å, which is much larger
128
than the length of a Pt–H bond (1.7 Å). 24 There is no chemical bond between Pt and H 7
ACS Paragon Plus Environment
ACS Applied Materials & Interfaces
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
129
Page 8 of 21
atoms.
130
Taking the first registry as an example, the changes of the structure with the height of
131
the supercell c and the resulting σ are shown in FIG. 3. When the stress increases to 24.0
132
GPa, it can been seen that one H atom is pushed close to a Pt atom, as shown in FIG. 4.
c = 16.46 0.00 GPa
c = 14.66 = 0.91 GPa
c = 12.26 = 21.8 GPa
c = 12.08 = 24.0 GPa
c = 14.06 = 2.60 GPa
c = 11.78 = 22.9 GPa
c = 13.46 = 8.40 GPa
c = 12.86 = 15.2 GPa
c = 11.66 = 29.4 GPa
c = 11.06 = 47.7 GPa
Figure 3: Ball–and–stick model of the entire compression process of the supercell with benzene on the Pt(111) surface. c is the length of the supercell along the z–direction. σ is the stress of the supercell at each c value. At c = 11.78 Å and σ = 22.9 GPa, the biphenyl–like structure forms.
133
As we compress the system further, a C–H bond breaks, and the H atom is attached
134
to the Pt substrate. The adsorption position of the dissociated H atom is the hollow site,
135
which is the most stable one. 25 Two C atoms which lose their H atoms form a C–C bond,
136
and connecting the carbon rings builds a biphenyl–like structure. The adsorbed benzene
137
molecules are on the way to more dramatic change: (1) H atoms can fall off due to further
138
compression. The detached atoms may form H2 and escape into the atmosphere, which
139
indicates that the dehydrogenation is a non–reversible process; (2) C atoms which lose their H
140
atoms become less saturated and more reactive. As these processes continue, more hydrogen
141
atoms may fall off, and reactive C atoms link with others or even pile up, forming a two 8
ACS Paragon Plus Environment
Page 9 of 21
(a)
(e)
(b)
120
(c)
Pt-C-H Angle q (°)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Applied Materials & Interfaces
(d)
Reaction occurs
110
100
90
H atom Dissociation
Formation of C-C bond 80
12
14
16
Height of the Supercell (Å)
Figure 4: Ball–and–stick model of detailed benzene bonding geometry and reactivity on Pt (111). (a) Benzene molecule adsorbed at the hollow site. The lengths of the Pt–C bonds are 2.1 Å, and the C–C bonds are 1.54 Å; (b) Part of the structure of the Pt–benzene system with c = 12.08 Å and σ = 24.0 GPa. One H atom is pressed close to a Pt atom; (c) The H atom dissociates from the adsorbed benzene; (d) A biphenyl–like structure forms; (e) The plot of supercell height vs. the Pt–C–H angle. As the compression goes on, the angle is reduced, indicating that the H atoms are pressed toward the Pt contact.
9
ACS Paragon Plus Environment
ACS Applied Materials & Interfaces
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 10 of 21
142
dimensional or three dimensional C–rich network. The biphenyl–like linked rings (Fig. 4)
143
can be viewed as the primary step of tribopolymer formation.
144
To further investigate the dehydrogenation and polymerization processes, we use the
145
NEB 26 method to calculate the energy barrier at different unit cell heights, as shown in Fig.
146
5. Initial structures with different supercell heights provide different stress/mechanical loads.
147
The calculated unit cell height dependence of reaction energies is shown in FIG. 6 (a). For
148
stress below 24.22 GPa, there is a linear relation between stress and reaction barrier, which
149
is consistent with a Brønsted–Evans–Polanyi (BEP) relationship 27–29 as expressed
∆Eact = α∆H = αF · ∆d
(2)
150
where ∆Eact is the change of activation energy Eact , α (0 ≤ α ≤ 1) is a constant which
151
characterizes the position of the transition state, H is the enthalpy of the reaction, F is the
152
applied force, which is proportional to the stress, and ∆d is the distance that the force F
153
moves though along the reaction coordinate. As we increase the stress and the dimension of
154
the system shrinks by ∆d, a work W = F · ∆d
(3)
155
is applied to the system. This energy is added to the adsorbates, promoting the reaction, as
156
shown in Fig. 5 (b). If the stress increases further, the energy barrier drops to zero within
157
just a 0.2 GPa period, as in Fig. 6 (a) and (b). This breakdown of the BEP relationship is
158
attributed to the catalysis effect of Pt. Compression not only adds energy to the adsorbates,
159
but also presses hydrogen atoms close to the Pt substrate, changing the reaction pathway.
160
Pt, as a catalyst, can weaken the C–H bond and reduce the reaction energy barrier. 30–32 The
161
stronger the Pt–H interaction is, the lower the energy barrier is.
162
The reaction occurs when the stress is 24.42 GPa. In the following part, we demonstrate
163
that the threshold stress can be reduced by including the electrochemical effect. Above 24.22
164
GPa, the relationship of the supercell height and energy barrier deviates from linearity. We 10
ACS Paragon Plus Environment
Page 11 of 21
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Applied Materials & Interfaces
165
deduce that from this point, Pt begins to affect the dehydrogenation as a catalyst. The
166
chemical reaction can be expressed as
C6 H6 → C6 H5 · + H+ + e− .
167
(4)
And the reaction energy change due to applied potential U is 33
∗ ∆Eact = ∆Eact − n · e · U
(5)
168
where e is the electron charge and n is the number of electrons involved in the reaction. At
169
24.22 GPa, the energy barrier is 0.051 eV, and in experiments, the applied voltage (5.0 V)
170
during device switching is larger than 0.051 V. Therefore, due to the electrochemical effect,
171
the H atom can get one electron and leave the benzene as soon as it makes electrical contact
172
with the electrode.
173
The process of adsorbed benzene molecules losing H atoms and linking together by C
174
bonds can be summarized as: (1) benzene molecules are adsorbed on the Pt surfaces; (2)
175
due to the closure of the Pt contacts, stress is applied to the adsorbate, flattening it and
176
making C–H bonds easier to break. We emphasize that while the H atoms are not close
177
to the Pt substrate, the fragility of adsorbed benzene mainly results from mechanical load
178
(rather than catalytic effect), and the activation energy–stress relationship obeys the BEP
179
relationship; (3) As the stress increases further, the benzene adsorbates become flat enough
180
and the H atoms begin to be attracted to the Pt atoms. Pt weakens the C–H bond and
181
accelerates the dehydrogenation process, as shown in Fig. 5 (b). This effect can also be
182
enhanced by electrochemistry.
183
In the following section, computational compression experiments of a Au–benzene system
184
are analyzed as a contrast with Pt, highlighting the importance of adsorption and catalysis in
185
polymerization. In FIG. 7, the changes of the Au–benzene structure during the compression
186
are shown. At the outset, benzene molecules are not strongly adsorbed on the contact 11
ACS Paragon Plus Environment
ACS Applied Materials & Interfaces
30
No Polymer Formed Polymer Formed
25 20
Energy (eV)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
15 10 5 0 10
12
14
16
Height of the Supercell (Å)
Figure 5: Energetic behavior of the compression process. The energies are relative energies, with the sum of the energies of a clean Pt slab and two free benzene molecules as zero. (a) Total energy vs. height of the supercell during the compression process. Energies of the structures without polymer formation are shown in red, whereas the black points correspond to the polymerized structures. The subplot shows the detailed stress effect on the polymer 1 evolves to , 2 which is the formation process. With mechanical load applied, structure 3 indicating structure one step before the formation of polymer. The energy drops to , the polymer formation. Continued compression of the cell does not lead to additional new 4 is an expansion step back up to the same supercell volume as . 2 species formation. 5 (b) The energy barrier to form the biphenyl Further expansion changes the structure to ; 1 and 5 are tribopolymer, calculated by NEB method. The supercell configurations at taken as the initial and final structures of the reaction. Schematic model demonstrating the change of the energy barrier is also shown. At the beginning, mechanical load raises the energy levels of the initial, final and transition states (red arrows), making the curve flatter and the reduction of the activation energy. The stress promotes the chemical reaction and the reduce of the energy barrier is proportional to the stress, as expressed in the BEP relationship ∆Eact = αF · ∆d. Then, the height of the barrier drops (blue arrow), due to the catalytic interaction of Pt and H. The three dimensional plot demonstrating the relationship of reaction path, mechanical load and energy barrier is shown in Fig. 6 (b).
12
ACS Paragon Plus Environment
Page 12 of 21
Page 13 of 21
(a)
(b) Height (Å) 12.20 12.16
12.11
12.05
Energy (eV)
7.5 7 6.5 6 5.5 5 4.5 4
0.12 8
Energy Barrier (eV)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Applied Materials & Interfaces
6
0.09
4 2
0.06
12 10 8
0.03
0 22
23.2
23
24
Stress (GPa)
25
23.4
6 23.6
23.8
26
Reaction Path
4 24
24.2
Stress (GPa)
24.4
2
Figure 6: (a) Energy barriers for benzene on Pt at different stresses. Above 24.22 GPa, the relationship of energy barrier and stress is no longer linear, demonstrating that the catalysis effect begins to act. The dashed blue line is extrapolated from the linear part. This line and the x–axis intersect at 26 GPa, indicating the estimated polymerization stress without catalysis; (b) The change of reaction paths with applied stress. The reaction barrier decreases sharply as the stress reaches 24.4 GPa.
c = 19.56 = -0.38 GPa
c = 13.56 = 7.33 GPa
c = 17.16 = 0.14 GPa
c = 12.96 = 10.1 GPa
c = 16.56 = 0.65 GPa
c = 12.36 = 12.1 GPa
c = 15.96 = 1.32 GPa
c = 15.36 = 2.40 GPa
c = 14.76 = 3.77 GPa
c = 11.76 = 22.9 GPa
c = 11.16 = 40.2 GPa
c = 10.56 = 54.6 GPa
Figure 7: Ball–and–stick model of the compression of benzene on the Au(111) surface. c is the supercell height along the z–axis, and σ is the stress of the supercell. Different from the benzene on Pt(111) case, no polymer formation is seen during the compression.
13
ACS Paragon Plus Environment
ACS Applied Materials & Interfaces
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
187
surfaces, with a -0.1 eV adsorption energy, much weaker than on Pt (-1.9 eV). After normal
188
load is applied, benzene molecules spontaneously glide and avoid each other, due to the lack
189
of strong adsorption or fixed location. As the compression increases, the benzene molecules
190
begin to distort, because of the reduction of free space and applied stress. Different from
191
the benzene on Pt contacts case, there is no C–H bonds broken. Even though the benzene
192
ring is broken at the high stress cases and Au–C bonds are able to form, no polymerization
193
is observed, even when the height of the supercell is reduced to 10.56 Å with stress of 54.6
194
GPa, which is far beyond the requirement of benzene polymerization on Pt contacts (≈ 24
195
GPa).
196
As we know, Au is the noblest metal, possessing a completely filled d–band, and the d–
197
band center is far from the Fermi level. 34 As a result, the bonds that form between Au and
198
adsorbates have anti–bonding orbitals that are partially or mostly filled, indicating that the
199
bonding itself is not stable. Therefore, Au contacts cannot adsorb molecules tightly, so there
200
is nothing to stop intermolecular gliding, and Au does not catalyze the dehydrogenation
201
reaction.
202
If the concentration of benzene is low, the molecules tend to escape from the center
203
space between the two Au contacts, making tribopolymer formation less likely. For the
204
high concentration case, if lots of benzene molecules are trapped in a certain region during
205
the switching, the mechanical load can eventually induce the polymerization. Previous
206
theoretical study demonstrates that benzene molecules can transform to a polymeric phase
207
due solely to pressure, 35 but the required polymerization pressure is relatively high, above
208
80 GPa.
209
The main differences between the Pt and Au cases are the adsorption and dehydrogena-
210
tion. Adsorption is important, because it not only prevents the repulsion and avoidance of
211
molecules, but also breaks the strong and stable delocalized π orbital. For the Au contacts
212
with little adsorption, the delocalized π system is crushed only by the mechanical load, and
213
there is a stricter requirement for this process, such as a higher molecule concentration or
14
ACS Paragon Plus Environment
Page 14 of 21
Page 15 of 21
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Applied Materials & Interfaces
214
larger mechanical load. However, for the polymerization on Pt case, the strong delocalized
215
π system is broken by the adsorption onto the substrate. Polymerization is induced by the
216
breaking of the C–H bonds, and this process is much easier.
217
For Pt contacts, computational compression experiments demonstrate that polymeriza-
218
tion can occur under stress of 24.22 GPa. Compared with the pure benzene polymerization
219
stress 80 GPa, it is much closer to the experimental situation, for which the stress between
220
contacts is estimated at several GPa. Besides, only uniaxial stress is investigated in this
221
study. Shear may also assist the dehydrogenation, accelerate the polymerization, and hence
222
lower the polymerization threshold stress. Moreover, in the study, a flat Pt surface is as-
223
sumed. Steps and vacancies can also enhance the catalysis effect and reduce the threshold
224
polymerization stress.
225
Conclusion
226
In this study, the mechanism of tribopolymer formation has been studied by computational
227
compression experiments. Benzene molecules adsorb on the surfaces of contacts, and then
228
original chemical bonds can be broken due to mechanical load. Dehydrogenation makes
229
adsorbed benzene molecules less saturated, so they bond with others to form tribopolymer.
230
In a certain stress range, the reaction energy decreases with mechanical load, and the linear
231
relationship between Eact and P follows the BEP relationship. When the stress is large
232
enough to press H atoms close to the Pt substrate, a catalytic effect takes over and the
233
energy barrier decreases to zero more rapidly as stress increases. Our study provides a
234
detailed analysis of the process and mechanism of the initial stage of tribopolymer formation
235
on contacts in MEMS and NEMS.
15
ACS Paragon Plus Environment
ACS Applied Materials & Interfaces
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 16 of 21
236
ACKNOWLEDGMENTS
237
Y.Q. was supported by the U.S. National Science Foundation, under grant DMR1124696.
238
J.Y. was supported by the U.S. National Science Foundation, under grant CMMI1334241.
239
A.M.R. was supported by the U.S. Department of Energy, under grant DE–FG02–07ER15920.
240
Computational support was provided by the National Energy Research Scientific Computing
241
Center of the Department of Energy.
242
References
243
244
245
246
247
248
249
250
(1) Chen, J.; Chan, T.; Chen, I.; Ko, P.; Hu, C. Subbreakdown Drain Leakage Current in MOSFET. IEEE Electron Device Lett. 1987, 8, 515–517. (2) Taur, Y.; Ning, T. H. Fundamentals of modern VLSI devices, 2nd ed.; Cambridge University Press: Cambridge, United Kingdom, 1998. (3) Salahuddin, S.; Datta, S. Use of Negative Capacitance to Provide Voltage Amplification for Low Power Nanoscale Devices. Nano Lett. 2008, 8, 405–410. (4) Loh, O. Y.; Espinosa, H. D. Nanoelectromechanical Contact Switches. Nature Nanotechnol. 2012, 7, 283–295.
251
(5) Spencer, M.; Chen, F.; Wang, C. C.; Nathanael, R.; Fariborzi, H.; Gupta, A.; Kam, H.;
252
Pott, V.; Jeon, J.; Liu, T.–J. K.; Markovic, D.; Alon, E.; Stojanovic, V. Demonstration
253
of Integrated Micro–Electro–Mechanical Relay Circuits for VLSI Applications. IEEE
254
J. Solid–State Circuits 2011, 46, 308–320.
255
(6) Lee, J. O.; Song, Y.–H.; Kim, M.–W.; Kang, M.–H.; Oh, J.–S.; Yang, H.–H.; Yoon, J.–
256
B. A Sub–1–Volt Nanoelectromechanical Switching Device. Nature Nanotechnol. 2013,
257
8, 36–40.
16
ACS Paragon Plus Environment
Page 17 of 21
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Applied Materials & Interfaces
258
(7) Sinha, N.; Jones, T. S.; Guo, Z.; Piazza, G. Body–Biased Complementary Logic Imple-
259
mented Using AlN Piezoelectric MEMS Switches. J. Microelectromech. Syst. 2012, 21,
260
484–496.
261
(8) Czaplewski, D. A.; Patrizi, G. A.; Kraus, G. M.; Wendt, J. R.; Nordquist, C. D.;
262
Wolfley, S. L.; Baker, M. S.; De Boer, M. P. A Nanomechanical Switch for Integration
263
with CMOS Logic. J. of Micromech. and Microeng. 2009, 19, 085003.
264
(9) Dickrell, D.; Dugger, M. T. Electrical Contact Resistance Degradation of a Hot–
265
Switched Simulated Metal MEMS Contact. IEEE Trans. Compon. Packag. Technol.
266
2007, 30, 75–80.
267
(10) Jensen, B. D.; Chow, L. L.–W.; Huang, K.; Saitou, K.; Volakis, J. L.; Kurabayashi, K.
268
Effect of Nanoscale Heating on Electrical Transport in RF MEMS Switch Contacts. J.
269
Microelectromech. Syst. 2005, 14, 935–946.
270
271
272
273
274
275
276
277
278
279
280
281
(11) Brand, V.; Baker, M. S.; de Boer, M. P. Impact of Contact Materials and Operating Conditions on Stability of Micromechanical Switches. Tribol. Lett. 2013, 51, 341–356. (12) Brand, V.; Baker, M. S.; de Boer, M. P. Contamination Thresholds of Pt– and RuO2 – Coated Ohmic Switches. J. Microelectromech. Syst. 2013, 22, 1248–1250. (13) Hermance, H.; Egan, T. Organic Deposits on Precious Metal Contacts. Bell Syst. Tech. J. 1958, 37, 739–776. (14) van Spengen, W. M. MEMS Reliability From a Failure Mechanisms Perspective. Microelectron. Reliab. 2003, 43, 1049–1060. (15) Hyman, D.; Mehregany, M. Contact Physics of Gold Microcontacts for MEMS Switches. IEEE Trans. Compon. Packag. Technol. 1999, 22, 357–364. (16) Patton, S.; Zabinski, J. Fundamental Studies of Au Contacts in MEMS RF Switches. Tribol. Lett. 2005, 18, 215–230. 17
ACS Paragon Plus Environment
ACS Applied Materials & Interfaces
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
282
(17) Yang, Z.; Lichtenwalner, D. J.; Morris, A. S.; Krim, J.; Kingon, A. I. Comparison of
283
Au and Au–Ni Alloys as Contact Materials for MEMS Switches. J. Microelectromech.
284
Syst. 2009, 18, 287–295.
285
286
(18) Tringe, J. W.; Uhlman, T. A.; Oliver, A. C.; Houston, J. E. A Single Asperity Study of Au/Au Electrical Contacts. J. Appl. Phys. 2003, 93, 4661–4669.
287
(19) Wabiszewski, G. E. Interrogation of Single Asperity Electrical Contacts Using Atomic
288
Force Microscopy with Application to Nems Logic Switches. Diploma Thesis 2013,
289
1–156.
290
(20) Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.;
291
Chiarotti, G. L.; Cococcioni, M.; Dabo, I.; Corso, A. D.; de Gironcoli, S.; Fabris, S.;
292
Fratesi, G.; Gebauer, R.; Gerstmann, U.; Gougoussis, C.; Kokalj, A.; Lazzeri, M.;
293
Martin–Samos, L.; Marzari, N.; Mauri, F.; Mazzarello, R.; Paolini, S.; Pasquarello, A.;
294
Paulatto, L.; Sbraccia, C.; Scandolo, S.; Sclauzero, G.; Seitsonen, A. P.; Smogunov, A.;
295
Umari, P.; Wentzcovitch, R. M. Quantum ESPRESSO: A Modular and Open–Source
296
Software Project for Quantum Simulations of Materials. J. Phys.: Condens. Matter
297
2009, 21, 395502–395520.
298
(21) Monkhorst, H. J.; Pack, J. D. Special Points for Brillouin–Zone Integrations. Phys. Rev.
299
B 1976, 13, 5188–5192.
300
(22) http://opium.sourceforge.net.
301
(23) Saeys, M.; Reyniers, M.–F.; Marin, G. B.; Neurock, M. Density Functional Study of
302
303
304
Benzene Adsorption on Pt (111). J. Phys. Chem. B 2002, 106, 7489–7498. (24) Walter, E. J.; Rappe, A. M. Coadsorption of Methyl Radicals and Oxygen on Rh(111). Surf. Sci. 2004, 549, 265–72.
18
ACS Paragon Plus Environment
Page 18 of 21
Page 19 of 21
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Applied Materials & Interfaces
305
(25) Watson, G. W.; Wells, R. P.; Willock, D. J.; Hutchings, G. J. A Comparison of the
306
Adsorption and Diffusion of Hydrogen on the {111} Surfaces of Ni, Pd, and Pt from
307
Density Functional Theory Calculations. J. Phys. Chem. B 2001, 105, 4889–4894.
308
(26) Henkelman, G.; Uberuaga, B. P.; Jónsson, H. A Climbing Image Nudged Elastic Band
309
Method for Finding Saddle Points and Minimum Energy Paths. J. Chem. Phys. 2000,
310
113, 9901–9904.
311
312
313
314
(27) Black, A. L.; Lenhardt, J. M.; Craig, S. L. From Molecular Mechanochemistry to Stress– Responsive Materials. J. Mater. Chem. 2011, 21, 1655–1663. (28) Bell, G. I. Models for the Specific Adhesion of Cells to Cells. Science 1978, 200, 618– 627.
315
(29) Konda, S. S. M.; Brantley, J. N.; Bielawski, C. W.; Makarov, D. E. Chemical Reactions
316
Modulated by Mechanical Stress: Extended Bell Theory. J. Chem. Phys. 2011, 135,
317
164103.
318
(30) Saeys, M.; Reyniers, M.–F.; Neurock, M.; Marin, G. B. Density Functional Theory
319
Analysis of Benzene (De)hydrogenation on Pt (111): Addition and Removal of the
320
First Two H–Atoms. J. Phys. Chem. B 2003, 107, 3844–3855.
321
(31) Trevor, D. J.; Whetten, R. L.; Cox, D. M.; Kaldor, A. Gas Phase Platinum Cluster
322
Reactions with Benzene and Several Hexanes: Evidence of Extensive Dehydrogenation
323
and Size–Dependent Chemisorption. J. Am. Chem. Soc. 1985, 107, 518–519.
324
(32) Tsai, M. C.; Friend, C. M.; Muetterties, E. L. Dehydrogenation Processes on Nickel
325
and Platinum Surfaces. Conversion of Cyclohexane, Cyclohexene, and Cyclohexadiene
326
to Benzene. J. Am. Chem. Soc. 1982, 104, 2539–2543.
327
(33) Martirez, J. M. P.; Kim, S.; Morales, E. H.; Diroll, B. T.; Cargnello, M.; Gordon, T. R.;
328
Murray, C. B.; Bonnell, D. A.; Rappe, A. M. Synergistic Oxygen Evolving Activity of a 19
ACS Paragon Plus Environment
ACS Applied Materials & Interfaces
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
329
TiO2 –Rich Reconstructed SrTiO3 (001) Surface. J. Am. Chem. Soc. 2015, 137, 2939–
330
2947.
331
332
(34) Hammer, B.; Norskov, J. Why gold is the noblest of all the metals. Nature 1995, 376, 238–240.
333
(35) Wen, X.–D.; Hoffmann, R.; Ashcroft, N. Benzene under High Pressure: a Story of
334
Molecular Crystals Transforming to Saturated Networks, with a Possible Intermediate
335
Metallic Phase. J. Am. Chem. Soc. 2011, 133, 9023–9035.
20
ACS Paragon Plus Environment
Page 20 of 21
Page 21 of 21
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Applied Materials & Interfaces
Stress
Stress
G id & A oi
Polymerize
C6H6 on Pt
C6H6 on Au
Table of Contents Graphics.
21
ACS Paragon Plus Environment