---
abstract: 'We investigate the effect of $A$-site size differences in the double perovskites BiScO$_3$-$M$NbO$_3$ ($M$$=$Na, K and Rb) using first-principles calculations. We find that the polarization of these materials is 70$\sim$90 $\mu$C/cm$^2$ along the rhombohedral direction. The main contribution to the high polarization comes from large off-centerings of Bi ions, which are strongly enhanced by the suppression of octahedral tilts as the $M$ ion size increases. A high Born effective charge of Nb also contributes to the polarization and this contribution is also enhanced by increasing the $M$ ion size.'
author:
- 'Shigeyuki Takagi$^{1,2}$'
- 'Alaska Subedi$^{1,2}$'
- 'Valentino R. Cooper$^1$'
- 'David J. Singh$^1$'
bibliography:
- 'perovskite2.bib'
title: ' Effect of $A$-site size difference on polar behavior in $M$BiScNbO$_6$, ($M$$=$Na, K and Rb): Density functional calculations '
---

INTRODUCTION
============

Ferroelectric perovskite oxides have attracted attention for many years because of their many applications in technology, including memory devices, electronic components, piezoelectric sensors and actuators, and pyroelectric devices. [@scott] Besides the technological importance, there is also a great scientific interest in understanding the underlying physics. This is stimulated by recent discoveries of entirely new ferroelectrics particularly in thin film form. [@ahn; @lee; @dawber] The key to these materials is competition between different interactions, leading to structural instability. Ferroelectricity in perovskite $AB\text{O}_3$ oxides depends on the delicate balance of the long-range Coulomb energy, which favors a ferroelectric distortion, and the short-range ionic repulsion, which stabilizes the cubic phase.[@cohen1990prb; @cohen1992nature] More generally, instabilities in perovskites are often discussed in terms of competition between preferred $A$-O and $B$-O bond lengths, which cannot both be accommodated in the ideal cubic structure for arbitrary $A$ and $B$ ions. This is characterized by a so-called tolerance factor, $t$, that describes the degree of non-satisfaction of ideal bond lengths in the cubic perovskite structure. The balance between the long-range Coulomb energy and short-range ionic repulsion can be modified by smaller effects related to the $A$-O hybridization, in particular the so-called lone pair physics associated with stereochemically active ions, particularly Pb$^{2+}$ and Bi$^{3+}$, as in the technologically important Pb(Zr,Ti)O$_3$ system and BiFeO$_3$. In thin films and superlattices the constraint imposed by lattice matching imposes strains that push bond lengths away from their ideal values, leading to highly enhanced ferroelectric instabilities in many cases. [@lee; @dawber; @cooper]

The piezoelectric performance in perovskite ferroelectrics is associated with polarization rotation near morphotropic phase boundaries (MPB). [@park1997jap; @noheda1999apl; @bellaiche2000prl; @fu2000nature; @guo2000prl; @noheda2001prl] That is, in materials where two different ferroelectric states with different lattice strains are close in energy, an electric field may rotate the polarization and therefore the ferroelectric state from the zero field equilibrium configuration towards the other nearby configuration, with resulting coupling to the lattice strain. Strong ferroelectricity is highly desirable for piezoelectricity, but epitaxial constraints on the lattice strain as in most thin films are not. It is of interest to examine strategies for producing structural frustration and perhaps enhance ferroelectricity in bulk materials.

