The interaction of ?Al(CH
3
)
2
with ?OH on a fully OH?terminated Si (001) surface was studied using density functional theory. Two sites for Al(CH
3
)
3
to react with the ?OH on the surface were identified. The ?Al(CH
3
)
2
product energetically favored the dimer?row site rather than the inter?row site because the Al atom of ?Al(CH
3
)
2
at the dimer?row site was attracted by the lone pair electrons of the O atom in the neighboring ?OH. The energy barrier for the transfer of the ?Al(CH3)2 between the two sites was 0.11 eV, and therefore, the ?Al(CH
3
)
2
at the inter?row site can easily transfer to the dimer?row site at room temperature.
1. INTRODUCTION
Functional one?dimensional metallic oxides have attracted much at gate oxides for metal oxide semiconductor field effect transistors (MOSFETs) have been fabricated using atomic layer deposition (ALD) because this technique allows the uniformity and thickness of the gate oxide film to be delicately controlled
[1]
,
[2]
. Among the potential gate oxide materials, aluminum oxide (Al
2
O
3
) has superior physical and electronic properties including high band gap (∼9 eV), high dielectric constant (∼9), and high breakdown field (5?10 MV/cm) [3]?[5]. Various aluminum precursors including AlCl
3
, AlMe
3
, AlBr
3
, AlMe
2
Cl, AlEt
3
, Al(OEt)
3
, and Al(O
i
Pr)
3
have been employed with O
2
, O
3
, H
2
O, and H
2
O
2
reactants for the growth of high quality Al2O3 thin films using the ALD process
[6]
?
[14]
.
Recently, many studies on tri?methylaluminum (TMA, AlMe
3
) have been performed using density functional theory (DFT)
[15]
?
[19]
. Widjaja and Musgrave
[15]
, Heyman and Musgrave
[16]
carried out ALD half reactions with TMA, AlCl
3
, and H
2
O. In addition, Lee et al.
[17]
were able to prepare a hydroxyl (OH)?terminated Si (100) 2 × 1 surface. Halls and Raghavachari
[18]
calculated the growth reaction of Al
2
O
3
with TMA on a reconstructed H?terminated Si (100) surface. Recently, Ghosh and Choi
[19]
researched the initial growth mechanism of Al
2
O
3
using TMA and H
2
O on a fully OH?terminated Si (100) 2 × 1 surface. They found that the fully OH?terminated Si surface allowed a ring?closing reaction, yielding a unique fivemembered ring structure, essentially saturating the surface bonding sites between the Si surface and the Al
2
O
3
junction region. However, the main focus of that work was on the TMA reaction with ?OH along the path of the dimer?row site, and the effect of the lone pair electrons of the ?OH near the di?methylaluminum (DMA) group (?Al(CH
3
)
2
) was only briefly mentioned.
In this study, we present two reaction sites of TMA with the ?OH on a fully OH?terminated Si (001) surface, and investigate the effect of the lone pair of ?OH electrons on the stability of ?Al(CH
3
)
2
. In addition, the transfer energy barrier of ?Al(CH
3
)
2
is calculated using the nudged elastic band (NEB) tool
[20]
.
Two sites of the tri?methylaluminum (TMA) transfer when TMA approached and reacted with ?OH on the fully OH?terminated Si (001) surface. The white, dark grey, red, orange, and yellow spheres are the H, C, O, Al, and Si atoms, respectively.
2. EXPERIMENTS
DFT calculations were performed using the Vienna ab?initio simulation package (VASP) code with the projector augmented wave (PAW) potentials, an all?electron method that combines the accuracy of augmented plane waves, the cost? effective pseudopotential implemented in VASP, and the generalized gradient approximation (GGA)
[21]
?
[24]
.
The residual minimization method of direct inversion in the iterative subspace (RMM?DIIS) was used for calculating the ground state of electrons
[25]
,
[26]
. More precisely, the Al35s
2
and 3p
1
states, O 2s
2
and 2p
4
states, C 2s
2
and 2p
2
states, H 1s
1
state, and Si 3s
2
and 3p
2
states were treated as valence wave functions. The cutoff energy was 500 eV, and the Monk? horst pack k?points mesh was 2×2×1 for the slab structure.
The most energetically favorable Si (001) 4×4 surface that was fully OH?terminated was first determined by calculation. H atoms of the ?OH groups that were attached to the Si dimers of the surface were aligned along the dimer?row direction in zigzags. There were two sites for the TMA transfer from the vacuum to the surface, the inter?row site and dimerrow site, when TMA approached and reacted with ?OH on the fully OH?terminated Si (001) surface, as shown in
Fig. 1
. The inter?row and the dimer?row sites indicate the sites of TMA with respect to ?OH attached to a Si atom of the Si dimer. The reaction between TMA and ?OH generated methane (CH
4
) as a by?product and the remaining TMA fragment, ?Al(CH
3
)
2
, was attached to the ?O radical. Si surfaces with ?Al(CH
3
)
2
located at the two sites were calculated by DFT and the resulting electron densities of ?Al(CH
3
)
2
groups were analyzed using the visualization for electronic and structural analysis (VESTA) tool
[27]
. In addition, the energy barrier was calculated when ?Al(CH
3
)
2
transferred from the inter?row site to the dimer?row site using the NEB tool.
3. RESULTS AND DISCUSSION
The chemical formula of the TMA reaction with ?OH is given by Si?OH + Al(CH
3
)
3
→ Si?O?Al(CH
3
)
2
+ CH
4
[19]
.
Fig. 2
The charge density of ?Al(CH3)2 at (a) the inter?row site and (c) the dimer?row site along the [110] direction, and the planar view of ?Al(CH3)2 at (b) the inter?row site and (d) the dimer?row site along the [001] direction. The blue, orange, and red arrows represent the transfer directions of the H and O atoms.
shows the charge density distribution and the transfer of atoms when ?Al(CH
3
)
2
is located at the inter?row site or the dimer?row site of the fully OH?terminated Si (001) surface along the [110] (
Figs. 2
(a) and (c)) and
[001]
(
Figs. 2
(b) and (d)) directions. An isosurface level is the amount of electron charge per a
0
3
(a
0
= Bohr radius); for this case its value in
Fig. 2
is 0.063. O′ and O′′ in
Fig. 2
indicate the oxygen atom reacting with the Al atom of ?Al(CH
3
)
2
and the oxygen atom of ?OH attached to the other Si atom of the same dimer, respectively. The O′ atom of ?O′H was bonded with the Al atom of ?Al(CH
3
)
2
when TMA transferred to the inter?row site and reacted with ?O′H, as shown in
Figs. 2
(a) and (b). H atoms of the ?OH groups located near the two CH
3
groups of ?Al(CH
3
)
2
were repelled (indicated by the blue arrows for the H atoms). The O′?Si?Si angle and the distance between O′ and O′′
The angles and distance between two atoms around ?Al(CH3)2.
The angles and distance between two atoms around ?Al(CH3)2.
atoms were decreased by 2.08% and 3.26%, respectively, because the O′ atom transferred under the influence of the ?Al(CH
3
)
2
, as shown in
Table 1
. Si atoms were covalently bonded to Si atoms, whereas O atoms were mainly ionically bonded to Si atoms. The O′ atom was bonded with the Al atom of ?Al(CH
3
)
2
by the p orbital electron of the O′ atom. When the ?Al(CH
3
)
2
was located at the inter?row site, where there was enough space for TMA to react with the OH group, the energy was ?1.76 eV, which was 0.12 eV higher than the value of Halls and Raghavachari
[18]
. The difference may be due to the differently terminated Si (001) surface since they used a ½ OH? and ½ H?terminated Si (001) surface. When ?Al(CH
3
)
2
was located at the dimer?row site between the O′ atom and ?O′′H, as shown in
Figs. 2
(c) and (d), the H atom of ?O′′H located near the ?Al(CH
3
)
2
transferred farther away from ?Al(CH
3
)
2
than in the case of the inter?row site, while the O atoms (O′ and O′′) transferred closer to each other (indicated by red arrows near the O atoms). In addition, the O′′?H bond length was decreased to 0.76 Å from 0.97 Å (i.e., by 21.65% as the O?H bond length located in the other sites was 0.97 Å), and the O′′?H bond rotated clockwise by approximately 55
°
. This indicates that the lone pair electrons of the O′′ atom transfer toward the Al atom of ?Al(CH
3
)
2
to create a bond between the O′′ and Al atom. The H atom located on the other dimer?row site was transferred (indicated by the orange arrow) because it was influenced by the H atom of ?O′′H due to repetition of the 4×4 size of the Si surface in this study. In this situation, the O′?Si?Si angle and the distance between O′ and O′ atoms were further decreased by 9.0% and 23.1%, respectively, relative to the case of the inter?row site because of the lone pair electrons of the O′′ atom. ?Al(CH
3
)
2
at the dimer?row site was energetically more favorable than that at the inter?row site by 0.43 eV, as shown in
Fig. 2
. The O′?Al bond length was increased compared to that of the inter?row site, as shown in
Table 1
, indicating that the Al atom was attracted to the O′′ atom.
Figure 3
shows the line profiles of the electron densities of the O′?Al and O′′?Al, which are the numerical values calculated by the linear interpolation of the volumetric data. The imaginary Al atoms (Ali, a green dotted circle in the inserted picture) were produced for comparison. When ?Al(CH
3
)
2
was at the inter?row site, the electron distributions around the non?bonded (green line, as shown in
Fig. 3
) and the bonded (blue line, as shown in
Fig. 3
) O atoms were different due to the Al atom, as shown in
Fig. 3
(a). We found that there was excess electron density at the O′?Al bond compared to that at the O?Ali bond because of the electron redistribution of the O′ atom to the O′?Al bond. When ?Al(CH
3
)
2
was at the dimer?row site, the O′′ atom transferred toward the Al atom
The line profiles of electron densities of (a) O′?Al with ?Al(CH3)2 on the inter?row site and (b) O′?Al?O′′ with ?Al(CH3)2 on the dimer?row site. The blue solid line and green dotted line indicate the real (Al) and imaginary Al atom (Ali), respectively.
Energy variation as a function of the ?Al(CH3)2 transfer fromthe inter?row site to the dimer?row site.
to form a bond by using the lone pair electrons of the O′′ atom, as shown in
Fig. 3
(b). The green dotted line was the shifted line from the green line to be compared with the blue line directly because the distance of the O?Al
i
?O was longer than that of the O′?Al?O′′. In
Fig. 3
(b), we can clearly see that the lone pair electrons of the O′′ atom are located between the O′′ and the Al atoms as the excess electron density of the O′′?Al bond. The excess electrons of the O′′ atom were two times richer than the excess electrons of the O′ atom because those of the O′′ atom had the lone pair electrons. Some electrons of the H atom of the O′′?H also transferred to the O′′?Al bond, which was indicated by the reduction of the O′′?H bond length, as shown in
Fig. 2
(d).
Figure 4
shows the calculated energy barrier of the ?Al(CH
3
)
2
transfer from the inter?row site to the dimer?row site, i.e., 0.11 eV. Since the energy barrier is quite small, the ?Al(CH
3
)
2
can easily transfer from the inter?row site to the dimer?row site. Therefore, the ?Al(CH
3
)
2
would primarily be located at the dimer?row site at room temperature.
4. CONCLUSIONS
We studied the interaction of ?Al(CH
3
)
2
with ?OH on a fully OH?terminated Si (001) surface in order to determine the mechanism of the initial Al
2
O
3
growth using DFT. There were two sites when TMA approached and reacted with ?OH on the fully OH?terminated Si (001) surface. ?Al(CH
3
)
2
favored the dimer?row site rather than the inter?row site, because the lone pair electrons of the O atom of ?OH located next to the ?Al(CH
3
)
2
on the same Si dimer attracted the Al atom of ? Al(CH
3
)
2
and reduced the overall total energy. Based on the low energy barrier from the inter?row site to the dimer?row site, the ?Al(CH
3
)
2
should mostly remain at the dimer?row site at room temperature.
Acknowledgements
This work was supported by the Korea Science and Engineering Foundation (KOSEF) through the Basic Science Programs for General Research Grants (No. 2009?0076318).
Ritala Mikko
,
Kukli Kaupo
,
Rahtu Antti
,
Raisanen Petri I.
,
Leskela Markku
,
Sajavaara Timo
,
Keinonen Juhani
2000
Atomic Layer Deposition of Oxide Thin Films with Metal Alkoxides as Oxygen Sources
Science
288
(5464)
319 -
DOI : 10.1126/science.288.5464.319
George S.M.
,
Ott A.W.
,
Klaus J.W.
1996
Surface Chemistry for Atomic Layer Growth
The Journal of Physical Chemistry
100
(31)
13121 -
DOI : 10.1021/jp9536763
Ye P.D.
,
Wilk G.D.
,
Kwo J.
,
Yang B.
,
Gossmann H.-J.L.
,
Frei M.
,
Chu S.N.G.
,
Mannaerts J.P.
,
Sergent M.
,
Hong M.
,
Ng K.K.
,
Bude J.
2003
GaAs MOSFET with oxide gate dielectric grown by atomic layer deposition
Electron Device Letters IEEE
24
(4)
209 -
211
DOI : 10.1109/LED.2003.812144
Klein T. M.
,
Niu D.
,
Epling W. S.
,
Li W.
,
Maher D. M.
,
Hobbs C. C.
,
Hegde R. I.
,
Baumvol I. J. R.
,
Parsons G. N.
1999
Evidence of aluminum silicate formation during chemical vapor deposition of amorphous AL2O3thin films on Si(100)
Applied Physics Letters
75
(25)
4001 -
4003
DOI : 10.1063/1.125519
Manchanda L.
,
Morris M.D.
,
Green M.L.
,
van Dover R.B.
,
Klemens F.
,
Sorsch T.W.
,
Silverman P.J.
,
Wilk G.
,
Busch B.
,
Aravamudhan S.
2001
Multi-component high-K gate dielectrics for the silicon industry
Microelectronic Engineering
59
(1-4)
351 -
359
DOI : 10.1016/S0167-9317(01)00668-2
Grigoras K
,
Sainiemi L
,
Tiilikainen J
,
Saynatjoki A
,
Airaksinen V-M
,
Franssila S
2007
Application of ultra-thin aluminum oxide etch mask made by atomic layer deposition technique
Journal of Physics: Conference Series
61
369 -
DOI : 10.1088/1742-6596/61/1/074
Ritala Mikko
,
Saloniemi Heini
,
Leskela Markku
,
Prohaska Thomas
,
Friedbacher Gernot
,
Grasserbauer Manfred
1996
Studies on the morphology of AL2O3thin films grown by atomic layer epitaxy
Thin Solid Films
266
(1-2)
54 -
58
DOI : 10.1016/s0040-6090(95)08524-6
Yates David J.C.
,
Dembinski G.W.
,
Kroll W.R.
,
Elliott James John
1969
Infrared studies of the reactions between silica and trimethylaluminum
The Journal of Physical Chemistry B
73
(4)
911?921 -
DOI : 10.1021/j100724a026
Brazeau Allison L.
,
Barry Sean T.
2008
Atomic Layer Deposition of Aluminum Oxide Thin Films from a Heteroleptic Amidinate-Containing Precursor
Chem Mater
20
(23)
7287 -
7291
DOI : 10.1021/cm802195b
Anikeev G.V.
,
Ezhovskii Y.K.
,
Koltsov S.I.
1988
Inorg Mater
24
514 -
Kukli Kaupo
,
Ritala Mikko
,
Leskela Markku
,
Jokinen Janne
1997
Atomic layer epitaxy growth of aluminum oxide thin films from a novel Al(CH3Cl precursor and H_2O
J Vac Sci Technol A
15
(4)
2214 -
2218
DOI : 10.1116/1.580536
Katamreddy Rajesh
,
Inman Ronald
,
Jursich Gregory
,
Soulet Axel
,
Takoudis Christos
2006
ALD and Characterization of Aluminum Oxide Deposited on Si(100) using Tris(diethylamino) Aluminum and Water Vapor
J Electrochem Soc
153
(10)
C701 -
C706
DOI : 10.1149/1.2239258
Huang Ron
,
Kitai Adrian H.
1993
Preparation and characterization of thin films of MgO AL2O3and MgAl_2O_4by atomic layer deposition
Journal of Electronic Materials
22
(2)
215 -
220
DOI : 10.1007/BF02665029
Hiltunen L.
,
Kattelus H.
,
Leskela M.
,
Makela M.
,
Niinisto L.
,
Nykanen E.
,
Soininen P.
,
Tiittad M.
1991
Growth and characterization of aluminium oxide thin films deposited from various source materials by atomic layer epitaxy and chemical vapor deposition processes
Materials Chemistry and Physics
28
(4)
379 -
388
DOI : 10.1016/0254-0584(91)90073-4
Widjaja Yuniarto
,
Musgrave Charles B.
2002
Quantum chemical study of the mechanism of aluminum oxide atomic layer deposition
Applied Physics Letters
80
(18)
3304 -
DOI : 10.1063/1.1473237
Heyman A.
,
Musgrave C. B.
2004
A Quantum Chemical Study of the Atomic Layer Deposition of Al2O3Using AlCl3and H2O as Precursors
The Journal of Physical Chemistry B
108
(18)
5718 -
DOI : 10.1021/jp049762x
Lee Sun S.
,
Baik Jae Y.
,
An Ki-Seok
,
Suh Yung D.
,
Oh Jin-Ho
,
Kim Yunsoo
2004
Reduction of Incubation Period by Employing OH-Terminated Si(001) Substrates in the Atomic Layer Deposition of AL2O3
J Phys Chem B
108
(39)
15128 -
15132
DOI : 10.1021/jp048038b
M D. Halls
,
K Raghavachari
2004
Atomic Layer Deposition Growth Reactions of AL2O3on Si(100)-2×1
J Phys Chem B
108
(13)
4058 -
4062
DOI : 10.1021/jp0378079
Ghosh Manik Kumer
,
Choi Cheol Ho
2006
The initial mechanisms of AL2O3atomic layer deposition on OH/Si(1 0 0)-2×1 surface by tri-methylaluminum and water
Chemical Physics Letters
426
(4-6)
365 -
369
DOI : 10.1016/j.cplett.2006.05.126
Sheppard Daniel
,
Terrell Rye
,
Henkelman Graeme
2008
Optimization Methods for Finding Minimum Energy Paths
The Journal of Chemical Physics
128
(13)
134106 -
134115
DOI : 10.1063/1.2841941
Kresse G.
,
Hafner J.
1994
Ab initio molecular-dynamics simulation of the liquid-metal?amorphous-semiconductor transition in germanium
Phys Rev B
49
(20)
14251 -
14269
DOI : 10.1103/PhysRevB.49.14251
Kresse G.
,
Furthmuller J.
1996
Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set
Computational Materials Science
6
(1)
15 -
50
DOI : 10.1016/0927-0256(96)00008-0
Kresse G.
,
Furthmuller J.
1996
Efficient Iterative Schemes for Ab initio Total-Energy Calculations Using a Plane-Wave Bases Set
Phys Rev B
54
(16)
11169 -
11186
DOI : 10.1103/PhysRevB.54.11169
Kresse G.
,
Joubert D.
1999
From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method
Phys Rev B
59
(3)
1758 -
1775
DOI : 10.1103/PhysRevB.59.1758
Momma K.
,
Izumi F.
2008
VESTA: a three-dimensional visualization system for electronic and structural analysis
Journal of Applied Crystallography
41
(3)
653 -
658
DOI : 10.1107/S0021889808012016