Molecular Orbital Calculation of Lead-Free Perovskite Compounds for Efficient Use of Alkaline and Alkaline Earth Metals

Takesue, Naohisa; Saito, Jun-ichi

doi:10.3390/cryst10110956

Open AccessArticle

Molecular Orbital Calculation of Lead-Free Perovskite Compounds for Efficient Use of Alkaline and Alkaline Earth Metals

by

Naohisa Takesue

^1,*

and

Jun-ichi Saito

²

¹

Department of Applied Physics, Graduate School of Science, Fukuoka University, Fukuoka 814-0180, Japan

²

Japan Atomic Energy Agency, Tsuruga 919-1279, Japan

^*

Author to whom correspondence should be addressed.

Crystals 2020, 10(11), 956; https://doi.org/10.3390/cryst10110956

Submission received: 28 August 2020 / Revised: 20 October 2020 / Accepted: 20 October 2020 / Published: 22 October 2020

(This article belongs to the Special Issue Multifunctional Crystal Materials with Non-centrosymmetry: Piezoelectric and Nonlinear Optical Properties and Their Applications)

Download

Browse Figures

Versions Notes

Abstract

:

The effective ionic charges of lead-free perovskite dielectric complex compounds were investigated with molecular orbital calculation. The base model was a double perovskite cluster that consisted of octahedral oxygen cages with a transition metal ion of titanium, niobium, or zirconium located at each of their centers, and alkali and/or alkaline earth metal ions located at the body center, corners, edge centers, or face centers of the cluster. The results showed significant covalent bonds between the transition metals and the oxygens, and the alkali metals, especially sodium and oxygen. On the other hand, the alkaline earth metals have weak covalency. Calculation was also performed with the replacement of some of the oxygens with chlorine or fluorine; such replacement enhances the covalency of the transition metals. These trends provide good guidelines for the design properties of lead-free perovskite piezoelectrics based on ubiquitous sodium use.

Keywords:

perovskite; dielectrics; piezoelectrics; ferroelectrics; molecular orbital calculation

1. Introduction

Lead (Pb)-based piezoelectrics have exclusively high electric polarization and piezoelectricity [1,2,3,4,5]. The typical substance is lead titanate, designated as PbTiO₃ (PT) [6,7]. Its crystal structure is of a typical simple perovskite, which is denoted as ABO₃ in Figure 1, where A, B, and O are ionic species; A and B each corresponds to cations of lead and titanium, and O is an anion of oxygen that forms an octahedral cage with the five other O ions. The chemical bonding of lead titanate was studied with first-principles calculation by R. E. Cohen [8]. The results indicated that the lead (Pb) ion is significantly covalent with the oxygen (O) and titanium (Ti) ions. This feature suggests that the covalency forms a molecular orbital (i.e., electronic channels) through which valence electrons are transferred back and forth between the cations and anions in response to external fields. The response may be soft or hard, depending on the significance of the covalency. Thus, covalency plays an important role in causing the piezoelectric strain, which is a reason PT and other Pb-based systems have excellent properties. Therefore, other ionic species that cause covalency at the site of Pb can be substitutes for harmless, lead-free piezoelectrics and their promising application. One of the authors, N. Takesue [9], confirmed with molecular orbital calculation that (Bi_1/2Na_1/2)TiO₃ (BNT) [10,11] could be such a substitute. This substance is now widely recognized as a potential substitute and has been studied through many types of fabrication [11,12,13], long before N. Takesue’s confirmation [9]. In addition to BNT, many excellent Pb-free piezoelectrics have been discovered [14,15,16].

For piezoelectric perovskites, alkaline and alkaline-earth metals are traditionally regarded as representative ionic species occupying the site of Pb or A [17]. Although there are more modern systems such as BNT [10,11,12,13] and BNT-based solid solutions [18], we focused on traditional systems, namely, BiFeO₃ [19] and GaFeO₃ [20], since our work was on the efficient use of ubiquitous alkaline and alkaline-earth metals. Both species of traditional systems are ionic; their formal ionic charges are 1+ and 2+, respectively. However, this idea indicates only the valence charge of neutral atoms. A portion of the charge is transferred to neighboring atoms in substances consisting of more than one atom. Therefore, atoms—even those of alkaline and alkaline-earth metals—are not neutral in general. Rather, they are partially covalent and ionic, especially for metals with high-electronegativity atoms such as oxygen atoms as neighbors. This idea suggests that a table of effective ionic charges enables the control of the properties of lead-free piezoelectrics [8,21] and is considered compensative to tolerance factors [22,23], which is a well-established criterion for the stability of perovskite structures.

Among the representative species, alkaline metals (sodium [Na] and potassium [K]) and alkaline-earth metals (calcium [Ca] and barium [Ba]) are ubiquitous; in particular, sodium contained in marine resources is abundant. These are constituents of perovskite piezoelectrics such as NaNbO₃ (Nb is niobium), KNbO₃, CaTiO₃, and BaTiO₃, and they occupy the Pb site. We included BaZrO₃ (zirconium [Zr] is chemically similar to Ti) in this category despite its cubic crystal structure (i.e., nonpolar) because it is frequently used as a solute component for this kind of piezoelectrics for their solid solutions.