In this regard, recent first principles calculations suggest that $A$-site off-centering can be enhanced by the use of ions with different ionic radii at the $A$-site.[@singh2008prl] There is experimental support for this in the (Ba,Ca)TiO$_3$ solid solution, which shows a strongly stabilized ferroelectric tetragonal phase over a wide composition range as compared to what would be expected from a classical averaging of the end-point properties.[@desheng2008prl] However, while this is shown experimentally to be effective in enhancing ferroelectricity in (Ba,Ca)TiO$_3$ it is not effective in the closely related (Ba,Ca)ZrO$_3$ system. [@bcz] We recently reported results of a first principles study of the double perovskites, BiPbZnNbO$_6$ and BiSrZnNbO$_6$, which have a large size difference on the $A$-site. [@takagi2010prb] We found high polarization values of $\sim$80 $\mu$C/cm$^2$ in these compounds, which is comparable to BiFeO$_3$ (90-100 $\mu$C/cm$^2$) [@neaton2005prb] even though these double perovskites have only half the $A$-sites occupied by the stereochemically active Bi$^{3+}$ ion. The high polarization in these materials results comes from large off-centerings of Bi ions, associated in part with the $A$-site size mismatch, cooperating with moderate off-centerings of Nb ions that have Born effective charges highly enhanced by covalency.

Here we use density functional calculations to study the double perovskites $M$BiScNbO$_6$ ($M$ $=$ Na, K and Rb) in order to further investigate the effect of $A$-site size differences on the polar behavior of perovskites. In relation to BiSrZnNbO$_6$, the compounds that we study here have lower ionic charges at the $A$-site and higher ionic charges at the $B$-site. The $M$ ions have increasing ionic radii ($r_{\textrm{Na}^+} < r_{\textrm{K}+} <
  r_{\textrm{Rb}^+}$), but are otherwise chemically similar, which is helpful in understanding how increasing $A$-site size difference affects the ferroelectricity. Finally, we note that these materials could potentially be of practical use. It is known that the end-point compound BiScO$_3$ (BS) strongly favors a rhombohedral state over a tetragonal state, [@eitel; @iniguez; @halilov2004prb] as does BiSrZnNbO$_6$. As such these systems, when alloyed with a tetragonal ferroelectric can show MPBs and possibly form useful piezoelectrics. [@eitel]

APPROACH
========

We investigated the lattice instabilities of three perovskite oxides (Na,Bi)(Sc,Nb)O$_6$, (K,Bi)(Sc,Nb)O$_6$ and (Rb,Bi)(Sc,Nb)O$_6$ using density functional theory (DFT) within the local-density approximation (LDA). All calculations, including the structure relaxations, polarizations and Born effective charges, were done using a plane wave basis as implemented in Quantum Espresso package.[@espresso] We used ultrasoft pseudopotentials[@uspp] and converged our results with respect to the cut-off energy and Brillouin zone sampling, similar to our previous study for BiPbZnNbO$_6$ and BiSrZnNbO$_6$.[@takagi2010prb] The polarization was calculated using the Berry’s phase method.

There is a large charge difference between Sc$^{3+}$ and Nb$^{5+}$, and this is expected lead to a strong ordering tendency on the $B$-site. This is actually the case in the relaxor compound PbSc$_{1/2}$Nb$_{1/2}$O$_3$ (PSN), where Sc and Nb have rock-salt ordering.[@chu1995jap; @malibert1997jpcond] Here we consider material with this double perovskite type $B$-site ordering. In addition, while ordering at the $A$-site is less clear, there is a large charge difference at the $A$-site, which may lead to an ordering tendency at least short range, as well. There have been reports of short-range cation ordering of Na$^+$ or K$^+$ and Bi$^{3+}$ in Na$_{1/2}$Bi$_{1/2}$TiO$_3$ (NBT) and K$_{1/2}$Bi$_{1/2}$TiO$_3$ (KBT),[@park1994jamceram; @jones2002powder; @jones2002acta; @petzelt2004jphys] which support this view. Here we consider the ordered system, both because it is simpler, which facilitates the calculations and analysis, and because it helps to avoid computational artifacts, such as polarization because of the particular symmetry breaking chemical ordering in a given supercell.

