2016 International Conference on Advanced Manufacture Technology and Industrial Application (AMTIA 2016) ISBN: 978-1-60595-387-8 The Study of Coordination Adsorption Effect that CO Adsorption on 4H-SiC (001) Surface about Different Placement Angles Yong-sheng WANG, Zheng-xin YAN * and Shao-hua ZHANG College of Science, Xi an University of Science and Technology, Xi an *Corresponding author Keywords: 4H-SiC, Coordination adsorption, Placement angle, First principles. Abstract. The bond structure, electron density difference, density of state, orbits of CO adsorption on 4H-SiC surface of different placement angles are calculated using the plane wave ultra-soft pseudo-potential based on density functional theory. The calculation results show that placement angles have no effect on adsorption process and the optimization results are basically same. It clear reveals that the main interaction in the process of adsorption comes from the 2p orbit of carbon atom and 3s and 3p orbit of silicon atom. CO would get electron not just from the silicon atom which bond with it but other silicon atoms which are near it in the process of adsorption. Introduction SiC, as the new type third generation semiconductor materials, has more superior performance than Si and GaAs which represent the first and second generation semiconductor materials, respectively. The adsorption properties of SiC surface have been studied by many scholars in recent years [1-3]. Jin [4] has studied the surfaces property from the β-sic (111) derived structure and its water molecule adsorption by first-principles. James Westover [5] has studied metallization of the β-sic (100)-3 2 surface by a DFT investigation, which revealed that metallization could be achieved via hydrogen atoms occupying the second silicon layer. Du [6] has studied the mechanism of hydrogen production via water splitting on 3C-SiC s different surfaces, which indicated that the adsorption behavior of water molecule could take place on 3C-SiC s different surfaces and it leads to the surface reconstruction. Although there are most study on SiC surface, the study of coordination adsorption effect that CO adsorption on 4H-SiC (001) surface about different placement angles are fewer. The purpose of this paper is to study the coordination adsorption effect of CO adsorption on 4H-SiC (001) surface. According to DFT (density functional theory), CO adsorption on 4H-SiC (001) surface models are set up by CASTEP (Cambridge serial total energy package). By geometry optimization and quantum calculation, we will obtain the bond structure, electronic orbital and density of states of models to analyze the coordination adsorption effect. Simulation Theories In this work, CASTEP [7-8] (Cambridge serial total energy package) package in Materials Studio software is used for calculation. The adsorption models are contrastively calculated through the first principles and DFT (density functional theory). To obtain better adsorption state, we selected GGA(Generalized Gradient Approximation) PBE (Perdew, Burke and Emzerhof) method [9-14] of ultrasoft pseudopotentials [15-17] in calculation. For an economical calculation, simulation was performed using 400eV cut off energy and a 4 3 1 k-point sampling grid. The convergence precision in iterative process was 2.0 10-6eV/atom. A (2 3) surface unit cell with a slab of one layer thickness was selected as the model. This slab was repeated periodically with a 10Åof vacuum region between the slabs. Four different angles that CO adsorption on 4H-SiC(001) surface are set to study coordination adsorption effect, adding atoms fixed direction when we made models to analyze its angled effects, as shown in Figure 1. From Figure 1(a),we can see carbon of CO have bonded two neighbor silicon atoms.
1.13 Å 1.13 Å 45 60 1.13 Å 75 1.13 Å 90 (a)45 (b)60 (c)75 (d)90 Figure 1. Four different angles between CO and 4H-SiC (001) surface. The Calculation Results and Discussion Bond Structure Change After geometry optimization, we can obtain the models as shown in Figure 2.Comparing Figure 2 with Figure 1, the bond structure of models has obvious change. Some data of bond structure are measured, as listed in Table 1. (a)45 (b)60 (c)75 (d)90 Figure 2. Four different angles between CO and 4H-SiC (001) surface after geometry optimization. We can see the site of CO adsorption on 4H-SiC, as shown in Figure 2, have obvious change. According to the data that list in Table 1, bond length of C-O and C-Si just have little increase comparing with the bond length we set CO molecule before. However, angles have clear variation after geometry optimization. The amusing is that four different angles between CO and 4H-SiC (001) surface have the same angles in the stable adsorption states. CO is mainly perpendicular to 4H-SiC surface, the reason we guess is that the coordination adsorption effect take some actions in this process. Table 1. Data of bond structure. C-O bond(å) C-Si bond(å) A 1 ( ) a A 2 ( ) b 45( ) 1.201 2.072 48.572 85.977 60( ) 1.200 2.076 48.910 85.100 75( ) 1.201 2.104 49.905 85.913 90( ) 1.201 2.095 48.831 87.188 a A 1 is the angle between C-Si bond and 4H-SiC(001) surface b A 2 is the angle between O-C bond and 4H-SiC(001) surface Electron Density Difference Analysis On the basis of the calculated results, we obtained the electron density difference that shows in Figure 3 and it definite amounts are listed in Table 2. Red and blue represent getting electron and losing electron respectively. The atoms of O and C get electron from silicon atom of 4H-SiC (001). An interesting characteristic in electron transfer is that the electron change is mainly identical of four different placement angles, as shown in Table 2. When CO adsorb on 4H-SiC (001) surface, we can reach, that the electron transfer condition never mind with placement angle. Comparing with atomic orbital charge population before and after adsorption, it clear reveals that the main
interaction in the process of adsorption comes from the 2p orbit of carbon atom and 3s and 3p orbit of silicon atom. Look through, the electron that CO got are more than the electron that silicon atom lost which bond with carbon atom. So, CO would get electron from other silicon atom which is near it in the process of adsorption. clean 45( ) 60( ) 75( ) 90( ) (a)45 (b)60 (c)75 (d)90 Figure 3. Electron density difference of four different angles. Table 2. Atomic orbital charge population of four different angles. Atoms S P Charge(e) O 1.84 4.58 C 1.68 1.91 Si 1.14 1.82 1.05 O 1.83 4.56-0.39 C 1.29 2.72-0.02 Si 1.06 1.79 1.15 O 1.83 4.56-0.39 C 1.29 2.73-0.02 Si 1.06 1.79 1.15 O 1.83 4.55-0.39 C 1.30 2.71 0.00 Si 1.08 1.79 1.14 O 1.83 4.68-0.39 C 1.30 2.71-0.01 Si 1.07 1.78 1.14 The electrons of 2s from O in 4 complexes keep a constant, 1.83, which show that the 2 orbital takes few actions in the bonding process. Contrary, 2 (px, py, pz) orbitals taking accepter which obtain mainly some electrons from C (and partly from Si) take positively part in this process. C-O bond length increases in complex 4 evidences this point of view. Density of State Analysis. To get more study the property of CO adsorption on 4H-SiC (001) surface in different placement angles, we have obtained the density of state after calculation, as shown Figure 4 and Figure 5. DOS (Density of states), including PDOS (partial DOS), can provide a fundamental understanding of the adsorbate-substrate interaction. Figure 4. Density of states of oxygen molecule before adsorption and that at four different placement angles after adsorption respectively (a) CO;(b) 45 ; (c) 60 ; (d) 75 ; (e) 90.
As shown in Figure 4, we can see that it is mainly contributed of P orbit of CO molecule near Fermi energy. In the Figure 4 (a), electron peak shows a strong activity at 0 ev. Compared Figure 4 (a) with (b)-(e), there are new waves produced and the waves move to low energy location, which shows there are new chemical bonds produced and the system is more stable after adsorption than before adsorption. Figure 5. The DOS of CO adsorption on 4H-SiC (001) surface of four different angles. (a) Clean surface (b) 45 ; (c) 60 ; (d) 75 ; (e) 90. Figure 5 (a) shows the density of states of clean 4H-SiC (001) surface. (b)-(e) represent the density of states of the whole system in the process of CO adsorption on 4H-SiC (001) surface at four different placement angles respectively. The first electron peak in Figure 5 (b)-(e) has move to high energy location in varying degrees. It is the result of the influence of distance that the oxygen atom leaves 4H-SiC (001) surface. At -8eV, the waves become smoothness especial in Figure 5(c). Orbits Analysis. Known by us, the electronic configuration of CO is (1 )2(2 )2(3 )2(4 )2(1 )4(5 )2(2 )0 (6 )0. HOMO LUMO (a)45 (b)60 (c)75 (d)90 Figure 6. The HOMO and LUMO of CO adsorption on 4H-SiC (001) surface of four different angles. From the chemical activity for the different initial angles, the highest occupied orbitals and the lowest unoccupied orbital were analyzed in detail. The hybrid orbital of s,p of SiC interacts with pz of C to change to bonding orbital. The lowest unoccupied orbital py of O interacts with py in C to become to exhibit a repulse effect to explain the increasing LOCS. The interactions of orbital of Figure 6(b) and (c) are similar to the orbital above. C atom gets electron from Si when its P orbital interacts with corresponding orbital in Si, meanwhile, O atom also gives some electrons to C atom to increase its negative electron, meanwhile, anti-bonding in C-Si is weaken duo a long distance
between them. In Figure 6(d), C-O keeps an initial perpendicular pose for the interactions from C counteracting with interaction from O atom. O lose some electrons in this process, which pz orbital obtains some electrons from pz of C atom, C gets some electrons from Si and takes as receiver, the transaction in Figure 6(b) and (c) is similar, interestingly, O gets some electrons when it loses some electrons. But it loses some electrons in fact. In Figure 6(d), O and C will lose some electrons to Si bulk which leads a more active Si surface. Common trait in four complex is that coordination covalent bond was strengthened in CO adsorbed process which would need a bigger space. This is main reason for the perpendicular adsorption structure. Conclusions The calculation results show that placement angles have no effect on adsorption process and the optimization results are basically same. It clear reveals that the main interaction in the process of adsorption comes from the 2p orbit of carbon atom and 3s and 3p orbit of silicon atom. CO would get electron not just from the silicon atom which bond with it but other silicon atoms which are near it in the process of adsorption. Acknowledgment The work was partly supported by National natural science foundation of China (No. 11405127 and No11404257.) and Department of Shanxi province natural science foundation of China (2015JM1027). References [1] J.C. Zhou, X.Q. Zheng, F. Liu, The Latest Research Progress of SiC Films and Devices Membrane, Materials Review,21(2007)112-114. [2] S. Zhao, X.G. Zhou, J.S. Yu, Research and Development in Fabrication and Properties of SiC/SiC Composites, Materials Review A, 21(2013)66-70. [3] D.J. Yang, L.L. Li, Development of SiC Semiconductor Materials Industry, Development and Application of Materials, 6(2015)110-114. [4] H.Y. Jin, First-principles study of surfaces property from the β-sic(111) derived structure and its water molecule adsorption, Yanshan University, 2010. [5] J. Westover, H. Enriquez, H Oughaddou, A. Kara, Metallization of the β-sic(100) 3 2 surface: A DFT investigation, Surface Science,606(2012)1471-1474. [6] J.L. Du, B. Wen, R. Melnik, Mechanism of hydrogen production via water splitting on 3C-SiC s different surfaces: A first-principles study, Computational materials Science, 95(2014) 451-455. [7] M.D. Segall, P.J.D. Lindan, M.J. Probert, C.J. Pickard, P.J. Hasnip, S.J. Clark, M.C. Payne, First-principles simulations: ideas, illustrations and the CASTEP code, Journal of Physics: Condens Matter, 14(2002)1727-1735. [8] S.J. Clark, M.D. Segall, C.J. Pickard, P.J. Hasnip, M.I.J. Probert, K. Refson, M.C. Payne, First Principles Methods Using CASTEP, Z. Kristallogr, 220(2005)567-572. [9] G. Kresse, J. Hafner, Ab initio molecule dynamics for liquid metals, Phys. Rev. B. 47 (1993) 558-564. [10] G. Kresse, J. Hafner, Ab initio molecule-dynamics simulation of the liquid-metal amorphous semiconductor transition in germanium, Phys. Rev. B.49(1994)14251-14257.
[11] G. Kresse, J. Furthmüller, Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set, Phys. Rev. B.54(1996)11169-11174. [12] G. Kresse, J. Furthmüller, Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set, Computational Materials Science,6(1996)15-20. [13] J.P. Perdew, J.A. Chevary, S.H. Vosko, K.A. Jackson, M.R. Pedersen, D.J. Singh, C. Fiolhais, Atoms, Molecules, Solids, and Surfaces: Applications of the Generalized Gradient Approximation for Exchange and Correlation, Phys. Rev. B. 46(1992)6671-6680. [14] R.O. Jones, O. Gunnarsson, The Density Functional Formalis, its applications and prospects, Rev. Mod. Phys, 61(1989)689-796. [15] A.D. Boese, J.M.L. Martin, Development of Density Functions for Thermochemical Kinetics, J. Chem. Phys, 121(2004)3405-3411. [16] C. Adamo, V. Barone, Toward Reliable Functions Methods Without Adjust Parameters. the PBEO Model, J. Chem. Phys, 110(1999)6158-6165. [17] F.A. Hamprecht, A.J. Cohen, D.J. Tozer, N.C. Handy, Development and assessment of new exchange-correlation functions, J. Chem. Phys, 109(1998)6264-6270.