Acessibilidade / Reportar erro

Study of the M-doped Effect (M=Al, Ag, W) on the Dissociation of O2 on Cu Surface Using a First Principles Method

Abstract

The O2 adsorption and dissociation on M-doped (M=Al, Ag, W) Cu (111) surface were studied by density functional theory. The adsorption energy of adsorbate, the average binding energy and surface energy of Cu surface, and the doping energy of doping atom were calculated. All the doped atoms can be stably combined with Cu atoms and improve the surface activity of Cu surface, while the Al and W atoms would strengthen the bonding effect between the atoms on the Cu (111) slab except Ag doping atom. Due to the different electronegativity of three metals of Al, Ag, W, these doping atoms can resist the dissociation of O2 by DOS analysis. The potential energy surface was computed, and the result showed that the dissociation reaction of O2 on the surfaces not only reflected in the barrier energy, but also the reaction energy. Ag atoms had the best resistance effect on the O2 dissociation comparing with Al and W atoms because of the large barrier energy and reaction energy.

Keywords:
O2 adsorption; Metal doped; Cu (111); Dissociation


1. Introduction

The reaction of gas molecule and metal is under study both experimentally and computationally for metal Cu, as it is related to many processes such as catalysis, deposition, coating, corrosion or oxidation11 Zhao J, Tao F. Dissociation mechanism of H2O on clean and oxygen-covered Cu (112) surfaces, A theoretical study. Vacuum. 2016;128:252-8.,22 Akca A, Genc AE, Kutlu B. BH4 dissociation on various metal (111) surfaces: A DFT study. Appl Surf Sci. 2019;473:681-92.. Among them, oxidation of Cu and Cu-based alloy surfaces is both common and important phenomenon because the oxide film formed in air protects the surface against further reaction and corrosion. The steps in the oxidation process involve dissociative adsorption of O2 on Cu-based alloy surfaces, migration of O into the surfaces, and then oxide formation. However, the initial reaction mechanism has not been well understood, and the energy barrier of the initial oxidation are still under debate.

The adsorption and reaction of O2 on Cu surface has been studied for years by experimental technologies. The interaction of Ni/Fe-covered Cu (100) surfaces with O2 were studied by x-ray photoelectron spectroscopy (XPS). The authors pointed that the composition and coverage of the oxides of Fe and Ni are related to the concentration of O2 and the exposure time33 Kishi K, Nishioka J. Interaction of Fe/Cu (100), Fe-Ni/Cu (100) and Ni/Fe/Cu (100) surfaces with O2 studied by XPS. Surf Sci. 1990;227:97-106.. Then A. R. Balkenende and collaborators found the exposure of O2 to Cu (100) and Cu (110) surfaces leads to dissociative adsorption and the dissociation on Cu (110) surface is about an order of magnitude faster than on Cu (100) surface44 Balkenende AR, Daas HD, Huisman M, Gijzeman OLJ, Geus JW. The interaction of NO, O2 and CO with Cu (100) and Cu (110). Appl Surf Sci. 1991;47:341-53.. Dissociative chemisorption of O2 on Cu (100), S/Cu (100) and Ag/Cu (100) surface alloy were investigated by Auger electron spectroscopy (AES). It is concluded that at very low Ag coverages, the reduced reactivity of Ag/Cu (100) towards O2 dissociation is primarily due to the steric blocking of the surface defects and that any electronic effects are only secondary and present only at higher Ag coverages55 Hirsimäki M, Chorkendorff I. Effects of steps and defects on O2 dissociation on clean and modified Cu (100). Surf Sci. 2003;538:233-9.. An energy selective molecular beam surface scattering experiments revealed a defect induced low-barrier dissociation pathway leading to enhanced dissociation of O2 molecules on Cu surface with translational energy up to 60 meV66 Ahonen M, Hirsimäki M, Puisto A, Auvinen S, Valden M, Alatalo M. Adsorption dynamics of O2 on Cu (100): the role of vacancies, steps and adatims in dissociative chemisorption of O2. Chem Phys Lett. 2008;456:211-4.. These qualitative and quantitative results are mainly due to the analysis of energy oscillation spectroscopy, without directly observing the microscopic process of the initial stage of O2 adsorption and dissociation on Cu surface.