There are various factors at play that may drive ferroelectricity in NaBiScNbO$_6$ (NBSN), KBiScNbO$_6$ (KBSN) and RbBiScNbO$_6$ (RBSN). These materials have Nb$^{5+}$ at the $B$-site which favors ferroelectricity, as in KNbO$_3$. The end-point perovskite BiScO$_3$ has a very small tolerance factor $t=0.85$ and a stereochemically active Bi$^{3+}$ at the $A$-site, two factors that are considered to be important in $A$-site driven materials, leading to strong lattice instabilities. Furthermore, the large size difference at $A$-site may also enhance the off-centering of the $A$-site ions through frustration of tilt modes. However, in addition to $A$-site off-centering, low tolerance factors also strongly favor rotation of the $B$O$_6$ octahedra. These rotations compete with ferroelectricity in general. In order to sort this out it is important to use supercells that have even numbers of units along the \[100\], \[010\] and \[001\] directions to accommodate the various Glazer tilt patterns.[@glazer1972acta] We used 40 atom 2$\times$2$\times$2 supercells and assumed both $A$- and $B$-site rocksalt chemical ordering. In addition to allowing various Glazer tilt patterns, these supercells are cubic and therefore non-polar as regards to chemical ordering. Thus, any off-centering is a consequence of lattice instabilities and not an artifact due to cation ordering in the supercell.

RESULTS AND DISCUSSION
======================

To our knowledge, syntheses of NaBiScNbO$_6$, KBiScNbO$_6$ and RbBiScNbO$_6$ have not been reported. We therefore used theoretical LDA lattice parameters. We note that in oxides with heavy elements the LDA typically leads to a small underestimation of the lattice parameter by $\sim$ 1-2%. If this is the case here, it will work against ferroelectricity and in favor of competing tilt instabilities. [@fornari-tilt] We assume a pseudocubic structure to evaluate the lattice parameters and relaxed all atomic coordinates with no symmetry constraints. This was done as a function of the lattice parameter to find the equilibrium volumes for these materials. The calculated LDA lattice parameters for the minimum energy were 3.98 Å, 4.03 Å and 4.07 Å for NaBiScNbO$_6$, KBiScNbO$_6$ and RbBiScNbO$_6$, respectively. These show the expected trend with increase in the $A$-site ion size.

[cc]{} ![\[structure\] (color online) Projections onto a (010) plane of the structures of the relaxed 40 atom pseudocubic (a) NaBiScNbO$_6$, (b) KBiScNbO$_6$ and (c) RbBiScNbO$_6$ supercells. The $B$-O bonds are shown. The Bi ions are denoted by light purple spheres, the Sc ions by blue, the Nb ions by light green, and the O ions by red. The Nb, K and Rb ions are shown by dark yellow spheres in each panel. The O octahedra for Sc are shown by (semi-transparent) blue boxes and those for Nb by (semi-transparent) green. The solid lines show the unit cell boundaries of each material. ](BNSN-pos.eps "fig:"){width="4.0cm"} & ![\[structure\] (color online) Projections onto a (010) plane of the structures of the relaxed 40 atom pseudocubic (a) NaBiScNbO$_6$, (b) KBiScNbO$_6$ and (c) RbBiScNbO$_6$ supercells. The $B$-O bonds are shown. The Bi ions are denoted by light purple spheres, the Sc ions by blue, the Nb ions by light green, and the O ions by red. The Nb, K and Rb ions are shown by dark yellow spheres in each panel. The O octahedra for Sc are shown by (semi-transparent) blue boxes and those for Nb by (semi-transparent) green. The solid lines show the unit cell boundaries of each material. ](BKSN-pos.eps "fig:"){width="4.0cm"}\
(a) NaBiScNbO$_6$ & (b) KBiScNbO$_6$\
\
\

