Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Recent Progress in Lectin-Based Biosensors
Next Article in Special Issue
Development of Porous Piezoceramics for Medical and Sensor Applications
Previous Article in Journal
The Corrosion Characteristics and Tensile Behavior of Reinforcement under Coupled Carbonation and Static Loading
Previous Article in Special Issue
Bismuth Sodium Titanate Based Materials for Piezoelectric Actuators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Piezoelectric, Mechanical and Acoustic Properties of KNaNbOF5 from First-Principles Calculations

1
Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China
2
School of Physics and Engineering, Henan University of Science and Technology, Luoyang 471003, China
3
Department of Electrical Engineering, University of California, Riverside, CA 92521-0204, USA
*
Authors to whom correspondence should be addressed.
Materials 2015, 8(12), 8578-8589; https://doi.org/10.3390/ma8125477
Submission received: 24 September 2015 / Revised: 13 October 2015 / Accepted: 12 November 2015 / Published: 9 December 2015
(This article belongs to the Special Issue Piezoelectric Materials)

Abstract

:
Recently, a noncentrosymmetric crystal, KNaNbOF5, has attracted attention due to its potential to present piezoelectric properties. Although α- and β-KNaNbOF5 are similar in their stoichiometries, their structural frameworks, and their synthetic routes, the two phases exhibit very different properties. This paper presents, from first-principles calculations, comparative studies of the structural, electronic, piezoelectric, and elastic properties of the α and the β phase of the material. Based on the Christoffel equation, the slowness surface of the acoustic waves is obtained to describe its acoustic prosperities. These results may benefit further applications of KNaNbOF5.

Graphical Abstract

1. Introduction

Global innovations in electronic devices are driving the demand for new piezoelectric materials with improved properties. These materials are of great importance for a variety of technological applications including the production and detection of sound, the generation of high voltages, positioning in scanning probe microscopes, and the implementation of second harmonic generation (SHG) devices [1]. Moreover, piezoelectric materials are also known as the “smart” components in radiation environments, especially in the nuclear industry. During the past several years, piezoelectric technology has been revolutionizing sensing in reactor and waste vessel environments for the control and safety protection of nuclear reactors [2,3]. Piezoelectric sensing is important to the measurement of many control variables, such as temperature, pressure, flow, and neutron flux. The importance of piezoelectricity in these areas demands the search for novel piezoelectric materials of different varieties with better performances.
A noncentrosymmetric crystal, KNaNbOF5, has been proposed as a candidate piezoelectric material. It has exhibited many interesting and intriguing properties. Firstly, oxy-fluoride perovskite is not a common piezoelectric structure [4] compared to the well-known perovskite materials with the general formula ABO3 [5,6,7,8]. The piezoelectric mechanism in KNaNbOF5 is very special. The oxygen and fluorine atoms can co-exist in one framework, [NbOF5]2−. In this framework, each Nb is surrounded by five fluorine anions and one oxygen anion. These atoms form a distorted octahedron, where the distortion is caused by the different bond lengths to the centric Nb5+ ion. Thus, a net dipole is generated by the Nb5+ ion moving from the center of the O-F cage to the oxygen side. Secondly, KNaNbOF5 is also a triboluminescent material. It was reported that strong triboluminescence can be visible to the naked eye under normal lighting conditions [9]. Last but not least, KNaNbOF5 has attracted great interest due to its polymorphism [10,11]. For the same identities of the cations and anions, KNaNbOF5 has two different phases that exhibit very different behaviors. One phase, the α-KNaNbOF5, is a polar noncentrosymmetric polymorph [12]. The crystalline KNaNbOF5 was first synthesized by Antokhina et al. [13]. However, the crystal structure was not verified at the time. The details of the atomistic structure in the α-phase were reported by Poeppelmeier et al. several years later [12]. The other phase, β-KNaNbOF5, is centrosymmetric [10,14]. This phase was introduced for the first time by Vasiliev et al. [10]. Experimentally, both polymorphic phases can be made from similar combinations of starting materials. The key difference comes from the different K:Na ratio in the synthetic route. The α-phase occurs when the ratio of K:Na is greater than 1:1, while a smaller ratio results in the β-phase. Although the α- and β-phases of KNaNbOF5 have the same stoichiometry, the same polar structure of the framework [NbOF5]2−, and very similar synthetic routes, they have very different properties.
In order to gain a comprehensive understanding of the physical properties of KNaNbOF5, detailed investigations are desirable. In this paper, employing the density functional theory (DFT), detailed computational studies of the electronic, mechanical, piezoelectric, and acoustic properties of KNaNbOF5 are carried out for both the α and β polymorphic phases. The property differences between the two phases are compared and discussed in detail.

2. Theoretical Method

The DFT calculations are performed with the Vienna Ab initio Simulation Package (VASP) [15]. The electron-ion interactions are represented by the projector augmented wave (PAW) method [16]. The electronic exchange correlation energy is treated using the generalized gradient approximation suggested by Perdew, Burke, and Ernzerhof (GGA-PBE) [17]. The wave functions are expanded on a plane wave basis with kinetic energy cutoff set to 600 eV. The Brillouin zone is sampled on a mesh of 6 × 10 × 8 k-points within the Monkhorst-Pack scheme [18] for α-KNaNbOF5, and a 10 × 10 × 8 k-point mesh for β-KNaNbOF5, respectively. The results of the total energy and the Hellmann-Feynman forces are convergent within 10−4 meV and 1 meV/Å, respectively.
It has been reported that the [NbOF5]2− polar structure results in the special piezoelectric properties of KNaNbOF5 [12]. In order to describe this special structure, some parameters are used to describe the distorted octahedral [NbOF5]2− [19]. The Baur’s distortion index D [20], based on bond lengths, is defined as:
D = 1 n i = 1 n | l i l a v | l a v
where li is the distance from the central atom to the ith coordinating atom, and lav is the average bond length. The quadratic elongation, <λ> [21], is defined as:
< λ > = 1 n i = 1 n ( l i l 0 ) 2
where l0 is the center-to-vertex distance of a regular polyhedron of the same volume. <λ> is a dimensionless quantity that provides a quantitative measurement of the polyhedral distortion in a crystal, independent of the polyhedral size. The bond angle variance, σ2 [21], is calculated by:
σ 2 = 1 m 1 i = 1 m ( ϕ i ϕ 0 ) 2
where m is (number of faces in the polyhedron) × 3/2, ϕi is the ith bond angle, and ϕ0 is the ideal bond angle for a regular polyhedron.

