search
for
 About Bioline  All Journals  Testimonials  Membership  News


Journal of Applied Sciences and Environmental Management
World Bank assisted National Agricultural Research Project (NARP) - University of Port Harcourt
ISSN: 1119-8362
Vol. 10, Num. 1, 2006, pp. 37-42

Journal of Applied Sciences & Environmental Management, Vol. 10, No. 1, March, 2006, pp. 37-42

Surface relaxation and surface energy of face-centered Cubic metals

1AGHEMENLO HE; *2IYAYI, S E; 3AVWIRI,G O

1, 3 Department of Physics, Ambrose Alli University, Ekpoma, Nigeria
2 Department of Physics, University of Benin, Benin City, Nigeria
3 Department of Physics, University of Port Harcourt, PH, Nigeria: E-mail: Sunnyiyayi@yahoo.com

*Corresponding author: E-mail: Sunnyiyayi@yahoo.com

Code Number: ja06006 

ABSTRACT: Surface relaxation and surface energy are calculated for Rh, Ir, Pb, Ca, Sr, and Th, for three low-index surfaces, using the equivalent crystal theory (ECT) method. The surfaces are the (100), (110) and (111) faces. In this study, we only considered perpendicular relaxation where the relaxation was determined by the minimization of the total energy. The results obtained are analyzed with good agreement to experiment and first- principles calculations where available. .@JASEM       

The phenomena of reconstruction and relaxation of metal surfaces has received considerable attention. Considerable experimental data indicate multilayer oscillatory relaxation of surfaces. Two methods, the Embedded Atom Method (EAM), based on density functional theory, and Equivalent Crystal Theory, ECT (Smith et al 1991), based on perturbation theory, have been extensively used to describe the energetics of defects in metals. Both methods have also been applied to the study of surface relaxation (Smith et al 1991). Surface energy is another surface property of great interest. The energy of a free surface plays an important role in several physical and chemical processes such as fracture, catalysis, etc. Experimental measurements of the surface energy are usually at high temperatures and are subject to errors due to surface – active contaminants and thus have a degree of uncertainty. Early theoretical calculations were based on perturbation theory or non-perturbative variational methods. Over the years, there has been an increasing effort on first-principles calculations as well as in the area of Semi-empirical methods: EAM and ECT have been applied to this and other surface properties (Smith et al 1991), with the latter method providing excellent agreement with available experimental data and first-principles calculations. Finally the work of Metheffessel, Hennig and Scheffler (1992) discusses trends in surface energies as well as different models that relate the surface energies to other crystal properties.

In this paper, we extend our previous application of ECT to surface energy calculation (Aghemenloh and Idiodi 1998), and in order to improve on our previous calculated results; we here include the effect of relaxation in the calculated surface energies. Several studies of surface phenomena have been done with the ECT on Al, Ni, Pd, Cu, Pt, Ag and Au (Smith et al 1991; Rodriguez et al.1993) obtaining remarkable agreement with experimental results. At the present time none has been performed to the author’s knowledge for Rh, Ir, Pb, Ca, Sr and Th metal surfaces. To fill this gap, we present in this paper the first ECT calculation of surface energy and surface relaxation for the three low-index faces (100), (110) and (111) of these metal surfaces. In this study, we ignore parallel relaxation as well as surface reconstruction and dealt only with rigid perpendicular translations of the near-surface planes. With these restrictions, we apply ECT and discuss the different contributions to the surface energy. 

Equivalent crystal theory: Equivalent crystal theory is based on an exact relationship between the total energy and atomic locations and applies to surfaces and defects in both simple and transition metals as well as in covalent solids. Lattice defects and surface energies are determined via perturbation theory on a fictitious, equivalent single crystal whose lattice constant is chosen to minimize the perturbation. The energy of the equivalent crystal, as a function of its lattice constant is given by a Universal binding energy relation (Smith et al 1991). The method has been applied previously to calculate surface energies (Rodriguez et al 1993; Aghemeloh and Idiodi 1998), surface reconstructions and bulk distortions of metals and semiconductors (Smith et al 1991). Since the details of the ECT method may be found in other previous works (Smith et al 1991; Rodriguez et al 1993; Aghemenloh and Idiodi 1998), we shall here restrict ourselves to the calculation of relaxed surface energies only. As given in section three of this work. 