In contrast, theoretical simulation work is more suitable for directly obtaining microscopic information such as the structure, energy, transition state, and potential energy surface of O2 adsorption and dissociation on Cu surface. Because O atoms exhibit insignificant selectivity with respect to positions and coverages on Cu surface. For examples, based on density functional theory (DFT) calculations, people found that on Cu (211) surface the step microfacets are very reactive and the dissociation of the O2 molecule at room temperature occurs mostly on those sites77 Lahti M, Nivalainen N, Puisto A, Alatalo M. O2 dissociation on Pd (211) and Cu (211) surfaces. Surf Sci. 2007;601:3774-7.. First principles simulations based on DFT were used to determine absorption energies and pathways to dissociation of O2 at finite temperature on the Cu (110) surface. The results suggested that adsorption kinetics play an important role in determining the mechanism of dissociation88 Liem SY, Clarke JHR, Kresse G. Pathways to dissociation of O2 on Cu (110) surface: first principles simulations. Surf Sci. 2000;459:104-14.. Another work showed that the expansive strain parallel to the surface plane could enhance the binding of atomic and molecular oxygen on Cu (111) surface as well as to decrease the transition state energy of O2 dissociation99 Xu Y, Mavrikakis M. Adsorption and dissociation of O2 on Cu (111): Thermochemistry, reaction barrier and the effect of strain. Surf Sci. 2001;494:131-44.. However, the theoretical sticking probability of O2 dissociation on Cu (100) is found non-activated and topology of the potential energy surface is very open in entrance channels1010 Martin-Gondre L, Crespos C, Larrégaray P. Dynamics of dissociative chemisorption of O2 on Cu (100) surface: A theoretical study. Surf Sci. 2019;688:45-50.. The previous work mentioned above show that the micro-mechanism of dissociation and adsorption of O2 on Cu surface is a very important research topic.

Dissociation of O2 on Cu (111) has been described as both a direct process and one mediated by a molecular precursor. The steering effect plays an important role in the oxygen molecule dissociative process. The potential-energy hypersurface (PES) for O2 on Cu (111) is of key importance to assess dissociative phenomena. However, the surface doped elements are also an important factor affecting the dissociation of adsorbates. The doping atom could penetrate the substrate actively, change the surface activity and affect the dissociation of adsorbate on the surface1111 Yang B, Shi KY, Li HY, Jiang LW, Zhang CH. Water dissociative adsorption on the precovered Fe (100) surface from DFT computation. Indian J Phys. 2019;93:1019-29.

12 Zhang CH, Chen B, Sun DB. A DFT study of H2O dissociation on metal-precovered Fe (100) Surface. Surf Interface Anal. 2018;50:420-9.
-1313 Yu C, Liu JY, Lu H, Chen JM. Study of the effects of an adatom Sn on the Cu surface electromigration using a first principles method. Appl Surf Sci. 2007;253:8652-6.. Some metals, such as Al, Ag, W, are common elements in electrical engineering materials, and their presence has an important effect on the performance of Cu wires. Therefore, studying the microscopic process of O2 dissociation by the doping atoms on the Cu surface will provide a direct help in understanding the oxidation and protection mechanisms of Cu-based products.

2. Computational Details

First principles calculations based on density functional theory (DFT)1414 Hohenberg P, Kohn W. s-wave K-N scattering by the N/D method. Phys Rev. 1964;136:B864.,1515 Kohn W, Sham LJ. Self-consistent equations including exchange and correlation effects. Phys Rev. 1965;140:A1133. were used to investigate O2 adsorption and dissociation on the Al, Ag and W doped Cu (111) surface in the CASTEP code1616 Segall MD, Lindan PJD, Probert MJ. First-principles simulation: ideas, illustrations and the CASTEP code. J Phys Condens Matter. 2002;14:2717-44.. The ultrasoft pseudopotential1717 Blöchl PE. Projector augmented-wave method. Phys Rev B Condens Matter. 1994;50:17953-79.,1818 Kresse G, Joubert D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys Rev B Condens Matter Mater Phys. 1999;59:1758-75. was used to describe electron-ion interaction for all atoms. The generalized gradient approximation-Perdew-Wang 1991 (GGA-PW91)1919 Perdew JP, Burke K, Ernzerhof M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996;77:3865-68.,2020 Perdew JP, Wang Y. Accurate and simple analytic representation of the electron-gas correlation Energy. Phys Rev B Condens Matter. 1992;45:13244-9. was used to describe the exchange and correlation functional. In our calculations, the cutoff energy of 500 eV was employed for plane wave expansions. The Brillouin zone is sampled by using a 10 × 10 × 10 Monkhorst-Pack k-point grid2121 Monkhorst HJ, Pack JD. Special points for brillouin-zone integrations. Phys. Rev. B. 1976;13:5188. for the Cu bulk lattice constant. The energy difference, the residual force, the maximum stress and the displacement of each atom were set as 10-5 eV/atom, 0.03 eV/Å, 0.05 GPa and 0.001 Å during the optimization to ensure the convergence of the self-consistent field calculations2222 Vanderbilt D. Soft self-consistent pseudopotentials in a generalized eigenvalue formalism. Phys Rev B Condens Matter. 1990;41:7892-5.. The climbing image nudged elastic band (CI-NEB) Method2323 Laasonen K, Pasqurello A, Car R, Lee C, Vanderbilt D. Car-Parrinello molecular dynamics with Vanderbilt ultrasoft pseudopotentials. Phys Rev B Condens Matter. 1993;47:10142-53. was used to estimate transition state and final state structures.

