Research Article  Open Access
Falko Schmidt, Marc Landmann, Eva Rauls, Nicola Argiolas, Simone Sanna, W. G. Schmidt, Arno Schindlmayr, "Consistent Atomic Geometries and Electronic Structure of Five Phases of Potassium Niobate from DensityFunctional Theory", Advances in Materials Science and Engineering, vol. 2017, Article ID 3981317, 13 pages, 2017. https://doi.org/10.1155/2017/3981317
Consistent Atomic Geometries and Electronic Structure of Five Phases of Potassium Niobate from DensityFunctional Theory
Abstract
We perform a comprehensive theoretical study of the structural and electronic properties of potassium niobate () in the cubic, tetragonal, orthorhombic, monoclinic, and rhombohedral phase, based on densityfunctional theory. The influence of different parametrizations of the exchangecorrelation functional on the investigated properties is analyzed in detail, and the results are compared to available experimental data. We argue that the PBEsol and AM05 generalized gradient approximations as well as the RTPSS metageneralized gradient approximation yield consistently accurate structural data for both the external and internal degrees of freedom and are overall superior to the localdensity approximation or other conventional generalized gradient approximations for the structural characterization of . Bandstructure calculations using a HSEtype hybrid functional further indicate significant near degeneracies of bandedge states in all phases which are expected to be relevant for the optical response of the material.
1. Introduction
The study of perovskites, oxidic compounds with the chemical formula and an ideal or distorted cubic crystal structure, is of both academic and technological interest. On one hand, perovskite minerals are the most abundant components in the Earth’s mantle [1]. On the other hand, they exhibit a large variety of interesting physical properties [2, 3], which can be modified by compositioninduced structural modifications.
Among the humanmade perovskites, potassium niobate () belongs to the important class of perovskitestructure ferroelectrics [4]. It crystallizes at about 1325 K in a paraelectric cubic phase (space group ) and, similar to [5], undergoes three successive ferroelectric phase transitions [6]: to a tetragonal phase (space group ) at 691 K, from the tetragonal to an orthorhombic phase (space group ) at 498 K, and from the orthorhombic to a rhombohedral phase (space group ) at 263 K. In addition, a new monoclinic phase (space group ) was very recently discovered in nanomaterials synthesized via a hydrothermal method [7]. It has long been known that the roomtemperature orthorhombic phase, in particular, is characterized by unusual dielectric properties and an exceptionally large nonlinear optical coefficient. For this reason, is widely used in nonlinear optics and holographic applications, where it represents a viable alternative to the more diffused .
Due to its technological importance, pure and defective has been subject to many experimental [8–12] and theoretical [13–19] studies. The previous theoretical investigations were typically aimed at the description of the groundstate structure and the identification of defect centers relevant for the optical properties. For the most part, they used densityfunctional theory [20, 21] with classical local or semilocal exchangecorrelation functionals or even simpler semiempirical models. While these approaches are appropriate for atomic geometries and the mechanical properties of a material if not too high accuracy is required, they are much less accurate with regard to the electronic and optical properties, because important factors like quasiparticle renormalization of the electronic band gap, thermal effects due to the electronphonon coupling, or excitonic effects in the linear and nonlinear optical absorption are missing. Furthermore, most of the published work so far concentrated on the structurally simple cubic and tetragonal phases, while less attention has been given to the remaining polymorphs. As a consequence, a thorough quantitative ab initio study of the electronic and optical properties of , especially of the more complex but technologically relevant roomtemperature phase, with stateoftheart methods that go beyond densityfunctional theory and take the abovementioned effects into account, as recently carried out for its sister material [22–25], is still missing. The accurate theoretical description of the atomic structure of the relevant phases remains a prerequisite for all such studies, however.
In order to pave the way for future investigations along these lines, here we perform a comprehensive study of the structural and electronic properties of all the five known phases of , including the newly discovered monoclinic phase, based on densityfunctional theory. The influence of different parametrizations of the exchangecorrelation functional on the investigated properties is discussed in detail, and the results are compared to available experimental data. Such a systematic approach is especially important for because the energetic separation and the structural differences between the various phases are small and therefore potentially sensitive to details of the computational procedure. Indeed, our results show that several recently introduced exchangecorrelation functionals give, overall, a better description of the atomic geometries throughout the various phases than common standard functionals as used in earlier studies. Bandstructure calculations are additionally performed with a nonlocal hybrid functional to correct the wellknown underestimation of the fundamental band gap incurred with conventional local or semilocal functionals, which stems from the absence of a derivative discontinuity with respect to the particle number as in the exact functional [26, 27].
This paper is organized as follows. In Section 2, we give an overview of our computational method. In Section 3, we present our results for the optimized atomic geometries of the five considered phases of and compare the performances of different exchangecorrelation functionals in relation to experimental data. Furthermore, we discuss the effects of the differently broken symmetries on the electronic band structure through the successive phase transitions. We finally summarize our conclusions in Section 4.
2. Computational Method
All calculations of this study are performed with the Vienna Ab initio Simulation Package (VASP) [28, 29], a planewave implementation of densityfunctional theory (DFT). The electronion interaction is described by the projectoraugmentedwave (PAW) scheme [29, 30], where the K and the Nb electrons are treated explicitly as valence states. For the exchange and correlation between the electrons, we test several different schemes in order to compare their performance, namely, the localdensity approximation (LDA), semilocal generalized gradient approximations (GGA), a semilocal metaGGA, and a nonlocal hybrid functional. In detail, we use the LDA as parametrized by Perdew and Zunger [31]. The conventional PerdewBurkeErnzerhof (PBE) parametrization [32] and the more recently introduced PBEsol [33] and AM05 [34] functionals are considered on the GGA side; the latter two are chosen because they are known to yield accurate lattice parameters comparable to hybridfunctional DFT calculations at a moderate computational cost [35]. Furthermore, the revised TaoPerdewStaroverovScuseria (RTPSS) metaGGA [36] is employed. Bandstructure calculations are additionally performed for the relaxed PBEsol geometries with the HeydScuseriaErnzerhof (HSE) screenedCoulombpotential hybrid functional [37, 38]. The fraction of the exact exchange from HartreeFock theory, 25% in the conventional HSE06 functional, is increased to 30% to match the experimentally observed band gap of the cubic phase. Nonlocal HSEtype hybrid functionals offer a reliable way to approximate more rigorous bandstructure calculation schemes within the framework of quasiparticle theory and thus effectively overcome the band gap problem of conventional DFT calculations with local or semilocal exchangecorrelation functionals at a lower numerical cost [39, 40].
The electronic wave functions are expanded into plane waves up to a kinetic energy of 600 eV for the latticeparameter optimization. All structural degrees of freedom are relaxed until the forces on each atom are below 0.001 meV/Å. The cutoff energy is reduced to 400 eV in the hybridfunctional bandstructure calculations. The Brillouinzone integrations are performed using a shifted MonkhorstPack 6 × 6 × 6 point mesh for all phases. In calculations with the HSE hybrid functional, a regular centered 6 × 6 × 6 point mesh is used instead.
3. Results and Discussion
Figure 1 illustrates the five different phases of considered here with exaggerated deformations from the cubic configuration (Figure 1(a)). Starting from the cubic phase, the tetragonal phase is reached by elongation along axis and compression along the other two axes (Figure 1(b)). Further distortions of the lattice parameters transform the tetragonal phase into an orthorhombic polymorph with distinct lattice parameters , , and (Figure 1(c), left). The crystal structure of orthorhombic , whose rectangular unit cell contains two formula units, that is, ten atoms, can alternatively be characterized by a smaller rhombic unit cell with lattice parameters and which contains only one formula unit (Figure 1(c), right); it is this nonrectangular primitive unit cell that we employ in all structure optimizations and electronicstructure calculations. The relation between the orthorhombic and the rhombic unit cells, including the ion displacements, is further illustrated in Figure 2. The orthorhombic phase transforms into the monoclinic phase by choosing mutually distinct lattice parameters (Figure 1(d)). Alternatively, the rhombohedral phase is obtained by choosing equal lattice parameters together with equal but acute angles (Figure 1(e)).
(a)
(b)
(c)
(d)
(e)
3.1. Structure Optimization
In order to study the influence of different approximations for exchange and correlation on the DFT results for , we perform complete structure optimizations of the lattice parameters and the atomic displacements for all five phases using the LDA, PBE, PBEsol, AM05, and RTPSS functionals. We start with the hightemperature cubic phase and then follow the order in which the phases would occur for decreasing temperature.
Table 1 lists the calculated lattice parameter and the optimized unitcell volume for the cubic phase obtained with the different exchangecorrelation functionals as well as experimental data for comparison. Consistent with the cubic symmetry, all lattice parameters have precisely the same value, and all lattice vectors are perpendicular to each other. There are no internal degrees of freedom and hence no atomic displacements in the cubic phase. At first glance, the PBE results appear to agree best with the experimental lattice constant, but it must be borne in mind that the measurement was performed at 708 K, whereas the calculations correspond to 0 K. The good quantitative agreement is in fact fortuitous and due to an error cancellation between the neglect of thermal expansion and the wellknown tendency of PBE (and other conventional GGAs) to overestimate lattice constants at zero temperature. The LDA predicts the smallest value for , while the results from other functionals are in between; the same applies to the equilibrium unitcell volume. In concordance with literature data on various other materials, we find that PBEsol, AM05, and RTPSS provide accurate and quite similar lattice parameters [36, 42–46].
 