Surface energy calculation: In a previous study (Aghemenloh and Idiodi 1998), we considered a rigid surface where no interlayer relaxation was allowed, hence all bond lengths and angles, retained their bulk equilibrium values, thus, F*(a 2 *) =F*(a 3 *) = F*(a 4 *) = 0. The surface energy was therefore obtained by solving for the ‘volume’ term represented by F*(a 1 *) only. If we now consider the simplest surface relaxation, that is, a rigid displacement of the surface layer towards the bulk, the contribution from the higher order anisotropic terms must be included. Letting the planes close to the surface relax turns on the other three contributions to the energy. Some bonds are compressed, contributing to F*(a 2 *), the bond angles near the surface as well as the difference between face diagonal are distorted, generating an increase of energy via F*(a 3 *) and F*(a 4 *). Also, allowing the atoms in the surface region to relax introduces the additional complexity of including in the calculation more nonequivalent atoms, located in deeper layers. Many – atom effects, which are represented in the ECT by the inclusion of the three – atom bond – angle anisotropy and the four – atom face diagonal distortion terms (see the third and fourth term in the RHS of Eq. (2) of Rodriguez et al 1993) are necessary but, in the case of surface energy calculations of metals, of very little relevance. As have been shown in other previous ECT works, they introduce a small correction, usually of the order of 1% of the leading term in Eq. (2) of Rodriguez et al (1993). This is however not the case for semiconductors, where angular anisotropies are a significant contribution to the surface energy.

Calculation of relaxed surface energy for (100) surface: This subsection provides an application of the ECT method to an fcc (100) metallic surface, where only the surface plane is allowed to relax. All atoms in a given plane are identical; therefore we only need to evaluate the contribution of a single nonequivalent atom per plane. Moreover, only two to three planes have to be considered in this calculation: atoms in the fourth plane and below find themselves in equilibrium, bulk-like environment. The data by Bozzolo et al (1993) shows the (100) surface of Fcc metals. Let assume that the top plane (j = 1) is lightly contracted by a distance x towards the first plane below the surface (j = 2), but both the j = 2 and j = 3 planes as well as any other plane below are at equilibrium positions. The rigid inter-planar spacing is d =a/2 (a is the equilibrium lattice parameter); therefore, the distance between j = 1 and j = 2 is a/2-x. Following the labeling of Bozzolo et al 1993, we write the ECT equations for the atoms A1, B1 and C1 in terms of the distances between atoms of nearest and next nearest neighbors:

Equations (1-3) are then solved, for each value of χ , for the equivalent crystal nearest neighbor distance Rec.

Finally, the ‘volume’ contribution to the surface energy is

A is the surface area and F* is given by Eq. (3) of Rodriguez et al. (1993).

Next, we compute the bond length anisotropy contribution to the surface energy (for atoms A1 and B1) associated with this defect:

The equation for atom B1, for this particular case is identical to Eq. (10). Eq. (10) is then solved with respect to Rec and the energy contribution is then

where a2* is given by Eq. (7) of Rodriguez et al (1993)

The third term in the ECT expansion deals with bond – angle anisotropies. Contributions to σ3 comes from atoms for which the angle between the nearest neighbours departs from its equilibrium value θoo = 90o for fcc metals). Following the conversion described by Smith et al (1991), i.e., if an atom is missing one or more nearest neighbours then F*(a3*(i,j,k))=0 for that atom and the atoms in the top layer do not contribute to σ3 .The only contribution then arises from the atoms in the first layer below the surface (j = 2), for which we solve the transcendental equation

                    

where

                                                                                                                  

Equation (12) is solved with respect to Rec and the bond – angle contribution is then

      

Summarizing, the surface energy of an fcc (100) face, where only the top layer is allowed to relax, is given by

The calculation of the contribution to the defect energy from a4* has been ignored in this study. Equations similar to Eq. (1-15) can also be written down for the (111) - (110) – faces. Such details are ignored here.

RESULTS AND DISCUSSIONS

In this section we report surface energies for six fcc metals. The ECT method requires as experimental input, the three bulk elastic constants, the bulk cohesive energy, the bulk lattice constants and the vacancy formation energy. These are collected together in Table 1 above and they have been employed to compute the ECT constants in Table 2. The different contributions of surface energy arising from the ECT expansion are displayed in Table 3, the surface energy are both for the rigid and relaxed case. From Table 3, it can be seen that the contributions from the bond – angle term s3, are very small when compared to the ‘volume’ and bond – compression term. It has been shown else where (Rodriguez et al. 1993; Bozzolo et al. 1993), that the bond-angle and the face – diagonal term of the ECT represents only 1 to 2% of the total energy, hence they are not of great relevance for the calculation of surface energies.     