The chemical bonding of BaTiO₃ was also studied by R. E. Cohen [8]. The results showed that an effective charge of Ba was about 2+ (i.e., ionic) and a strong covalency of Ti with O, similar to the case of PT, was observed. However, the other elements including those in the alkaline and alkaline-earth columns of the periodic table, have not been considerably investigated. Na and Ca may be covalent because of their higher electronegativity than that of Ba. If their effective charges show an evident systematic variation along the two columns, then the piezoelectricity of these kinds of substances can be controlled through the synthesis of solid solutions and/or heterogeneous structuralization such as coherent interfaces; the way(s) may be realized by matching the one base piezoelectric with the other piezoelectric component(s). The successful case leads us to certain achievements, after which lead-free substance(s) with properties are as good as, or better than, those of Pb-based systems will become obtainable. This scheme motivated us to investigate the effective charges of perovskites with alkaline and/or alkaline-earth metals located at the Pb site including the charges of PT as the reference data through molecular orbital calculation. Then, we performed a calculation including the partial replacement of O with a halogen such as chlorine (Cl) and fluorine (F), or nitrogen (N) on the other side to investigate their electronegative effects on the effective charges of the cations.

2. Calculation

A discrete-variational (DV) Xα method of molecular orbital calculation [24,25] was used throughout this study by using the SCAT code for calculating the electronic states. First, the covalency and ionicity of alkaline metals lithium (Li), Na, and K (designated as M_A) were investigated based on a fundamental compound cluster model (M_A)₂O with an antifluorite structure, as illustrated in Figure 2a. Previously reported lattice parameters [26] were used for these calculations. Second, the alkaline-earth metals magnesium (Mg), Ca, strontium (Sr), and Ba (designated as M_AE) were investigated based on a cluster model (M_AE)O with a sodium chloride structure, as illustrated in Figure 2b. The same lattice parameters [26] were used for these calculations.

In the third step, the perovskite was investigated based on double perovskite cluster models for PT, M_ANbO₃, M_AETiO₃, BaZrO₃, (Ba, Na)(Ti, Nb)O₃, (Ba, Na)Ti(O, Cl)₃, Na(Ti, Nb)(O, Cl)₃, (Ba, Na)Ti(O, F)₃, Na(Ti, Nb)(O, F)₃, (Ba, Na)Nb(O, N)₃, and Ba(Ti, Nb)(O, N)₃, as illustrated in Figure 2c. The sites of ions in each model are indicated by colored spheres, and the correspondence of the colors to the ions is summarized in Table 1. Each model is a rectangular cluster and was given a nonpolar symmetry D₂ (Schoenflies notation), common to all perovskite models employed; the reason is that the preliminary calculation with ferroelectric strains with a smaller symmetry than D₂ indicated no essential effect on the bonds [9]. Furthermore, this simplified modeling reduced the time needed to complete the calculations while providing the essential chemical features. A typical value of the primitive cubic perovskite, 0.4 nm [7,17], was employed for the lattice parameters, so each of the three edge lengths of each model was given a value of 0.8 nm. The calculation on only PT was followed by a calculation based on its model, which incrementally provided an isotropic compression and expansion of 95%–105%; therefore, the above-mentioned effect of the electronic channels can be investigated [8,21].

All the above calculations were performed with the Madelung potential applied to each model, and Mulliken population analysis was conducted to determine the effective charges [24,25,27,28]. The lattice parameters used for modeling are listed in Table 2.

3. Results

The effective charges of the cations of each (M_A)₂O is listed in Table 3 in a downward sequence of the M_A column in the periodic table. The charges are shown as a histogram in Figure 3a. Each effective charge is the average of the eight M_A values in each model. The table and the figure indicate that the effective charges were noticeably below the formal charge of +1; in particular, the difference for sodium was significant. This evident charge transfer is considered to be caused by a large difference in electronegativity between neighboring M_A and O; this transfer can be used to control the covalency and ionicity of the perovskite A site.

The effective charges of the cations of each (M_AE)O are listed in Table 4 in the same manner as those of (M_A)₂O and are drawn as a histogram in Figure 3b. Two types of sodium chloride cluster models, namely, M_AE- and O-body-centered, for each (M_AE)O were employed. For the former model, the effective charge of M_AE at the body center, bonded with the six nearest neighbor oxygens, was employed.

For the latter model, each of the eight M_AE species at the corners was bonded with the three nearest neighbors, and each of the six at the face centers with the five; thus, the effective charges at these two non-body-centered locations were multiplied by 2 and 6/5, respectively, and the weighted values were averaged over each model.

The table and the figure indicate that the effective charges of M_AE increased downward along the column toward the formal charge of +2. This finding differed from those of M_A and was also related to the possible electronegativity of M_AE and O. Therefore, the results of M_AE may provide another aspect to control the chemical bonds of the perovskite A site.

The effective charges of the cations of the double perovskites of PT, M_ANbO₃, M_AETiO₃, and BaZrO₃ are listed in Table 5 and drawn as a histogram in Figure 4. The effective charge of each of the A-site ions Pb, M_A, and M_AE, each located at the body center, was employed; that of each of the B-site ions in the cages was averaged over each model. The table and the figure showed variations in the effective charges of M_A and M_AE for the perovskite models, similar to those along the columns of the periodic table for (M_A)₂O and (M_AE)O, respectively. For the effective charges of the B site, Nb showed a value of almost +3, regardless of the value of M_A in Table 5. Apart from PT, the values of Ti increased downward along the periodic column of M_AE. This variation of Ti is similar to that of M_AE and considered relevant to its electronegativity. The values of Zr of BaZrO₃, which lies immediately below Ti along the periodic column, were higher than those of Ti. This variation was the same as that of M_AE.