In order to check the validity of our methodology, the lattice parameter of bulk Cu which has a three-dimension primitive cell with bcc symmetry was computed. Our calculated value (3.491 Å) has a difference less than 3.5% compared to the experimental value of 3.615 Å2424 Dinsdale AT. Sgte data for pure elements. Calphad. 1991;15:317-425. and theoretical value of 3.607 Å2525 Zhang CH, Han JJ, Huang S, Shen J. Chen’s Lattice inversion embedded-atom method for FCC Metal. Adv Mat Res. 2011;320:415-20.). The surface size was modelled by using super approach with periodic boundary conditions. A vacuum layer of 20 Å was found to be sufficient to prevent the interactions between the periodic images. We used the slabs consisting of 7-layer, and relaxed the first 3 layers including doping atoms and the bottom 4 layers. By using a 3 × 3 × 1 Monkhorst-Pack k-point grid, the surface structural relaxation and the total energy calculation were performed in our study. Besides, the p(3x3) super cell was used to study the properties of O2 adsorption and dissociation on the pristine/doped Cu (111) seeing in Figure 1.

Figure 1
The initial structures of O2 adsorption pristine/doped Cu (111) surface. (a) is the O2 adsorption pristine surface. (b) is the M-doped surface. (c) is O2 adsorption on TOP-type doped surface. (d) is the O2 adsorption on 1NN-type doped surface. The green, bronzing and blue particles denote the first, second, third layers of Cu atoms. The oxygen and doping atoms are shown in red and yellow color, respectively.

The adsorption energy (Eads) of adsorbate, the average binding energy (Eb) and surface energy (γ) of Cu surface, and the doping energy (Eim) of doping atom are calculated. The four kinds of energies are calculated in following equations 1-4.

E a d s = ( E t o t a l E s l a b E O 2 ) / 2 (1)
E b = ( E s l a b n m E C u m E m ) / n + m (2)
γ = ( E s l a b n m E b C u / 4 m E b m / 4 ) / 2 A (3)
E i m = E s l a b E s l a b C u E m (4)

where Etotal is the total energy of the surface system, Eslab is the energy of pristine/doped Cu surface, EslabCu is the energy of Cu surface with the doping atom removed, and n and m denote the number of Cu and metal atoms, respectively. The ECu and Em denote the atomic energy of Cu and metal atoms, respectively. EO2 is the energy of the isolated O2 in gas phase, while EbCu and Ebm are the bulk energy of fcc-type crystal of Cu and metals, respectively.

The energy barrier (EB) and reaction energy (ER) are defined in equation 5 and equation 6.

E B = E T S E I S (5)
E R = E F S E I S (6)

The EIS, ETS, and EFS are the energies of the initial state, transition state, and the final state of the configuration, respectively.

3. Results and Discussions

3. 1. Al, Ag, W doped Cu (111) surface