[41]; measured at 708 K. 
For the tetragonal phase, the data listed in Table 2 is obtained with equal lattice parameters and and perpendicular lattice vectors. In contrast to the cubic phase, the displacements do not vanish in this case. If we keep the Nb atom fixed at its ideal position in the center of the unit cell, we observe that all other atoms shift towards positive values along axis; displacements along the other axes are forbidden due to symmetry constraints. Since the experimental lattice parameters shown for comparison refer to 543 K and include effects of thermal expansion, the ratio is also listed. The lattice parameter is best matched by PBEsol, while the PBE value for is closest to the experiment. The latter is fortuitous, however, as indicated by the large overestimation of which exceeds 3% and constitutes the largest deviation among all tested functionals. The deficiency of the PBE functional is further reflected in the ratio of 1.0498, which shows the strongest deviation from the experimental value of 1.0165, while LDA (1.0112) slightly underestimates the experimental value, PBEsol (1.0224) and RTPSS (1.0233) slightly overestimate it by a similar amount, and AM05 (1.0328) yields a somewhat larger ratio. Despite the reasonably accurate LDA latticeparameter ratio, the associated equilibrium volume of the unit cell of 62.08 shows the poorest agreement with the experimental volume of 64.91 . The equilibrium volume is excellently reproduced by the RTPSS value of 65.10 , followed by the AM05 value of 64.52 Å. As for the relative atomic displacements, all functionals except the LDA, which predicts systematically too small shifts, are in reasonable agreement with the experimental data. However, RTPSS clearly shows the best performance in the description of these internal degrees of freedom, aside from the accidentally close PBE displacement of the K atom. Overall, we thus conclude that PBEsol and RTPSS are well suited to describe the tetragonal phase of , although RTPSS provides the most balanced description of all structural degrees of freedom. This highlights the importance of incorporating all external as well as internal degrees of freedom in the evaluation of functional performance for structural properties.
 