3. Results and Discussion

The noncentrosymmetric α-KNaNbOF5 belongs to the orthorhombic system with the space group Pna21, while the centrosymmetric β-KNaNbOF5 is a tetragonal system with a space group of P4/nmm. The lattice parameters in [10,12] are: a = 11.8653(11) Å, b = 5.8826(6) Å, c = 8.1258(8) Å for α-KNaNbOF5; a = 5.9352(2) Å, c = 8.5487(5) Å for β-KNaNbOF5, respectively. Starting from these experimental parameters, the structures in our calculations are obtained by minimizing the total energy.

3.1. Atomic Structure and Bonding Properties

Orthorhombic (α-) and tetragonal (β-) KNaNbOF5 are investigated both within the local density approximation (LDA) and the generalized gradient approximation (GGA-PBE). Our estimations of the lattice parameters are listed in Table 1, with LDA, PBE, and experimental results listed for comparison. It is well known that LDA lever calculations usually underestimate and GGA-lever calculations overestimate the lattice constants, which also occurs in our calculations. Since the PBE lattice constants of the β phase are closer to the experimental values (the error is around 1.47%–1.77%), the PBE function is used for further calculations.
Table 1. The calculated lattice constants and the unit-cell volume in α- and β-KNaNbOF5. Δa, Δb and Δc are the relative errors of the calculated lattice constants compared to the experimental values. LDA = Local density approximation; PBE = Perdew, Burke, and Ernzerhof approximation; Exp. = Experimental values.
Table 1. The calculated lattice constants and the unit-cell volume in α- and β-KNaNbOF5. Δa, Δb and Δc are the relative errors of the calculated lattice constants compared to the experimental values. LDA = Local density approximation; PBE = Perdew, Burke, and Ernzerhof approximation; Exp. = Experimental values.
PhaseMethoda (Å)b (Å)c (Å)Volume (Å3)Δa (%)Δb (%)Δc (%)
α-phaseExp. [12]11.8655.8838.126567.171---
LDA11.6435.7797.994537.880−1.87−1.76−1.63
PBE12.0985.9828.300600.777+1.96+1.70+2.15
β-phaseExp. [10]5.9355.9358.549301.142---
LDA5.8165.8168.385283.629−2.01−2.01−1.91
PBE6.0406.0408.674316.502+1.77+1.77+1.47
It can be found that α- and β-KNaNbOF5 are very similar to each other, as shown in Figure 1. Both crystal structures consist of [NbOF5] and [NaOF5] octahedron sharing vertices to form three-dimensional networks, with the K atoms occupying the interstitial sites. In the framework of [NaOF5] and [NbOF5], the Na and Nb atoms are coordinated by five F and one O atoms in a distorted octahedral arrangement. There are also obvious differences between them in crystal structure. Firstly, in α-KNaNbOF5, the alternating NbOF5 octahedra share both vertices and F-F edges, while in β-KNaNbOF5, all the NbOF5 octahedra share only vertices. Secondly, the K atoms are coordinated by six F atoms and one O atom in α-KNaNbOF5, while the K atoms are coordinated by 8 F atoms in the β phase, without any surrounding O atoms in β-KNaNbOF5. Thirdly, in α-KNaNbOF5, all the atoms occupy the Wyckoff positions 4a, while in β-KNaNbOF5, the K atoms occupy the Wyckoff position 2a; the Na, Nb and O atoms occupy positions 2c; and the remaining inequivalent fluorides occupy two kinds of Wyckoff positions (2c and 8j).
In order to make a detailed analysis, the Baur’s distortion index D, the quadratic elongation <λ>, and the bond angle variance σ2 can be used to evaluate the distortion of the structure, which are listed in Table 2. The effective charge for each atom (charge difference after bonding) is determined using Bader charge analysis [22,23], which is given in Table 3 with the corresponding calculated atomic positions. It can be found that the calculated values agree well with the experimental data. From these values, several conclusions can be made: (1) The Na-O and Na-F bond lengths are larger than that of Nb-O and Nb-F. Thus, the [NaOF5] occupies a larger volume than [NbOF5]; (2) The bond angle variances, σ2, are larger in [NaOF5] than [NbOF5] octahedra, both in α- and β-KNaNbOF5; (3) The Nb-O bond length in the α-phase (1.776 Å) is noticeably longer than that in the β-phase (1.757 Å); (4) The ionic formula of both α- and β-KNaNbOF5 can be defined as K+0.90Na+0.90[NbOF]5−1.80; (5) The effective charges on fluorides are different due to their different surroundings (−0.791–−0.698 for α-phase, and −0.813–−0.708 for β-phase).
To analyze the electronic structure of KNaNbOF5, the electronic density of states (DOS), is illustrated in Figure 2. Although there are many differences in the structures between α- and β-KNaNbOF5, their density of states is similar to each other, indicating that their bonding characters are similar. The projected density of states (PDOS) of F, O and Nb are localized, indicating their ionic characteristics. Both α- and β-KNaNbOF5 are insulators with band gaps of about 4 eV. Although the band gap error is generally not definite due to the PBE approximation, the calculation method usually underestimates the value in most cases. Thus, it is very probable that the actual band gap is greater than 4 eV. Further experimental investigations might be necessary to confirm this result. From the PDOS, we find the peak near the Fermi surface (the top of the valence band) mainly consists of O-2p, F-2p and Nb-4d states. The −0.4 eV peak is mainly contributed by the F-2p state. Compared with the PDOS of F, O, and Nb, the PDOS of Na and K are negligible between −5 and 6 eV, as shown in Figure 2d,e. The bottom of the conduction band is mainly contributed by the Nb-4d states, which splits into two sub-bands. The peaks at 5.4 eV for the α-phase and at 5.5 eV for the β-phase should be the anti-bonding states contributed by O-2p and Nb-4d. This reflects a little covalent characteristic in the O-Nb interaction.
In order to further reveal the bonding properties of KNaNbOF5, the charge differences between the crystal charge and the atomic charge are depicted in Figure 3. The positive values mainly locate on the F and the O atoms, while the negative values mainly locate on the Nb atoms. This charge transfer indicates that the crystal is held together mainly by the ionic interactions between the Nb cations and the F/O anions. Moreover, the distribution of electrons gained by O-2p is obviously asymmetric. The electron charge densities on the O-Nb bond are higher, illustrating a little covalent bonding characteristic between the Nb and the O. All these results are consistent with our previous DOS analysis.
Figure 1. The crystal structures of KNaNbOF5 in different phases. The purple, golden, green, red and white spheres are K, Na, Nb, O and F atoms, respectively. The Na- and Nb-centered octahedra are also shown. (a) The orthorhombic noncentrosymmetric phase of KNaNbOF5 (α-phase, space group Pna21); (b) The tetragonal centrosymmetric phase of KNaNbOF5 (β-phase, space group P4/nmm). The crystal axes are also shown for clarity.
Figure 1. The crystal structures of KNaNbOF5 in different phases. The purple, golden, green, red and white spheres are K, Na, Nb, O and F atoms, respectively. The Na- and Nb-centered octahedra are also shown. (a) The orthorhombic noncentrosymmetric phase of KNaNbOF5 (α-phase, space group Pna21); (b) The tetragonal centrosymmetric phase of KNaNbOF5 (β-phase, space group P4/nmm). The crystal axes are also shown for clarity.
Materials 08 05477 g001
Figure 2. The calculated Density of States (DOS) of α- and β-KNaNbOF5 are plotted in (a), while the Projected Density of States (PDOSs) of Nb, O and F are plotted in (b,c). The negligible PDOSs of K and Na are also shown in (d,e) for comparison. The Fermi level is set to zero. The PDOSs of different atoms are plotted by the same colors as in Figure 1.
Figure 2. The calculated Density of States (DOS) of α- and β-KNaNbOF5 are plotted in (a), while the Projected Density of States (PDOSs) of Nb, O and F are plotted in (b,c). The negligible PDOSs of K and Na are also shown in (d,e) for comparison. The Fermi level is set to zero. The PDOSs of different atoms are plotted by the same colors as in Figure 1.
Materials 08 05477 g002
Figure 3. Deformation charge density (difference between the crystal charge and the atomic charge distribution) of KNaNbOF5 in (a) the α-phase; and (b) the β-phase. The yellow and blue isosurfaces (1.5 × 10−2 electrons/Bohr3) correspond to the electron increase and the depletion zone, respectively.
Figure 3. Deformation charge density (difference between the crystal charge and the atomic charge distribution) of KNaNbOF5 in (a) the α-phase; and (b) the β-phase. The yellow and blue isosurfaces (1.5 × 10−2 electrons/Bohr3) correspond to the electron increase and the depletion zone, respectively.
Materials 08 05477 g003
Table 2. The geometric properties of the [NbOF5] and [NaOF5] octahedra of α-NaNbOF5 and β-NaNbOF5, including the average bond length lav (in Å), the polyhedral volume V (in Å3), the distortion index D, the quadratic elongation <λ>, and the bond angle variance σ2 (in degree2). The values calculated from experimental data in [10,12] are listed in brackets.
Table 2. The geometric properties of the [NbOF5] and [NaOF5] octahedra of α-NaNbOF5 and β-NaNbOF5, including the average bond length lav (in Å), the polyhedral volume V (in Å3), the distortion index D, the quadratic elongation <λ>, and the bond angle variance σ2 (in degree2). The values calculated from experimental data in [10,12] are listed in brackets.
PhaseOctahedralavVD<λ>σ2
α-phase[NbOF5]1.98510.2150.0361.01748.784
(1.942)(9.592)(0.034)(1.016)(44.311)
[NaOF5]2.34216.7000.0201.01858.982
(2.300)(15.786)(0.027)(1.019)(63.760)
β-phase[NbOF5]1.99010.2300.0391.02365.121
(1.944)(9.546)(0.040)(1.022)(63.543)
[NaOF5]2.33716.6630.0101.01551.677
(2.310)(16.118)(0.011)(1.014)(48.468)
Table 3. The calculated atomic positions for KNaNbOF5 with Pna21 symmetry for the α-phase and P4/nmm for the β-phase, respectively. The calculated lattice parameters are listed in Table 1. WP refers to Wyckoff position, MP refers to multiplicity, and Q refers to the atomic effective charge. The experimental data in [10,12] are listed in brackets.
Table 3. The calculated atomic positions for KNaNbOF5 with Pna21 symmetry for the α-phase and P4/nmm for the β-phase, respectively. The calculated lattice parameters are listed in Table 1. WP refers to Wyckoff position, MP refers to multiplicity, and Q refers to the atomic effective charge. The experimental data in [10,12] are listed in brackets.
PhaseAtomWPMPxyzQ (e)
α-phaseKa40.0398(0.0385)0.4685(0.4680)0.0955(0.0927)0.903
Naa40.1373(0.1389)0.9509(0.9533)0.3401(0.3414)0.896
Nba40.8480(0.8492)0.0407(0.0420)0.3328(0.3334)2.796
Oa40.7219(0.7225)0.8998(0.9016)0.2919(0.2940)−0.977
F1a40.0023(0.0029)0.2167(0.2169)0.3701(0.3700)−0.791
F2a40.8731(0.8725)0.1485(0.1472)0.1093(0.1090)−0.702
F3a40.9529(0.9527)0.7849(0.7853)0.2979(0.2976)−0.720
F4a40.8595(0.8609)0.9625(0.9650)0.5655(0.5650)−0.708
F5a40.7874(0.7879)0.3353(0.3356)0.3956(0.3956)−0.698
β-phaseKa20.5000(0.5000)0.5000(0.5000)0.0000(0.0000)0.903
Nac20.0000(0.0000)0.5000(0.5000)0.2668(0.2660)0.891
Nbc20.0000(0.0000)0.5000(0.5000)0.7348(0.7351)2.768
Oc20.0000(0.0000)0.5000(0.5000)0.5322(0.5351)−0.915
F1c20.0000(0.0000)0.5000(0.5000)0.9891(0.9865)−0.813
F2j80.2304(0.2295)0.2696(0.2705)0.7723(0.7719)−0.708