The effective charges of Pb of PT were approximately in the middle of the M_A and M_AE levels, and this Ti showed a smaller value than those of Ti and Zr in the other models. If the value of Pb (1.302) and Ti (2.369) were both magic numbers for the excellent performance of Pb-based substances, piezoelectricity may be designed or enhanced as a function of both the effective charges of the A and B sites averaged with respect to any arbitrarily modeled perovskite solid solution.

The effective charges of PT for each ratio of the isotropic deformation given for the model are listed in Table 6 and drawn as a histogram in Figure 5a. The effective charges were averaged as described above. The results show that the effective charges of Pb and Ti decreased with increasing expansion. Therefore, a compressed electronic state of the vacant 6s of Pb²⁺ and fully occupied 2p of O²⁻ benefits from their ionic nature. The reason is that these oppositely charged ions, which are closely located, decrease the electrostatic energy [21]. In contrast, a covalent nature is favorable for the expanded state. The reason is that they reduce energy through their orbital hybridization of the low energy states [21]. The whole trend is illustrated in Figure 5b and is similar to the case of 3d⁰ of Ti⁴⁺ [21]. Therefore, the effective charges can be utilized as a measure of covalency, similar to the Born effective charge [29].

An example of averaging the effective charges is given below. The effective charges of our solid solution model (Ba, Na)(Ti, Nb)O₃ are listed in Table 7 and drawn as a histogram in Figure 6. This case led us to assume two types of solid solution models, namely, Ba- and Na-body-centered. For the former model, the effective charge of Ba at the body center, bonded with the 12 nearest neighbor oxygens, was employed. For Na, each of the eight Na ions at the corners was bonded with the three nearest neighbors, and each of the six at the face centers with the eight. Thus, their charge transfers were multiplied by 4 and 12/8, respectively, to compute the effective charges at these two non-body-centered locations. These weighted transfers were each averaged over their equivalent locations of this model, and these averaged values were each subtracted from each formal charge to obtain each effective charge. Furthermore, the average effective charges of Ba and Na were further averaged in a manner we called the cluster chemical mean field. The further averaged value is also given in Table 7. In the same way, the effective charges of the latter model were obtained by reversing Ba and Na in the former model (reverse Na and Ba, respectively).

For the former and latter models, the effective charges of Ti were simply averaged over the four species, and the same was done for Nb. The averaged effective charges of Ti and Nb were further averaged, and this value is also given in Table 7.

4. Discussion and Guidelines for Designing Lead-Free Piezoelectrics

We checked the uniqueness of the above-mentioned magic numbers of PT and other Pb-based systems such as relaxor Pb(Mg_1/3Nb_2/3)O₃, called PMN [30], which is a good component for high piezoelectric performance. Since there is no adequate symmetry for modeling the PMN cluster, we believed that the calculation of PMN would consume a long time. Thus, the calculation was not conducted; instead, the charges of Pb and Nb in Table 5 and that of Mg in Table 4 were employed to obtain the charges. The obtained values of A and B were approximately 1.3 and about 2.4, respectively, and these values for PMN were close to those for PT. Therefore, the magic numbers are likely true, although they must be verified through the fabrication of Pb-free substances, at least indirectly.

We also checked the relaxor PbZn_1/3Nb_2/3O₃, called PZN [31], which is also a good component. The Zn contained in this substance did not belong to the related columns of the periodic table, so the effective charge was not calculated. However, the Zn of its oxide is expected to have a charge that is considerably less than the formal charge of +2 because Zn is a relatively highly electronegative metal. Let us assume that the effective charge of Zn in its oxide is +2, which is the maximum value. An estimated value of the effective charge of this B site, designated as B_avereged, given by averaging the formal charge of Zn with the charge of Nb of about +3 in Table 5, was about 2.7, which was evidently beyond the proposed magic number, but the inherent electronegativity of Zn causes B_avereged to be reduced in reality.

We also examined the well-known, excellent piezoelectric Pb(Zr_xTi_1−x)O₃ (x~0.5) called PZT [32]. This B_avereged averaged with the charges of Ti and Zr in Table 5 was approximately 2.6, which was also considerably higher than the proposed magic number. Thus, the good property is considered purely as the morphotropic effect of PT with PbZrO₃, which is the case with BT–BZ solid solutions [14,15,16].

A check of the uniqueness of the magic numbers did not reveal any morphotropic effect, piezoelectric coefficient, or electromechanical coupling coefficient, but it confirmed whether the effective charges of the modeled perovskites were close to the magic numbers. It also suggests the fabrication target systems that may have Pb–O-type covalency.

We set up two orthogonal axes so that they both lay on one plane. We assumed that one of the two axes indicated the averaged effective charges of either A or B ions for each model; the other axis, the other charge of B or A. In this way, both charges for each model give a coordinate in the plane. We marked, with a dot, each coordinate determined from the charges in Table 5, Table 6 and Table 7 including those of the BNT of the Na-centered model employed in a previous study [9]. The dots are given in Figure 7, which shows the following features.

