Andreas
Klein
* and
Denis
Sudarikov
Technical University of Darmstadt, Electronic Structure of Materials, Darmstadt, Germany. E-mail: aklein@esm.tu-darmstadt.de
First published on 12th March 2025
Defect chemistry is the classical approach to evaluate point-defect concentrations in solids depending on the chemical activity of the (n − 1) of n constituents by evaluating the mass action laws of a number of defect reactions conserving species, lattice sites, and charge. In an alternative approach, formation energies of individual defects can be calculated to determine the dependence on the Fermi level and on the chemical potentials of the reservoirs. This contribution provides the quantitative relationship between the two approaches, offering the opportunity to compare calculated defect formation energies with experimentally determined quantities. As an example, the application of the two approaches to the comparison of electronic and ionic compensation of doping and the influence of the band edge energies on it is given. This example demonstrates that the Gibbs energy of reduction and oxidation are essentially aligning the energy axis of ionic defects relative to that of electronic defects. In conjunction with the dependence of the valence band maximum and conduction band minimum energies on material composition, this offers the opportunity to tune the preference for electronic, ionic or mixed compensation of doping by two independent quantities.
This article will first provide a brief introduction to the two approaches for describing defect properties, restricted to the simple case of dilute, non-interacting defects and to a situation with only four different native defect species, an ionic acceptor and an ionic donor species and free electrons and holes in the conduction and valence bands. The defect reactions and their relationship to defect formation energies will be exemplified using the example of Schottky disorder (cation and anion vacancies as native defect species), but are equally valid for any other kind of disorder. This simplified scenario is sufficient to provide a quantitative relationship between equilibrium constants of defect reactions and calculated defect formation energies, which remain valid if the complexity is increased, for example by including valence changes of the ionic defects, trapped electronic carriers, and charged dopants. The influence of the latter is specifically elaborated for the important competition between ionic and electronic compensation of doping, which determines the nature of conductivity (ionic, electronic, or mixed). The scenario will also be used to highlight differences and synergies of the two descriptions of defect properties.
μi = μ0i + kBT![]() ![]() | (1) |
![]() | (2) |
Together with the tabulated values for the standard chemical of molecular oxygen, μ0(O2) (equivalent to the standard molar Gibbs energy of oxygen, ΔG0m(O2), see e.g. ref. 22), the relationship between the chemical potential and the partial pressure of oxygen is obtained from eqn (2). Data for 300, 600 and 1000 K are shown in Fig. 1. The significant influence of temperature is evident from the calculation.
![]() | ||
Fig. 1 Relationship between the oxygen chemical potential and the oxygen partial pressure as calculated from eqn (2) using data for the standard Gibbs energy of molecular oxygen22 (100 kJ mol−1 = 1.0364 eV ≈ 1 eV). The insert depicts the standard molar Gibbs energy of molecular oxygen, G0m(O2) = μ0 (data taken from ref. 22). |
Defect concentrations are obtained by evaluating the defect reactions specific to the material. When only charged intrinsic defects are assumed and defect association is neglected, the relevant ionic defects in a binary compound MXz are the metal and anion vacancies, vM and vX, and interstitials, Mi and Xi, respectively (we use lowercase “v” and “i” for vacancies and interstitials for clarity1). In most cases, only two ionic defects are important. The resulting different disorders are of Schottky (S: vM & vX), Frenkel (F: vM & Mi), anti-Schottky (: Mi and Xi), or anti-Frenkel (
: vX & Xi) type. The different types of disorder are independent but connected to each other by F +
− S =
. This contribution will restrict the comparison between the chemical and physical approaches for describing properties to the example of Schottky disorder, but it is equally valid for any kind of ionic disorder. For an oxide with composition MO, the ionic equilibrium for Schottky disorder is determined by transferring a metal and an oxygen species to the surface of the material, which is described using the following reaction:
![]() | (3) |
with the mass action law:
![]() | (4) |
In addition to the ionic defects, electronic disorder is related to the excitation of electrons from the valence band into the conduction band, thereby generating an electron in the conduction band and a hole in the valence band according to:
nil ⇌ e′ + h˙ | (5) |
with the mass action law:
![]() | (6) |
Here, the band gap of the material, Eg, is the free energy of formation of the electron–hole pair, n and p are the concentrations of free (non-trapped) electrons and holes, and NC and NV are the effective density of states in the conduction and valence band for non-degenerate semiconductors with parabolic band dispersion. The link connecting electronic and ionic disorder to each other (see the section on ionic vs. electronic compensation for an explanation) and to the oxygen partial pressure is the reduction (or oxidation) reaction of the sample. For example, the reduction of an oxide may proceed according to:
![]() | (7) |
The mass action law related to this reduction is
![]() | (8) |
Instead of the reduction reaction, it is sometimes more convenient to consider the oxidation reaction. To be consistent with the discussion below, we express the oxidation reaction for Schottky disorder as:
![]() | (9) |
![]() | (10) |
The significance of the distinction between the reduction and oxidation reaction will become evident during the discussion of the influence of the band edge energies.
Eqn (3)–(8) include only doubly charged oxygen and metal vacancies. In principle, also singly charged and neutral vacancies may be present. For every added defect species, an additional reaction is required. For differently charged ionic defects, these would be the ionization reactions,3 which are omitted here for clarity. Eventually, the defect concentrations for a compound MO with Schottky and electronic disorder are obtained from the Gibbs energies of Schottky disorder and of the reduction reaction by solving the set of the three mass action laws, eqn (4), (6) and (8) (or (10)), in conjunction with the charge neutrality condition. The latter is given for Schottky disorder with the concentrations (per volume) of doubly charged metal and oxygen vacancies, , by:
![]() | (11) |
For the charge neutrality condition, concentrations of defects have to be used instead of the site fractions as the number of potential defect sites can differ for defects, for example, the concentration of potential defect sites for oxygen vacancies is the concentration of oxygen lattice sites, while that of the electronic defects is the effective density of states, NC and NV, which is about three orders of magnitude lower than that of the lattice sites.
One possible practical solution to the set of equations is obtained by calculating the concentrations of all defects depending on the concentrations of, for example, oxygen vacancies for a given set of oxygen partial pressure and equilibrium constants K of the three reactions. For the given example it is convenient to start with eqn (4) to obtain the metal vacancy concentration as a function of oxygen vacancy concentration. Eqn (8) then provides the electron concentration, and eqn (6) the hole concentration as a function of . With the dependence of all defect concentrations on
in the form
, the equilibrium oxygen vacancy concentration,
, is obtained by directly solving the charge neutrality equation, eqn (11). Eventually, all equilibrium defect concentrations are obtained from the equilibrium oxygen vacancy concentration,
. Doing the calculation with a variation of the oxygen partial pressure finally reveals the Brouwer diagram, i.e. the concentrations of all defects depending on oxygen partial pressure. An example of this kind of diagram is presented in Fig. 2(a).
The charge state of the defect is accounted for by adding either positive or negative background charges. As for atoms, addition or removal of an electron from the supercell changes its energy by the chemical potential of the electron μe. In summary, the defect formation energy ΔGqD of a defect in charge state q can be written as:17
![]() | (12) |
![]() | (13) |
It is mentioned that the DFT calculations reveal formation energies of the supercells at zero Kelvin, HD(0 K). The temperature dependence of the enthalpy and entropy contribution to the Gibbs formation energy of defects are not included. It is also mentioned that additional posteriori corrections to eqn (12) are required, for example, to account for effects related to the small size of the supercells. These corrections do not affect the dependence of ΔGqD on the chemical potentials and the Fermi level. Details for the corrections can be found, for example, in ref. 18.
According to eqn (12), the Gibbs formation energies of charged defects depend linearly on the Fermi energy. This can be rationalized as the electrons added or removed from a defect have an energy corresponding to the (electro)chemical potential of the electrons, which is called the Fermi level.20 Following eqn (12), it is sufficient to calculate the defect formation energy in a given charge state at any value of the Fermi level and chemical potential. To display the result of the calculation, the defect formation energies are typically plotted as a function of the Fermi level for the most oxidizing and most reducing conditions. An example for a simple hypothetical binary metal oxide MO with Schottky disorder is given in Fig. 2(b). Only doubly charged states of the metal and oxygen vacancies are considered for clarity.
Practically, the equilibrium defect concentrations are obtained by calculating the defect concentrations as a function of the Fermi level using the defect formation energies for a given chemical potential of oxygen. The equilibrium Fermi level is then determined by solving the charge neutrality eqn (11) and used to obtain the equilibrium defect concentrations. Repeating the calculations for the allowed range of oxygen chemical potentials finally results in the Brouwer diagram. The formation energies given in Fig. 2(b) reveal defect concentrations at 1000 °C identical to those displayed in Fig. 2(a).
![]() | (14) |
Although the Gibbs formation energies of the metal and oxygen vacancies both depend on the chemical potentials and on the Fermi level, the sum of the two does not. Whenever the Gibbs formation energy of one of the two defects is changing with EF or μO, the formation energy of the other defect is changing by the same amount in the opposite direction. The first is the case due to the opposite charge state of the metal and oxygen vacancy in eqn (12) and the second is the case as the chemical potential of the oxygen and the metal in the last term of eqn (12) are connected by:
μ(O) + μ(M) = μ(MO) | (15) |
The relationship for the free reaction energy of the reduction can be obtained via the mass action law for the reduction reaction (8). In this equation, the site fraction for oxygen vacancies, , which is a function of the Fermi level and the oxygen chemical potential, is replaced by
, that for electrons, [e′], which is a function of the Fermi level, by exp(−(ECB − EF)/kT), and the oxygen partial pressure by using the expression given in eqn (2). This leads to:
![]() | (16) |
Taking the logarithm of (16) and multiplying the result by kBT reveals:
![]() | (17) |
This equation is valid for any combination of EF and μO. It is convenient to evaluate eqn (17) for μO = 0 eV (oxygen-rich limit) and EF = ECB (Fermi level at the conduction band minimum), resulting in:
![]() | (18) |
While ΔGred can be obtained from experiment (for BaTiO3 see e.g., ref. 7 and references therein), is only accessible from computations (for BaTiO3 see, e.g., ref. 11, 15, and 23). Therefore, eqn (18) directly offers the opportunity to compare calculated defect formation energies with experimentally derived quantities, which is not possible otherwise. The difference between ΔGred and
is given by
, which is a function of temperature. The oxygen chemical potential can be obtained from thermodynamic data22 and can be up to several electronvolts as illustrated by the insert in Fig. 1. For given formation energies of the cation and anion vacancies, eqn (14) and (18) provide the quantitative relationship to calculate the free energies and equilibrium reaction constants for Schottky disorder and the reduction reaction. For the example given in Fig. 2, the sum of the formation energies of the metal and oxygen vacancy is 2.5 eV at any value of the Fermi energy and oxygen chemical potential. The formation energy of the oxygen vacancy at the conduction band minimum (EF = ECB) under oxidizing conditions (μO = 0 eV is 6.2 eV, extrapolation of the red dashed line to E = ECB in Fig. 2(b)). To obtain the free reduction energy, the standard molar free formation energy of oxygen molecules, μ0(O2) at 1273 K, is required. This is obtained from thermodynamic reference data22 as −287.67 kJ mol−1 = −2.98 eV, resulting in ΔGred = 6.2 − 1.5 = 4.7 eV, which is used to calculate the Brouwer diagram in Fig. 2(a). The value of μ0(O2) can also be obtained by extrapolating the respective curve in Fig. 1 to μO = 0 and then evaluating eqn (2).
The connection between the formation energy of the metal vacancy and the Gibbs energy of the oxidation can be obtained in the same manner as that for the reduction (eqn (16)) as:
![]() | (19) |
The free reaction energy of the oxidation, ΔGox, is connected to that of the reduction, ΔGred, the Schottky pair, ΔGS, and the electron–hole pair, ΔGe–h = Eg, via:
ΔGox = ΔGS + 2·Eg − ΔGred | (20) |
The relationship between Brouwer diagrams and Gibbs energies of the formation of metal and oxygen vacancies of a binary metal oxide MO with Schottky disorder for three different values of the Gibbs energy of reduction is illustrated in Fig. 3. Three scenarios with identical Gibbs energy for Schottky disorder (ΔGS = 4.5 eV) and identical band gaps (Eg = 2.7 eV) but different Gibbs energies of reduction are compared. The scenario in the middle column of Fig. 3 is identical to the one depicted in Fig. 2. In the left and right columns, the Gibbs energy of reduction is changed by ±1.3 eV. Lowering (left column) and raising (right column) the Gibbs energy of reduction shifts the Brouwer diagram towards higher and lower oxygen activity, respectively.
The Gibbs formation energies of metal and cation vacancies representing the three scenarios are given below the respective Brouwer diagrams. As expressed by eqn (18), lowering the Gibbs energy of reduction lowers the Gibbs formation energy of the oxygen vacancy by the same amount for both oxidizing and reducing conditions (dashed and solid red lines). As the Gibbs energy for Schottky disorder is the same in all three scenarios, the changes of the Gibbs formation energy of the metal vacancy are opposite to those of the oxygen vacancy. This leads to an apparent shift of the formation energies on the energy axis, which amounts to half of the change in the Gibbs energy of reduction due to the double charge of the defects considered.
The difference between ionic and electronic compensation is directly related to the formation energies of the ionic defects.20,27,28 Donor or acceptor doping adds or removes electrons from the host material, which raises or lowers the Fermi level, respectively. As illustrated in Fig. 2 and in the lower row of Fig. 3, raising the Fermi level will lower the Gibbs formation energy of the intrinsic acceptor (metal vacancy), while lowering the Fermi level will reduce the formation energy of the intrinsic donor (oxygen vacancy). Eventually, the formation energies may become zero or even negative. Negative formation energies correspond to a thermodynamic instability. If the Gibbs formation energies of defects are negative, the material can gain energy unlimited by the formation of defects with an equilibrium defect concentration being higher than the concentration of potential defect sites ([…] > 1). If not kinetically limited, the decrease in the formation energy towards zero will result in the formation of a high concentration of ionic defects. These prevent the Fermi level from shifting further. Using this picture, ionic compensation will occur if the Gibbs energy of the formation of the ionic defect becomes zero before the Fermi level reaches either the conduction band minimum upon donor doping or the valence band maximum upon acceptor doping. Electronic compensation corresponds to the formation of electrons or holes, which appears if the Fermi level approaches the conduction or valence band, respectively. The limits of the Fermi level for the different conditions are indicated in the bottom row of Fig. 3 by red arrows for acceptor doping and by blue arrows for donor doping. The limits are depending on the chemical boundary conditions, i.e. the oxygen partial pressure in the present scenario. Based on this description, the occurrence of ionic compensation can be clearly identified. For ΔGred = 3.4 eV, for example, oxidizing conditions lead to ionic compensation for both donor and acceptor doping, while reducing conditions lead to ionic compensation for acceptor doping and to electronic compensation for donor doping. The respective compensation mechanisms are indicated in the graphs in the upper row of Fig. 3.
The relationship between the Gibbs energy of reduction, i.e. the Gibbs formation energy of the oxygen vacancy, and the charge compensation for the case of donor doping is illustrated in Fig. 4. For the calculation, the concentration of singly positively charged donors is included in the charge neutrality eqn (11) when determining the equilibrium concentrations. For simplicity, a shallow donor with a valence independent of the Fermi level is used. For donor doping, charge compensation can be either by electrons or, for the case of Schottky disorder, the metal vacancy. To illustrate the effect, Gibbs energies for Schottky disorder of ΔGS = −3.2 eV and for reduction of ΔGred = −5.4 eV and a band gap of Eg = 3.0 eV are used for the calculation of defect concentrations as a function of donor concentration. The Brouwer diagram for the undoped material is again symmetric, while adding 1018 cm−3 donors leads to an increase of the concentration of metal vacancies for oxygen chemical potentials μO > −3 eV. The dependence of defect concentrations on doping concentrations for μO > −3 eV (case (1) in Fig. 4) are displayed in the lower left graph of Fig. 4. For all doping concentrations, the metal vacancy remains the dominant compensating defect.
The situation is changing with decreasing oxygen activity. For μO = −3.7 eV, case (2) in Fig. 4, electrons are the dominant compensating defect for donor concentrations up to ND ≈ 3 × 1018 cm−3. For higher ND, metal vacancies dominate. In this situation, which corresponds to a mixed compensation, the formation energy of the metal vacancy vanishes exactly at the conduction band minimum (see the upper right graph in Fig. 4). For even more reducing conditions, case (3) in Fig. 4, electronic compensation is dominating for the whole range of donor concentration. Fig. 4 illustrates the change in compensation mechanism for donor doping and its relationship to the formation energy of the metal vacancies and the conduction band minimum. In none of the cases illustrated in Fig. 4 is saturation of the electron carrier concentration with increasing doping concentration observed. This kind of saturation is expected as the Fermi level cannot increase beyond the energy, at which the formation energy of the cation vacancy vanishes. The saturation will occur but at much higher doping concentrations. These are not included in Fig. 4 as the approximation of dilute defect concentrations is not valid in this case. The case of acceptor doping and its relationship to the formation energy of the oxygen vacancy and the valence band maximum can be demonstrated analogously for oxygen chemical potentials of μO = −1.5…0 eV.
The examples illustrate that the preference for electronic or ionic compensation is determined by the defect formation energy of the ionic defects. Comparing the magnitudes of the formation energy of the ionic defects, e.g. Schottky disorder, and the electronic defects is not sufficient. In order to reveal the compensation mechanism, it is necessary to plot the formation energy of the ionic defects on the same energy axis as the formation energies of the electronic defects as in Fig. 2(b). Keeping the formation energy of the electronic defects fixed by fixing the energies of the band edges, it is the energy of the reduction (or oxidation) reaction, which shifts the formation energy of the ionic defects relative to the electronic ones as illustrated by the top graphs in Fig. 3. As the plots with the formation energy directly reveal the preferred compensation mechanism, the representation of the defect properties in the physical picture is more intuitive in this case.
As discussed in the previous section, the preference for electronic or ionic compensation is determined by the energies, at which the formation energies of the respective defects vanish. The preference can hence be changed either by moving the formation energies of the ionic defects or by moving the formation energies of the electronic defects on the common energy axis. In principle, the two effects can be considered independently. The discussion will start with the variation of the ionic defects.
According to eqn (18) and (19), the Gibbs energies of reduction and oxidation, ΔGred and ΔGox, affect the formation energy of the oxygen and metal vacancies as illustrated at the right in Fig. 5. The effect of ΔGred and ΔGox on the formation energies has to be added to the difference between oxidizing (high p(O2)) and reducing environments (low p(O2)). For a fixed formation energy of Schottky disorder, as exemplified in Fig. 2, 3 and 5, the formation energy of the metal and oxygen vacancies are inversely affected by changes of ΔGred and ΔGox if the metal has a valency of +II. For materials with different stoichiometry, such as M2O, M2O3, etc., the change of the formation energy of the cation vacancy has a different magnitude than that of the oxygen vacancy.
If a material is more resistive to reduction, the Gibbs energy of reduction increases and that of the formation of oxygen vacancies (blue lines in Fig. 5) increases accordingly. Consequently, the lower limit of the Fermi level, labelled Eredlo and Eoxlo in Fig. 5 moves downwards on the energy axis. If a material is more resistive to oxidation, the upper limits of the Fermi level, labelled Eredhi and Eoxhi in Fig. 5, move upwards on the energy axis. With the labels for the upper and lower limits of the Fermi level in Fig. 5, the condition for electronic and ionic compensation of acceptor and donor doping, which are the essence of the practical doping principles of semiconductors,33 can be expressed as:
![]() | (21) |
The distinction between electronic and ionic compensation according to eqn (21) also includes the variation between oxidizing and reducing conditions. For example, donor-doped SrTiO3 exhibits electronic compensation under reducing but ionic compensation under oxidizing conditions.24 In the context of Fig. 5, this can be expressed by Eredhi > ECB > Eoxhi.
In addition to changing the preference for electronic or ionic compensation by changing the Gibbs energies of reduction and oxidation, the preference is also changed when the formation energies of the ionic defects is kept constant on the common energy axis but the formation energies of the electronic defects are shifted in energy as illustrated on the left side of Fig. 5. The latter are determined by the position of the band edges. By definition, the formation energy of electrons will be zero at ECB and that of holes will be zero at EVB. Therefore, moving the position of the band edges up or down as indicated on the left side of Fig. 5 can be used to change the preference for ionic or electronic compensation. This dependence of the charge compensation mechanism on the band edge energies is well known in semiconductor physics and explains the different dopability of materials, i.e., the achievable electronic carrier concentration upon donor or acceptor doping.25–27,33–39
It is also well-known that the energies of the valence band maximum and conduction band minimum can vary quite substantially between materials. The valence band maximum energy, for example, is affected by the choice of the anion due to the energy of the p-orbitals,37,40 but also by cationic contributions to the valence band density of states. Such cationic contributions are present for transition metal compounds with shallow d orbitals, e.g. Cu, or partially, e.g. Fe, but also for compounds containing heavier ions with filled s-orbitals, such as Pb2+, Sn2+ or Bi3+.36,37,41–43 These contributions result in an upward shift in the valence band maximum, which is favorable for electronic compensation by holes. In fact, room temperature p-type electronic conduction by holes is exclusively found in materials with a high valence band maximum energy originating from such cation orbital contributions to the valence band.27,44–48 Further discussion of the factors determining the energy band alignment of materials can be found, for example, in ref. 36, 37, and 43.
The competition between ionic and electronic compensation decides if a material becomes an ionic conductor (or an insulator if the ionic defects are not mobile), an electronic (semi)conductor, or a mixed ionic-electronic conductor. The prevailing compensation mechanism can directly be extracted from the dependence of the defect formation energies on the Fermi level in the physical approach. In contrast, the equilibrium constants of the defect reactions do not provide a direct connection to the compensation mechanism. A simple comparison of the Gibbs energies of ionic and electronic disorder is not sufficient. This is not only the case because the charge states of ionic and electronic defects are different. More importantly, it is necessary to connect the Fermi level dependence of the formation energies of electronic and ionic defects on a common energy axis in order to extract the prevailing compensation mechanism for donor and acceptor doping. It has been demonstrated that the connection of the energy axes of ionic and electronic defects is directly provided by the Gibbs energy of reduction and oxidation. While the latter two are connected in a specific material by its ionic disorder and the band gap (eqn (20)), they can, in principle, be tuned independently by material composition. In addition to tuning the Gibbs energies of reduction and oxidation, the prevalent charge compensation mechanism is affected by the valence band maximum and conduction band minimum energies. A high valence band maximum favors electronic compensation for acceptor doping and a low conduction band minimum favors electronic compensation for donor doping. The energies of the band edges also depend on material composition.
As the band edge energies and the Gibbs energies of reduction and oxidation both depend on the material composition, it would, perspectively, be interesting to figure out if the two dependencies are connected or independent. To our knowledge, no direct connection between the band edge energies and the Gibbs energies of reduction/oxidation has been established. It is unclear whether this kind of relationship is physically reasonable at all but it is interesting to note that Cu-containing materials, which exhibit a significantly higher valence band maximum than comparable compounds without shallow occupied d-states,37,43 typically show a high concentration of Cu vacancies,49,50 indicating that they are easy to oxidize. This indicates that it may not be possible to tune the electronic and ionic disorder independently, which would considerably reduce the degree of freedom in the design of materials with the desired electrical properties. It would therefore be important to evaluate whether such a fundamental connection exists or not.
This journal is © the Owner Societies 2025 |