3.2. Piezoelectricity and Acoustic Properties

The crystallographic symmetry of materials plays an important role in the piezoelectric phenomena. According to the definition of the piezoelectric effect, the piezoelectric tensor of α-KNaNbOF5 is:
[ 0 0 0 0 e x 5 0 0 0 0 e y 4 0 0 e z 1 e z 1 e z 3 0 0 0 ]
The calculated piezoelectric components of α-KNaNbOF5 are listed in Table 4, in which the data of ZnO are also listed. It is well known that ZnO has good piezoelectric properties, and has been widely used in filters for incoming television signals [24,25]. ZnO ceramics are also widely used as varistors for surge protection [26]. Compared with ZnO, the piezoelectric constants of α-KNaNbOF5 are much lower. In order to make a comparison to the experimental piezoelectric response measurement reported in [9], the calculated piezoelectric stress matrix [e] is converted into the piezoelectric strain matrix [d] by the relationship [e] = [c][d], where [c] is the elastic matrix. Based on Table 4 and Table 5, the calculated dz3 is 1.5 pC/N, which is lower than the experimental value 6.7 pC/N [9]. Since our calculated lattice parameters, atomic positions, and relative permittivity agree well with experiments, the errors in the calculated dz3 may originate from three aspects: (1) The experiment was carried out at room temperature—since finite temperatures are not included in the calculation, the calculated value of dz3 is smaller than that observed in experiments; (2) The value of dz3 is sensitive to the crystal structure—the calculated lattice constants deviate from the experimental values by ~1%, which affects the evaluation of dz3; (3) Compared to other commonly used piezoelectric materials, such as BaTiO3 (~190 pC/N) and K0.5Na0.5NbO3 (~160 pC/N) [9], the dz3 value of KNaNbOF5 is much lower. Thus, although the relative error of the calculated dz3 is large, the absolute error is only 4.8 pC/N. From these values, we can conclude that the piezoelectricity of α-KNaNbOF5 is not very strong. This is consistent with the fact that the piezoelectricity of α-KNaNbOF5 arises from the competition between primary and secondary distortions [11,12]. In the centrosymmetric structure of β-KNaNbOF5, the calculated piezoelectric components are all zero as expected from the symmetry of the crystal.
Table 4. The calculated piezoelectric tensor elements eij of α-KNaNbOF5. The data of ZnO are listed for comparison. The piezoelectric components of β-KNaNbOF5 are all zero.
Table 4. The calculated piezoelectric tensor elements eij of α-KNaNbOF5. The data of ZnO are listed for comparison. The piezoelectric components of β-KNaNbOF5 are all zero.
eij (C/m2)ex5ey4ez1ez2ez3
α-KNaNbOF50.110.12−0.07−0.080.05
ZnO [27]-−0.59−0.611.14-
Table 5. The calculated elastic constants of the two phases of KNaNbOF5 (in GPa).
Table 5. The calculated elastic constants of the two phases of KNaNbOF5 (in GPa).
PhaseC11C12C13C22C23C33C44C55C66
α-phase50.926.813.851.219.764.211.017.120.2
β-phase51.125.622.5--63.010.3-26.7
In experiments, the ultrasonic wave velocities of certain directions can be measured to obtain the elastic constant by the Christoffel equation [27,28]. The Christoffel equation for β-KNaNbOF5 is written as:
( l i K C K L l L j δ j k ρ V 2 ) ( α j ) = 0   ( j ,   k = 1 , 2 , 3 )
where CKL is the elastic constant (from Table 5), l is the propagation matrix, δjk is the Kronercker sign, ϱ is the density of crystal, αj is the eigenvector, and V is the velocity of acoustic wave.
Since α-KNaNbOF5 has piezoelectric property, the quasistatic approximation is used to transform the conventional Christoffel equation to the stiffened Christoffel equation. In this approximation, the effect of the quasistatic electric field is saved. The stiffened Christoffel equation is written as:
k 2 { l i K [ C K L + ( e K j l j ) ( l i e i L ) l i ( ε i j e i I S I J e J i ) l j ] l L j } α j = ρ ω 2 α i
where CKL is the elastic constant (from Table 5), SIJ corresponds to the components of the elastic compliance matrix, ϱ is the density of crystal (from Table 1), ω is the phase velocity of the acoustic wave, α is the eigenvector, e is the piezoelectric stress tensor (from Table 4), ε is the dielectric constant (from Table 6), and l is the propagation matrix. All of these values are obtained from the first-principles calculations.
Table 6. Values of electronic and ionic contributions and total value of the relative dielectric permittivity (εij) for α- and β-KNaNbOF5.
Table 6. Values of electronic and ionic contributions and total value of the relative dielectric permittivity (εij) for α- and β-KNaNbOF5.
PhaseεijElectronicIonicTotal
α-phaseεxx4.052.246.29
εyy3.702.165.86
εzz3.682.135.81
β-phaseεxx4.742.036.77
εyy4.742.036.77
εzz2.512.194.70
The elastic stiffness (CKL) defines the resistance of the material to undergo strain under the action of a mechanical stress within the elastic regime. The orthorhombic (α-KNaNbOF5) and the tetragonal crystal (β-KNaNbOF5) have nine and six independent elastic constants, respectively. The calculated elastic constants are listed in Table 5. The relative dielectric permittivities for α-KNaNbOF5 are εxx = 6.29, εyy = 5.86 and εzz = 5.81. For the β phase, they are εxx = εyy = 6.77 and εzz = 4.70, as shown in Table 6. These calculated results are close to the experimental measured relative permittivity (ε = 7.1) for α-KNaNbOF5 [9].
For a given direction, three velocities are determined by solving the Christoffel and the stiffened Christoffel equation. By changing the propagation direction, we obtain the velocity as a function of propagation direction. Usually, in practice, it is more convenient to use the slowness surface (the inverse of the velocity). We calculate the slowness curves of KNaNbOF5 for both the α- and β-phases and illustrate them in Figure 4 (for the tetragonal phase, the (010) plane is the same as the (001) plane). It is clear that there are three acoustic waves in each direction: one quasi-longitudinal mode and two quasi-transverse modes. The velocities of quasi-longitudinal modes are much larger than those of the other two in both phases of KNaNbOF5. These acoustic properties of KNaNbOF5 significantly change when a phase transition occurs: (1) For α-KNaNbOF5, the largest velocity is the quasi-longitudinal mode along the (001) direction, indicating that the velocity reaches maximum along this direction, while for β-KNaNbOF5, the velocity is at a maximum along (110) and (001); (2) For α-KNaNbOF5, the maximum values of the slowness surface are along (100) and (001), while for β-KNaNbOF5, all the directions in the (001) plane have the same maximum values of slowness; (3) All the acoustic waves of α-KNaNbOF5 are anisotropic, however, in the (001) plane of β-KNaNbOF5, the propagation of the outer acoustic wave is obviously isotropic. This mode is a pure transverse mode, in which the polarization is normal to the direction of propagation. (4) The slowness curves are very different on the (100) plane between α- and β-phases.
Figure 4. Slowness curves of acoustic waves for α- and β-KNaNbOF5. The slowness curves of α-KNaNbOF5 in (001), (010) and (100) planes are (a), (b) and (c), while those of β-KNaNbOF5 in (001) and (010)/(100) planes are (d) and (e), respectively. The inner blue circle represents the quasi-longitudinal wave and the middle green and outer red circles represent the other two quasi-transverse waves. The units of the slowness values are 10−4 s/m.
Figure 4. Slowness curves of acoustic waves for α- and β-KNaNbOF5. The slowness curves of α-KNaNbOF5 in (001), (010) and (100) planes are (a), (b) and (c), while those of β-KNaNbOF5 in (001) and (010)/(100) planes are (d) and (e), respectively. The inner blue circle represents the quasi-longitudinal wave and the middle green and outer red circles represent the other two quasi-transverse waves. The units of the slowness values are 10−4 s/m.
Materials 08 05477 g004