The one circled dot is at the coordinate of the unstrained PT.
The three dots of M_ANbO₃ are in the region where the effective charges of A and B are relatively low and high, respectively.
The three dots of M_AETiO₃ are in the region where the charges of A and B are high and low, respectively.
The two dots of (Ba, Na)(Ti, Nb)O₃ are about in the middle of the regions of M_ANbO₃ and M_AETiO₃.
The one dot of BaZrO₃ is where the charges of A and B are both high.
The one dot of BNT is in the vicinity of PT.

Figure 7. Map of the values in Table 5, Table 6 and Table 7 including those of BNT (Bi_1/2Na_1/2TiO₃) of the Na-centered model. LN, NN, KN are LiNbO₃, NaNbO₃, KNbO₃. CT, ST, BT are CaTiO₃, SrTiO₃, and BaTiO₃.

From the above features, the solid solutions of M_ANbO₃ and M_AETiO₃ formed a narrow band from one region to the other, as shown in the figure. This band includes the locations of (Ba, Na)(Ti, Nb)O₃, but not PT, so the lead-free property improvement through this synthesis of solid solutions is chemically limited. Moreover, BaZrO₃ is not effective in the property improvement according to this idea. Nonetheless, this substance solely has a morphotropic effect with BT [14,15,16].

BNT is close to PT; BNT, which has been noteworthy for decades, is a successful case of solid solutions [10,11,12,13]. However, the bismuth (Bi) contained in BNT is volatile at high temperatures for sintering [33,34,35], which often deteriorates the quality during fabrication. Therefore, Bi may be an unwanted element from a fabrication point of view. Furthermore, the hardness of BNT is high compared with those of Pb-based systems. Therefore, an ultimate enhancement of the piezoelectricity may be naturally difficult, although BNT is a good piezoelectric, and solid solutions are promising because they can be fabricated at relatively low sintering temperatures [18], thereby suppressing vaporization.

Figure 7 shows that reducing the effective charge of the B site of (Ba, Na)(Ti, Nb)O₃ improved its properties. The reduction can be done by replacing O with halogens such as F and Cl, which are more electronegative than O. The replacement is performed to enhance the covalency and reduce the charge of B, thus effectively mitigating the limitation of the property improvement. According to this idea, the replacement with N is the opposite case because N is less electronegative than O.

The results of the replacement are listed in Table 8 and drawn as a histogram in Figure 8. The four oxygens marked in Figure 2c were replaced with Cl, F, and N. Since their formal charges are different from that of O, the simple replacement causes the total charges, as the double perovskite unit cell, to be nonzero. Therefore, the zero charge was maintained by further replacing Ba with Na or the opposite, and Ti with Nb or the reverse, depending on each model given by replacing with Cl, F, and N. We marked at each coordinate with a dot as determined from the charges in Table 8. The table and figure indicate an expected reduction in the effective charges. The dots of the coordinates of the charges in the table are given in Figure 9. The dots of PT and (Ba, Na)(Ti, Nb)O₃ in Figure 7 were also given to confirm any effect of the replacement. Figure 9 shows the two dots of (Ba, Na)Ti(O, Cl)₃ and the two dots of (Ba, Na)Ti(O, F)₃. Their locations were very close to PT, which means that the replacement with halogens is effective. Furthermore, Nb is not necessary in these substances; otherwise, the opposite effect is expected, since the charge of Nb is relatively high compared with that of Ti. The figure also shows the data of the other halogen-replaced models, namely, Na(Ti, Nb)(O, Cl)₃, and Na(Ti, Nb)(O, F)₃. The data indicate that the charges of A were excessively low because there was no Ba included, and B was also excessively high because Nb was included. Therefore, Ba is also effective for the design of piezoelectrics as well as halogens. Finally, as expected, the calculation based on the N-replaced models, namely, (Ba, Na)Nb(O, N)₃ and Ba(Ti, Nb)(O, N)₃, showed the reverse effect.

5. Conclusions

The effective ionic charges of perovskite dielectric complex compounds were investigated with molecular orbital calculation. Aside from PT, the results showed remarkable covalency for almost all cations and systematic variations in the effective charges of the A-site ions along their column in the periodic table. This variation corresponded to their electronegativity, and this was also the case with Ti and Zr, which belong to the same column. However, the charge of Pb for PT, which is an important component of good piezoelectrics, was a stand-alone value. Therefore, both charges of Pb and Ti may be magic numbers. In fact, the calculation results for the good piezoelectrics BNT and PMN were close to the proposed numbers. Hence, they may provide us with good guidelines for the design properties of lead-free piezoelectrics that will lead to excellent performance. Through the cluster chemical mean field, the calculation of (Ba, Na)(Ti, Nb)O₃ indicates values close to those of PT, but this improvement was limited by the excessively high charge of B. Thus, the O in the model was replaced with the halogen Cl or F, which was more electronegative. The obtained charges were very close to that of PT, and Ba, Na, Ti, O, Cl, or F was necessary for the simplified synthesis of the perovskite substances. This idea does not explain the morphotropic effect, but may support it in terms of obtaining good Pb-free substances from a chemical point of view. This finding may also be very useful for the efficient use of ubiquitous alkaline metals and alkaline-earth metals, especially sodium.