The surface energies obtained from the three low – index faces of the six fcc metals are summarized in Table 4, and are compared with the results from first – principles calculations (Methfessel et al 1992; Skriver and Rosengaard 1992), the Tight – binding (TB) method (Mehl and Papconstantopoulos 1996), the modified embedded atom method (MEAM) (Baskes 1992), and experiment (de Boer et al 1988). In all cases we find σ111< σ100110. Thus close packed surfaces are the most stable for the fcc metals. The ECT surface energies are uniformly larger and closer to experiment, than those obtained by TB. Our surface energies are generally closer to experiment than those obtained by the MEAM and in good agreement with first – principles calculations for the surface energy of Ca (111) and Sr (Skriver and Rosengaard 1992). Our results for the top-layer relaxations are presented in Table 5, and are compared with the results of other theoretical calculations and experiment of which the authors are aware. Table 5 shows that our ECT values are in good agreement with the first– principles calculation of Methfessel et al (1992) for the (100) and (111) faces of Rh. The ECT value of the (110) face of Rh are also in better agreement with experiment (Begley et al 1993) than the theoretical calculations of Sinnott et al. (1991) and Methfessel et al (1992).

Conclusion: We have in this study extended the surface energy results of Aghemenloh and Idiodi (1998), by including the effect of relaxation on the calculated surface energy of Rh, Ir, Pb, Ca, Sr and Th. Very generally, the surface energies are found to be in good agreement with the results from first principles calculation and experiment. Our work on relaxation is good, as it gives the correct sign and magnitude of the top-layer relaxations where experimental and first – principles calculations are available.               

Acknowledgement: We gratefully acknowledge enlightening communication with G. Bozzolo of Ohio Aerospace Institute, Cleveland, Ohio and J. Ferrante of Cleveland State University, Ohio, U.S.A.

REFERENCES

  • Aghemenloh, E; Idiodi, J O A (1998). Equivalent crystal theory of Fcc metal surfaces. J. Nig. Ass. Math. Phys. Vol. 2, 271-284.
  • Baskes, M I (1992). Modified embedded-atom potentials for cubic materials and impurities. Phys. Rev. B46, 2727-2742
  • Begley, A.M; Kim, S K; Jona, F; Marcus, P M (1993). Surface relaxation of Rh (001). Phys. Rev. B48, 12326-12329.
  • Bozzolo, G; Ferrante, J; Rodriguez, A M (1993). Modelling of surfaces.1. Monatomic metallic surfaces using equivalent crystal theory. J. Comput – Aided Mater. Design. 1, 285-304 
  • De’ Boer, F R; Boom, R; Mattens, W C M; Miedema, A R; Niessen, A K (1988).  Cohesion in Metals (Eds. F.R. deBoer and D.G. Pettifor) Vol.1, North – Holland Amsterdam.
  • Mehl, M J; Papaconstantopoulos, D A (1996).  Application of a tight-binding total-energy method for transition and noble metals: Elastic constants, vacancies, and surfaces of monatomic metals. Phys. Rev. B54, 4519-4530
  • Methfessel, M; Hennig, D; Scheffler, M (1992). Trends of the surface relaxations, surface energies, and work functions of the 4d transition metals.  Phys. Rev. B46, 4816-4829
  • Rodrguez, A M; Bozzolo, G; Ferrante, J (1993). Multilayer relaxation and surface energy of fcc and bcc metals. Surf. Sci. 289, 100-126
  • Smith, J R; Perry, T; Banerjea, A; Ferrante , J; Bozolo, G (1991). Equivalent crystal theory of metal and semiconductors surfaces and defects. Phys. Rev. B44, 6444.-6465
  • Sinnott, S B; Stave, M S; Raeker, T J;  Depristo, A E (1991). Corrected effective-medium study of metal-surface relaxation. Phys. Rev. B44, 8927-8941
  • Skriver, H L; Rosengaard, M N (1992). Surface energy and work function of elemental metals. Phys. Rev. B46, 7157-7168

Copyright 2006 - Journal of Applied Sciences & Environmental Management


The following images related to this document are available:

Photo images

[ja06006t5.jpg] [ja06006t2.jpg] [ja06006t4.jpg] [ja06006t3.jpg] [ja06006t1.jpg]
Home Faq Resources Email Bioline
© Bioline International, 1989 - 2024, Site last up-dated on 01-Sep-2022.
Site created and maintained by the Reference Center on Environmental Information, CRIA, Brazil
System hosted by the Google Cloud Platform, GCP, Brazil