3.3. Mechanical Properties

From the elastic constants (CKL) listed in Table 5, more elastic-relevant properties can be obtained. The mechanical stability of KNaNbOF5 can be determined by the Born-Huang criterion [29]. For the α-phase, they are:
C 11 > 0 , C 22 > 0 , C 33 > 0 , C 44 > 0 , C 55 > 0 , C 66 > 0 , [ C 11 + C 22 + C 33 + 2 ( C 12 + C 13 + C 23 ) ] > 0 , ( C 11 + C 22 2 C 12 ) > 0 , ( C 11 + C 33 2 C 13 ) > 0 , ( C 22 + C 33 2 C 23 ) > 0
And for the β-phase, they are:
C 11 > 0 , C 33 > 0 , C 44 > 0 , C 66 > 0 , C 11 C 12 > 0 , C 11 + C 33 2 C 13 > 0 , 2 ( C 11 + C 12 ) + C 33 + 4 C 13 > 0
The calculated elastic constants listed in Table 5 satisfy the mechanical stability criterion. Using the calculated elastic constants, we further calculate the bulk modulus, shear modulus, Young’s modulus, and Poisson ratio to obtain a complete description of the mechanical behavior. The Voigt and the Reuss approximations [30] are used to estimate the bulk modulus and the shear modulus.
Since α-KNaNbOF5 belongs to the orthorhombic structure, the bulk and the shear moduli (B and G) for orthorhombic structures can be obtained from the general expression of the bulk and the shear modulus in Voigt-Reuss-Hill approximation. For the Voigt approximation, they are:
B V = [ ( C 11 + C 22 + C 33 ) + 2 ( C 12 + C 13 + C 23 ) ] / 9
G V = [ ( C 11 + C 22 + C 33 ) ( C 12 + C 13 + C 23 ) + 3 ( C 44 + C 55 + C 66 ) ] / 15
While by the Reuss approximation, they are:
B R = Δ / [ C 11 ( C 22 + C 33 2 C 23 ) + C 22 ( C 33 2 C 13 ) 2 C 33 C 12 + C 12 ( 2 C 23 C 12 ) + C 13 ( 2 C 12 C 13 ) + C 23 ( 2 C 13 C 23 ) ]
G R = 15 / { 4 [ C 11 ( C 22 + C 33 + C 23 ) + C 22 ( C 33 + C 13 ) + C 33 C 12 C 12 ( C 23 + C 12 ) C 13 ( C 12 + C 13 ) C 23 ( C 13 + C 23 ) ] / Δ + 3 ( 1 / C 44 + 1 / C 55 + 1 / C 66 ) }
Δ = C 13 ( C 12 C 23 C 13 C 22 ) + C 23 ( C 12 C 13 C 23 C 11 ) + C 33 ( C 11 C 22 C 12 2 )
where B and G are the bulk and the shear moduli obtained by Voigt or Reuss approximations, respectively (with subscripts V or R).
The β-KNaNbOF5 belongs to the tetragonal structure. Therefore, the formulas become:
B V = [ ( C 11 + C 22 + C 33 ) + 2 ( C 12 + C 13 + C 23 ) ] / 9
G V = [ ( C 11 + C 22 + C 33 ) ( C 12 + C 13 + C 23 ) + 3 ( C 44 + C 55 + C 66 ) ] / 15
B R = Δ / [ C 11 ( C 22 + C 33 2 C 23 ) + C 22 ( C 33 2 C 13 ) 2 C 33 C 12 + C 12 ( 2 C 23 C 12 ) + C 13 ( 2 C 12 C 13 ) + C 23 ( 2 C 13 C 23 ) ]
G R = 15 / { 4 [ C 11 ( C 22 + C 33 + C 23 ) + C 22 ( C 33 + C 13 ) + C 33 C 12 C 12 ( C 23 + C 12 ) C 13 ( C 12 + C 13 ) C 23 ( C 13 + C 23 ) ] / Δ + 3 ( 1 / C 44 + 1 / C 55 + 1 / C 66 ) }
Δ = C 13 ( C 12 C 23 C 13 C 22 ) + C 23 ( C 12 C 13 C 23 C 11 ) + C 33 ( C 11 C 22 C 12 2 )
Within the Voigt-Reuss-Hill approximation [30], the bulk modulus B and the shear modulus G are the average of the values by Voigt and Reuss approximations. The Young’s modulus (E), the Poisson ratio (v), the velocity of the transverse (Vt), the longitudinal (Vl) acoustic wave, the average velocity (Va), and the Debye temperature (ΘD) can be obtained by:
E = 9 B G / ( 3 B + G ) , v = ( 3 B 2 G ) / ( 6 B + 2 G )
V l = ( 3 B + 4 G ) / 3 ρ , V t = G / ρ , V a = [ 1 3 ( 2 V t 3 + 1 V l 3 ) ] 1 / 3
Θ D = h k B [ 3 n 4 π ( N A ρ M ) ] 1 / 3 V a
where h is the Planck constant, kB is the Boltzmann constant, n is the number of atoms in the formula unit, NA is the Avogadro number, ϱ is the density of the crystal, and M is the molecular weight. The constants describing the mechanical behavior of KNaNbOF5 in the elastic regime are given in Table 7.
The bulk modulus B of a material determines the resistance to compression under a given hydrostatic pressure [31]. The shear modulus G describes the resistance of a material to deform under a shear stress. The calculated bulk modulus of KNaNbOF5 is only 31.9–34.0 GPa, which is much lower than steel (approximately 160 GPa) and comparable to glass (35–55 GPa) [32]. Therefore, neither α- nor β-KNaNbOF5 are hard materials. The bulk modulus B is twice the shear modulus, indicating that the parameter limiting the material’s stability is the shear modulus.
For the application of KNaNbOF5, the brittle or ductile behavior is of great importance. The ratio of bulk modulus to shear modulus (B/G) is frequently used to discriminate the brittle properties of materials [33]. According to the criterion given by Pugh [34], a material is brittle if the B/G ratio is less than 1.75. Otherwise, it behaves in a ductile manner. In the case of KNaNbOF5, the B/G ratio is 1.99–2.31, indicating that KNaNbOF5 is predominantly ductile. Since the B/G ratio of the β-phase is larger than that of the α-phase, we predict that β-KNaNbOF5 is more ductile.
Table 7. The calculated elasticity-relevant properties of three different phases of KNaNbOF5. The bulk modulus B and the shear modulus G within the Voigt and Reuss approximation (with subscripts V and R, respectively) are listed. Based on B and G, the Young’s modulus E, the Poisson ration v, the acoustic wave velocities (Vl, Vt, Va), and the Debye temperature ΘD are calculated (referring to Equations (9)–(21)).
Table 7. The calculated elasticity-relevant properties of three different phases of KNaNbOF5. The bulk modulus B and the shear modulus G within the Voigt and Reuss approximation (with subscripts V and R, respectively) are listed. Based on B and G, the Young’s modulus E, the Poisson ration v, the acoustic wave velocities (Vl, Vt, Va), and the Debye temperature ΘD are calculated (referring to Equations (9)–(21)).
Propertyα-Phaseβ-Phase
BV GPa)31.934.0
BR (GPa)31.833.8
B (GPa)31.833.9
GV (GPa)16.715.7
GR (GPa)15.413.7
G (GPa)16.014.7
E (GPa)41.238.6
v0.2840.310
Vl (m/s)4253.14381.4
Vt (m/s)2335.52297.3
Va (m/s)2603.32569.3
ΘD (K)303.3294.2
At low temperatures, the vibrational excitations arise solely from the acoustic modes of the phonon spectrum. Hence, the Debye temperature calculated using the elastic constants can be compared to the experimentally measured values. From Table 7, we find that the α-phase has a higher velocity of the longitudinal acoustic wave and a lower velocity of the transverse acoustic wave compared to the β-phase. The calculated velocity of the longitudinal acoustic wave is 4253–4381 m/s and the transverse acoustic wave is 2297–2336 m/s. The calculated velocities using the classic Debye model compare well with the velocities calculated from the Christoffel equation discussed above. The Debye temperature from the calculated acoustic velocities is 303 K for α-KNaNbOF5 and 294 K for β-KNaNbOF5. If the Debye temperature is low, the thermal conductivity is also low for insulators in general [31]. Heat is transported by two mechanisms in solids: lattice vibrations (phonons) and free electrons. Since KNaNbOF5 is an insulator, the scattering becomes quite large above the Debye temperature, making the solid a poor thermal conductor. The Debye temperature of KNaNbOF5 (294–303 K) is much lower than that of diamond (approximately 2230 K), indicating a poor thermal conductivity of the material.