Fig. \[structure\] shows the structures of the relaxed 40 atom pseudocubic NaBiScNbO$_6$, KBiScNbO$_6$ and RbBiScNbO$_6$ supercells at the equilibrium volume. As expected, the octahedral tilting is reduced with larger $A$-site ion size difference (note that the ionic radii of $A$-site ions are $r_{\text{Bi}^{3+}}=1.31$ Å, $r_{\text{Na}^{+}}=1.53$ Å, $r_{\text{K}^{+}}=1.78$ Å and $r_{\text{Rb}^{+}}=1.86$ Å, respectively). In particular, there is relatively little tilt in RbBiScNbO$_6$, which has the largest $A$-site size difference. Fig. \[structure\] also shows large displacements of the Bi ions from the ideal positions at the centers of the O cages.

![\[offcenter\_NBSN\] (color online) Cation off-centerings along the Cartesian directions and their magnitude with respect to the center of their O cages (the 12 nearest O ions for the $A$-site and the 6 nearest for the $B$-site) in the relaxed 40 atom NaBiScNbO$_6$ supercell. ](BNSN-disp.eps){width="8cm"}

![\[offcenter\_BMSN\] (color online) Average magnitude of cation off-centerings with respect to the center of their O cages in the relaxed 40 atom NaBiScNbO$_6$, KBiScNbO$_6$ and RbBiScNbO$_6$ supercells. ](BMSN-disp.eps){width="8cm"}

The cation off-centerings in NaBiScNbO$_6$ from the center of their O cages is shown in Fig. \[offcenter\_NBSN\]. All cations off-center along \[111\] direction and the magnitude of displacement is almost the same for a given cation. The Bi ions have a large off-centering, which is a consequence of the difference in size compared to the other $A$-site ion and also its stereochemical activity. In addition to the off-centering of $A$-site cations, there are displacements of $B$-site cations as well. For NaBiScNbO$_6$, the average magnitudes of the cation off-centerings are 0.62 Å for Bi$^{3+}$, 0.28 Å for Na$^{+}$, 0.18 Å for Sc$^{3+}$ and 0.24 Å for Nb$^{5+}$, respectively. This cooperation of all cations in the ferroelectric distortion is not unusual and is also seen in KNbO$_3$, BaTiO$_3$, PbTiO$_3$, BiPbZnNbO$_6$ and BiSrZnNbO$_6$.[@singh1992ferro; @ghita2005prb; @takagi2010prb] The cation off-centerings in KBiScNbO$_6$ and RbBiScNbO$_6$ are also collinear along \[111\] direction.

As mentioned, the LDA lattice parameters of these supercells may be underestimated by 1-2$\%$. To study the robustness of ferroelectric displacements, we performed relaxation calculations for larger values of lattice parameters that lie between the equilibrium lattice parameters (3.98 Å, for NaBiScNbO$_6$, 4.03 Å, for KBiScNbO$_6$ and 4.07 Å, for NaBiScNbO$_6$) and the higher value of 4.23 Å. We found ferroelectric cation off-centerings for each such calculation, confirming the robustness of ferroelectric distortions against LDA volume errors.

Fig. \[offcenter\_BMSN\] shows the average magnitudes of the cation off-centerings in NaBiScNbO$_6$, KBiScNbO$_6$ and RbBiScNbO$_6$. The displacements of Bi and Nb increase as the $A$-site size difference is increased. The reason for this increase in displacement is two-fold. First, the larger size difference decreases the tilting of $B$O$_6$ octahedra and provides more volume for displacement of Bi and Nb ions. Secondly, the lattice parameter increases as the size of $M$ ion increases. The ferroelectricity becomes stronger with increasing volume, as is usual in perovskite oxides. [@samara] Competing with this is the fact that, since polarization is dipole per unit volume, larger unit cell volume works against high polarization. It is interesting to note that in contrast to Bi and Nb ions, the off-centerings become smaller for the $M$ ions, while there is no appreciable change in the displacement of Sc ion. As discussed below, there is little covalency between the unoccupied states of the $M$ ions with the occupied O 2$p$ states in these materials. This works against off-centering. Furthermore, $r_{\textrm{Na}^+}
  < r_{\textrm{K}+} < r_{\textrm{Rb}^+}$, and as the radius of an atom increases, there is less volume available for off-centering. This explains the smaller displacement of Rb ions as compared to Na ions.