[41]; measured at 543 K. 
Table 3 contains our results for the orthorhombic phase. We point out again that the primitive rhombic unit cell is not rectangular, as only the double unit cell exhibits the actual orthorhombic structure. The arrows in Figure 2 indicate the direction of the displacements in the plane perpendicular to axis; there are no detectable displacements parallel to . As can be seen, the atoms shift mainly along direction of the double unit cell. In the rotated basis system of the primitive single unit cell, the displacements have roughly equal components along and directions. As before, the LDA lattice parameters and volume substantially underestimate the experimental data. In addition, LDA tends strongly towards cubic symmetry, as indicated by the almost right angle of 89.91°, ratio of 0.9922 that is close to unity, and the significantly underestimated atomic displacements. The PBE parametrization shows the smallest (for ) and largest (for ) deviations from the experimental values among all considered functionals, indicating a similar imbalance in the description of lattice parameters as observed for the tetragonal phase. In contrast, the functionals PBEsol, AM05, and RTPSS provide reasonable lattice parameters with comparably small deviations from the experimental values for both and . This more balanced description of lattice parameters is underlined by improved ratios with PBEsol (0.9851) and RTPSS (0.9842) which show virtually identical but opposing minimal deviations from the experimental value (0.9846). While all functionals predict a very similar small reduction of angle , the experimental value of 89.73° is most accurately reproduced by AM05 (89.76°). AM05 also yields the equilibrium volume (64.70) closest to experiment (64.68). The three functionals PBEsol, AM05, and RTPSS not only yield results in reasonable quantitative agreement with the experimental data but also give a good account of the atomic displacements, although the best description of the internal degrees of freedom is provided by the AM05 parametrization. Overall, AM05 provides the best balanced description of all structural degrees of freedom in orthorhombic .
 
[41]; measured at 295 K. 
The results for the newly discovered monoclinic phase depend sensitively on the starting parameters for the atomic displacements, as the structure tends to relax to the orthorhombic phase for poorly chosen initial displacements. The only previous theoretical study of in the monoclinic phase was reported by Zhang et al. [48]. In our simulations, we take the displacements given in [48] as the initial configuration. In this way, we eventually obtain a properly relaxed monoclinic structure, which corresponds to a local minimum of the totalenergy landscape, indicating that the monoclinic phase is metastable not only for nanomaterials but also for the bulk. The results listed in Table 4 show that the performance of the different functionals is the same as that for the related orthorhombic phase, which occurs in the same temperature range. The universal trends of underestimation and overestimation of lattice parameters by the LDA and PBE functionals are clearly reflected in equilibrium volumes of 62.13 and 66.99 which fall below or exceed the experimental volume of 64.86 , respectively. Surprisingly, LDA performs best for angle as well as for the ratios and . A notable improvement over the LDA and PBE values is observable for the other functionals, but no clear trend emerges. PBEsol, AM05, and RTPSS provide virtually the same accuracy for the lattice parameters. The atomic displacements have again roughly equal components along and directions and no detectable component along . Compared to PBEsol and RTPSS, AM05 shows a slight tendency towards larger atomic displacements. Although there are no directly comparable experimental values, the calculated displacements are compatible with those of Zhang et al. [48] based on densityfunctional theory in combination with the otherwise rarely used WuCohen functional [49].
 
[47]; measured at 298 K. 
Finally, for the rhombohedral phase, we set the starting values for the displacements of all atoms except Nb to two percent along each lattice vector, because the structure tends to relax to the cubic phase for too small initial displacements. The resulting optimized values are listed in Table 5. The lattice parameter and angle obtained from PBEsol, AM05, and RTPSS are overall in good agreement with the experimental data. The performance of LDA and PBE follows that for the other phases and shows characteristic underestimation and overestimation of the lattice parameters, respectively. Our calculated atomic displacements are notably larger than those reported by Fontana et al. in 1984 [41], especially for the oxygen atoms. The reason for this discrepancy remains unclear. On one hand, the theoretical results are very consistent for all functionals considered here and thus appear robust. On the other hand, the lowtemperature rhombohedral phase has received relatively little attention in the literature so far, so that there are no other, more recent experimental characterizations of this structure which could be used for further comparison. In our simulations, all displacements are predominantly along the spatial diagonal of the unit cell with roughly equal components along all three lattice vectors.
 
[41]; measured at 230 K. 
3.2. Phase Stability
For the purpose of comparison, Table 6 summarizes our results for the totalenergy differences per formula unit for the various phases of , which correspond to the formation energies relative to the cubic polymorph. Although the numerical differences are generally small, we find that the rhombohedral structure is energetically favored by all exchangecorrelation functionals considered in this work. As the calculations refer to 0 K, this is gratifyingly consistent with the experimental fact that the rhombohedral phase is stable in the lowtemperature limit. It is followed, in ascending order, by the monoclinic, the orthorhombic, the tetragonal, and finally the cubic structure, consistent among all functionals. The energetic ordering of the different phases thus follows the sequence in which they become stable for increasing temperature. A proper ab initio description of the phase transitions requires a numerical evaluation of the temperaturedependent free energies with the full incorporation of phonon contributions [54] which goes beyond the present study, however.