4. Conclusions

In this work, we used density functional theory to study the electronic, piezoelectric, mechanical, and acoustic properties of the recently synthesized KNaNbOF5 crystal. The bonding and mechanical properties of α- and β-KNaNbOF5 are very similar. KNaNbOF5 is found to present a strong ionic characteristic with weak covalent bonding. The ionic formula can be defined as K+0.90Na+0.90[NbOF]5−1.80. Both α- and β-KNaNbOF5 are insulators with an energy band gap of about 4 eV. The top of the valence band mainly consists of O-2p, F-2p and Nb-4d states, while the bottom of the conduction band is mainly contributed by the Nb-4d states, which split into two sub-bands. The B/G ratio of KNaNbOF5 is 1.99–2.31, which is higher than the critical value of 1.75, indicating that it is predominantly ductile. The calculated bulk modulus of KNaNbOF5 is only 31.8–33.9 GPa, which is comparable to glass. In contrast, the calculated piezoelectric and acoustic properties of the two phases are very different. The small values of the piezoelectric components of α-KNaNbOF5 illustrate its poor piezoelectricity compared with ZnO, while the values of β-KNaNbOF5 are all zero, as expected from the symmetry of the crystal. Based on the quasistatic approximation, the slowness surface of the acoustic waves is calculated so as to describe its acoustic properties. A pure transverse mode is found in the (001) plane of β-KNaNbOF5, in which the propagation of the outer acoustic wave is obviously isotropic. These calculated results are an initial step towards characterizing the properties of KNaNbOF5.