![\[dos\_NBSN\] (color online) Electronic DOS of the pseudocubic relaxed NaBiScNbO$_6$ and the Bi $p$, Sc $d$, Nb $d$ and O $p$ projections. The valence band maximum is taken as the energy zero. ](BNSN-dos.eps){width="9cm"}

We now consider the electronic structure of these materials as it pertains to ferroelectricity. Fig. \[dos\_NBSN\] shows the electronic density of states (DOS) of the pseudocubic relaxed NaBiScNbO$_6$ and the projections onto Bi $p$, Sc $d$, Nb $d$ and O $p$ pseudo-wavefunctions (the DOS of KBiScNbO$_6$ and RbBiScNbO$_6$ are almost identical and are not shown here). The O 2$p$ states provide the DOS from -5 eV to the valence band maximum and there is a gap that separates the O 2$p$ derived states from the unoccupied metal states. The calculated LDA band gap is 3.2 eV, 3.1 eV and 3.1 eV for NaBiScNbO$_6$, KBiScNbO$_6$ and RbBiScNbO$_6$, respectively. The actual band gaps are likely to be larger than these calculated ones due to LDA band gap underestimation. In any case, based on the calculated gaps, these materials are expected to be good insulators at ambient temperature. The O 2$p$ derived valence bands contain significant admixtures of the nominally unoccupied metal states, which are mainly of Bi 6$p$, Sc 3$d$ and Nb 4$d$ character. This mixing of occupied O 2$p$ states with unoccupied $d$ states of metals implies chemical bonding and is responsible for enhanced Born effective charges and ferroelectric off-centering, as in PbTiO$_3$ and BaTiO$_3$.[@cohen1992nature; @zhong1994prl] Nb 4$d$ states are dominant near the conduction band minimum and this low energy favors hybridization with O 2$p$ states, as discussed in the case of PbMg$_{1/3}$Nb$_{2/3}$O$_3$ and KNbO$_3$ in comparison with KTaO$_3$.[@singh1996prb; @suewattana2006prb] However, the Sc 3$d$ states lie further above the conduction band minimum and this results in a somewhat weaker hybridization with O 2$p$ states. We do not find significant mixing of unoccupied Na/K/Rb $s$ states with the occupied O 2$p$ states as may be expected considering the very large electronegativity difference between these elements and O. In NBSN, the influence of covalency is reflected in the enhanced Born effective charges of 4.38 for Bi, 4.73 for Sc and 6.42 for Nb as compared to the very modest enhancement of (z$^*$=1.14) for Na. Similar numbers are obtained for KBSN and RBSN (see Table \[tab:comparison\]).