Author Contributions

Conceptualization, N.T.; Methodology, N.T.; Validation, N.T. and J.-i.S.; Investigation, N.T. and J.-i.S.; Writing—original draft preparation, N.T.; Writing—review and editing, N.T. and J.-i.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the research team of science and technology of Fukuoka University.

Acknowledgments

N. T. thanks H. Zhang and S. Park for extending this work to the next fabrication stage. Figure 1 and Figure 2 in this text were partially drawn with graphic software VENUS developed by R. A. Dilanian and F. Izumi.

Conflicts of Interest

The authors declare no conflict of interest.

References

Kuwata, J.; Uchino, K.; Nomura, S. Dielectric and piezoelectric properties of 0.91 Pb (Zn_1/3Nb_2/3) O₃-0.09 PbTiO₃ single crystals. Jpn. J. Appl. Phys. 1982, 21, 1298–1302. [Google Scholar] [CrossRef]
Shrout, T.R.; Chang, Z.P.; Kim, N.; Markgraf, S. Dielectric behavior of single crystals near the (1−x) Pb(Mg_1/3Nb_2/3)O₃-(x) PbTiO₃ morphotropic phase boundary. Ferroelectr. Lett. 1990, 12, 63–69. [Google Scholar] [CrossRef]
Yamashita, Y.; Harada, K. Crystal growth and electrical properties of lead scandium niobate-lead titanate binary single crystals. Jpn. J. Appl. Phys. 1997, 36, 6039–6042. [Google Scholar] [CrossRef]
Yasuda, N.; Ohwa, H.; Kume, M.; Yamashita, Y. Piezoelectric properties of a high curie temperature Pb(In_1/2Nb_1/2)O₃-PbTiO₃ binary system single crystal near a morphotropic phase boundary. Jpn. J. Appl. Phys. 2000, 39, 66–68. [Google Scholar] [CrossRef]
Sui, H.; Sun, H.; Liu, X.; Zhou, D.; Xu, R. Ferroelectric and dielectric behaviors of sol-gel derived perovskite PMN-PT/PZT heterostructures via compositional development: An interface-dependent study. J. Eur. Ceram. Soc. 2018, 38, 5382–5387. [Google Scholar] [CrossRef]
Cole, S.S.; Espenschied, H. Lead titanate: Crystal structure, temperature of formation, and specific gravity data. J. Phys. Chem. 1937, 41, 445–451. [Google Scholar] [CrossRef]
Megaw, H.D. Crystal structure of double oxides of the perovskite type. Proc. Phys. Soc. 1946, 58, 133–152. [Google Scholar] [CrossRef]
Cohen, R.E. Origin of ferroelectricity in perovskite oxides. Nature 1992, 358, 136–138. [Google Scholar] [CrossRef]
Takesue, N.; Ishibashi, K.; Asakura, K. Evaluation of covalency of ions in lead-free perovskite-type dielectric oxides. AIP Adv. 2017, 7, 105016. [Google Scholar] [CrossRef] [Green Version]
Nagata, H.; Shinya, T.; Hiruma, Y.; Takenada, T.; Sakaguchi, I.; Haneda, H. Piezoelectric properties of bismuth sodium titanate ceramics. Ceram. Trans. 2005, 167, 213–221. [Google Scholar]
Fan, P.; Zhang, Y.; Zhang, Q.; Xie, B.; Zhu, Y.; Mawat, M.A.; Ma, W.; Liu, K.; Xiao, J.; Zhang, H. Large strain with low hysteresis in Bi₄Ti₃O₁₂ modified Bi_1/2(Na_0.82K_0.18)_1/2TiO₃ lead-free piezoceramics. J. Eur. Ceram. Soc. 2018, 38, 4404–4413. [Google Scholar] [CrossRef]
Kim, M.; Ito, R.; Kim, S.; Khanal, G.P.; Fujii, I.; Suzuki, T.S.; Uchikoshi, T.; Moriyoshi, C.; Kuroiwa, Y.; Wada, S. Fabrication of lead-free piezoelectric (Bi_0.5Na_0.5)TiO₃-BaTiO₃ ceramics using electrophoretic deposition. J. Mater. Sci. 2018, 53, 2396–2404. [Google Scholar] [CrossRef]
Zhou, X.; Qi, H.; Yan, Z.; Xue, G.; Luo, H.; Zhang, D. Large energy density with excellent stability in fine-grained (Bi_0.5Na_0.5)TiO₃-based lead-free ceramics. J. Eur. Ceram. Soc. 2019, 39, 4053–4059. [Google Scholar] [CrossRef]
Saito, Y.; Takao, H.; Tani, T.; Nonoyama, T.; Takatori, K.; Homma, T.; Nagaya, T.; Nakamura, M. Lead-free piezoceramics. Nature 2004, 432, 84–87. [Google Scholar] [CrossRef] [PubMed]
Liu, W.; Ren, X. Large piezoelectric effect in Pb-free ceramics. Phys. Rev. Lett. 2009, 103, 257602. [Google Scholar] [CrossRef] [Green Version]
Zhang, H.; Zhu, Y.; Fan, P.; Marwat, M.A.; Ma, W.; Liu, K.; Liu, H.; Xie, B.; Wang, K.; Koruza, J. Temperature-insensitive electric-field-induced strain and enhanced piezoelectric properties of <001> textured (K,Na)NbO₃-based lead-free piezoceramics. Acta Mater. 2018, 156, 389–398. [Google Scholar] [CrossRef]
Megaw, H.D. Ferroelectricity in Crystals; Methuen & Co. Ltd.: London, UK, 1957. [Google Scholar]
Kujirai, T.; Takagi, Y.; Nagata, H.; Takenaka, T. Low sintering temperature and electrial properties on cu-doped (Bi_1/2Na_1/2)TiO₃-based solid solution ceramics. In Proceedings of the IEEE International Symposium on Applications of Ferroelectrics (ISAF) 2019, Lausanne, Switzerland, 14–19 July 2019. [Google Scholar] [CrossRef]
Tuluk, A.Y.; Mahon, T.R.; Van Der Zwaag, S.; Groen, P. BiFeO₃ synthesis by conventional solid-state reaction. In Proceedings of the IEEE International Symposium on Applications of Ferroelectrics (ISAF) 2019, Lausanne, Switzerland, 14–19 July 2019. [Google Scholar] [CrossRef]
Tachiyama, K.; Yasui, S.; Rao, B.N.; Dazai, T.; Usami, T.; Taniyama, T.; Katayama, T.; Hamasaki, Y.; Yu, J.; He, H.; et al. Magnetic properties of single crystal GaFeO₃. MRS Adv. 2019, 4, 61–66. [Google Scholar] [CrossRef]
Hill, N.A. Why are there so few magnetic ferroelectrics? J. Phys. Chem. 2000, 104, 6694–6709. [Google Scholar] [CrossRef]
Goldschmidt, V.M. Die gesetze der krystallochemie. Naturwissenschaften 1926, 21, 477–485. [Google Scholar] [CrossRef]
Bartel, C.J.; Sutton, C.; Goldsmith, B.R.; Ouyang, R.; Musgrave, C.B.; Ghiringhelli, L.M.; Scheffler, M. New tolerance factor to predict the stability of perovskite oxides and halides. Sci. Adv. 2019, 5, 1–9. [Google Scholar] [CrossRef] [Green Version]
Slater, J.C. Quantum, Theory of Molecules and Solids; McGraw Hill: New York, NY, USA, 1974; Volume 4. [Google Scholar]
Adachi, H.; Tsukada, M.; Satoko, C. Discrete variational Xα cluster calculations. I. Application to metal clusters. J. Phys. Soc. Jpn. 1978, 45, 875–883. [Google Scholar] [CrossRef]
Editorial Commitees of Encyclopedic Dictionary of Chemistry. Encyclopedic Dictionary of Chemistry Maruzen eBook Library; Kyoritsu Shuppan Co., Ltd.: Tokyo, Japan, 1963. [Google Scholar]
Mulliken, R.S. Electronic population analysis on LCAO-MO molecular wave functions. I. J. Chem. Phys. 1955, 23, 1833–1840. [Google Scholar] [CrossRef] [Green Version]
Csizmadia, I.G. Theory and Practice of MO Calculations on Organic Molecules; Elsevier: Amsterdam, The Netherlands, 1976. [Google Scholar]
Rabe, K.; Ahn, C.H.; Tricsone, J.M. Physics of Ferroelectrics-Modern Perspective; Springer: Berlin/Heidelberg, Germany, 2007. [Google Scholar]
Cowley, R.A.; Gvasaliya, S.N.; Lushnikov, S.; Roessli, B.; Rotaru, G.-M. Relaxing with relaxors: A review of relaxor ferroelectrics. Adv. Phys. 2011, 60, 229–327. [Google Scholar] [CrossRef]
Kuwata, J.; Uchino, K.; Nomura, S. Phase transitions in the Pb (Zn_1/3Nb_2/3)O₃-PbTiO₃ system. Ferroelectrics 1981, 37, 579–582. [Google Scholar] [CrossRef]
Rouquette, J.; Haines, J.; Bornand, V.; Pintard, M.; Papet, P.; Bousquet, C.; Konczewicz, L.; Gorelli, F.A.; Hull, S. Pressure tuning of the morphotropic phase boundary in piezoelectric lead zirconate titanate. Phys. Rev. 2004, 70, 014108. [Google Scholar] [CrossRef]
König, J.; Spreitzer, M.; Jancar, B.; Suvorov, D.; Samardzija, Z.; Popović, A. The thermal decomposition of K_0.5Bi_0.5TiO₃ ceramics. J. Eur. Ceram. Soc. 2009, 29, 1695–1701. [Google Scholar] [CrossRef]
Kainz, T.; Naderer, M.; Schutz, D.; Fruhwirth, O.; Mautner, F.; Reichmann, K. Solid state synthesis and sintering of solid solutions of BNT–xBKT. J. Eur. Ceram. Soc. 2014, 34, 3685–3697. [Google Scholar] [CrossRef]
Taghaddos, E.; Charalambous, H.; Tsakalakos, T.; Safari, A. Electromechanical properties of flash sintered BNT-based piezoelectric ceramic. J. Eur. Ceram. Soc. 2019, 39, 2882–2888. [Google Scholar] [CrossRef]

