1. Introduction
Ferroelectric materials are of interest for such diverse applications as non-volatile Ferroelectric Random Access Memory (FRAM),1) electro-optic modulators2) and frequency converters.3) Especially, LiNbO3 has excellent piezoand pyro-electricity properties, is photo-refractive, and displays nonlinear optical properties. It also has a high spontaneous polarization, 70 μC/cm2, and high Curie temperature, ~1480 K,4,5) so it is suitable for high temperature applications as well. Generally, congruent LiNbO3, which has an R3c structure, is easy to grow under Li deficient conditions. 6,7) However, recent successes using the vapor-transport equilibration (VTE) method8) have enabled the growth of stoichiometric LiNbO3. Although various defect models including those for pseudo-Schottky defects, Schottky defects, and Frenkel defects have been suggested9) to explain energetically preferred defects and defect clusters in congruent LiNbO3, only an Nb antisite, compensated for by four Li ion vacancies9-11) and five Nb antisites compensated for by four Nb vacancies11) were found to successfully explain the stoichiometry of congruent LiNbO3. Atomistic modeling by Donnerberg et al.12) showed that compensation by Nb vacancies is energetically less favorable than compensation by Li vacancies. Recent DFT calculations13) showed that Nb antisite compensated for by Li vacancies, and the Li Frenkel pair, are the most energetically favorable defect clusters in LiNbO3 under Li deficient (congruent) and Li rich (stoichiometric) conditions, respectively. An analysis on various Nb antisite-Li vacancy cluster models determined that several arrangements of Li vacancies around an Nb antisite are nearly energetically equivalent.14)
Interfaces between two opposite polarization domains (Domain walls) are also studied using both experiments and simulations. Gopalan et al.15) suggested two different domain walls, the Y- and X-walls, based on crystallographic orientation; the Y-wall lies parallel to the (112̄0) plane, while the X-wall lies parallel to the (101̄0) plane. In subsequent works,16) the X-wall was identified as a mixed anion-cation plane, while the Y-wall was determined to be alternating planes of cations only and anions only. Energetic studies on the two domain walls17) show that the Y-wall is energetically favored but less mobile. It was also found that non-uniaxial polarization components near the domain walls can lead to Bloch-like characteristics of the X-wall and both Bloch-like and Néel-like characteristics of the Y-walls. 18) Previous studies have demonstrated that defect/domain wall interactions can lead to switching in the preferred orientation of the domain wall from the Y-wall to the X-wall19) or to a reduction of the electric polarization.20) In extending the previous work on defects and domain walls, this study mainly focuses on the interaction between defects and domain walls in LiNbO3 and the effect of those defects on the domain dynamics. Various aspects of defects, including type and charge state, are investigated to understand the effect of domain walls on defect energetics.
2. Simulation Methods
All the calculations are performed within the density functional theory (DFT)21) level with generalized gradient approximation22) of exchange and correlation potential. Monkhorst-Pack23) k-point sampling with a grid of 1 × 4 × 2 is used. The plane wave DFT calculations are performed with the Vienna Ab Initio Simulation Package (VASP)24,25) using the projected augmented wave (PAW)26) pseudopotential. The generalized gradient approximation (GGA)22) is used to evaluate the exchange and correlation interactions. The plane waves are expanded up to 400 eV of kinetic energy.27) Ionic relaxation is performed until the maximum atomic force is below 0.01 eV/Å. The self-consistent solution of the Kohn-Sham functional is obtained using the residual minimization method direct inversion in the iterative subspace (RMM-DIIS) algorithm,28) which optimizes several individual energy bands at the same time. The pseudopotential and methodologies used here are the same as those used in our previous studies of intrinsic defects13) and domain walls.17) The climbing-image Nudged Elastic Band (NEB) method29,30) with 10 images is employed to predict the energy barrier of defect migration.
3. Results and Discussion
Although the ferroelectric fatigue and optical instability of ferroelectric LiNbO3 are known to result from defects/domain wall interactions, these drawbacks still limit the usage of ferroelectric materials. To overcome the current application limits, thorough understanding of the interaction mechanism is necessary. In the current study, the interaction energetics are determined by putting three major point defects
V L i ′ , L i i · and
N b L i · · · · at various distances from both the Y-walls and the X-wall. Three major point defects and two distinctive domain walls are identified from previous studies. From the full geometry relaxation, the interaction energetics are determined by looking at the energy difference between the defects near the wall and the defects far away from the wall. Although the size limitation constrains the maximum distance of our study at 6.46 Å for the Y-wall and at 7.75 Å for the X-wall, empirical study with a larger system proves that the defects/domain-wall interactions mainly occur within 5 Å.31) Because the interaction between the defects and the domain wall means that their energies cannot be uniquely separated from each other, we here ascribe the change in the energy of the system to a change in the defect formation energy; we could, with equal merit, ascribe this energy change to a change in the domain wall energy.
Figure 1 provides a schematic view of the Y-wall in LiNbO 3 and the spatial location of each point defect near the Y-wall. As can be seen, the Y-wall sits between two atom planes; point defects are located at different distances from the domain walls. For all cases, the defect/domain wall interactions decreased the energy, indicating that defects are preferentially positioned near the domain wall. For the Y-wall, the quadrupley charged
N b L i · · · · shows a stronger interaction (~ −0.26 eV), while the singly charged
V L i ′ or
L i i · shows weaker interaction (~ −0.12 eV). The singly charged defects are preferentially located at the nearest plane from the Y-wall, while the
N b L i · · · · prefer to be at the second nearest plane from the wall. For comparison, the energetics of the neutral defects are also considered. Interestingly, three neutral defects show similar energetics regardless of their kinds and the values are quite close to the case of the singly charged defects (
V L i x = ~ - 0.11 eV , L i i x = ~ - 0.12 eV and
N b L i x = ~ - 0.14 eV ). Thus, the charge states of defects play an important role in determining the energetics for the defect/Y-wall interactions.
A similar study has been also performed for the X-walls. Because our previous study identified two inequivalent X-walls related with the oxygen sublattice, current studies have calculated the interaction energies of intrinsic defects with both the XI-wall and the XII-wall. Although maximum interactions are observed within the second nearest plane from both the X-walls regardless of the charge states of the intrinsic defects, the energetics for their interactions are different. Fig. 2 shows a comparison of the maximum defects/X-wall interaction energy within the three nearest planes. As can be seen in Fig. 2, the majority of defects show stronger interaction with the XI-wall than with the XII-wall. The most prominent difference in the interaction energetics is observed for the quadrupley charged
N b L i · · · · . The DFT predicts that the interaction of
N b L i · · · · with the XI-wall (~ −0.55 eV) will be much stronger than that with the XII-wall (~ − 0.05 eV). For single charged defects, the maximum interaction energies of
L i i · are −0.32 eV for the XI-wall and −0.16 eV for the XII-wall; the maximum interaction energies of
V L i ′ are −0.14 eV for the XI-wall and −0.13 eV for the XII-wall. Similar to the fully charged defects, the XI-wall also shows stronger interaction with neutral defects; the value of −0.17 eV of the
N b L i x / X I -wall interaction energy is stronger than the value of −0.03 eV of the
N b L i x / X I I -wall interaction energy. Previous study identified that two inequivalent X-walls have essentially identical values of domain wall energy in a perfect system. With the existence of defects, however, the current study observed that the XI-wall could be energetically more stable than the XII-wall. Thus, the defects/X-walls interactions explain the experimental preference of the XI-wall over the XII-wall.16,19)
Previous study has shown that the domain wall in PbTiO332) acts as a trap site for oxygen vacancies. As an extension of that study, the current work shows that the domain wall does not act as a trap site for only oxygen vacancies, but for all other intrinsic defects regardless of their charge states or defect types. For the Y-wall, the quadrupely charged
N b L i · · · · will more strongly interact with the wall than the singly charged
V L i ′ and
L i i · . On the other hand,
L i i · . and
N b L i · · · · should interact more strongly with the XI-wall than with the XII-wall. Because the strong interaction of
N b L i · · · · with the domain walls increases the stability of the domain wall, the domain wall motion will be more difficult with the existence of
N b L i · · · · . Thus, among the three intrinsic defects,
N b L i · · · · will show the strongest pinning effect.
Because the domain walls act as trap sites for defects, defects might form a cluster near the domain walls. In order to verify the energetics of single defect/domain wall interactions, we have also considered defect clusters. For this study, we have considered two different defect pairs,
L i i · + V L i ′ and
N b L i · · · · + 4 V L i ′ . By assuming that defects do not interact with each other and are equally distant from the wall, we were able to determine the maximum interaction energy of the two defect clusters with one Y-wall and two X-walls; the defect cluster/domain wall interaction energies are shown in the change of DFE from the reference state and listed in Table 1. Consistent with the results of the single defect case, three different types of domain walls act as trap sites for both defect pairs. Also, the XI-wall interacts more strongly with both defect clusters.
We have looked at the energetics of the defects interacting with domain walls in LiNbO3. Now we will discuss how the defects can migrate near the domain wall. In order to study the migration of various point defects, we have employed the NEB method along the migration pathway. Because the previous study17) has shown that the Y-wall, lying parallel to the (112̄0) plane, is energetically much more favorable than the two X-walls, our study focused on the migration behavior of defects around the Y-wall. Especially, the current study has focused on the
N b L i · · · · - V L i ′ pair, because
N b L i · · · · is the most highly interacting type of defect and
V L i ′ is necessary for mitigating the highly positive charge state of
N b L i · · · · . Fig. 3 shows the potential energy profile along the migration pathway of the
N b L i · · · · - V L i ′ pair across the Y-wall. The bottom figures show the atomic configuration of the
N b L i · · · · - V L i ′ pair across the Y-wall in the initial state (IS), transition state (TS) and final state (FS). As can be seen,
N b L i · · · · and
V L i ′ sit in two oxygen octahedral cages across the Y-wall that are next to each other. Thus, it is reasonable to assume that
N b L i · · · · will migrate through the shortest pathway to
V L i ′ , as shown by the cyan arrow in the left bottom of Fig. 3. The migration barrier for
N b L i · · · · moving toward
V L i ′ the site is 3.44 eV at the transition state (TS). This shows that significant energy is required to move
N b L i · · · · across the Y-wall. Therefore, it is reasonable to think that
N b L i · · · · , which can be easily captured by the Y-wall, will have difficulty migrating across the Y-wall and thus will pin the Y-wall.
We now consider how the migration behavior of the
N b L i · · · · - V L i ′ pair changes under surplus electrons. Fig. 4 shows the potential energy profile along the migration pathway of the
N b L i · · · · - V L i ′ pair under one excess electron.
The addition of an extra electron decreases the migration barrier of the
N b L i · · · · - V L i ′ pair from 3.44 to 3.16 eV. Thus, the extra electron will enhance the migration of
N b L i ⋯ into
V L i ′ . In order to understand the decrease in the migration barrier of
N b L i · · · · into
V L i ′ , the electronic density of state (DOS) is analyzed. The inset of Fig. 4 shows the site projected DOS at the TS during the migration of
N b L i · · · · into
V L i ′ . At the TS,
N b L i · · · · sits exactly on the Y-wall. (See Fig. 3) As a result, a defect energy level is created near the conduction band minimum (CBM) of bulk LiNbO3. This energy level created by
N b L i · · · · is shown in red color in the DOS. We term this energy level a gap state because it is located within the bandgap of bulk LiNbO3. This gap state plays a major role in the change in energy barrier during the migration of
N b L i · · · · . At the IS,
N b L i · · · · does not yield any defect energy level within the bandgap. On the other hand,
N b L i · · · · , which sits in the Y-wall, will create an additional gap state at the TS. Because this gap state is only available at the TS, extra electrons will go to the bottom of the conduction band at the IS. In contrast, the peeling off of a gap state from the bottom of the conduction band reduces the energy of the system at the TS of the
N b L i · · · · - V L i ′ migration. Thus, the migration barrier of the
N b L i · · · · - V L i ′ pair is reduced by the addition of an electron. Because
N b L i · · · · can pin the motion of the Y-wall, the reduction of the migration barrier of
N b L i · · · · may enhance the movement of the Y-wall by weakening the pinning effect.
4. Conclusions
Current works have performed a systematic study of the energetics of defects/domain wall interaction in LiNbO3. Our energetic study predicts that the formation energy of intrinsic defects can be lowered near the domain wall. So, the interactions of intrinsic point defects with domain walls decrease the DFE and trap the defects in the domain wall. In the following study, we find that the migration barrier of
N b L i · · · · is significantly high near the Y-wall. Thus, the trapped defects increase the coercive energy and threshold field for the domain wall motion, in what is known as the pinning effect. In addition, trapped defects can cause a local field and this might cause optical degradation. The degradation in both ferroelectricity and the optical property is known as ferroelectric fatigue. By identifying the defects and defect pairs leading to the strongest interaction with the domain walls, we can potentially prevent ferroelectric fatigue or optical degradation by constraining the formation of specific types of defect. From our systematic study of defect/domain interactions, the current works show that
N b L i · · · · has the strongest interaction with both one Y-wall and two X-walls. Therefore, reducing or replacing
N b L i · · · · with low interacting defects may be one way to overcome the known limitations of ferroelectric LiNbO3.