The coefficients of thermal expansion of thin metal film with face-centered cubic structure at zero pressure are investigated using the statistical moment method (SMM), including the anharmonic effects of thermal lattice vibration. The Helmholtz free energy, mean-square atomic displacement and linear thermal expansion coefficients are derived in closed analytic forms in terms of the power moments of the atomic displacement.
JOURNAL OF SCIENCE OF HNUE Mathematical and Physical Sci., 2013, Vol 58, No 7, pp 109-116 This paper is available online at http://stdb.hnue.edu.vn COEFFICIENTS OF THERMAL EXPANSION OF THIN METAL FILMS INVESTIGATED USING THE STATISTICAL MOMENT METHOD Duong Dai Phuong1 , Vu Van Hung2 and Nguyen Thi Hoa3 Tank Armour Officers Training School, Tam Duong, Vinh Phuc Viet Nam Education Publishing House, Hanoi Hanoi University of Transport and Communications Abstract The coefficients of thermal expansion of thin metal film with face-centered cubic structure at zero pressure are investigated using the statistical moment method (SMM), including the anharmonic effects of thermal lattice vibration The Helmholtz free energy, mean-square atomic displacement and linear thermal expansion coefficients are derived in closed analytic forms in terms of the power moments of the atomic displacement Numerical calculations for Al, Au and Ag thin films are found to be in good and reasonable agreement with other theoretical results and with the experimental data Keywords: Coefficients of thermal expansion, statistical moment method, anharmonicity, thin film Introduction Materials in the form of thin film have come under the spotlight in recent years They have been found to show physical, chemical and mechanical properties that differ from corresponding bulk materials [3, 4] Properties, limitations and advantages obtained on thin film geometry have been widely studied and they were found to depend on different factors: structure, size, film thickness and different substrates [3, 5] Thin film is used in a vast range of applications which include microelectromechanical and nanoelectrome-chanical systems, sensors and electronic textiles [1] Metallic thin film displays a common geometry and presents enormous scientific interest, mainly due to their attractive and novel properties for technological applications [1] In the early twentieth century, the theory of size-dependent effects in thin metal layers was developed by a number of scientists [5] For ultrathin metal layers, both surface effects and quantum size effects must be considered [8] A large number of experimental and Received September 20, 2013 Accepted October 30, 2013 Contact Vu Van Hung, e-mail address: bangvu57@yahoo.com 109 Duong Dai phuong, Vu Van Hung and Nguyen Thi Hoa theoretical studies have been carried out on the thermal behavior of metal and nonmetal thin films Most of the previous theoretical studies, however, are concerned with the material properties of thin metal film at low temperature, and the temperature dependence of thermodynamic quantities has not been studied in detail The purpose of the present article is to investigate the temperature dependence and the thickness dependence of the coefficients of thermal expansion of thin metal film with face-centered cubic structure using the analytic statistical moment method (SMM) [6, 10].The coefficients for thermal expansion are derived from the Helmholtz free energy, and the explicit expression of the thermal expansion coefficient is presented taking into account the anharmonic effects of the thermal lattice vibrations In the present study, the influence of surface and size effects on the coefficients of thermal expansion have also been studied We compared the results of the present calculations with those of the previous theoretical calculations as well as with the available experimental results Content 2.1 Theory 2.1.1 The anharmonic oscillations of thin metal film Investigating a metallic free standing thin film includes n layers with film thickness d Suppose of the thin film consists of two surface layers of atoms, the next two deeper layers of atoms and the (n-4) internal layers atoms (Fig 1) Figure The thin metal film 110 Coefficients of thermal expansion of thin metal films investigated First of all, we present the statistical moment method (SMM [6, 10]) formulation > as a solution of for the displacement of the surface layer atoms of the thin film < usur i the equation γsur θ2 d2 i dp2 d i dp + 3γsur θ < usur > i +γsur θ ksur + γsur < usur >3 +ksur < usur > i i (xsur coth xsur − 1) < usur > −p = 0, i (2.1) ωsur ; θ = kB T, 2θ (2.2) where xsur = ksur γ1sur γsur 1∑ = i ( ∂ φin io ∂u2iα ) , ≡ m0 ωsur (2.3) eq ( ) ∑ ∑ ∂ φsur io , γ = = 2sur 48 i ∂u4iα eq 48 i ( ∂ φsur io ∂uiβ ∂u2iγ ) ( ) ( sur ) sur ∑ ∂ φio ∂ φio = (γ1sur + γ2sur ) ; +6 = 12 i ∂uiα eq ∂u2iβ ∂u2iγ (2.4) eq α, β, γ = x, y, z eq (2.5) Here, KB is the Boltzmann constant; T is the absolute temperature; m0 is the mass of the atoms at the lattice node; ωsur is the frequency of lattice vibration of surface layer atoms; ksur , γ1sur , γ2sur and γsur are the parameters of crystals depending on the structure of crystal lattice and the interaction potential between atoms at the node; is the effective interatomic potential between 0th and ith surface layer atoms; φsur i0 , uiy and uiz are the displacement of ith atom from the equilibrium position on direction x, y, z respectively; and the subscript eq indicates the evaluation at equilibrium To determine the atomic displacement < usur >, the symmetry property i appropriate for cubic crystals is used sur sur sur > < usur ix >=< uiy >=< uiz >=< ui (2.6) Then, the solutions of the nonlinear differential equation of Eq (2.1) can be expanded in the power series of the supplemental force p as < usur >=< usur >0 +A1 p + A2 p2 i i (2.7) 111 Duong Dai phuong, Vu Van Hung and Nguyen Thi Hoa Here, < usur >0 is the average atomic displacement in the limit of zero of i supplemental force p Substituting the above solution of Eq (2.7) into the original differential Eq (2.1), one can get the coupled equations on the coefficients A1 and A2 , from which the solution of < usur >0 is given as [6, 10] i √ < usur >0 ≈ i 2γsur θ2 Asur , 3ksur (2.8) where Asur = asur + γsur θ2 sur γsur θ3 sur γsur θ4 sur γsur θ5 sur γsur θ6 sur a + a + a + a + a6 , (2.9) 10 12 ksur ksur ksur ksur ksur with asur (η = 1, 2, , 6) being the values of parameters of crystals depending on the η structure of the crystal lattice [6, 10] A similar derivation can be also done for atoms at layers just below the surface and internal layers of atoms of thin film, their displacement noted as < usur1 >0 and < uin i i >0 , respectively 2.1.2 Free energy of the thin metal film In the present research, the influence of the size effect on the coefficients of thermal expansion of the thin metal film is studied, introducing the surface energy contribution in the free energy of the system atoms Following the general formula in the SMM formulation [6, 10], one can get the free energy of the surface layers atoms as { ( )} θ2 Ψsur ≈ {U0sur + 3Nsur θ[xsur + ln(1 − e−2xsur )]} + 3Nksur γ2sur Xsur + Xsur − 2γ1sur sur { } θ3 Xsur Xsur + 6Nksur γ )X − (γ )(1 + X ) , (1 + + 2γ γ ) (1 + sur sur 1sur 2sur 2sur 1sur 2 sur (2.10) here, Xsur = xsur cthxsur The free energy of the layers of atoms just below the surface and internal layers of atoms have fomulae similar to Eq (2.10) and can be noted as Ψsur1 , Ψin , respectively The free energy in the harmonic approximation for surface layers, layers of atoms just below the surface and internal layers of atoms of thin film are given as } { −2xsur )] ,} Ψsur ≈ 3Nsur { 61 usur + θ[xsur + ln(1 − e Ψsur1 ≈ 3Nsur1 {16 usur1 + θ[xsur1 + ln(1 − e−2xsur1 } )] ; in Ψin ≈ 3Nin u0 + θ[xin + ln(1 − e−2xin )] , where ≡ usur ∑ i 112 sur1 ≡ φsur i0 , u0 ∑ i in φsur1 i0 , u0 ≡ ∑ i φin i0 (2.11) (2.12) Coefficients of thermal expansion of thin metal films investigated Let us consider a system where N atoms consist of n layers, the number of atoms on each layer are the same and and NL , ψsur denotes the free energy of the two surface layers of atoms, ψsur1 denotes the free energy of the two layers of atoms just below the surface and ψin denotes the free energy of the (n-4) internal layers atoms The number of internal layer atoms, atoms in layers just below the surface, and surface layer atoms, respectively ( ) N Nin = (n − 4) NL = − NL = N − 4NL , (2.13) NL Nsur1 = 2NL = N − (n − 2)NL and Nsur = 2NL = N − (n − 2)NL (2.14) Free energy of the system and of an atom, respectively, is given by ψ = Nin ψin + Nsur1 ψsur1 + Nsur ψsur = (N − 4NL ) ψin + 2NL ψsur1 + 2NL ψsur (2.15) [ ] ψ 2 = 1− ψin + ψsur1 + ψsur N n n n (2.16) Using a ¯ as the average nearest-neighbor distance and ¯b as the average thickness of two-layers and a ¯c as the average lattice constant, respectively, we can easily calculate the relation √ a ¯ b = √ and a a (2.17) ¯c = 2¯b = 2¯ The thickness d of thin film can be given by √ n−5 a ¯ 2asur + 2asur1 + √ ain = (n−1) √ 2 (2.18) In Eq (2.18), the average nearest-neighbor of surface layers atoms, atoms of layers just below the surface and internal layer atoms in thin metal film at a given temperature T can be determined as d = 2bsur +2bsur1 +(n − 5) bin = (n − 1) b = √ asur = asur (0) + < usur >0 , asur1 = asur1 (0) + < usur1 >0 , ain = ain (0) + < uin i i i >0 , (2.19) where asur (0), asur1 (0) and ain (0) denote the nearest-neighbor at zero temperature of surface layers atoms, atoms of layers just below the surface and internal layer atoms, respectively, which can be determined experimently or derived considering the minimum condition of the potential energy of the system ∂U0sur ωsur ∂ksur + = 0; ∂asur 4ksur ∂asur ∂U0sur1 ωsur1 ∂ksur1 + = 0; ∂asur1 4ksur1 ∂asur1 ∂U0in ωin ∂kin + = 0, ∂ain 4kin ∂ain (2.20) 113 Duong Dai phuong, Vu Van Hung and Nguyen Thi Hoa The average nearest-neighbor of thin film is determined as a ¯= 2asur + 2asur1 + (n − 5)ain n−1 (2.21) Finally, we obtain the expression of the free energy through the number of layers as follows √ ψ d − 3¯ a 2¯ a 2¯ a = √ ψin + √ ψsur + √ ψsur1 (2.22) N d 2+a ¯ d 2+a ¯ d 2+a ¯ 2.1.3 The coefficients of thermal expansion The average thermal expansion coefficient of metallic thin film can be calculated as follows KB d¯ a dsur αsur + dsur1 αsur1 + (d − dsur − dsur1 ) αin = , a ¯ dθ d (2.23) kB ∂y0in (T ) kB ∂y0sur (T ) kB ∂y0sur1 (T ) ; αsur = ; αsur1 = ain ∂θ asur ∂θ asur1 ∂θ (2.24) α= with αin = with dsur and dsur1 being the thickness of the surface layers and the layers just below the surface, respectively 2.2 Numerical results and discussion In order to apply the theoretical calculation to the Au, Ag and Al of metallic thin film, we use the pair interaction potential in the following form with the potential parameters D, m, n and r0 in Table φ(r) = [ ( r )n ( r )m ] D 0 m −n (n − m) r r (2.25) Using the expressions obtained in Section 2, we calculate the values of the average lattice constant, a ¯c , the thickness, d, and the linear thermal expansion coefficient, α for Au, Ag and Al thin metal thin film at various temperatures The calculated results are presented in Figures 2-5 Table Lennars-Jones potential parameters for Au, Ag and Al of metallic thin film [7] Metal Au Ag Al 114 n 10.5 11.5 12.5 m 5.5 5.5 5.5 r0 , A0 2.8751 2.8760 2.8541 D/kB , K 4683.0 3325.6 2995.6 Coefficients of thermal expansion of thin metal films investigated In Figures and 3, we present the thickness dependence of the average lattice constant of Au thin metal film as a function of the thickness, and it has been compared to gold splattered on glass presented in [3] We obtained the results of an average lattice constant for an Ag and Al decrease, also with increasing thickness In Figure 4, we present the temperature dependence of the thermal expansion coefficients of thin metal film for Al as a function of the temperature The thermal expansion coefficients increase with absolute temperature T, and our calculated results have been compared with the results presented in [2] for Al thin film, and Al bulk [9] are shown to be in good agreement In Figure 5, we also show the thickness dependence of the thermal expansion coefficients of Au and Al thin film, and the results calculated show that the thermal expansion coefficients for Au and Al decrease with increasing thickness When the thickness increases, the thermal expansion coefficients approach the bulk value in agreement with the results presented in [9] 115 Duong Dai phuong, Vu Van Hung and Nguyen Thi Hoa Conclusion We have presented the SMM formalism using Lennard-Jones interaction potentials and investigated the thermal expansion coefficients of Au, Ag and Al thin metal film The SMM is simple and physically transparent, and thermal expansion coefficients of thin metal film with FCC structures can be expressed in closed forms within the fourth order moment approximation of the atomic displacements In general, we have obtained good agreement in the thermal expansion coefficients between our theoretical calculations, and other theoretical results and experimental values Acknowledgement: This research is funded by Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant number 103.01-2011.16 REFERENCES [1] Bonderover E and Wagner S, 2004 A woven inverter circuit for e-textile applications, JEEE Elektron Dev Lett., 25:295 [2] Cornella G et al., 1998 Determination of temperature dependent unstressed lattice spacings crystalline thin films on substrates MRS online proc Lib., Vol 505, pp 527-532 [3] Kolska Z et al., 2010 Lattice parameter and expected density of Au nano-structures sputtered on glass Materials Letters 64, pp 1160-1162 [4] Liang L H and Li B., 2006 Size-dependent thermal conductivity of nanoscale semiconducting systems Physical Review B73 (15), 153303 [5] Liang L H et al., 2002 Size-dependent elastic modulus of Cu and Au thin films Solid State Communications, 121 (8), pp 453-455 [6] Masuda-Jindo K et al., 2003 Thermodynamic quantities of metals investigated by an analytic statistical moment method Phys Rev B67, 094301 [7] Mazomedov M N., 1987 J Fiz Khimic, 61, 1003 [8] Roduner E, 2006 Size dependent chemistry: properties of nanocrystals Chem Soc Rev., 35:583 [9] Simmons R.O and Balluffi R.W, 1960 Phys Rev 117, 52 [10] Tang N., Hung V V, 1988 Investigation of the thermodynamic properties of anhamonic crystal by the moment method I General Results for Face-Centred Cubic Crystals Phys.stat sol(b), 149, pp 511-519 116 ... Figure The thin metal film 110 Coefficients of thermal expansion of thin metal films investigated First of all, we present the statistical moment method (SMM [6, 10]) formulation > as a solution of. .. and the thickness dependence of the coefficients of thermal expansion of thin metal film with face-centered cubic structure using the analytic statistical moment method (SMM) [6, 10] .The coefficients. .. Coefficients of thermal expansion of thin metal films investigated In Figures and 3, we present the thickness dependence of the average lattice constant of Au thin metal film as a function of