The effect of doping atoms (Al, Ag, W) on the surface binding and stability are studied firstly. We replace one Cu atom of the first layer with the metal atom to get a M-doped (M= Al, Ag, W) Cu (111) surface. We optimize the stoichiometric geometries of M-doped Cu (111) surface until the external and internal degree of freedom can relax to the force and stress vanished in the bulk calculations. The lattice constants change after metal atoms have been doped on the Cu (111) surface. The relaxed structure of the M-doped Cu (111) surfaces are shown in Figure 2. All the doping atoms can be stably combined with Cu atoms, while being slightly embedded in the surface to a certain depth. The new lattice constants for M-doped Cu surfaces are 7.617 Å, 7.632 Å, 7.638 Å, respectively, which are all slightly larger than those of pristine Cu (111) surface (7.611 Å) because of larger radius of the doping atoms. The average bond length between Cu-Cu atoms in the first layer of pristine surface is 2.537 Å. While the average bond length between Al-Cu, Ag-Cu, W-Cu atoms in the first layer of doped surfaces are 2.543 Å, 2.582 Å, 2.554 Å, respectively, which means the doping atoms make the surrounding Al atoms slightly excluding to them. This repulsion effect enhances the bonding between surface atoms and improves the energy stability of the doped Cu surfaces.

Figure 2
The top view of optimized structure of the M-doped Cu (111) surfaces. The purple, dark brown and white particles denote the Al, Ag and W atoms.

The doping energy (Eim), binding energy (Eb) and surface energy (γ) of pristine/doped Cu surface are listed in Table 1. It is clear seen that all the doping atoms can be stably combined with Cu atoms due to the negative values of Eim. The W-doped surface has the lowest Eim with -7.614 eV than Al-doped and Ag-doped surfaces because W atom has more valence electrons that can be transferred to and from Cu atoms. Then we compared and analyzed the change of Eb and γ between the doped surface and the pristine surface. In this paper, the Eb reflects the strength of the bonding between atoms in the slab. The Eb of pristine Al surface is -3.628 eV and lower than the Eb of Al-doped and W-doped surfaces, so the binding energy data confirms that the Al and W atoms will strengthen the bonding effect between the atoms on the Cu (111) slab. But on the Ag-doped surface, the result is exactly the opposite. In addition, the γ of pristine Al surface (3.6451 J/m2) is lower than Al-doped and W-doped surfaces. That means Al and W atoms improves the surface activity of Cu surface, which may be more favorable for small molecules. In the same situation, Ag-doped surface has lower γ than the pristine Cu surface. Based on above results, Ag atom reduces the surface activity of Cu and inhibit the reaction between Cu surface and oxygen molecules.

Table 1
The calculated the doping energy (Eim), binding energy (Eb), adsorption energy (Eads) and surface energy (γ) of pristine/doped and O2 adsorption Cu surface.

3.2. O2 adsorption on pristine/doped Cu (111) surface

It is known that O2 prefer to adsorb at the top site of metal surface. In this work, there are one site for O2 adsorption on pristine Cu (111) surface and two sites for O2 adsorption on doped surface. Figure 3a shows the optimized O2 adsorption on pristine Cu surface, and Figure 3b are the optimized O2 adsorption on the top site of doped Cu surfaces (TOP-type) which named Al_TOP, Ag_TOP, W_TOP, respectively, and Figure 3c are the optimized O2 adsorption on the 1st near-neighbor top site of doped Cu surface (1NN-type) which named Al_1NN, Ag_1NN, W_1NN, respectively. For all kinds of doped surfaces, it should be pointed out that the O2 molecules prefer to adsorb in a parallel to the TOP-type surfaces, and they adsorb with a little tilt angle related to the 1NN-type surfaces. This difference is mainly caused by different atomic radii and different surface flatness under the adsorbate.

Figure 3
The top view of optimized structure of O2 adsorption on the pristine/doped Cu surfaces. The white, maroon and brown particles denote the Al, Ag and W atoms.

