Chemical bonding and elastic properties of quaternary arsenide oxides YZnAsO and LaZnAsO investigated by first principles
SHI Yi-ming1, YE Shao-long2
1. Department of Basic Teaching, Hunan Institute of Technology, Hengyang 421000, China;
2. School of Metallurgical Science and Engineering, Central South University, Changsha 410083, China
Received 9 March 2011; accepted 8 June 2011
Abstract:
The structural parameters, chemical bonding and elastic properties of the tetragonal phase quaternary arsenide oxides YZnAsO and LaZnAsO were investigated by using density-functional theory (DFT) within generalized gradient approximation (GGA). The GGA calculated structural parameters are in agreement with the experimental results. Population analysis suggests that the chemical bonding in YZnAsO and LaZnAsO can be classified as a mixture of ionic and covalent characteristic. Single-crystal elastic constants were calculated and the polycrystalline elastic modules were estimated according to Voigt, Reuss and Hill’s approximations (VRH). The result shows that both YZnAsO and LaZnAsO are relatively soft materials exhibiting ductile behavior. The calculated polycrystalline elastic anisotropy result shows that LaZnAsO is more anisotropy in compressibility and YZnAsO is more anisotropy in shear.
Key words:
YZnAsO; LaZnAsO; chemical bonding; elastic properties; density-functional theory; generalized gradient approximation; Voigt, Reuss and Hill’s approximations;
1 Introduction
The quaternary oxypnictide layered compounds REZnPO (RE=Y, La-Nd, Sm, Gd, Dy, Ho) and REZnAsO (RE=Y, La-Nd, Sm, Gd-Er) with the ZrCuSiAs-type structure have received considerable attention [1-6] because of the discovery of superconductivity in LaOFeP, LaNiOP and LaFeAsO(F)[7-11]. They have similar layered structures with alternating [ReO] and [ZnPn] (Pn=As, P) layers, and exhibit a wide variety of interesting electrical, optical, and magnetic properties. YZnAsO and LaZnAsO, which are isostructural to the layered rare-earth oxypnictides REZnAsO (RE=Y, La-Nd, Sm, Gd, Er), are semiconductor compounds with band gap of 1.5-2.3 eV in experimental investigation [2, 5] or 0.65-1.3 eV in theoretical calculation [6]. KAYANUMA et al [2] have successfully grown LaZnAsO thin films on MgO substrates with an epitaxial relationship of (001) LaZnAsO || (001) MgO and [110] LaZnAsO || [110] MgO. The electrical conductivity of the undoped LaZnAsO epitaxial films was 0.2 S/cm at room temperature. The positive Seebeck coefficients of the epitaxial film indicated a p-type semiconductor. The optical bandgap of LaZnAsO was estimated to be about 1.5 eV. It is expected that LaZnAsO may be a promising matrix for diluted magnetic semiconductors if Zn2+ is partially replaced by a divalent transition metal, such as Mn2+. LINCKE et al [5] have synthesized well crystallized form REZnAsO (RE = Y, La) from the ingots of rare earth elements, arsenic, and ZnO in NaCl/KCl fluxes in sealed silica ampoules. The REZnAsO structures consist of alternate stacks of (RE3+O2-) and (Zn2+As3-) layers with covalent RE—O and Zn—As intralayer and weak ionic interlayer bonding. LaZnAsO and YZnAsO show an optical bandgap of Eg≈2.3 eV and Eg≈1.9 eV respectively. BANNIKOV et al [6] investigated the electronic properties and chemical bonding picture of LaZnAsO and YZnAsO. The band structures result showed that both of them are semiconductors with the band gaps of about 0.65 eV for LaZnAsO and 1.3 eV for YZnAsO. The bonding can be classified as a mixture of ionic and covalent contributions. To the best of our knowledge, there are no reports on the elastic properties of YZnAsO and LaZnAsO. However, elastic properties are very important for materials because they provide information on interatomic potentials and relate to various fundamental solid state phenomena such as interatomic bonding, equations of state, phonon spectra as well as specific heat capacity, thermal expansion, Debye temperature [12]. Most importantly, elastic properties are essential for many practical applications related to the mechanical properties of solid: load deflection, internal strain, thermoelastic stress, sound velocities and fracture toughness [12].
In this study, the chemical bonding and elastic properties of the tetragonal phase of YZnAsO and LaZnAsO are investigated using first principles calculations based on DFT.
2 Calculation details and models
The first principle calculations were based on DFT using a plan-wave expansion of the wave function implanted in CASTEP code [13]. The exchange correlation energy was calculated by the GGA with the Perdew-Burke-Ernzerhof (PBE) function [14]. The ionic cores were represented by norm-conserving pseudo- potentials [15-17] for Y, La, O, Zn and As atoms. The Y 4d15s2 electrons, La 5d16s2 electrons, O 2s22p4 electrons, Zn 3d104s2 electrons, and As 4s24p3 electrons were explicitly treated as valence electrons. The Monkhorst and Pack scheme of k-point sampling was used for integration over the first Brillouin zone [18]. The cutoff energy was chosen to be 900 eV, and the Brillouin-zone sampling k-point set mesh parameters were 12×12×6. This set of parameters assure the total energy convergence of 5.0×10-6 eV/atom, the maximum force of 0.01 eV/?, the maximum stress of 0.02 GPa and the maximum displacement of 5.0×10-4 ?.
3 Results and discussion
3.1 Geometry and structure properties
The crystal structure of the tetragonal phase of YZnAsO and LaZnAsO belongs to the space group P4/nmm, Z=2 (ZrCuSiAs type), where blocks [YO/LaO] are sandwiched with [ZnAs] blocks as depicted in Fig. 1. There are four inequivalent atomic positions: Y/La at 2c site (1/4, 1/4, zY/La), O at 2a site (3/4, 1/4, 0), Zn at 2b site (3/4, 1/4, 1/2), and As at 2c site (1/4, 1/4, zAs) [2, 5-6], in which zY/La and zAs are the internal coordinates of Y/La and As, respectively. The experimental lattice parameters are a=b=0.394 2 nm and c=0.884 0 nm for YZnAsO, and a=b=0.409 6 nm and c=0.907 0 nm for LaZnAsO; and the internal coordinates of Y/La and As are reported as zY/La=0.121 7 and zAs=0.682 5 for YZnAsO and zY/La=0.133 7 and zAs=0.671 1 for LaZnAsO [2, 5-6], respectively.
Fig. 1 Crystal structure sketch of YZnAsO
At the first stage the full structural optimization of this phase was performed both over the lattice parameters and the atomic positions including the internal coordinate zY/La and zAs. The calculated optimization lattice parameters a, c, V and atomic coordinates compared with available experimental data [2, 5] for YZnAsO and LaZnAsO are summarized in Table 1, which show that the calculated structural parameters are in agreement with the experimental data.
3.2 Chemical bonding
The Mulliken bond population is useful for evaluating the bonding character in a material. A high value of the bond population indicates a covalent bond, and a low value indicates an ionic bond. Positive and negative values indicate bonding and anti-bonding states, respectively [17, 19]. The Mulliken atomic population analyses of YZnAsO and LaZnAsO listed in Tables 2 and 3 show a charge transfer from [YO/LaO] block to [ZnAs] block, indicating that the internal of [YO/LaO] block and the [ZnAs] block is ionic character; and the bond population shows the intra-block of Y—O/ La—O and Zn—As are covalent character. The chemical bondings in YZnAsO and LaZnAsO have a mixture of ionic and covalent characteristic, which is in agreement with the results of Refs. [5-6].
Table 1 Calculated lattice parameters and atomic internal coordinate compared with experimental data [2, 5] for YZnAsO and LaZnAsO
Table 2 Mulliken atomic population analysis of YZnAsO and LaZnAsO
Table 3 Mulliken bond population analysis of YZnAsO and LaZnAsO
3.3 Elastic properties
Elastic constants are defined by means of a Taylor expansion of the total energy, namely, the derivative of the energy as a function of a lattice strain [13]. The tetragonal phase YZnAsO and LaZnAsO crystals have six independent elastic constants C11, C33, C44, C66, C12 and C13 [20-21] which are presented in Table 4. For the tetragonal crystal, its mechanical stability requires that its independent elastic constants should satisfy the Born’s stability criteria [20-21]:
Cii>0 (i=1, 3, 4, 6),
C11-C12>0, C11+C33-2C13>0,
2(C11+C12)+C33+4C13>0 (1)
Table 4 Calculated single crystal elastic constants Cij, bulk modulus B and compressibility coefficient β of YZnAsO and LaZnAsO
From Table 4, it can be seen that these criteria are all satisfied, which indicates that YZnAsO and LaZnAsO are mechanically stable. The single crystal bulk moduli B are 78.38 and 87.00 GPa for YZnAsO and LaZnAsO, respectively, which means that YZnAsO and LaZnAsO are soft materials. The polycrystal bulk modulus B and shear modulus G are estimated using the VRH approach in the following forms [20-23]:
,
,
,
,
,
,
,
(2)
The elastic modulus E and Poisson ratio v are estimated by [20-21]
,
(3)
All the calculated results are presented in Table 5. It can be seen that the value of BH/GH ratio for YZnAsO is 2.57 and for LaZnAsO is 1.87, which are larger than the critical value (1.75) separating ductile and brittle materials [20-21], indicating that both YZnAsO and LaZnAsO behave in a ductile manner, and YZnAsO is more ductile than LaZnAsO. It is known that the values of the Poisson ratio (ν) are minimal for covalent materials (ν=0.1), and increase for ionic systems [20-21]. In our case, the calculated Poisson ratios (ν) are 0.327 7 for YZnAsO and 0.272 9 for LaZnAsO, which means a sizable ionic contribution in intra-bonding.
The elastic anisotropy in compressibility (AB) and shear (AG) using the model in Ref. [24] for polycrystalline materials is
,
(4)
Table 5 Calculated polycrystalline elastic constants, elastic modulus E, Poisson ratio ν and BH/GH of YZnAsO and LaZnAsO
A value of zero represents elastic isotropy and a value of 100% is the largest possible anisotropy. The calculated polycrystalline elastic anisotropy in bulk modulus AB and shear modulus AG are listed in Table 6. It shows that the shear modulus anisotropy AG is larger than the bulk modulus anisotropy AB for both YZnAsO and LaZnAsO. And LaZnAsO is more anisotropy in compressibility (4.35%) and YZnAsO is more anisotropy in shear (22.04%). The larger anisotropy may produce cracking propagation during sample preparing and using. Utilization of very small YZnAsO and LaZnAsO particles would be a way of minimizing cracking.
Table 6 Calculated polycrystalline elastic anisotropy for bulk modulus AB and shear modulus AG of YZnAsO and LaZnAsO
4 Conclusions
1) The structural parameters, chemical bonding and elastic properties of YZnAsO and LaZnAsO were investigated by first principles. The GGA calculated structural parameters are in agreement with the available experimental data.
2) The chemical bonding was analyzed. The results show that YZnAsO and LaZnAsO are a mixture of ionic and covalent characteristic.
3) The single crystal elastic constants were calculated, and the polycrystal bulks and shear modulus, elastic modulus, Poisson ratios and BH/GH were derived. The results show that both YZnAsO and LaZnAsO are mechanically stable and relatively soft materials behaving in a ductile manner.
4) The calculated polycrystalline elastic anisotropy in bulk modulus and anisotropy in shear modulus show that LaZnAsO is more anisotropy in compressibility and YZnAsO is more anisotropy in shear.
References
[1] KAYANUMA K, HIRAMATSU H, HIRANO M, KAWAMURA R, YANAGI H, KAMIYA T, HOSONO H. Apparent bipolarity and seebeck sign inversion in layered semiconductor: LaZnOP [J]. Phys Rev B, 2007, 76: 195325.
[2] KAYANUMA K, KAWAMURA R, HIRAMATSU H, YANAGI H, HIRANO M, KAMIYA T, HOSONO H. Heteroepitaxial growth of layered semiconductors, LaZnOPn (Pn=P and As) [J]. Thin Solid Films, 2008, 516: 5800-5804.
[3] TAKANO Y, KOMATSUZAKI S, KOMASAKI H, WATANABE T, TAKAHASHI Y, TAKASE K. Electrical and magnetic properties of LnOZnPn (Ln = rare earths; Pn = P, As, Sb) [J]. J Alloys Compd, 2008, 451: 467-469.
[4] LINCKE H, GLAUM R, DITTRICH V, TEGEL M, JOHRENDT D, HERMES W, M?LLER M H, NILGES T, P?TTGEN R. Magnetic, optical, and electronic properties of the phosphide oxides REZnPO (RE = Y, La-Nd, Sm, Gd, Dy, Ho) [J]. Z Anorg Allg Chem, 2008, 634, 1339-1348.
[5] LINCKE H, GLAUM R, DITTRICH V, M?LLER M H, P?TTGEN R. Structure and optical properties of the arsenide oxides REZnAsO (RE=Y, La-Nd, Sm, Gd-Er) [J]. Z Anorg Allg Chem, 2009, 635: 936-941.
[6] BANNIKOV V V, SHEIN I R, IVANOVSKII A L. Electronic properties and chemical bonding in quaternary arsenide oxides LaZnAsO and YZnAsO [J]. Mater Chem Phys, 2009, 116: 129-133.
[7] KAMIHARA Y, HIRAMATSU H, HIRANO M, KAWAMURA R, YANAGI H, KAMIYA T, HOSONO H. Iron-based layered superconductor: LaOFeP [J]. J Am Chem Soc, 2006, 128: 10012-10013.
[8] WATANABE T, YANAGI H, KAMIYA T, KAMIHARA Y, HIRAMATSU H, HIRANO M, HOSONO H. Nickel-based oxyphosphide superconductor with a layered crystal structure, LaNiOP [J]. Inorganic Chemistry, 2007, 46(19): 7719-7722.
[9] KAMIHARA Y, WATANABE T, HIRANO M, HOSONO H. Iron-based layered superconductor La[O1-xFx]FeAs (x= 0.05-0.12) with Tc=26 K [J]. J Am Chem Soc, 2008, 130: 3296-3297.
[10] WEN H H. Developments and perspectives of iron-based high-temperature superconductors [J]. Adv Mater, 2008, 20(19): 3764-3775.
[11] REN Z A, ZHAO Z X. Research and prospects of iron-based superconductors [J]. Adv Mater, 2009, 21(45): 4584-4593.
[12] PENG F, CHEN D, YANG X D. First-principles calculations on elasticity of OsN2 under pressure [J]. Solid State Commun, 2009, 149: 2135-2138.
[13] SEGALL M D, LINDAN PHILIP J D, PROBERT M J, PICKARD C J, HASNIP P J, CLARK S J, PAYNE M C. First-principles simulation: Ideas, illustrations and the CASTEP code [J]. J Phys: Condens Matter, 2002, 14(11): 2717-2744.
[14] PERDEW J P, BURKE K, ERNZERHOF M. Generalized gradient approximation made simple [J]. Phys Rev Lett, 1996, 77: 3865-3868.
[15] HAMANN D R, SCHLUTER M, CHIANG C. Norm-conserving pseudopotentials [J]. Phys Rev Lett, 1979, 43: 1494-1497.
[16] ZHANG X Y, CHEN Z W, ZHANG S L, LIU R P, ZONG H T, JING Q, LI G, MA M Z, WANG W K. Electronic and optical properties of rock-salt aluminum nitride obtained from first principles [J]. J Phys: Condens Matter, 2007, 19: 425231.
[17] XIA Q L, YI J H, LI Y F, PENG Y D, WANG H Z, ZHOU C S. First-principles investigations of the band structure and optical properties of γ-boron [J]. Solid State Commun, 2010, 150: 605-608.
[18] MONKHORST H J, PACK J D. Special points for brillouin-zone integrations [J]. Phys Rev B, 1976, 13: 5188-5192.
[19] HERMET P, GOUMRI-SAID S, KANOUN M B, HENRARD L. First-principles investigation of the physical properties of magnesium nitridoboride [J]. J Phys Chem C, 2009, 113: 4997-5003.
[20] SHEIN I R, IVANOVSKII A L. Elastic properties of quaternary oxyarsenide LaOFeAs and LaOFeP as basic phases for new 26-52 K superconducting materials from first principles [J]. Scripta Materialia, 2008, 59: 1099-1102.
[21] SHEIN I R, IVANOVSKII A L. Elastic properties and chemical bonding in ternary arsenide SrFe2As2 and quaternary oxyarsenide LaFeAsO—Basic phases for new 38-55 K superconductors from first principles [J]. Physica C, 2009, 469: 15-19.
[22] HILL R. The elastic behaviour of a crystalline aggregate [J]. Proc Phys Soc London A, 1952, 65: 349-355.
[23] ANDERSON O L. A simplified method for calculating the Debye temperature from elastic constants [J]. J Phys Chem Solids, 1963, 24(7): 909-917.
[24] CHUNG D H, BUESSEM W R. Anisotropy in single crystal refractory compounds, edited by vahldiek F W and mersol S A [M]. Vol.2. New York: Plenum Press, 1968: 217.
四元砷氧化物YZnAsO和LaZnAsO的
化学键和弹性性质的第一性原理研究
施毅敏1, 叶绍龙2
1. 湖南工学院 基础教学部，衡阳 421000；
2. 中南大学 冶金科学与工程学院，长沙 410083
摘 要：利用基于密度泛函理论(DFT)的广义梯度近似(GGA)研究四方相四元砷氧化物YZnAsO和LaZnAsO的结构参数、化学键和弹性性质。结果表明，GGA计算的结构参数与实验值吻合较好。布居分析显示YZnAsO和LaZnAsO中的化学键具有混合的离子-共价特征。计算得到了单晶的弹性常数，由Voigt, Reuss和Hill (VRH)近似导出了多晶弹性模量，结果表明YZnAsO和LaZnAsO为具有一定延展性、相对较软的材料。多晶的弹性各向异性结果显示，LaZnAsO的体模量各向异性相对较大，而YZnAsO的剪切模量各向异性相对较大。
关键词：YZnAsO；LaZnAsO；化学键；弹性性质；密度泛函理论；广义梯度近似；VRH近似
(Edited by YUAN Sai-qian)
Foundation item: Project (50474051) supported by the National Natural Science Foundation of China
Corresponding author: SHI Yi-min; Tel: +86-15111157635; E-mail: hn_syb@126.com
DOI: 10.1016/S1003-6326(11)60869-8