Figure 1. Simple perovskite structure expressed as ABO₃.

Figure 2. Cluster models for the calculation. (a) is an antifluorite structure (M_A)₂O with oxygens colored red. (b) is a sodium chloride structure (M_AE)O with M_AE body-centered; M_AE and O are each colored green and blue, and exchanging the colors provides the oxygen body-centered structure. (c) is a double perovskite structure. The green spheres and the blue spheres are M_A, M_AE, Pb, Na, or Ba. The gold spheres and purple spheres are Ti, Nb, or Zr. The black spheres are O, Cl, F, or N. The red spheres are oxygens. Chemical species are located so that symmetry D₂ is maintained. The colors of ions in each model are summarized in Table 1.

Figure 3. Histograms of the effective charges of the cations of (a) (M_A)₂O and (b) (M_AE)O. The values in Table 3 and Table 4 were used.

Figure 4. Histograms of the effective charges of the cations of double simple perovskites ABO₃. The values in Table 5 were used.

Figure 5. Case similar to 3d⁰ of Ti⁴⁺: (a) histogram of the effective charges of the cations of strained and unstrained PT, (b) physical and chemical image of energy gain and loss with respect to strain. (a) was drawn by using the values in Table 6. r, r_c, r_e in (b) are the distance between the two ions in general meaning, the distance in the compressed state, that in the expanded state, respectively.