We calculated the adsorption energy (Eads) of O2 on the pristine/doped surfaces to obtain the most stable adsorption configurations. As is shown in Table 1, for the M-doped surfaces, we found that most TOP-type surfaces are more stable than pristine surface for O2 adsorption. The Eads of Al_TOP, W_TOP and W_1NN surfaces are -1.011 eV, -2.175 eV and -1.609 eV, which are still much lower than pristine surface (-0.589 eV). The reason is that the distance between adsorbates and top-site Cu is smaller than these of top-site Al and W atoms. In another word, Al and W atoms embedded in the surface increase the internal binding of the system, and weaken the adsorption effect with the adsorbate. In addition, the O atom prefers to adsorb on the top of the doping atom rather than on the top site of Cu atom, because the 1NN-type surfaces have higher values of Eads than TOP-type surfaces. However, the results of Ag-doped surfaces are opposite to those of other surfaces. Because of the inhibition effect of Ag atoms, the Eads on both types of Ag-doped surfaces is lower than that of pristine surface. Besides that, we also computed the Eim, Eb and γ of O2 adsorption on the M-doped Cu surfaces and shown the data in Table 1. After O2 adsorption, the values of Eim and Eb of pristine/doped surfaces are lower than before except Ag-doped surfaces, but the values of γ of pristine/doped surfaces are higher than the clean surfaces except W-doped surfaces. The reason for our analysis is that the adsorption of O2 causes the doping atoms to be further embedded in the surface, which increases the doping energy. This conclusion is consistent with the result of adsorption energy we mentioned above. As the doping atoms continue to penetrate into the Cu surface, the space between the Cu atoms is reduced, thereby increasing the binding energy inside the material. The reason for the increase in surface energy is also due to the presence of adsorbates. However, since the radius of Ag atom is much larger than that of Cu atom, the obvious repulsion effect causes the weak bonding between surface atoms to no longer follow the above rules.

To further explore the interaction between the substrate, doping atom and adsorbate, we introduced the electron density of states (DOS) for O2 adsorbed pristine/doped Cu surfaces. As shown in Figure 4, we plotted the DOS of one doping atoms, two oxygen atoms and two kinds of Cu atoms, one is the Cu atom just below the adsorbed O atoms and the other one is the nearest neighbor to both the doping atom and the adsorbed O atom. We can see that the contribution of Cu-d elctorns to the DOS is obviously large, although the Cu-s and Cu-p electrons have hybridization and little contributions in the energy range of -10 eV to + 10 eV. Among them, the number of d-electrons of Cu atoms on the pristine surface is higher than that of O adsorbed surfaces. In another word, the Cu-d orbital electrons in the former are concentrated in the high energy region (-5 eV to 0 eV), while O-p orbitals in the latter are concentrated in wider energy level range of -7 eV to +2 eV. At the same time, the O-s and O-p orbital electrons have obvious peaks in the energy range of -7 eV to + 10eV, which indicates that after O adsorption, an interaction occurs between O-Cu, and some electrons are transferred from Cu atom to O atom, which in turn causes the remaining electrons of Cu atom to move to a lower energy level and improves the bonding between the inner atoms in the Cu surface, and finally shows the effect on the value of the binding energy. However, O atoms and doping atoms contribute significantly to the DOS, especially the O-p orbital electrons and the d orbital electrons of doping atoms. Due to the different electronegativity of three metals Al, Ag, W, the hybridization rules of the electron orbitals on the three doped surfaces are different. The Al-p orbital peaks and O-p orbital peaks have overlap in a wide energy range, but there are nearly not hybridization peaks between Cu-d and Al-p orbits. The hybridization of the peaks of the Ag-d orbit and O-p orbit is affected by the O adsorbed sites and appears to overlap at high and low levels, respectively. But Ag-d orbit has the highest peak at low energy level of all the three doping atoms. The W-d orbital and O-p orbital have a broad peak hybridization only at high energy levels, but the number of electrons of these overlapping peaks is the smallest of the three doped metals. The analysis of the electronic structure shows that the difference in energy calculated in Table 1 is not only influenced by the orbital hybridization between the adsorbate and the surface atoms, but also related to the electronegativity and the orbital energy level of the doping atoms. In one sense, our results provide fundamental insights into O2 impossible dissociative adsorption on these M-doped Cu surfaces.

Figure 4
The DOS of O2 adsorbed the pristine/doped Cu surfaces.

3.3. O2 dissociation on pristine/doped Cu (111) surface

In this section, we mainly investigate the O2 dissociative processes (O2 → O +O) on the pristine/doped Cu surface. We computed the initial state (IS) and final state (FS) in order to obtain the potential energy surface (PES). Then we obtained the intermediate states using a linear interpolation. At last, we used the same convergence criteria of force and energy as the whole simulation process. The energy barrier (EB) and reaction energy (ER) and structure parameters of O2 are listed in Figure 5 and Table 2.

Figure 5
The potential energy surface for O2 dissociation on the pristine/doped Cu surfaces. The initial state (IS), transition state (TS) and final state (FS) are shown as solid lines.
Table 2
The bond length of O atoms (LO-O), the energy barrier (EB) and reaction energy (ER) of O2 dissociation on the pristine/doped Cu surfaces.