Figure 3 illustrates the correlation between the unitcell volume and the relative formation energies. The remarkably narrow spread of the LDA results and the much larger spread of the PBE results for the different phases are clearly visible. This is in accordance with our earlier observation that the LDA tends to underestimate and PBE tends to overestimate symmetrybreaking deformations relative to the cubic structure. The PBEsol and AM05 results lie between the LDA and PBE values, indicating a more balanced description of the structure and energetics. The RTPSS results span a similar energy range as PBEsol and AM05 but are shifted to slightly larger volumes. In addition, Figure 3 shows an almost linear dependence of the formation energies on the unitcell volume, suggesting a dominant volume effect on the phase stability of the polymorphs. Similar trends in the correlation between formation energies and equilibrium volume have also been observed for other oxide materials that exhibit extended structure polymorphism [55, 56].
3.3. Functional Performance
The underestimation and overestimation of lattice parameters by LDA and by conventional generalized gradient approximations like PBE, respectively, observed consistently throughout this work for all investigated polymorphs, are well known and documented for a large variety of materials [45, 57, 58]. In particular, these trends were confirmed in a DFT study of various ferroelectric perovskites by Tinte et al. [59]. Nevertheless, the reasonable overall agreement with experimental data together with the wellunderstood performance in numerical simulations makes these established functionals a sound reference for the benchmarking of other exchangecorrelation functionals. The improved prediction of structural data by PBEsol, AM05, and RTPSS observed here for is likewise in accordance with recent studies for simple metallic and semiconducting systems [36, 42–46, 58], but our results demonstrate that this extends to oxides with various stable and metastable polymorphs of complex symmetry, which exhibit manifold internal degrees of freedom. We note in this context that it is also possible to improve conventional semilocal functionals through restriction to specific classes of materials by introducing systemdependent parameters that modify the localization of the exchangecorrelation hole [59], an idea that is not strictly ab initio but presently has a strong influence on the construction and use of hybrid functionals.
The energetic ordering (Table 6) of the five polymorphs is identical for all considered functionals and, most importantly, in accordance with the observed experimental sequence. Hence, even though the LDA energy differences are close to the overall numerical accuracy of current implementations of planewave DFT for solidstate systems, this does not disqualify any of the considered exchangecorrelation functionals from applications to . This happenstance of correct phase ordering independent of specific functionals, especially among structurally complex oxides characterized by polymorphrich phase diagrams, is in fact rather fortuitous. In the case of , for example, common local (LDA) and semilocal (PBE and PBEsol) as well as hybrid functionals fail to reproduce the experimentally observed energy ordering among the three main crystallographic phases rutile, anatase, and brookite, which structurally differ solely by distortions and rearrangements of the octahedral building blocks [60–63]. Furthermore, the choice of exchangecorrelation functionals in general, as well as the fraction of exact exchange in various hybrid functionals, has been found to alter the phase stability of oxide systems [63, 64].
The body of data for the atomic structure and phase energetics suggests that there is no good justification to stick to the LDA or to conventional GGAs when dealing with complex oxide materials. More recent parametrizations like PBEsol or AM05, which specifically aim at an accurate description of solidstate systems, as well as metaGGAs like RTPSS, constructed to satisfy a maximum of exact constraints for the exchangecorrelation functional without empirical parameters, offer a reliable and numerically efficient way to improve the accuracy of stateoftheart ab initio simulations. All three functionals provide a sound basis for subsequent calculations. In the following, we largely restrict the analysis of the electronic structure to the PBEsol geometries summarized in Table 7, as we do not expect qualitatively different trends for geometries from AM05 or RTPSS. In view of future studies, these functionals seem to be well suited as a starting point not only for hybridfunctional DFT but also for higherlevel perturbative approaches within the framework of manybody perturbation theory [65]. The analysis of vibrational properties in will equally benefit from the qualitative and quantitative reliability provided by these functionals, as recently demonstrated for the related material [24].

3.4. Band Structures
We now proceed to analyzing the effects of the structural deformations in the various phases on the electronic band structures. In a preliminary step, we study the influence of the structural parameters and the exchangecorrelation functionals for the cubic phase, which serves as a test case. Table 8 shows results for the fundamental band gaps with different combinations of structural parameters (see Table 1) and functionals. In all cases, the band gap is predicted to be indirect and corresponds to a transition between the valenceband maximum at , the corner of the cubic Brillouin zone, and the conductionband minimum at , its center. Furthermore, we find that the band gap is not very sensitive to the lattice parameter and that, as expected, all (semi)local functionals like LDA, PBE, or PBEsol systematically underestimate the reported experimental values, which lie between 3.14 eV [50] and 3.24 eV [52]. According to the preceding section, we consistently use structural parameters from PBEsol (see Table 7) from now on, and we perform a series of additional bandstructure calculations with the nonlocal HSE functional. In our notation, HSE means that the HSE functional is used with percent exact exchange, which we treat as a tunable parameter. Depending on , the hybridfunctional band gap can be positioned between the lower limit of (semi)local DFT and the upper limit of HartreeFock theory. In Figure 4, the dependence of the fundamental band gap of cubic on the fraction of exact exchange is illustrated for . In the most relevant region , close to the experimentally reported band gap, an almost linear dependence on is observed. The LDA, PBE, and PBEsol band gaps are shown for comparison and exhibit an underestimation by a factor of one half. According to Table 8, a choice of 30% exact exchange leads to much better agreement with the experimental band gap than the standard setting of 25%, so that we use the HSE30 variant in the following. The improved description of band gaps and the electronic structure in general with an increased fraction of exact exchange in HSEtype hybrid functionals has previously been demonstrated for oxide semiconductors [61] and nitride semiconductor alloys [66, 67], for instance. Wang et al. [68] also noted that an increased fraction of exact exchange improves the band gap of . Although adjusting the fraction of exact exchange for a single polymorph might not yield perfect agreement between experiment and theory for other phases, it represents a reasonable approach to study physical trends among the various polymorphs.

Figure 5 shows the calculated band structures for the five different phases of . Furthermore, in Table 9, we compare the fundamental band gaps with experimental values. The calculated HSE30 band gaps show a broad spread of 0.66 eV between the fundamental band gaps of cubic (3.14 eV) and rhombohedral (3.80 eV) . It should be noted that our results refer to electronic band gaps, which can be measured, in principle, by photoemission spectroscopy. No such measurements have been reported for so far, however. Instead, the experimental values listed in the table are deduced from optical spectroscopies and thus refer to optical band gaps. The difference between the two equals the exciton binding energy, whose magnitude is unknown for , but detailed investigations of the related material [25] suggest that it may be significant. In addition, the experimental values include thermal effects due to thermal expansion and electronphonon coupling, particularly in the hightemperature phases, as well as size and geometry effects, because the experiments were typically not carried out for bulk samples but for nanostructures, such as microcubes in [52] or nanowires in [48]. Nanoparticles of , for instance, exhibit much smaller band gaps than in the bulk [69]. For all these reasons, the comparison in Table 9 can only be indicative. No experimental band gap is available for the rhombohedral phase, but our results are consistent with earlier calculations by Wang and Rappe [70] which yielded a fundamental band gap of 3.33 eV with the HSE25 functional and LDA structural data; our own HSE25 calculation with structural data from PBEsol for the rhombohedral phase yields 3.38 eV.
 