Acknowledgments

We acknowledge the financial support by the Program of International S&T Cooperation (Grant No. 2014DFG60230), National Natural Science Foundation of China (Nos. 91326105, U1404111, 11504089, 21501189, 21571185, 21306220), the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA02040104).

Author Contributions

Han Han and Ping Huai conceived and designed the whole study. Han Han and Cheng Cheng performed DFT calculations. Han Han, Xiao-Gen Xiong, Jing Su and Jian-Xing Dai analyzed the data. Han Han and Hui Wang wrote a program for acoustic properties analysis and co-wrote the manuscript. Ping Huai and Gen Yin helped with correction, modification, and revision of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ye, N.; Chen, Q.; Wu, B.; Chen, C. Searching for new nonlinear optical materials on the basis of the anionic group theory. J. Appl. Phys. 1998, 84, 555–558. [Google Scholar] [CrossRef]
  2. Zhang, S.J.; Yu, F.P. Piezoelectric materials for high temperature sensors. J. Am. Ceram. Soc. 2011, 94, 3153–3170. [Google Scholar] [CrossRef]
  3. Schlesinger, T.E.; Toney, J.E.; Yoon, H.; Lee, E.Y.; Brunett, B.A.; Franks, L.; James, R.B. Cadmium zinc telluride and its use as a nuclear radiation detector material. Mater. Sci. Eng. 2001, 32, 103–189. [Google Scholar] [CrossRef]
  4. Flerov, I.N.; Gorev, M.V.; Tressaud, A.; Laptash, N.M. Perovskite-like fluorides and oxyfluorides: Phase transitions and caloric effects. Crystallogr. Rep. 2011, 56, 9–17. [Google Scholar] [CrossRef]
  5. Stramare, S.; Thangadurai, V.; Weppner, W. Lithium lanthanum titanates: A review. Chem. Mater. 2003, 15, 3974–3990. [Google Scholar] [CrossRef]
  6. Rooksby, H.P. Compounds of the structural type of calcium titanate. Nature 1945, 155, 484. [Google Scholar] [CrossRef]
  7. Wang, H.; Che, J.G. A route to multiferroics by non-d0 cation B in magnetic perovskites. Eur. Phys. Lett. 2011, 96, 67012. [Google Scholar] [CrossRef]
  8. Han, H.; Che, J.G. Dependence of silicon oxidation channel on distribution of surface electrons at initial stage of oxide growth on Si(001). Appl. Phys. Lett. 2013, 103, 163113. [Google Scholar] [CrossRef]
  9. Chang, K.B.; Edwards, B.W.; Frazer, L.; Lenferink, E.J.; Stanev, T.K.; Stern, N.P.; Nino, J.C.; Poeppelmeier, K.R. Hydrothermal crystal growth, piezoelectricity, and triboluminescence of KNaNbOF5. J. Solid State Chem. 2015, in press. [Google Scholar] [CrossRef]
  10. Vasiliev, A.D.; Laptash, N.M. Polymorphism of KNaNbOF5 crystals. J. Struct. Chem. 2012, 53, 902–906. [Google Scholar] [CrossRef]
  11. Mishra, A.K.; Marvel, M.R.; Poeppelmeier, K.R.; Waghmare, U.V. Competing Cation-anion interactions and noncentrosymmetry in metal oxide-fluorides: A first-principles theoretical study. Cryst. Growth Des. 2014, 14, 131–139. [Google Scholar] [CrossRef]
  12. Marvel, M.R.; Lesage, J.; Baek, J.; Halasyamani, P.S.; Stern, C.L.; Poeppelmeier, K.R. Cation-anion interactions and polar structures in the solid state. J. Am. Chem. Soc. 2007, 129, 13963–13969. [Google Scholar] [CrossRef] [PubMed]
  13. Antokhina, T.F.; Ignat’eva, L.N.; Savchenko, N.N.; Tkachenko, I.A.; Kaidalova, T.A. Synthesis and physicochemical properties of new niobium fluoro complexes with mixed alkali-metal cations. Russ. J. Inorg. Chem. 2003, 48, 472–479. [Google Scholar]
  14. Pinlac, R.A.F.; Stern, C.L.; Poeppelmeier, K.R. New layered oxide-fluoride perovskites: KNaNbOF5 and KNaMO2F4 (M = Mo6+, W6+). Crystals 2011, 1, 3–14. [Google Scholar] [CrossRef]
  15. Kresse, G.; Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169–11186. [Google Scholar] [CrossRef]
  16. Kresse, G.; Joubert, J. Fromultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 1999, 59, 1758–1775. [Google Scholar] [CrossRef]
  17. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [PubMed]
  18. Monkhorst, H.J.; Pack, J.D. Special points for Brillouin-zone integrations. Phys. Rev. B 1976, 13, 5188–5192. [Google Scholar] [CrossRef]
  19. Han, H.; Yin, G.; Wickramaratne, D. A first-principles investigation of the electronic, elastic, piezoelectric and acoustic properties of K3B6O10C. Comput. Mater. Sci. 2013, 69, 81–86. [Google Scholar] [CrossRef]
  20. Baur, W.H. The geometry of polyhedral distortions. Predictive relationships for the phosphate group. Sect. B: Struct. Sci. 1974, 30, 1195–1215. [Google Scholar] [CrossRef]
  21. Robinson, K.; Gibbs, G.V.; Ri, P.H. Quadratic elongation: A quantitative measure of distortion in coordination polyhedra. Science 1971, 172, 567–570. [Google Scholar] [CrossRef] [PubMed]
  22. Bader, R.F.W. Atoms in Molecules: A Quantum Theory; Oxford University Press: Oxford, UK, 1990. [Google Scholar]
  23. Henkelman, G.; Arnaldsson, A.; Jónsson, H. A fast and robust algorithm for Bader decomposition of charge density. Comput. Mater. Sci. 2006, 36, 354–360. [Google Scholar] [CrossRef]
  24. Corso, A.D.; Posternak, M.; Resta, R.; Baldereschi, A. Ab initio study of piezoelectricity and spontaneous polarization in ZnO. Phys. Rev. B 1994, 50, 10715–10721. [Google Scholar] [CrossRef]
  25. Catti, M.; Noel, Y.; Dovesi, R. Full piezoelectric tensors of wurtzite and zinc blende ZnO and ZnS by first-principles calculations. J. Phys. Chem. Solids 2003, 64, 2183–2190. [Google Scholar] [CrossRef]
  26. Hutson, A.R. Piezoelectricity and conductivity in ZnO and CdS. Phys. Rev. Lett. 1960, 4, 505–507. [Google Scholar] [CrossRef]
  27. Auld, B.A. Acoustic Fields and Waves in Solids Vol. I; John Wiley & Sons Press: New York, NY, USA, 1973. [Google Scholar]
  28. Kim, K.Y.; Bretz, K.C.; Every, A.G.; Sachse, W.G. Ultrasonic imaging of the group velocity surface about the cubic axis in silicon. J. Appl. Phys 1996, 79, 1857–1863. [Google Scholar] [CrossRef]
  29. Born, M.; Huang, K. Dynamical Theory of Crystal Lattices; Clarendon Press: Oxford, UK, 1998. [Google Scholar]
  30. Hill, R. The Elastic behaviour of a crystalline aggregate. Proc. Phys. Soc. Lond. 1952, 65, 349–354. [Google Scholar] [CrossRef]
  31. Kittel, C. Introduction to Solid State Physics, 8th ed.; John Wiley & Sons Inc.: New York, NY, USA, 2005; pp. 47–128. [Google Scholar]
  32. Han, H. Density-functional theory study of the effect of pressure on the elastic properties of CaB6. Chin. Phys. B 2013, 22, 077101. [Google Scholar] [CrossRef]
  33. Frantsevich, I.N.; Voronov, F.F.; Bokuta, S.A. Elastic Constants and Elastic Moduli of Metals and Insulators Handbook; Naukova Dumka: Kiev, Ukraine, 1982. [Google Scholar]
  34. Pugh, S.F. Relations between the elastic moduli and the plastic properties of polycrystalline pure metals. Philos. Mag. 1954, 45, 823–843. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Han, H.; Cheng, C.; Xiong, X.-G.; Su, J.; Dai, J.-X.; Wang, H.; Yin, G.; Huai, P. Piezoelectric, Mechanical and Acoustic Properties of KNaNbOF5 from First-Principles Calculations. Materials 2015, 8, 8578-8589. https://doi.org/10.3390/ma8125477

AMA Style

Han H, Cheng C, Xiong X-G, Su J, Dai J-X, Wang H, Yin G, Huai P. Piezoelectric, Mechanical and Acoustic Properties of KNaNbOF5 from First-Principles Calculations. Materials. 2015; 8(12):8578-8589. https://doi.org/10.3390/ma8125477

Chicago/Turabian Style

Han, Han, Cheng Cheng, Xiao-Gen Xiong, Jing Su, Jian-Xing Dai, Hui Wang, Gen Yin, and Ping Huai. 2015. "Piezoelectric, Mechanical and Acoustic Properties of KNaNbOF5 from First-Principles Calculations" Materials 8, no. 12: 8578-8589. https://doi.org/10.3390/ma8125477

Article Metrics

Back to TopTop