It is seen that the degree of O2 dissociation is related to the doping atoms on Cu surfaces. This result is as similar as our previous calculations1111 Yang B, Shi KY, Li HY, Jiang LW, Zhang CH. Water dissociative adsorption on the precovered Fe (100) surface from DFT computation. Indian J Phys. 2019;93:1019-29.,1212 Zhang CH, Chen B, Sun DB. A DFT study of H2O dissociation on metal-precovered Fe (100) Surface. Surf Interface Anal. 2018;50:420-9. that pointed that the precovered-atoms affect the subsequent H2O dissociation. The EB is 0.405 eV on the pristine Cu surface. Except W_TOP surface, the value is lower than the EB on other M-doped surfaces. It shows that the doping atom can resist the dissociation of O2. Further studying the data, we found that the EB on three kinds of Al_Top, Ag_1NN and W_1NN surfaces are 0.505 eV, 0.905 eV and 0.814 eV, respectively, which are much larger than other surface models. This seems to indicate that the adsorbed O2 in a certain range near the doping atoms has good effect of accelerated dissociation. However, by analyzing the reaction energy, we found that the dissociation reaction of O2 on the surfaces is not only reflected in the barrier energy, but the reaction energy can be also an important index. The ER are -5.355 eV, -6.426 eV, -5.771 eV, -6.113 eV, -5.386 eV, -7.191 eV, -5.356 eV for pristine, Al_Top, Al_1NN, Ag_Top, Ag_1NN, W_Top, W_1NN surfaces, respectively. The results show that the O2 can be dissociated more easily on the M_1NN surfaces than that of M_Top surfaces. We found that the reaction energy of O2 on the pristine Cu surface is larger than that of all doped surfaces, and the ER of O2 at the M_TOP surface is also low. For the previous results, we believe that the doping Ag atoms actually resist the dissociation of O2. The results of bond length (LO-O) also prove our inference, because the value of LO-O for Ag doped surfaces are smaller than that on the pristine Cu surface. Therefore, Ag atoms have the best resistance effect on the O2 dissociation comparing with Al and W atoms because of the large barrier energy and reaction energy.

4. Conclusions

The O2 adsorption and dissociation on the pristine/doped Cu (111) surface were studied by density functional theory. The Al, Ag and W atoms were doped on the surface to study the interaction between O2 and Al. All the doped atoms could be stably combined with Cu atoms, while being slightly embedded in the surface to a certain depth. The Eb of pristine Al surface is lower than that of Al-doped and W-doped surfaces, so the binding energy data confirms that the Al and W atoms would strengthen the bonding effect between the atoms on the Cu (111) slab except Ag doping atom. The O atom prefers to adsorb on the top of the doping atom rather than on the top site of Cu atom. However, Ag atom reduces the surface activity of Cu and inhibit the reaction between Cu surface and oxygen molecules. Due to the different electronegativity of three metals Al, Ag, W, the hybridization rules of DOS for O2 adsorbed pristine/doped Cu surfaces are different. After O adsorption, an interaction occurs between O-Cu, and some electrons are transferred from Cu atom to O atom, but the doping atom can resist the dissociation of O2. Moreover, we computed the initial state and final state in order to obtain the PES. The result shows that the dissociation reaction of O2 on the surfaces not only reflected in the barrier energy, but the reaction energy can be also an important index. Ag atoms have the best resistance effect on the O2 dissociation comparing with Al and W atoms because of the large barrier energy and reaction energy.

5. Acknowledgements

The authors gratefully acknowledge the financial support of the Science and Technology Foundation of State Grid Corporation of China (Contract No. SGZJ0000KXJS1800302): Research on metallic material selection, manufacturing, installation process for power transmission and transformation equipment and nondestructive detection method for components.