[50]; optical gap from absorption spectrum. [51]; optical gap from UVvis diffuse reflectance spectrum. [52]; optical gap from absorption spectrum. [53]; optical gap from Xray photoelectron spectrum. [47]; optical gap from absorption spectrum. [48]; optical gap from diffuse reflectance UVvis spectroscopy. 
The graphical band structures in Figure 5 as well as the numerical data in Table 10, which lists energy values of the highest valence band and the lowest conduction band at selected highsymmetry points, confirm that in all phases there are states elsewhere in the Brillouin zone which are nearly degenerate with the valenceband maximum or the conductionband minimum. For example, in the orthorhombic phase that is stable at room temperature, the valenceband edge at lies only 0.04 eV below the valenceband maximum, which is located at . It is expected that such near degeneracies strongly influence the optical properties close to the absorption edge due to strong mixing of electronic transitions at different points in the optical response. The corresponding densities of states (DOS) in Figure 5 share many similar features, independent of the particular polymorph. The DOS contributions from the K atoms are negligible around the band edges in all phases. The topmost valence bands are consistently formed by oxygen states and give rise to two prominent peaks close to the valenceband edge: the first peak arises from the valenceband maximum and the nearly degenerate electronic states at other highsymmetry points, while the second peak originates from an almost degenerate continuum of states around ; its energetic distance from the conductionband edge equals the direct band gap at for all phases. The lowest conduction bands, in contrast, are mainly formed by niobium states.