Figure 6. Histograms of the effective charges of the cations of the solid solution model (Ba, Na)(Ti, Nb)O₃. A corresponds to (Ba, Na), and B (Ti, Nb). The values in Table 7 were used.

Figure 8. Histogram of the effective charges of A and B of cluster model (Ba, Na)(Ti, Nb)O₃ replaced with Cl, F, N, respectively. The values in Table 8 were used. The values of PT in Table 5 and (Ba, Na)(Ti, Nb)O₃ in Table 7 were also used for comparison.

Figure 9. Map of the values in Table 8. The dots of PT and (Ba, Na)(Ti, Nb)O₃ in Figure 7 are given for comparison.

Table 1. Colors of the ions in each perovskite model.

Model	Ion	Color
PT	Pb	green, blue
	Ti	gold, purple
	O	red, black
M_ANbO₃	M_A	green, blue
	Nb	gold, purple
	O	red, black
M_AETiO₃	M_AE	green, blue
	Ti	gold, purple
	O	red, black
BaZrO₃	Ba	green, blue
	Zr	gold, purple
	O	red, black
(Ba, Na)(Ti, Nb)O₃ Ba-centered	Ba	green
	Na	blue
	Ti	gold
	Nb	purple
	O	red, black
(Ba, Na)(Ti, Nb)O₃ Na-centered	Na	green
	Ba	blue
	Ti	gold
	Nb	purple
	O	red, black
(Ba, Na)Ti(O, Cl)₃ Ba-centered	Ba	green
	Na	blue
	Ti	gold, purple
	O	red
	Cl	black
(Ba, Na)Ti(O, Cl)₃ Na-centered	Na	green
	Ba	blue
	Ti	gold, purple
	O	red
	Cl	black
Na(Ti, Nb)(O, Cl)₃	Na	green, blue
	Ti	gold
	Nb	purple
	O	red
	Cl	black
(Ba, Na)Ti(O, F)₃ Ba-centered	Ba	green
	Na	blue
	Ti	gold, purple
	O	red
	F	black
(Ba, Na)Ti(O, F)₃ Na-centered	Na	green
	Ba	blue
	Ti	gold, purple
	O	red
	F	black
Na(Ti, Nb)(O, F)₃	Na	green, blue
	Ti	gold
	Nb	purple
	O	red
	F	black
(Ba, Na)Nb(O, N)₃ Ba-centered	Ba	green
	Na	blue
	Nb	gold, purple
	O	red
	N	black
(Ba, Na)Nb(O, N)₃ Na-centered	Na	green
	Ba	blue
	Nb	gold, purple
	O	red
	N	black
Ba(Ti, Nb)(O, N)₃	Ba	green, blue
	Ti	gold
	Nb	purple
	O	red
	N	black

Table 2. Lattice parameters used for modeling.

Model	Crystal System	Lattice Parameter (nm)
Li₂O	cubic	a = 0.461
Na₂O	cubic	a = 0.555
K₂O	cubic	a = 0.644
MgO	cubic	a = 0.421
CaO	cubic	a = 0.480
SrO	cubic	a = 0.5144
BaO	cubic	a = 0.55391
Double perovskites	a, b, c orthogonal, D₂	a = b = c = 0.4 $\times$ 2 = 0.8

Table 3. Effective charges of M_A of reverse fluoride structures.

Model	Effective Charge of M_A
Li₂O	0.601
Na₂O	0.409
K₂O	0.561

Table 4. Effective charge of M_AE of the sodium chloride structures.

Model	Effective Charge of M_AE
MgO Mg-centered	1.325
MgO O-centered	1.295
CaO Ca-centered	1.597
CaO O-centered	1.584
SrO Sr-centered	1.791
SrO O-centered	1.762
BaO Ba-centered	1.824
BaO O-centered	1.800

Table 5. Effective charges of the cations of the ABO₃ double perovskites.