6. References

  • 1
    Zhao J, Tao F. Dissociation mechanism of H2O on clean and oxygen-covered Cu (112) surfaces, A theoretical study. Vacuum. 2016;128:252-8.
  • 2
    Akca A, Genc AE, Kutlu B. BH4 dissociation on various metal (111) surfaces: A DFT study. Appl Surf Sci. 2019;473:681-92.
  • 3
    Kishi K, Nishioka J. Interaction of Fe/Cu (100), Fe-Ni/Cu (100) and Ni/Fe/Cu (100) surfaces with O2 studied by XPS. Surf Sci. 1990;227:97-106.
  • 4
    Balkenende AR, Daas HD, Huisman M, Gijzeman OLJ, Geus JW. The interaction of NO, O2 and CO with Cu (100) and Cu (110). Appl Surf Sci. 1991;47:341-53.
  • 5
    Hirsimäki M, Chorkendorff I. Effects of steps and defects on O2 dissociation on clean and modified Cu (100). Surf Sci. 2003;538:233-9.
  • 6
    Ahonen M, Hirsimäki M, Puisto A, Auvinen S, Valden M, Alatalo M. Adsorption dynamics of O2 on Cu (100): the role of vacancies, steps and adatims in dissociative chemisorption of O2 Chem Phys Lett. 2008;456:211-4.
  • 7
    Lahti M, Nivalainen N, Puisto A, Alatalo M. O2 dissociation on Pd (211) and Cu (211) surfaces. Surf Sci. 2007;601:3774-7.
  • 8
    Liem SY, Clarke JHR, Kresse G. Pathways to dissociation of O2 on Cu (110) surface: first principles simulations. Surf Sci. 2000;459:104-14.
  • 9
    Xu Y, Mavrikakis M. Adsorption and dissociation of O2 on Cu (111): Thermochemistry, reaction barrier and the effect of strain. Surf Sci. 2001;494:131-44.
  • 10
    Martin-Gondre L, Crespos C, Larrégaray P. Dynamics of dissociative chemisorption of O2 on Cu (100) surface: A theoretical study. Surf Sci. 2019;688:45-50.
  • 11
    Yang B, Shi KY, Li HY, Jiang LW, Zhang CH. Water dissociative adsorption on the precovered Fe (100) surface from DFT computation. Indian J Phys. 2019;93:1019-29.
  • 12
    Zhang CH, Chen B, Sun DB. A DFT study of H2O dissociation on metal-precovered Fe (100) Surface. Surf Interface Anal. 2018;50:420-9.
  • 13
    Yu C, Liu JY, Lu H, Chen JM. Study of the effects of an adatom Sn on the Cu surface electromigration using a first principles method. Appl Surf Sci. 2007;253:8652-6.
  • 14
    Hohenberg P, Kohn W. s-wave K-N scattering by the N/D method. Phys Rev. 1964;136:B864.
  • 15
    Kohn W, Sham LJ. Self-consistent equations including exchange and correlation effects. Phys Rev. 1965;140:A1133.
  • 16
    Segall MD, Lindan PJD, Probert MJ. First-principles simulation: ideas, illustrations and the CASTEP code. J Phys Condens Matter. 2002;14:2717-44.
  • 17
    Blöchl PE. Projector augmented-wave method. Phys Rev B Condens Matter. 1994;50:17953-79.
  • 18
    Kresse G, Joubert D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys Rev B Condens Matter Mater Phys. 1999;59:1758-75.
  • 19
    Perdew JP, Burke K, Ernzerhof M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996;77:3865-68.
  • 20
    Perdew JP, Wang Y. Accurate and simple analytic representation of the electron-gas correlation Energy. Phys Rev B Condens Matter. 1992;45:13244-9.
  • 21
    Monkhorst HJ, Pack JD. Special points for brillouin-zone integrations. Phys. Rev. B. 1976;13:5188.
  • 22
    Vanderbilt D. Soft self-consistent pseudopotentials in a generalized eigenvalue formalism. Phys Rev B Condens Matter. 1990;41:7892-5.
  • 23
    Laasonen K, Pasqurello A, Car R, Lee C, Vanderbilt D. Car-Parrinello molecular dynamics with Vanderbilt ultrasoft pseudopotentials. Phys Rev B Condens Matter. 1993;47:10142-53.
  • 24
    Dinsdale AT. Sgte data for pure elements. Calphad. 1991;15:317-425.
  • 25
    Zhang CH, Han JJ, Huang S, Shen J. Chen’s Lattice inversion embedded-atom method for FCC Metal. Adv Mat Res. 2011;320:415-20.

Publication Dates

  • Publication in this collection
    29 Mar 2021
  • Date of issue
    2021

History

  • Received
    30 Mar 2020
  • Reviewed
    25 May 2020
  • Accepted
    01 June 2020
ABM, ABC, ABPol UFSCar - Dep. de Engenharia de Materiais, Rod. Washington Luiz, km 235, 13565-905 - São Carlos - SP- Brasil. Tel (55 16) 3351-9487 - São Carlos - SP - Brazil
E-mail: pessan@ufscar.br