4. Conclusions
In summary, we have performed a comprehensive analysis of the structural and electronic properties of in its five known phases, including the only very recently discovered monoclinic phase, based on densityfunctional theory. As the energetic separation and the structural differences between the phases are small and thus potentially sensitive to details of the computational procedure, we have placed special emphasis on the influence of different approximations for the exchangecorrelation functional. To this end, we compared the LDA with a range of GGAs, a metaGGA, and, in the case of band structures, a HSEtype hybrid functional. As an important result, we find that all functionals predict the same energetic ordering of the phases, which, moreover, coincides with the order in which the phases become stable for increasing temperature. Our results further confirm that the monoclinic configuration, which has so far only been observed in synthesized nanostructures, is also metastable in the bulk. Among the considered functionals, PBEsol as well as AM05 and RTPSS yield structural data in consistently good quantitative agreement with experimental measurements for all phases; in contrast, LDA and PBE systematically underestimate and overestimate the lattice parameters, respectively, and are also less accurate with respect to the atomic displacements inside the unit cell. For this reason, we evenhandedly recommend PBEsol, AM05, or RTPSS as reasonable starting points for ab initio investigations with higherlevel perturbative methods that require accurate equilibrium geometries. In particular, phonons, which not only are crucial in the context of phase transitions but also play a significant role for the temperature dependence of the electronic and optical properties due to the electronphonon coupling, are sensitive to small changes in the atomic spacing.
Our bandstructure calculations performed with the HSE hybrid functional and 30% exact exchange confirm that the fundamental band gap is indirect for all phases. The variation in the size of the gap amounts to 0.66 eV through the successive phase transitions. In particular, in the roomtemperature orthorhombic phase, the band gap is appreciably larger than in the tetragonal phase, which is often used as a structurally simpler substitute in theoretical simulations. For all phases, our results further indicate a range of near degeneracies of bandedge states in different parts of the Brillouin zone which are expected to be relevant for the linear and nonlinear optical properties of . More elaborate ab initio bandstructure calculations with a proper treatment of quasiparticle and thermal effects, which go beyond the present study, are highly desirable, however.
Competing Interests
The authors declare that there are no competing interests regarding the publication of this paper.
Acknowledgments
The authors gratefully acknowledge financial support from the Deutsche Forschungsgemeinschaft (DFG) via Sonderforschungsbereich TRR 142.
References
 C.M. Fang and R. Ahuja, “Structures and stability of ABO_{3} orthorhombic perovskites at the Earth's mantle conditions from firstprinciples theory,” Physics of the Earth and Planetary Interiors, vol. 157, no. 12, pp. 1–7, 2006. View at: Publisher Site  Google Scholar
 R. H. Mitchell, Perovskites: Modern and Ancient, Almaz Press, Thunder Bay, Canada, 2002.
 R. J. D. Tilley, Perovskites: StructureProperty Relationships, John Wiley & Sons, Chichester, UK, 2016. View at: Publisher Site
 B. T. Matthias, “New ferroelectric crystals,” Physical Review, vol. 75, no. 11, p. 1771, 1949. View at: Publisher Site  Google Scholar
 S. Sanna, C. Thierfelder, S. Wippermann, T. P. Sinha, and W. G. Schmidt, “Barium titanate ground and excitedstate properties from firstprinciples calculations,” Physical Review B, vol. 83, no. 5, Article ID 054112, 2011. View at: Publisher Site  Google Scholar
 A. W. Hewat, “Cubictetragonalorthorhombicrhombohedral ferroelectric transitions in perovskite potassium niobate: neutron powder profile refinement of the structures,” Journal of Physics C: Solid State Physics, vol. 6, no. 16, pp. 2559–2572, 1973. View at: Publisher Site  Google Scholar
 S. Kim, J.H. Lee, J. Lee et al., “Synthesis of monoclinic potassium niobate nanowires that are stable at room temperature,” Journal of the American Chemical Society, vol. 135, no. 1, pp. 6–9, 2013. View at: Publisher Site  Google Scholar
 S. Simsek, H. Koc, V. A. Trepakov, A. M. Mamedov, and E. Ozbay, “Electron spectroscopy and the electronic structure of KNbO_{3}: first principle calculations,” Ferroelectrics, vol. 461, no. 1, pp. 99–105, 2014. View at: Publisher Site  Google Scholar
 L. Grigorjeva, D. K. Millers, V. Pankratov et al., “Experimental and theoretical studies of polaron optical properties in KNbO_{3} perovskite,” Solid State Communications, vol. 129, no. 11, pp. 691–696, 2004. View at: Publisher Site  Google Scholar
 E. A. Kotomin, R. I. Eglitis, G. Borstel, L. Grigorjeva, D. Millers, and V. Pankratov, “Theoretical and experimental study of primary radiation defects in KNbO_{3} perovskite crystals,” Nuclear Instruments and Methods in Physics Research B, vol. 166167, pp. 299–304, 2000. View at: Publisher Site  Google Scholar
 E. E. Krasovskii, O. V. Krasovska, and W. Schattke, “Ab initio calculation of the optical and photoelectron spectra of KNbO_{3} and KTaO_{3},” Journal of Electron Spectroscopy and Related Phenomena, vol. 83, no. 23, pp. 121–127, 1997. View at: Publisher Site  Google Scholar
 S. Torbrügge, M. Imlau, B. Schoke et al., “Optically generated small electron and hole polarons in nominally undoped and Fedoped KNbO_{3} investigated by transient absorption spectroscopy,” Physical Review B, vol. 78, no. 12, Article ID 125112, 2008. View at: Publisher Site  Google Scholar
 H. Donnerberg and M. Exner, “Derivation and application of ab initio Nb^{5+}–O^{2−} shortrange effective pair potentials in shellmodel simulations of KNbO_{3} and KTaO_{3},” Physical Review B, vol. 49, no. 6, pp. 3746–3754, 1994. View at: Google Scholar
 T. Neumann, G. Borstel, C. Scharfschwerdt, and M. Neumann, “Electronic structure of KNbO_{3} and KTaO_{3},” Physical Review B, vol. 46, no. 17, pp. 10623–10628, 1992. View at: Publisher Site  Google Scholar
 S. Cabuk, “Electronic structure and optical properties of KNbO_{3}: first principles study,” Optoelectronics and Advanced Materials—Rapid Communications, vol. 1, no. 3, pp. 100–107, 2007. View at: Google Scholar
 A. V. Postnikov, T. Neumann, G. Borstel, and M. Methfessel, “Ferroelectric structure of KNbO_{3} and KTaO_{3} from firstprinciples calculations,” Physical Review B, vol. 48, no. 9, pp. 5910–5918, 1993. View at: Publisher Site  Google Scholar
 R. I. Eglitis, A. V. Postnikov, and G. Borstel, “Semiempirical HartreeFock calculations for KNbO_{3},” Physical Review B, vol. 54, no. 4, pp. 2421–2427, 1996. View at: Publisher Site  Google Scholar
 E. A. Kotomin, R. I. Eglitis, A. V. Postnikov, G. Borstel, and N. E. Christensen, “Firstprinciples and semiempirical calculations for boundhole polarons in KNbO_{3},” Physical Review B, vol. 60, no. 1, pp. 1–5, 1999. View at: Publisher Site  Google Scholar
 N. E. Christensen, E. A. Kotomin, R. I. Eglitis et al., “Quantum mechanical modelling of pure and defective KNbO_{3} perovskites,” in Defects and SurfaceInduced Effects in Advanced Perovskites, G. Borstel, A. Krumins, and D. K. Millers, Eds., pp. 3–16, Kluwer, Dordrecht, The Netherlands, 2000. View at: Google Scholar
 P. Hohenberg and W. Kohn, “Inhomogeneous electron gas,” Physical Review, vol. 136, no. 3B, pp. B864–B871, 1964. View at: Google Scholar  MathSciNet
 W. Kohn and L. J. Sham, “Selfconsistent equations including exchange and correlation effects,” Physical Review, vol. 140, no. 4A, pp. A1133–A1138, 1965. View at: Google Scholar
 C. Thierfelder, S. Sanna, A. Schindlmayr, and W. G. Schmidt, “Do we know the band gap of lithium niobate?” Physica Status Solidi (C) Current Topics in Solid State Physics, vol. 7, no. 2, pp. 362–365, 2010. View at: Publisher Site  Google Scholar
 A. Riefer, S. Sanna, A. Schindlmayr, and W. G. Schmidt, “Optical response of stoichiometric and congruent lithium niobate from firstprinciples calculations,” Physical Review B, vol. 87, no. 19, Article ID 195208, 2013. View at: Publisher Site  Google Scholar
 M. Friedrich, A. Riefer, S. Sanna, W. G. Schmidt, and A. Schindlmayr, “Phonon dispersion and zeropoint renormalization of LiNbO_{3} from densityfunctional perturbation theory,” Journal of Physics Condensed Matter, vol. 27, no. 38, Article ID 385402, 2015. View at: Publisher Site  Google Scholar
 A. Riefer, M. Friedrich, S. Sanna, U. Gerstmann, A. Schindlmayr, and W. G. Schmidt, “LiNbO_{3} electronic structure: manybody interactions, spinorbit coupling, and thermal effects,” Physical Review B, vol. 93, no. 7, Article ID 075205, 2016. View at: Publisher Site  Google Scholar
 L. J. Sham and M. Schlüter, “Densityfunctional theory of the energy gap,” Physical Review Letters, vol. 51, no. 20, pp. 1888–1891, 1983. View at: Publisher Site  Google Scholar
 J. P. Perdew, “Density functional theory and the band gap problem,” International Journal of Quantum Chemistry, vol. 28, no. 19S, pp. 497–523, 1985. View at: Publisher Site  Google Scholar
 G. Kresse and J. Furthmüller, “Efficiency of abinitio total energy calculations for metals and semiconductors using a planewave basis set,” Computational Materials Science, vol. 6, no. 1, pp. 15–50, 1996. View at: Publisher Site  Google Scholar
 G. Kresse and D. Joubert, “From ultrasoft pseudopotentials to the projector augmentedwave method,” Physical Review B, vol. 59, no. 3, pp. 1758–1775, 1999. View at: Publisher Site  Google Scholar
 P. E. Blöchl, “Projector augmentedwave method,” Physical Review B, vol. 50, no. 24, pp. 17953–17979, 1994. View at: Publisher Site  Google Scholar
 J. P. Perdew and A. Zunger, “Selfinteraction correction to densityfunctional approximations for manyelectron systems,” Physical Review B, vol. 23, no. 10, pp. 5048–5079, 1981. View at: Publisher Site  Google Scholar
 J. P. Perdew, K. Burke, and Y. Wang, “Generalized gradient approximation for the exchangecorrelation hole of a manyelectron system,” Physical Review B, vol. 54, no. 23, pp. 16533–16539, 1996. View at: Publisher Site  Google Scholar
 J. P. Perdew, A. Ruzsinszky, G. I. Csonka et al., “Restoring the densitygradient expansion for exchange in solids and surfaces,” Physical Review Letters, vol. 100, no. 13, Article ID 136406, 2008. View at: Publisher Site  Google Scholar
 R. Armiento and A. E. Mattsson, “Functional designed to include surface effects in selfconsistent density functional theory,” Physical Review B, vol. 72, no. 8, Article ID 085108, 2005. View at: Publisher Site  Google Scholar
 A. E. Mattsson, R. Armiento, J. Paier, G. Kresse, J. M. Wills, and T. R. Mattsson, “The AM05 density functional applied to solids,” Journal of Chemical Physics, vol. 128, no. 8, Article ID 084714, 2008. View at: Publisher Site  Google Scholar
 J. Sun, M. Marsman, G. I. Csonka et al., “Selfconsistent metageneralized gradient approximation within the projectoraugmentedwave method,” Physical Review B, vol. 84, no. 3, Article ID 035117, 2011. View at: Publisher Site  Google Scholar
 J. Heyd, G. E. Scuseria, and M. Ernzerhof, “Hybrid functionals based on a screened Coulomb potential,” Journal of Chemical Physics, vol. 118, no. 18, pp. 8207–8215, 2003. View at: Publisher Site  Google Scholar
 J. Heyd, G. E. Scuseria, and M. Ernzerhof, “Erratum: hybrid functionals based on a screened Coulomb potential [Journal of Chemical Physics 118, 8207 (2003)],” Journal of Chemical Physics, vol. 124, no. 21, Article ID 219906, 2006. View at: Google Scholar
 J. Heyd and G. E. Scuseria, “Efficient hybrid density functional calculations in solids: assessment of the HeydScuseriaErnzerhof screened Coulomb hybrid functional,” Journal of Chemical Physics, vol. 121, no. 3, pp. 1187–1192, 2004. View at: Publisher Site  Google Scholar
 J. Heyd, J. E. Peralta, G. E. Scuseria, and R. L. Martin, “Energy band gaps and lattice parameters evaluated with the HeydScuseriaErnzerhof screened hybrid functional,” Journal of Chemical Physics, vol. 123, no. 17, Article ID 174101, 2005. View at: Publisher Site  Google Scholar
 M. D. Fontana, G. Metrat, J. L. Servoin, and F. Gervais, “Infrared spectroscopy in KNbO_{3} through the successive ferroelectric phase transitions,” Journal of Physics C: Solid State Physics, vol. 17, no. 3, pp. 483–514, 1984. View at: Publisher Site  Google Scholar
 G. I. Csonka, J. P. Perdew, A. Ruzsinszky et al., “Assessing the performance of recent density functionals for bulk solids,” Physical Review B, vol. 79, no. 15, Article ID 155107, 2009. View at: Publisher Site  Google Scholar
 J. P. Perdew, A. Ruzsinszky, G. I. Csonka, L. A. Constantin, and J. Sun, “Workhorse semilocal density functional for condensed matter physics and quantum chemistry,” Physical Review Letters, vol. 103, no. 2, Article ID 026403, 2009. View at: Google Scholar
 R. Peverati and D. G. Truhlar, “Performance of the M11L density functional for bandgaps and lattice constants of unary and binary semiconductors,” Journal of Chemical Physics, vol. 136, no. 13, Article ID 134704, 2012. View at: Publisher Site  Google Scholar
 L. He, F. Liu, G. Hautier et al., “Accuracy of generalized gradient approximation functionals for densityfunctional perturbation theory calculations,” Physical Review B, vol. 89, no. 6, Article ID 064305, 2014. View at: Publisher Site  Google Scholar
 C. Lin, Y. Zhao, and G. Yin, “Calculation of the lattice constant of solids with the use of valence electron structure parameters,” Computational Materials Science, vol. 97, pp. 86–93, 2015. View at: Publisher Site  Google Scholar
 P. G. Kang, B. K. Yun, S. Shin et al., “Possible origin of stabilized monoclinic structure of KNbO_{3} nanomaterials at room temperature,” Materials Science and Engineering B, vol. 210, pp. 19–23, 2016. View at: Publisher Site  Google Scholar
 T. Zhang, W. Lei, P. Liu et al., “Insights into the structurephotoreactivity relationships in welldefined perovskite ferroelectric KNbO_{3} nanowires,” Chemical Science, vol. 6, no. 7, pp. 4118–4123, 2015. View at: Publisher Site  Google Scholar
 Z. Wu and R. E. Cohen, “More accurate generalized gradient approximation for solids,” Physical Review B, vol. 73, no. 23, Article ID 235116, 2006. View at: Publisher Site  Google Scholar
 J. W. Liu, G. Chen, Z. H. Li, and Z. G. Zhang, “Hydrothermal synthesis and photocatalytic properties of ATaO_{3} and ANbO_{3} (A = Na and K),” International Journal of Hydrogen Energy, vol. 32, no. 13, pp. 2269–2272, 2007. View at: Google Scholar
 T. Su, H. Jiang, H. Gong, and Y. Zhai, “An alternative approach of solidstate reaction to prepare nanocrystalline KNbO_{3},” Journal of Materials Science, vol. 45, no. 14, pp. 3778–3783, 2010. View at: Publisher Site  Google Scholar
 T. Zhang, K. Zhao, J. Yu et al., “Photocatalytic water splitting for hydrogen generation on cubic, orthorhombic, and tetragonal KNbO_{3} microcubes,” Nanoscale, vol. 5, no. 18, pp. 8375–8383, 2013. View at: Publisher Site  Google Scholar
 Y. Shiozaki, E. Nakamura, and T. Mitsui, Ferroelectrics and Related Substances: Oxides Part 1: PerovskiteType Oxides and LiNbO_{3} Family, New Series III/36A1, LandoltBörnstein, 2001.
 M. Friedrich, A. Schindlmayr, W. G. Schmidt, and S. Sanna, “LiTaO_{3} phonon dispersion and ferroelectric transition calculated from first principles,” Physica Status Solidi B, vol. 253, no. 4, pp. 683–689, 2016. View at: Publisher Site  Google Scholar
 J. K. Dewhurst and J. E. Lowther, “Highpressure structural phases of titanium dioxide,” Physical Review B, vol. 54, no. 6, pp. R3673–R3675, 1996. View at: Publisher Site  Google Scholar
 A. Beltrán, L. Gracia, and J. Andrés, “Density functional theory study of the brookite surfaces and phase transitions between natural titania polymorphs,” The Journal of Physical Chemistry B, vol. 110, no. 46, pp. 23417–23423, 2006. View at: Publisher Site  Google Scholar
 A. van de Walle and G. Ceder, “Correcting overbinding in localdensityapproximation calculations,” Physical Review B, vol. 59, no. 23, pp. 14992–15001, 1999. View at: Publisher Site  Google Scholar
 P. Haas, F. Tran, and P. Blaha, “Calculation of the lattice constant of solids with semilocal functionals,” Physical Review B, vol. 79, no. 8, Article ID 085104, 2009. View at: Publisher Site  Google Scholar
 S. Tinte, M. G. Stachiotti, C. O. Rodriguez, D. L. Novikov, and N. E. Christensen, “Applications of the generalized gradient approximation to ferroelectric perovskites,” Physical Review B, vol. 58, no. 18, pp. 11959–11963, 1998. View at: Publisher Site  Google Scholar
 J. Muscat, V. Swamy, and N. M. Harrison, “Firstprinciples calculations of the phase stability of TiO_{2},” Physical Review B, vol. 65, no. 22, pp. 4112–4126, 2002. View at: Google Scholar
 M. Landmann, E. Rauls, and W. G. Schmidt, “The electronic structure and optical response of rutile, anatase and brookite TiO_{2},” Journal of Physics Condensed Matter, vol. 24, no. 19, Article ID 195503, 2012. View at: Publisher Site  Google Scholar
 F. Labat, P. Baranek, C. Domain, C. Minot, and C. Adamo, “Density functional theory analysis of the structural and electronic properties of TiO_{2} rutile and anatase polytypes: performances of different exchangecorrelation functionals,” Journal of Chemical Physics, vol. 126, no. 15, Article ID 154703, 2007. View at: Publisher Site  Google Scholar
 M. Gerosa, C. E. Bottani, L. Caramella, G. Onida, C. Di Valentin, and G. Pacchioni, “Electronic structure and phase stability of oxide semiconductors: performance of dielectricdependent hybrid functional DFT, benchmarked against GW band structure calculations and experiments,” Physical Review B, vol. 91, no. 15, Article ID 155201, 2015. View at: Publisher Site  Google Scholar
 M. T. Curnan and J. R. Kitchin, “Investigating the energetic ordering of stable and metastable TiO_{2} polymorphs using DFT+U and hybrid functionals,” Journal of Physical Chemistry C, vol. 119, no. 36, pp. 21060–21071, 2015. View at: Publisher Site  Google Scholar
 A. Schindlmayr, “The GW approximation for the electronic selfenergy,” in ManyElectron Approaches in Physics, Chemistry and Mathematics, V. Bach and L. Delle Site, Eds., vol. 29 of Mathematical Physics Studies, pp. 343–357, Springer International, Cham, Switzerland, 2014. View at: Publisher Site  Google Scholar
 M. Landmann, M. Feneberg, M. Röppischer et al., “Transition energies and directindirect band gap crossing in zincblende Al_{x}Ga_{1−x}N,” Physical Review B, vol. 87, no. 19, Article ID 195210, 2013. View at: Google Scholar
 C. Mietze, M. Landmann, E. Rauls et al., “Band offsets in cubic GaN/AlN superlattices,” Physical Review B, vol. 83, no. 19, Article ID 195301, 2011. View at: Publisher Site  Google Scholar
 G. Wang, H. Chen, Y. Li, A. Kuang, H. Yuan, and G. Wu, “A hybrid density functional study on the visible light photocatalytic activity of (Mo,Cr)–N codoped KNbO_{3},” Physical Chemistry Chemical Physics, vol. 17, no. 43, pp. 28743–28753, 2015. View at: Publisher Site  Google Scholar
 M. M. Khan, S. A. Ansari, D. Pradhan et al., “Band gap engineered TiO_{2} nanoparticles for visible light induced photoelectrochemical and photocatalytic studies,” Journal of Materials Chemistry A, vol. 2, no. 3, pp. 637–644, 2014. View at: Publisher Site  Google Scholar
 F. Wang and A. M. Rappe, “Firstprinciples calculation of the bulk photovoltaic effect in KNbO_{3} and (K,Ba)(Ni,Nb)O_{3−δ},” Physical Review B, vol. 91, no. 16, Article ID 165124, 2015. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Falko Schmidt et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.