In contrast to the trend in the cation off-centerings, the Born effective charges of Bi and Nb decrease, while those of $M$ and Sc ions increase when the $A$-site size difference becomes larger. These decreases in Born charges of Bi and Nb ions would mean lower values of polarization if the Born charges were independent of atom position and the off-centering were of the same magnitude in all compounds. In fact, however, the polarizations as calculated with the Berry phase method are 73 $\mu$C/cm$^2$ for NBSN, 77 $\mu$C/cm$^2$ for KBSN, and 87 $\mu$C/cm$^2$.

                         NaBiScNbO$_6$       KBiScNbO$_6$      RbBiScNbO$_6$
  -------------------- ------------------ ------------------ ------------------
  $r_{M^+}$(Na/K/Rb)         1.53 Å             1.78 Å             1.86 Å
  $a$(LDA)                   3.98 Å             4.03 Å             4.07 Å
  $\delta$(Bi)               0.62 Å             0.67 Å             0.77 Å
  $z^*$(Bi)                   4.38               4.21               3.73
  $\delta$(Na/K/Rb)          0.28 Å             0.22 Å             0.20 Å
  $z^*$(Na.K/Rb)              1.14               1.17               1.22
  $\delta$(Sc)               0.18 Å             0.16 Å             0.16 Å
  $z^*$(Sc)                   4.73               4.81               4.92
  $\delta$(Nb)               0.24 Å             0.27 Å             0.31 Å
  $z^*$(Nb)                   6.42               6.23               5.80
  $P$(111)              73 $\mu$C/cm$^2$   77 $\mu$C/cm$^2$   87 $\mu$C/cm$^2$

  : \[tab:comparison\] Comparison of properties of NaBiScNbO$_6$, KBiScNbO$_6$ and RbBiScNbO$_6$, where $r_{\text{M}^+}$ denotes the ion radii of $M$ ion, $a$ denotes the pseudocubic lattice parameter, $z^*$ denotes the average Born effective charge, $\delta$ denotes the average cation off-centering with respect to the O cages, and $P$ is the magnitude of polarization.

As mentioned, high piezoelectric performance is induced by the polarization rotation from the rhombohedral state to tetragonal state near the MPB. We studied the tendency towards tetragonality by calculating the total energies of NBSN, KBSN, and RBSN as a function of imposed tetragonal strain ($0.94<c/a<1.06$) while keeping the volumes equal to their equilibrium volume for the pseudocubic cells. In all materials, the lowest energy is at $c/a=1.00$ and the cation off-centerings are collinear and along \[111\] direction. The cations do not off center along \[001\] directions even with large $c/a$ and the off-centerings remain practically along the \[111\] direction over the range of $c/a$ values that we tried, implying robustness of the rhombohedral ferroelectricity. Therefore, these are rhombohedral ferroelectric materials, and as such it is of interest to alloy with tetragonal ferroelectrics in order to produce a potential MPB.

Related to this, a similar compound Na$_{0.5}$Bi$_{0.5}$TiO$_3$ (NBT) is also a rhombohedral ferroelectric[@buhrer62] with $T_C
  = 593.15$ K and polarization 38 $\mu$C/cm$^2$.[@takenaka91] Solid solutions of NBT with tetragonal BaTiO$_3$ (NBBT) show a useful MPB with polarization 20 $\mu$C/cm$^2$.[@takenaka91] While the material has not proved superior to other compounds such as Pb(Zr,Ti)O$_3$ in transducers, it does have a lower free permittivity and a high electromechanical coupling factor that make it an interesting lead-free actuator material. The materials in the present study are similar in that they are rhombohedral ferroelectrics and may be considered to be related to NBT with Ti replaced by Sc$_{0.5}$Nb$_{0.5}$. By analogy, solid solutions of these materials with BaTiO$_3$ may provide materials that are useful as lead-free piezoelectrics. Another possibility may be alloying with the super tetragonal material Bi(Zn,Ti)O$_3$. [@suchomel]

CONCLUSIONS
===========

We investigated the polar behavior of double perovskites NaBiScNbO$_6$, KBiScNbO$_6$ and RbBiScNbO$_6$ using density functional calculations on small supercells. We find high polarization along \[111\] directions. The main contribution to the polarization comes from large displacements of Bi and more modest displacements of Nb ions, which, however, have high Born charges. We find that the polarization increases as the $A$-site size difference becomes larger. These materials also have substantial band gaps and the rhombohedral ferroelectricity is robust. It will be of interest to investigate whether solid solutions of these materials with tetragonal ferroelectrics show MPBs possibly providing another avenue for development of lead-free piezoelectrics.

Acknowledgements {#acknowledgements .unnumbered}
================

This work was supported by the Materials Sciences and Engineering Division, Office of Basic Energy Sciences, U.S. Department of Energy (AS,VRC,DJS) and the Office of Naval Research (ST).