Model	M_A/M_AE/Pb	Effective Charge of A	Ti/Zr/Nb	Effective Charge of B
LiNbO₃	M_A	0.703	Nb	3.020
NaNbO₃	M_A	0.608	Nb	3.007
KNbO₃	M_A	0.791	Nb	3.032
CaTiO₃	M_AE	1.581	Ti	2.451
SrTiO₃	M_AE	1.838	Ti	2.480
BaTiO₃	M_AE	1.900	Ti	2.507
BaZrO₃	M_AE	1.910	Zr	2.823
PbTiO₃	Pb	1.302	Ti	2.369

Table 6. Effective charges of PT (PbTiO₃) for each ratio of the isotropic deformation.

Model	Pb	Ti
PbTiO₃ 95%	1.381	2.584
PbTiO₃ 97.5%	1.337	2.471
PbTiO₃ 100%	1.302	2.369
PbTiO₃ 102.5%	1.275	2.274
PbTiO₃ 105%	1.254	2.188

Table 7. Effective charges of the solid solution models (Ba, Na)(Ti, Nb)O₃. Each of the detailed chemical notations of the models is given in parenthesis in each cell on the left.

Model	Ba/Na	A Averaged	Ti/Nb	B Averaged
(Ba, Na)(Ti, Nb)O₃ ((Ba_1/2Na_1/2)(Ti_1/2Nb_1/2)O₃) Ba-centered	Ba 1.900 Na 0.659	1.279	Ti 2.592 Nb 2.832	2.712
(Ba, Na)(Ti, Nb)O₃ ((Ba_1/2Na_1/2)(Ti_1/2Nb_1/2)O₃) Na-centered	Ba 1.914 Na 0.607	1.261	Ti 2.583 Nb 2.841	2.712

Table 8. Effective charges of solid solution models (Ba, Na)(Ti, Nb)O₃; O ions were partially replaced with halogen Cl or F, or not halogen N. Each of the detailed chemical notations of the models is given in parenthesis in each cell on the left.

Model	A	A Averaged	B	B Averaged
(Ba, Na)Ti(O, Cl)₃ ((Ba₄Na₄)Ti₈(O₂₀Cl₄)) Ba-centered	Ba 1.920 Na 0.500	1.210	Ti 2.344	2.344
(Ba, Na)Ti(O, Cl)₃ ((Ba₄Na₄)Ti₈O₂₀Cl₄) Na-centered	Ba 1.873 Na 0.459	1.166	Ti 2.355	2.355
Na(Ti, Nb)(O, Cl)₃ (Na₈(Ti₄Nb₄)(O₂₀Cl₄))	Na 0.445	0.445	Ti 2.434 Nb 2.611	2.523
(Ba, Na)Ti(O, F)₃ ((Ba₄Na₄)Ti₈(O₂₀F₄)) Ba-centered	Ba 1.865 Na 0.523	1.194	Ti 2.422	2.422
(Ba, Na)Ti(O, F)₃ ((Ba₄Na₄)Ti₈(O₂₀F₄)) Na-centered	Ba 1.856 Na 0.528	1.192	Ti 2.423	2.423
Na(Ti, Nb)(O, F)₃ (Na₈(Ti₄Nb₄)(O₂₀F₄))	Na 0.545	0.545	Ti 2.542 Nb 2.761	2.651
(Ba, Na)Nb(O, N)₃ ((Ba₄Na₄)Nb₈(O₂₀N₄)) Ba-centered	Ba 1.912 Na 0.724	1.318	Nb 2.983	2.983
(Ba, Na)Nb(O, N)₃ ((Ba₄Na₄)Nb₈(O₂₀N₄)) Na-centered	Ba 1.943 Na 0.609	1.276	Nb 2.983	2.983
Ba(Ti, Nb)(O, N)₃ (Ba₈(Ti₄Nb₄)(O₂₀N₄))	Ba 1.899	1.899	Ti 2.589 Nb 2.817	2.703

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Takesue, N.; Saito, J.-i. Molecular Orbital Calculation of Lead-Free Perovskite Compounds for Efficient Use of Alkaline and Alkaline Earth Metals. Crystals 2020, 10, 956. https://doi.org/10.3390/cryst10110956

AMA Style

Takesue N, Saito J-i. Molecular Orbital Calculation of Lead-Free Perovskite Compounds for Efficient Use of Alkaline and Alkaline Earth Metals. Crystals. 2020; 10(11):956. https://doi.org/10.3390/cryst10110956

Chicago/Turabian Style

Takesue, Naohisa, and Jun-ichi Saito. 2020. "Molecular Orbital Calculation of Lead-Free Perovskite Compounds for Efficient Use of Alkaline and Alkaline Earth Metals" Crystals 10, no. 11: 956. https://doi.org/10.3390/cryst10110956

APA Style

Takesue, N., & Saito, J. -i. (2020). Molecular Orbital Calculation of Lead-Free Perovskite Compounds for Efficient Use of Alkaline and Alkaline Earth Metals. Crystals, 10(11), 956. https://doi.org/10.3390/cryst10110956

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Molecular Orbital Calculation of Lead-Free Perovskite Compounds for Efficient Use of Alkaline and Alkaline Earth Metals

Abstract

1. Introduction

2. Calculation

3. Results

4. Discussion and Guidelines for Designing Lead-Free Piezoelectrics

5. Conclusions

Author Contributions

Funding

Acknowledgments

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI