Study of some thermodynamic properties of ZrxCe1xO2/XeO2 systems by statistical moment method

Chia sẻ: Y Y | Ngày: | Loại File: PDF | Số trang:11

Thêm vào BST

Báo xấu

14
lượt xem 1
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Some thermodynamic quantities of ZrxCe1-xO2/CeO2 systems are investigated by using the moment method in statistical dynamics taking into account the anharmonicity effects of lattice vibrations. The analytic expression of the Gibbs free energy, Helmholtz free energy, and specific heats at the constant volume of ZrxCe1-xO2/CeO2 systems are obtained.

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Study of some thermodynamic properties of ZrxCe1xO2/XeO2 systems by statistical moment method

HNUE JOURNAL OF SCIENCE DOI: 10.18173/2354-1059.2018-0068 Natural Sciences 2018, Volume 63, Issue 11, pp. 41-51 This paper is available online at http://stdb.hnue.edu.vn STUDY OF SOME THERMODYNAMIC PROPERTIES OF ZrxCe1-xO2/CeO2 SYSTEMS BY STATISTICAL MOMENT METHOD Vu Van Hung1, Le Thi Thanh Huong2 and Dang Thanh Hai3 1 University of Education, Vietnam National University, Hanoi 2 Hai Phong University, 3Vietnam Education Publishing House Abstract. Some thermodynamic quantities of ZrxCe1-xO2/CeO2 systems are investigated by using the moment method in statistical dynamics taking into account the anharmonicity effects of lattice vibrations. The analytic expression of the Gibbs free energy, Helmholtz free energy, and specific heats at the constant volume of ZrxCe1-xO2/CeO2 systems are obtained. The lattice parameter and specific heats at the constant volume of the ZrxCe1-xO2/CeO2 systems are calculated as functions of the temperature, pressure, concentration of Zr and the thickness ratio d2/d1 of the ZrxCe1-xO2/CeO2 systems by using the Buckingham potential. The effects of temperature (for the temperature range T = 0K ÷ 2900K), pressure (P = 0; 5; 10; 20; 30; and 40GPa), concentration of Zr (x = 2%; 4%; 6%; 8%), and the thickness ratio d2/d1 (d2/d1 = 1÷ 20) on the lattice parameter and specific heats at the constant volume of the ZrxCe1- xO2/CeO2 systems will be discussed in details. The dopant dependence of the lattice constants of Ce1−xZrxO2 at T = 300K and zero pressure, and the pressure dependence of the lattice constants of CeO2 at T = 300K are compared with the available experimental and other theoretical results. Keywords: Statistical moment method, high pressure, oxide superlattice. 1. Introduction Study of the thermodynamic properties of semiconductor and oxide superlattices has become quite interesting in recent years since semiconductor and oxide super lattices play an important role in technological applications, especially in the electronic industry. The structural and thermo- mechanical properties of semiconductor superlattices were studied by using the lattice dynamics approach [1] and ab initio molecular dynamic simulations [2]. The structural characteristics of undoped CeO2/ZrO2 superlattices were investigated by Wang et al. [3]. This study revealed several interesting properties of the layered film structure beyond their improved ionic conduction. When the individual layer thickness in superlattices consisting of alternating layers of ceria and zirconia doped with gadolinia (Gd2O3) was reduced below approximately 50 nanometers, the ionic conduction of this superlattice began to increase greatly with respect to either bulk material [4]. The increase in conduction observed in the ceria/zirconia superlattice could prove to be a breakthrough in this regard. Received September 7, 2018. Revised October 15, 2018. Accepted October 22, 2018. Contact Dang Thanh Hai, e-mail address: dthai@nxbgd.vn 41
Vu Van Hung, Le Thi Thanh Huong and Dang Thanh Hai It is remarkable that despite the widespread use of CeO2/ZrO2 and doped ceria/zirconia superlattices like as ZrxCe1-xO2/CeO2 in oxygen sensors, and other related applications for several decades, the studies of thermodynamic quantities of doped ceria/zirconia superlattice are not fully understood in these materials and continue to be an important focus of experimental and theoretical studies. The purpose of present article is to study effects of the temperature, pressure and dopant concentration on the thermodynamic properties of the ZrxCe1-xO2/CeO2 superlattice using the statistical moment method (SMM) [5-6]. The lattice parameter and specific heats at the constant volume of ZrxCe1-xO2/CeO2 superlattice are calculated as functions of the temperature, pressure and doped concentration by using the Buckingham potential. 2. Content 2.1. Theory The ZrxCe1-xO2/CeO2 oxide superlattice consists of ZrxCe1-xO2 and CeO2 layers, where both the ceria (CeO2) and doped ceria (ZrxCe1-xO2) layers are assumed to be in the cubic fluorite structure, consistent with experimental observations [3]. In the fluorite structure, cations are located in the face-centered cubic positions, with eight oxygen ions in the positions. d1 d2 Zrx Ce1 x O2 CeO2 Figure 1. Structure of ZrxCe1-xO2/CeO2 superlattices In the ZrxCe1-xO2/CeO2 present oxide superlattice with very thin layers and small lattice mismatch, the difference in lattice parameter between the layers can be completely compensated for by strain. We assume that the thickness of ZrxCe1-xO2 layer is d1 and this oxide layer consists of N1 atoms with NZr atoms of Zr, NCe atoms of Ce and NO atoms of O (as shown in Fig.1) then N1  NZr  NCe  NO , (1) N Zr N x  Zr , N Zr  NCe N1 / 3 xN1 N N 2N (2) N Zr  , NCe  1  N Zr  1 (1  x), NO  1 . 3 3 3 3 Similarly, CeO2 ceria layer is supposed to have the thickness d2 and consist of N2 atoms with NCe atoms of Ce and N O* atoms of O then * 42
Study of some thermodynamic properties of ZrxCe1-xO2/CeO2 systems by statistical moment method N2  NCe *  NO* . (3) So N  ( N1  N2 )n, (4) where N is the number of atoms and n is the period of the ZrxCe1-xO2/CeO2 oxide superlattice. In order to investigate the thermodynamic properties of the ZrxCe1-xO2/CeO2 oxide superlattice, we firstly consider the change of Gibbs free energy of ZrxCe1-xO2 system when NZr atoms of Ce are replaced by Zr atoms in CeO2 crystal. The substitution of an atom Ce by an atom Zr causes the change of the free Gibbs energy g Zrf as g Zrf  u0Ce  Zr , (5) where  Zr is the free energy of an atom Zr in the ZrxCe1-xO2 system, and u0Ce is the cohensive energy associated with atom Ce of the cubic-fluorite CeO2. Because of the ZrxCe1-xO2 system is supposed to be built by substituting NZr atoms Zr into the positions of Ce atoms of cubic-fluorite CeO2 system then the Gibbs free energy of system has an approximate form as d1 GCeO2  N Zr g Zrf  TSc , xN1 f G  GCeO d1  g Zr  TSc , (6) 2 3 d1 where Sc is the configuration entropy of ZrxCe1-xO2 system, and GCeO2 is the Gibbs free energy of CeO2 system with N1 of Ce and O atoms ( N1  NCe  NO ). From Eqs. (5), and (6), we obtain the Gibbs free energy of ZrxCe1-xO2 system G  GCeO d1 2  N Zr (u0Ce  Zr )  TSc , xN1 xN G  GCeO d1   Zr  1 u0Ce  TSc , (7) 2 3 3 d1 where GCeO2 is the Gibbs free energy of the CeO2 system with N1 of Ce and O atoms d1 GCeO2   CeO d1 2  PV1 , (8) and P denotes the hydrostatic pressure, V1 is the volume of the ZrxCe1-xO2 system. Because of CeO2 ceria layer of oxide superlattice is supposed to have the thickness d2 and consist of N2 atoms with NCe * atoms of Ce and N O* atoms of O, then the Gibbs free energy is given by d2 GCeO2   CeO d2 2  PV2 . (9) From Eqs. (7), (8) and (9), it is easy obtain the Gibbs free energy of the ZrxCe1-xO2/CeO2 oxide superlattice Gsup  n(G  GCeO d2 2 ), 43
Vu Van Hung, Le Thi Thanh Huong and Dang Thanh Hai  d1 xN xN  Gsup  n  CeO  1  Zr  PV1  1 u0Ce  TSc   CeO d2  PV2 . (10)  2 3 3 2  The Helmholtz free energy  CeO d1 2 ,  CeO d2 2 , of ZrxCe1-xO2 and CeO2 layers is then written by taking into account the configurational entropies Scd1 , and Scd2 via the Boltzmann relation as [5,6]  CeO d1 2  CCe N1 Ce  CO N1 O  TScd1 , N1 2N  CeO d1   Ce  1  O  TScd1 , (11) 2 3 3  CeO d2 2  CCe N2 Ce  CO N2 O  TScd2 , N2 2N  CeO d2   Ce  2  O  TScd2 , (12) 2 3 3 where  Ce ,  O are the total Helmholtz free energies of an Ce and O atoms, CCe , CO denote concentrations of Ce, O atoms ( CCe = , ), respectively. Using SMM the quasi-harmonic contributions to the free energies of Ce, and O atoms in the CeO2 crystal are treated as [5] 1  Ce  u0Ce  0Ce   iCe0 (| ri |)  3  xCe  ln(1  e2 xCe ) , 2 i (13) 1 2   O  u0O  0O   O (| ri |)  3  xO  ln(1  e2 xO )  . i i0 (14) In Eqs. (13), (14) xR , xO … are given by Ce kCe / m 1   2iCe  xCe  2  2 , kCe    i 0 2   , (15) 2  ui  eq 1   2iO0   kO / m kO    2   mO2 , xO  O  i  2 2 . (16) 2  ui  eq 0 (or i 0 ) is the interaction potential between the 0-th and the where       x, y, or z, and iCe O i-th Ce (or O) atoms, and ui , ui are  ,  -Cartesian components of the displacement of i-th ion, and Ce (or O ) is the vibration frequency of Ce (or O) atoms, and   kBT (kB - the Boltzmann constant), and m is the average atomic mass of the system, m  CCe mCe  CZr mZr  CO mO . Using Eqs. (1), (2), (3), (4), (10), (11), (12), (13) and (14) we obtain the Gibbs free energy of ZrxCe1-xO2/CeO2 oxide superlattice Gsup  N 3  Ce  2N 3    O  N Zr u0Ce   Zr  PV  nTSc* , (17) Gsup  3  N Ce  u0  3  xCe  ln 1  e2 xCe      2N O 3  u0  3  xO  ln 1  e2 xO       44
Study of some thermodynamic properties of ZrxCe1-xO2/CeO2 systems by statistical moment method  N x 3  d2  u0Ce  N x 3  d2     u0Zr  3  xZr  ln 1  e2 xZr   PV  nTSc* .   1   1    d1   d1      G sup  N 3  1 x  d   u Ce  2 N u O  N  0 3 0 x 3  d2     u0Zr  N  xCe  ln 1  e 2 xCe    1  2   1     d1    d1  N x    d2      xZr  ln 1  e2 xZr   PV  nTSc* , (18) 1    d1  where V  n V1  V2  , Sc*  Scd1  scd2  Sc . (19) In order to calculate the average nearest-neighbor distance (NND) between two intermediate atoms a(P,T) in ZrxCe1-xO2/CeO2 oxide superlattice at temperature T and various pressure P, we can use the minimum condition of the Gibbs free energy of ZrxCe1-xO2/CeO2. G sup  0. (20) a From Eqs. (18) and (20), we obtain the following equation G SM   1 x   u0Ce  2  u0O  x  u0Zr    N            a  3 3(1  X )   a  3  a  3(1  X )  a     x coth xCe  kCe  x coth xO  kO   x xZr coth xZr  kZr   N  Ce    2 O        2kCe  a  2kO  a  (1  X ) 2kZr  a  V p , (21) a with X  d2 / d1. By numerically solving Eq.(21), one can determine the average nearest-neighbor distance (NND) between two intermediate atoms a(P,T) at temperature T and pressure P of ZrxCe1-xO2/CeO2. In the case of zero pressure (P = 0), using Eq. (17) we can find the Helmholtz free energy of ZrxCe1-xO2/CeO2 N 2N  sup   Ce   O  N Zr (u0Ce  Zr )  nTSc* , 3 3 N 2N N x   Ce  O  (u0Ce   Zr )  nTSc* , 3 3 3  d2  (22) 1    d1  45
Vu Van Hung, Le Thi Thanh Huong and Dang Thanh Hai N x  N CeO2  (u0Ce   Zr )  nTSc* , 3  d2  (23) 1    d1  where N 2N  CeO2   Ce   O  nT (Scd1  Scd2 ). (24) 3 3    Applying the Gibbs-Helmholtz relation, E       , and using Eq.(23), we obtain the    expression for the energy of ZrxCe1-xO2/CeO2 oxide superlattice 1 x N x   u0Ce  E sup  ECeO2  EZr     u0Ce   , 3  d2  3  d 2       (25) 1   1    d1   d1    CeO2    (26)   NT 2  , T ECeO2 T     Zr  EZr   NT 2  T  (27) . T Furthermore, the specific heat at constant volume CV of ZrxCe1-xO2/CeO2 oxide superlattice is determined as E sup  E sup  CV   kB  ,       u Ce     0  u0Ce   CV  CVCeO2  1 x CVZr k N  B x      . (28) 2  d2  3  d2   1   1    d1   d1  2.2. Results and discussion To calculate the thermodynamic quantities of ZrxCe1-xO2/CeO2, we will use the Buckingham potential [7] qi q j Cij ij   Aij exp(r / Bij )  , (29) r r6 where qi and qj are the charges of the i-th and the j-th ions, r is the distance between them and the parameters Aij, Bij and Cij are empirically determined by Ref.[8-11]. The potential parameters of ZrxCe1-xO2/CeO2 systems are given in Table 1. 46
Study of some thermodynamic properties of ZrxCe1-xO2/CeO2 systems by statistical moment method Table 1. The potential parameters Aij, Bij and Cij of ZrxCe1-xO2/CeO2 systems Interaction A(eV) B(Å) C(eV.Å6) O2- - O2- 9547.92 0.2192 32.00 Potential 1 Ce4+ - O2- 1809.68 0.3547 20.40 Zr4+- O2- 1502.11 0.3477 5.10 O2- - O2- 9547.92 0.2192 32.00 Potential 2 Ce4+ - O2- 2531.50 0.3350 20.40 Zr4+- O2- 1502.11 0.345 5.10 O2- - O2- 22764.3 0.1490 27.89 Butler Ce4+ - O2- 1986.83 0.3511 20.40 Zr4+- O2- 985.87 0.3760 0.00 In the case of the zero thickness d2 (d2  0), using Eq.(21) one can determine the average nearest-neighbor distance (NND) between two intermediate atoms a(P,T) of ZrxCe1-xO2 layer. The lattice constants of the ZrxCe1-xO2 system with the different dopant concentrations at temperature T = 300K and zero pressure are presented in Fig.2. One can see that the lattice constant of ZrxCe1-xO2 decreases with the increasing dopant concentration. The variation of the lattice parameter of ZrxCe1- xO2 as a function of Zr concentration in the present work by SMM (empty square), and our SMM results are in good agreement with the results of empirical equations [12], MD simulation [13], and experiments [14 -17]. Figure 2. The dopant dependence of the lattice constants of the cubic ﬂuorite Ce1−xZrxO2 system at T = 300K and zero pressure. The dark yellow circles are the results from Emprirical equation [12], the pink triangles with line and symbol are the results from the molecular dynamics (MD) simulations [13], red triangle, purple dimonds, blue stars and olive pentagons are the experimental results [14-17]. In Figures 3 and 4, we present the lattice parameter of the ZrxCe1-xO2/CeO2 systems as functions of thickness ratio d2/d1 at the different composition Zr (x = 2%; 4%; 6%; 8%) at room temperature, pressure P = 5GPa, and T = 900K, and P = 15GPa, respectively. As it can be seen from these two figures, the lattice parameters of the ZrxCe1-xO2/CeO2 systems are the increasing functions of thickness ratio d2/d1. Furthermore, when the thickness ratio d2/d1 of the ZrxCe1-xO2/CeO2 systems increases to 15, the superlattice constants have the same value as ceria's lattice constant. 47
Vu Van Hung, Le Thi Thanh Huong and Dang Thanh Hai Figure 3. Thickness dependence of lattice Figure 4. Thickness dependence of lattice constants of ZrxCe1-xO2/CeO2 with Zr different constants of ZrxCe1-xO2/CeO2 with Zr different concentrations (x = 2%; 4%; 6%, and 8%) at concentrations (x = 2%; 4%; 6%, and 8%) at T room temperature and pressure P = 5GPa = 900K, and P = 15GPa using potential 2. using potential 2. Figure 5. Pressure dependence of lattice Figure 6. Pressure dependence of lattice constants of ZrxCe1-xO2/CeO2 system with constants of ZrxCe1-xO2/CeO2 system with thickness d2 = d1, and Zr concentration x = thickness d2 = 20d1, and Zr concentration x = 2% at various temperature (T = 0K; 300K; 2% at room temperature using potentials 1,2 600K; 900K; 1200K; 1500K; 1800K; 2300K; and Butler potential and experimental results 2600K and 2900K) using potential 2 of CeO2 [18] The pressure dependence of the lattice constants of the ZrxCe1-xO2/CeO2 system with Zr concentration x = 2% and thickness d2 = d1, and d2 = 20d1 at various temperature and room temperature are presented in Figs.5 and 6, respectively. One can see that the lattice constant of ZrxCe1-xO2/CeO2 system decreases with the increasing pressure. In the case of ZrxCe1-xO2/CeO2 48
Study of some thermodynamic properties of ZrxCe1-xO2/CeO2 systems by statistical moment method system with a thickness ratio d2/d1 of 20 and a small dopant concentration (x = 2%), the thermodynamic properties of the ZrxCe1-xO2/CeO2 system are similar to those of the ceria (CeO2 layers with thickness d2). The variation of the lattice parameter of ZrxCe1-xO2/CeO2 as a function of pressure in the present work by SMM calculations using potentials 1, 2 and Butler potential (as presented in Figs. 6 and 7), and our SMM results are in good agreement with the experiments [18] and the ab initio calculations [19]. One can see in Fig. 6 that the lattice parameters calculated by using potentials 1 and 2 are very similar. The small difference between the two calculations simply comes from the difference in cerium-oxygen interaction potentials, since the ionic Coulomb contribution and the oxygen-oxygen potential are the same for potentials 1 and 2. In Fig.8 we show the temperature dependence of the lattice constants of ZrxCe1-xO2/CeO2 system with thickness d2 = d1, and Zr concentration x = 2% at various pressures calculated by using potential 2 for the temperature range T = 0K - 2900K. The lattice constant of ZrxCe1-xO2/CeO2 system increased smoothly with an increase of temperature due to the thermal expansion. Figure 7. Pressure dependence of lattice Figure 8. Temperature dependence of constants of ZrxCe1-xO2/CeO2 system with lattice constants of ZrxCe1-xO2/CeO2 thickness d2 = 20d1, and Zr concentration x = system with thickness d2 = d1, and Zr 2% at room temperature using potential 2 and concentration x = 2% at various pressures experimental results [18] and ab initio (P = 0; 5; 10; 20; 30, and 40GPa) using calculations of CeO2 [19] potential 2 The thickness ratio dependence of the specific heats at constant volume CV of the ZrxCe1-xO2/CeO2 system with Zr concentration x = 2% at zero pressure and room temperature are presented in Fig.9. As it can be seen from this figure, the specific heats at constant volume CV of the ZrxCe1-xO2/CeO2 systems are the increasing functions of thickness ratio d2/d1. But the value of the specific heats at constant volume CV of the ZrxCe1-xO2/CeO2 system with different thickness ratios of d2/d1 at the same temperature does not differ much. In Fig.10, we show the temperature dependence of the specific heats at constant volume CV of ZrxCe1-xO2/CeO2 system with thickness d2 = d1, and composition of Zr (x = 0.2) at zero pressure. The calculated specific heats CV increase briskly with temperature in the temperature range below 750K. As also shown in Fig.10, the specific heat CV depends weakly on the temperature for the temperature region from 750K up to the melting temperature. These properties of ZrxCe1-xO2/CeO2 system with thickness ratio d2/d1 = 1, and a small dopant concentration, are similar to those of ceria (CeO2 layers with thickness d2). 49
Vu Van Hung, Le Thi Thanh Huong and Dang Thanh Hai Figure 9. Thickness ratio dependence of the Figure 10. Temperature dependence of the specific heats at constant volume CV of ZrxCe1- specific heats at constant volume CV of xO2/CeO2 system with composition of Zr (x = ZrxCe1-xO2/CeO2 system with thickness d2 = 0.2) at zero pressure and room temperature d1, and composition of Zr (x = 0.2) at zero calculated using potential 2 (P2) pressure calculated using potentials 1, 2 (P1, P2) and Butler potential (B) 3. Conclusions In conclusion, the SMM calculations have been performed to investigate the thermodynamic properties of ZrxCe1-xO2/CeO2 systems. Our development is establishing and we have derived the analytical expressions of the Helmholtz free energy, specific heats at the constant volume of the cubic fluorite ZrxCe1-xO2/CeO2 systems. The lattice constants of ZrxCe1-xO2/CeO2 are calculated as a function of the dopant, temperature and pressure. The present calculated results are compared with the available experimental and theoretical results. REFERENCES [1] P.S. Branicio, J.P. Rino, R.K. Kalia, A. Nakano and P. Vashishta, 2003, J. Appl. Phys. 94, 3840. [2] S. Baron, S. de Giroconli, and P. Giannozzi, 1990, Phys. Rev. Lett. 65 84. [3] C.M. Wang, S. Azad, V. Shutthanandan, D.E. McCready, C.H.F. Peden, L. Saraf, S. Thevuthasan, 2005, Acta Materialia. 53, 1921. [4] S. Azad, O.A. Marina, C.M. Wang, L. Saraf, V. Shutthanandan, D.E. McCready, A. El-Azab, J.E. Jaffe, M.H. Engelhard, C.H.F. Peden, and S. Thevuthasan, 2005, Applied Physics Letters. 86, 131906. [5] V.V. Hung, J. Lee, and K. Masuda-Jindo, 2006, J. Phys. And Chem. Solids. 67, 682. [6] V.V. Hung, L.T.M. Thanh, K. Masuda-Jindo, 2010, Comp. Mater. Sci. 49, S355. [7] R. Devanathan, W.J. Weber, S.C. Singhal, J.D. Gale, 2006, Solid State Ionics. 177, 1251. 50
Study of some thermodynamic properties of ZrxCe1-xO2/CeO2 systems by statistical moment method [8] Beth Savoy, 2007. An investigation of ceria’s oxygen storage capacity by computer simulation. Department of Chemical and Biomolecular Engineering, Cornell University. MAE715 project. [9] M.O. Zacate, L. Minervini, D.J. Bradfield, R.W. Grimes, and K. E. Sickafus, 2000. Solid State Ionics. 128, 243. [10] P.K. Schelling and S.R. Phillpot, 2001. Journal of the American Ceramic Society. 84, 2997. [11] Shyam Vyas, Robin W. Grimes, David H. Gay and Andrew L. Rohl, 1998. J. Chem. Soc., Faraday Trans. 94, 427. [12] D.-J. Kim, 1989. J. Am. Ceram. Soc. 72, 1415. [13] G. Balducci, M.S. Islam, J.K Par, P. Fornasiero, M. Graziani, 2000, Chem. Mater. 12, 677. [14] Rodriguez, J.A. Hanson, J.C. Kim, J.Y. Liu, 2003. J. Phys. Chem. B107, 3535. [15] Cabanas et al, 2001, J. Mater. Chem. 11, 561. [16] T.A. Lee, C.R. Stanek, K.J. McClellan, J.N. Mitchell, Navrotsky, 2008. J. Mater. Res. 23,1105. [17] G. Dutta, U.V. Waghmare, T. Baidy, M.S. Hegde, K.R. Priolkar, P.R. Sarode, 2006. Catalysis Letters. 108, 165. [18] Z.W. Wang, Y. Zhao, D. Shiferl, C.S. Zha, and R. Downs. 2004. J Appl. Phys. Lett. 85, 124. [19] Li-Li Sun, Yan Cheng and Guang-Fu Ji, 2010. J. At. Mol. Sci. 1, 143. 51