Small Polaron Hopping in Fe : LiNbO 3 as a Function of Temperature and Composition

Small-polaron hopping involved in charge transport in Fe-doped congruent lithium niobate is investigated as a function of temperature and composition by means of light-induced transient absorption spectroscopy. The relaxation dynamics of the light-induced polaron population is characterized by individual activation energies within different temperature ranges. A numerical investigation carried out by Monte Carlo simulations reveals that these findings may be understood in terms of the varying abundance of the different types of hops that the polarons may perform among regular or defective lattice sites. The role of the temperature and of the sample composition on the distribution of the different hop types is thus explored for a wide range of parameters, allowing one to preview the charge transport properties for a given set of experimental conditions.


Introduction
The study of small-polaron charge transport is emerging as a hot topic for the important ferroelectric crystal lithium niobate (LiNbO 3 , LN).In the field of nonlinear, ultra-fast and integrated optics, it is the physical mechanism determining how photo-excited charges redistribute in the material, giving rise to internal space charge fields that affect, via the electro-optic effect, the interaction of the material with laser beams [1,2].This is the origin of the so-called optical damage that often constitutes a major problem for laser applications.On the other hand, in photorefractive optics and holography, those mechanisms allow for low-intensity optical interaction and for recording of high fidelity optical gratings and real-time holograms [3].In analogous materials, polaron formation and transport play a key role as well in a number of important applications such as ferroelectric photovoltaics [4][5][6] and in oxide catalysis [7][8][9].
As demonstrated in [10], by increasing the temperature and/or changing the sample composition, the conduction of LN increases by orders of magnitude, canceling out any photo-generated internal field.Analogous effects may be obtained by doping the crystal with the so-called Optical Damage-Resistant Ions (ODRIs) such as Mg, Zn or Zr.This is obviously of great importance for practical applications because it allows one to control the material response against light.However, to date, there has not been a quantitative interpretation of those phenomena in terms of polaron transport, so that it appears mandatory to check the dependence of the above-mentioned properties in this framework.
Polarons are quasiparticles made up of an electrical charge that, by the interaction with the polar crystalline environment, is able to distort the neighboring lattice, creating a local potential well.As a net result, the particle becomes self-localized.If the confinement effect is strong, the charge is concentrated on a single lattice site (small strong-coupling polaron) and moves by randomly hopping among different sites in response to lattice thermal motion.The hopping frequency of the polaron can be computed by the Marcus-Holstein theory [11][12][13].In an ideal, defect-free material, the polaron mobility can be then calculated using the standard Einstein relation [14]: where e is the elementary charge, d is the characteristic hopping length (here equal to the anionic distance), kT is the thermal energy and U is the hopping barrier, which in standard polaron theory is equal to one half of the polaron energy U = E/2.The term w 0 exp − U kT is thus the hopping frequency among equally-spaced identical sites.The pre-factor w 0 can be also obtained by the Marcus-Holstein theory, as will be explained in the following.
Unfortunately, in real materials, this situation is hardly realized: standard LN crystals of congruent composition contain a high density of point defects, such as the substitutional "antisite" defect Nb 5+ Li or the Li vacancy V − Li , bothextrinsic (dopants or impurities).Those defects constitute preferential sites for the formation of polarons, so that three types of small electron polarons are intrinsically recognized in LN [15]: the free polaron FPforming on regular Nb 5+  Nb sites; the bound polaron GPforming on antisite defects Nb 4+  Li ; the bipolaron BP, which is a combination of a free and of a bound polaron on neighboring sites Nb 4+ Nb : Nb 4+ Li .LN may also host O − holes, created when energetic photons capable of band-to-band excitation are involved.In this case, hole polarons HPare formed.Finally, impurities with two valence, states such as the prototypical case of Fe 2+/3+ Li , have the possibility to capture an electron with a potential that is a mixture of Coulomb attraction and polaronic lattice deformation.Those charged defects can be described in the polaron formalism, as well [15].The energy barriers for hopping processes involving those sites are quite different from each other, so that there now is a set of possible activation energies instead of a single one.Moreover, the w 0 term in Equation (1) depends on the distance between the hopping sites, and if direct hopping among defective sites is considered, this may entail an anomalous diffusion regime for which Einstein's relation is no longer valid [16].As a consequence, the relation (1) cannot be used in real crystals.Nevertheless, in recent works [17,18], it was shown that by an appropriate choice of a few microscopic parameters, the problem can be tackled numerically by resorting to Monte Carlo simulations incorporating the presence of defects in the lattice structure of LN.In line with those studies, one may expect that depending on the temperature and composition of the crystal, one type of hop may dominate the transport process.
The aim of this work is thus to quantify the number and type of hops among the different polaronic centers upon variation of the temperature and composition.We will focus on the case of Fe doping in LN, because here, we can assume that, under the experimental conditions used, all the charges are generated by photoexcitation of Fe 2+ centers.This allows us to neglect the presence of hole polarons (that would have been present in undoped material [19]) for which the polaronic parameters are not well known yet.At the same time, Fe 3+ can be assumed to be the sole deep trap state.The situation herein addressed is of practical importance because it corresponds to the typical situation encountered in photorefractive experiments.The problem is attacked both experimentally and numerically.Transient absorption data are collected as a function of temperature to provide direct evidence of a change in the activation energy of the transport, hinting at the activation of different hopping processes with temperature.The situation is then explored in greater detail with the help of a dedicated simulation software based on the Monte Carlo approach.

Samples
Three Fe:CLN crystal boules were Czochralski-grown at the University of Padova mixing together congruent LiNbO 3 (Sigma Aldrich, St. Louis, MO, USA) and Fe 2 O 3 powders, with Fe contents in the melt of 0.02%, 0.05% and 0.1% mol, respectively.Being grown in congruent conditions, the Li/Nb ratio of the final crystals was assumed to be 0.94 (neglecting the compensation of intrinsic defects by Fe impurities).We verified that small deviations from this value and yet below the typical compositional uncertainties of the samples produced in our laboratories had no effect on the results presented in the following paragraphs.Parallelepipeds 4 × 5 × 1 mm 3 cut from these crystals were optically polished.The sample with the highest doping level was chemically reduced under a mixture of Ar − H 2 atmosphere at 500 K for several hours.The absorption coefficient at 532 nm was measured for samples Fe:LN 0.02 mol % and Fe:LN 0.05 mol %, from which the Fe 2+ concentration was obtained [20].For sample Fe:LN 0.1 mol %, instead, due to the high absorption level, we preferred to use the datum at 1128 nm [21] from which the Fe 2+ concentration was also deduced.The total Fe content was assumed for simplicity to be the same as in the melt, and the Fe 3+ concentration was obtained by the difference.The samples compositions are given in Table 1.

Transient Absorption
The polaron decays were studied by light-induced Transient Absorption Spectroscopy (TAS) using the setup described elsewhere [22,23] based on a Q-switched, frequency-doubled Nd:YAG laser at a wavelength of λ = 532 nm with a pulse duration of τ FWHM 8 ns and a maximum pulse energy of E max = (290 ± 10) mJ.The pump intensities in the various experiments were adjusted to maximize the signal to noise ratio of the experimental decay: for Fe:LN 0.02 mol %, I ∼ 60 MW/cm 2 ; for Fe:LN 0.05 mol %, I ∼ 200 MW/cm 2 ; and for Fe:LN 0.1 mol %, I ∼ 40 MW/cm 2 .Four ordinarily-polarized continuous-wave probe lasers enabled a time-resolved detection of changes in the absorption simultaneously at different wavelengths (λ = 445 nm, 488 nm, 785 nm and 633 nm or 1310 nm).The sample temperature was set by a cryostat equipped with four plane-parallel optical windows for all-optical access.
In Figure 1, typical TA measurements obtained on sample Fe:LN 0.1 mol % are shown at T = 298 K and 198 K, probing simultaneously the transient absorption at different combinations of the wavelengths: 488, 785, 1310 nm.
All the data demonstrated the previously-reported behavior [22,23]: after the pulse, a marked increase of the absorption was detected for some of the probe wavelengths, which decay to their dark value after some time, with a dynamics that depends on the probe.The signal at 785 nm was attributed mainly to bound polaron absorption, in view of the spectral dependence of the polaron cross-section, and to a minor extent, to the absorption tail of the free polaron at this wavelength.In the following, we will therefore consider that this signal is essentially proportional to the instantaneous concentration of GP decaying towards Fe traps, with a possible correction for the presence of free polarons.Similar considerations indicated that the signal at 1310 nm can be attributed essentially to free polarons with a contribution due to the bound polarons' absorption.The initial amplitude of the TA signal in the red/infrared correlated directly with the concentration Fe 2+ present in the samples (see Table 1), suggesting that Fe 2+ was the main source for photoexcitation of the centers probed at 785 and 1310 nm [22].In the following, we will restrict our attention to these probes and to free and bound polarons.At room temperature, the absorption signals decayed to zero in the range of milliseconds, while at a low temperature, the processes were slowed down.The observed decay spectra can in general be described using an empirical stretched exponential function, named after Kohlrausch, Williams and Watts (KWW) [24,25].This so-called KWW function is given by: where α 0 TA is the absorption coefficient for t → 0 (not to be confused with the signal amplitude at the beginning of the time window of our data), β is the stretching factor 0 < β < 1, where β = 1 denotes the case of a mono-exponential relaxation, and τ K is the characteristic decay time.

Decay Model
The pump pulse generated an initial polaron population by direct excitation from Fe 2+ donors.Due to the moderate pump pulse intensity, the photo-excitation of electrons from the valence to the conduction band by the two-photon processes was disfavored [22], so that the presence of hole polarons was neglected here.Furthermore, by comparing the absorption coefficients with the peak photon fluence of our pump, we may consider that in samples Fe:LN 0.02 mol % and Fe:LN 0.05 mol %, all Fe 2+ ions were likely to be ionized after the pulse so that the initial concentration of deep Fe 3+ traps can be assumed to be equal to the total Fe concentration.In sample Fe:LN 0.1 mol %, instead, due to the high Fe 2+ content, the absorption of the pump beam along the sample thickness cannot be neglected so that the number of photo-excited carriers became depth-dependent.To deal with this, we considered in this case an effective initial population calculated from the amplitude of the absorption signal by comparison with the other samples.Finally, since the polaron formation time (<1 ps, see [1]) was much shorter than our pulse duration, we can assume that the polarons' position after the pulse was completely uncorrelated from the initial photogeneration site [26].The so-produced polaron population performed a random walk hopping either among defective Nb Li sites (bound polarons) or regular Nb Nb sites (free polarons) until they were trapped by a Fe 3+ ion.Both bound and free polarons contributed to the absorption at 785 nm [23].

Simulation Code
According to the Marcus-Holstein polaron hopping model [11][12][13], the non-adiabatic hopping frequency for a (i → f ) hop was: In Equation (3), r is the distance between initial and final sites and kT the absolute temperature (in energy units).λ i, f is the reorganization energy of Marcus' theory corresponding to the energy paid to rearrange the lattice, here equal to (E i + E f ), the sum of the elastic energies of the two polarons.a i, f is an orbital parameter describing the overlap between the electronic wavefunctions at site i and f .U i, f = U f ,i is the hopping barrier, given by [17,27]: with ε i and ε f the binding energies of the electron at zero deformation.When the hopping was between sites of the same type (i = f ), U i,i = E i /2, recovering the standard result that the hopping activation energy was one half of the polaron stabilization energy [1,14,15].The pre-exponential factor I i, f describes the intrinsic hopping rate between the two sites and was determined by the choice of the (i, f ) combination.Using the data reported in [15,18] the hopping barrier (4) for any process and the hopping rates (3) can be calculated.This result was used in a dedicated Monte Carlo code based on a classical Gillespie algorithm.The LN structure was generated at the beginning of the code with a randomly-placed set of Nb Li and Fe defects, the concentration of which was set according to Table 1, assuming that they were the same as in the melt.The Nb 4+ Li and Nb 4+ Nb population decays were computed by an MC code in a 80 × 80 × 80 super-cell of the LN structure with periodic boundary conditions.

Activation Energies Deduced via TAS
According to the discussion reported elsewhere [17,18,22,23] it can be considered that α TA (785 nm, t) was proportional to the weighted sum of the bound and free polarons populations.
Since the signal at 1310 nm remained always much smaller than the one at 785 nm and in general proportional to it, we may in a first approximation consider that the signal at 785 nm was proportional to the instantaneous bound polaron concentration and that the signal at 1310 was due essentially to the absorption tail of the bound polaron.From the measured KWW parameters, the average bound polaron lifetime can thus be computed: where Γ(x) is the gamma Euler function.This approach has the advantage of considering both the parameters τ K and β, allowing for a comparison between decays with different stretching exponents [28].The result of the τ calculation for the three samples is shown in Figure 2 as a function of the temperature in an Arrhenius plot.The results gave evidence of two distinct activation energies with a transition temperature around 180 K for all samples.The high temperature activation energy was the same for all three samples U(T > 200 K) = (0.42 ± 0.01) eV.At low T, it was possible to measure an activation energy only for the less doped sample, providing U(T < 200 K) = (0.08 ± 0.01) eV.Nevertheless, all the samples appeared to follow the same qualitative behavior.What was the origin of the observed slope change?The total mean lifetime defined in Equation ( 5) was determined mainly by the slowest processes leading to the bound polaron population disappearance.In this vision, the slope of the Arrhenius plot corresponded to the activation energy of the slowest hopping process.However, if more processes with similar hopping frequency occurred in a comparable amount, the observed activation energy would somehow lie in the middle.Using Equation ( 4) and the known values of the polaron energies (see [18] and the references therein), it is possible to compute the characteristic hopping barriers for the different processes, reported in Figure 3.As can be verified by inserting some typical values in Equation (3), the processes requiring the longest time to hop away from an antisite were the Nb Li − Nb Li and Nb Li − Nb Nb hops, characterized by an activation energy of 0.375 eV and 0.55 eV.The measured activation energy was found in between (0.42 eV), accordingly.This means that the polaron, initially stuck on an antisite defect, either hopped directly towards another antisite or became a free polaron (which can eventually move on the Nb Nb lattice or re-captured by an antisite, but since those processes were very fast, they did not significantly affect the polaron lifetime).4)) calculated using the values for the different polaron energies detailed in [18].FP, free polaron.

Fe
Those claims were quantitatively checked against the numerical simulations for the sample Fe:LN 0.02 mol %, for which two distinct slopes are clearly visible in the Arrhenius plot.An initial population of 5000 polarons, all forming on a Nb 4+  Nb site just after the pulse, were considered.The TAS signal at 785 nm was computed by considering the weighted sum of the FP and GP populations α LI (785 nm, t) ∝ σ GP (785 nm)N GP (t) + σ FP (785 nm)N FP (t), where σ GP (785 nm) = 7 × 10 −22 m and σ FP (785 nm) = 3.1 × 10 −22 m were the bound and free polaron cross-sections at 785 nm, respectively [23].The numerically-computed decays were fitted by a KWW function in the experimental range starting from 10 −8 s, and the average lifetime was calculated using Equation (5).The results are plotted in Figure 2. It can be seen that in the range between 180 and 400 K, the agreement between the simulated and experimental data was excellent, and the same situation was obtained for the other samples (not shown).Below 180 K, however, the agreement is less satisfactory: although the simulated data correctly reproduced the fact that the activation energy passed from U ∼ 0.42 eV to U ∼ 0.08 eV, the experimentally-measured transition appeared to be sharper, and the simulated decays were several times slower than the experimental ones.One possible explanation for this disagreement was that at low temperatures, the Marcus-Holstein hopping frequency became inaccurate [13,14].Another reason may stem from the fact that in our simulation, we disregarded the presence of holes, and this assumption may be not be fully adequate in sample Fe:LN 0.02 mol %, which contained the lowest concentration of Fe.A third possible explanation was that the absorption at 785 nm may have contained an additional contribution from some other center than GP and FP, decaying by an alternative relaxation channel.This problem was already apparent also in [18] in which the sample Fe:LN 0.02 mol % was analyzed from the point of view of the decay shape, obtaining a set of parameters somehow different from that in other samples, especially concerning the trapping parameter a GP−Fe = c.
Despite the partial agreement, our simulations reproduced qualitatively the observed activation energy change.An inspection of the results indicated that at low temperatures, the number of hops decreased significantly to a few units and, in many cases, even to a single trapping event.Therefore, the change in the activation energy can be associated with what we may call the "trapping" regime, in which the polaron performed a very limited hop number, so that the pertinent activation energy was somehow in between the energy barrier for the different trapping processes: Nb 4+  Li − Fe 3+ Li and Nb 4+  Nb − Fe 3+ Li (see Figure 3).

Hopping Regimes' Landscapes Inspected by Monte Carlo Simulations
One of the most intriguing question concerns the understanding of the role played by the different types of polarons, according to temperature, Nb Li antisites and Fe concentration.With the analysis of the activation energy of the mean lifetime, we obtained an indication that, depending on the temperature, different hopping processes can take place.It was only possible to conclude that at a high temperature, a mixed transport made by Nb 4+  Nb free and Nb 4+ Li bound polarons occurred, while at a low temperature, the iron had a greater impact on the results.These respective parts in light-induced transport are still an open question.The present section aims to study the interplay between these species and predict which type of transport should predominate, depending on temperature and trap concentrations, with the help of the Monte Carlo simulation.

Free Polaron Transport
To evaluate quantitatively the role of FP in charge transport, we can calculate, for a polaron performing N tot hops from its birth to its death, the fraction N FP−FP /N tot of Nb 4+  Nb → Nb 4+ Nb hops performed as a free polaron.This quantity was then averaged over a large number of realizations.The calculation was performed by simulating 5000 polarons using the microscopic parameters reported in [18].The experimental conditions (antisites and Fe traps concentration and temperature) were varied, to explore the typical range of parameters encountered in standard situations.
On the left side of Figure 4, we investigate the effect of the sample composition and of the temperature on the free polaron contribution to conduction.Here, the Fe concentration was set to zero.This situation therefore describes the impact of the stoichiometry or of doping with photorefractive-resistant ions (Mg, Zn, Zr) on the conduction.It can be seen that at room temperature, the reduction of the [Nb Li ] concentration initially had a small influence on the amount of FP-FP hops.It was necessary to eliminate about 90% of the antisites before observing a steep increase in the free polaron contribution to the conduction.This result explained quantitatively the origin of the so-called photorefractive threshold, which was necessary to surpass in order to observe a marked decrease of the photorefractive effect in undoped LN [29].The situation was even more extreme at low temperatures and, conversely, more gradual by heating the sample.On the right side of Figure 4, we consider now the impact of Fe doping from the point of view of free polaron conduction.We considered a stoichiometric sample to distinguish the effect of shallow Nb 5+  Li traps.The qualitative behavior was obviously similar.In this case, it is interesting to note that at room temperature, even for high Fe concentrations (2 × 10 25 m −3 ), the free polarons were always able to perform on average several hops before being trapped.This is in contrast, as we will show below, with the behavior of bound polarons.

Bound Polaron Transport
As has been shown in the above paragraph, when some antisites were present and for not too elevated temperatures (Figure 4 (left)), only a limited percentage of the hops that a polaron performed during its life were of the Nb 4+ Nb → Nb 4+ Nb type.In Figure 5 (left), the same analysis is carried out for Nb 4+  Li → Nb 4+ Li hops.Here, the situation was reversed: if the temperature was not too high and the antisite concentration was close to the congruent one, a large (on average) percentage of the total hops was constituted by direct hops among antisite defects.This led us to the important conclusion that in standard conditions, i.e., a congruent material at room temperature (and below), the dominating polaronic species was the bound polaron, hopping directly among antisite defects.This situation was completely different with respect to the free polaron case, because bound polarons hopped on a disordered lattice, for which standard diffusion equations may not be valid [16].Two important remarks are necessary: (i) the fact that a large percentage of hops was performed as a bound polaron does not necessarily mean that charge transport was ruled by bound polarons: free polarons were so much more mobile than bound ones that, even if they were scarce, they could carry charge in a very efficient way; (ii) even if some hops were performed as free polarons, their concentration remained small, because they performed their hops very quickly before being recaptured very soon either by a Nb 5+  Li shallow trap or by a Fe 3+ deep trap.Therefore, the rate at which they were produced was much slower than the rate at which they disappeared, so that their equilibrium concentration remained very limited.This was confirmed by our experimental TAS results showing that the signal at 1310 nm remained small or close to zero, especially at a low temperature.The influence of iron traps on Nb 4+ Li polarons is shown in Figure 5 (right).In this case, the contribution of free polaron can be disregarded because this analysis was performed at a low temperature for congruent composition, to reduce the computational cost of this simulation.Here, the situation was even more dramatic if compared to the free polarons case, because at a low temperature, even with a relatively low concentration of deep traps, the bound polaron could perform at most one single jump.The data shown in Figure 5 (right) may be taken as an overestimation because of the so-called "soft pair" effect.In other words, the Nb Li polaron can perform in some cases a large number of hops back and forth between two close-by antisite defects, without contributing effectively to the transport.In any case, this graph can be used to look for the conditions for which the number of hops becomes smaller than one.This corresponds to the situation for which at a low temperature, the bound polarons were directly trapped by a Fe 3+ , and the transport process was no longer based on Nb 4+  Li → Nb 4+ Li hops, but on Nb 4+ Li → Fe 3+ ones.This is in qualitative agreement with the results of Figure 2, which show the gradual change in the activation energy of the average hopping frequency by cooling down the sample.The measured activation energy (E A = (0.08 ± 0.01) eV) at a low temperature in sample Fe:LN 0.02 mol % corresponded to a mix of hopping processes of the type Nb 4+  Li → Nb 4+ Li and Nb 4+ Li → Fe 3+ , with an apparent activation energy in between the ones of the corresponding processes reported in Equation (4) (U GP−GP = 0.375 eV and U GP−Fe = 0.038 eV, respectively).

Conclusions
In this work, we investigated how the relative weight of the different hopping processes for small polarons in Fe:LN was influenced by sample composition and temperature.A change in the hopping processes was evidenced experimentally by measuring the activation energy of the photo-generated electron polarons' lifetime for different trap concentrations and different temperatures.The activation energy passes from about 0.42 eV at room temperature to 0.08 eV at a low temperature.By a numerical analysis based on the Marcus-Holstein theory, we were able to ascribe the transport to a mix of hop Nb 4+  Li → Nb 4+ Li (GP-GP) and Nb 4+ Nb → Nb 4+ Nb (FP-FP) types.At high temperatures and low antisite concentrations, the number of hops performed as free polarons was dominant, so that the effective mobility approached the ideal "free polaron" case described by Equation (1).In standard conditions, i.e., at room temperature and congruent composition, the "antisite" hopping GP-GP occurred more frequently than free polaron hopping FP-FP, and by cooling down the sample, it became preponderant.Still, the role of free polaron conduction cannot be neglected: because of their much higher mobility, even a scarce presence of free polarons may significantly affect the conductivity.Finally, at very low temperatures (below 150 K) and for the range of compositions herein considered, one entered another regime in which the polarons typically performed only one or two hops before being captured by a deep Fe 3+ trap.In these conditions, the activation energy of the transport process was close to the one of the trapping event, which was U GP−Fe = 0.038 eV.We point out that in recent works [17,18], a useful scheme, based on the idea of associating a finite volume with deep and shallow traps, has been proposed that allows for a semi-quantitative description of our findings without resorting to numerical simulations.
Our results confirmed that in real materials, the situation was much more intricate than in the ideal case, and that numerical studies can effectively be used to investigate the role of the different polaronic species.As a future perspective, we plan to apply this approach to the practically important topic of transport processes under continuous illumination.As shown here, the effective polaron mobility resulted in general from a mix between free and bound polarons contributions.While for the former, a standard theoretical description is available, for the latter, the situation is more complicated, because they can be considered as particles hopping on the disordered lattice constituted by antisite defects, for which standard diffusion equations may not be valid.This may ultimately lead to a revision of the theoretical description of polaron conduction in lithium niobate and related materials.

Figure 1 .
Figure 1.Typical experimental results of TA spectroscopy on the sample Fe:LN 0.1 mol %.(a) Room temperature; (b) Low temperature.

Figure 2 .
Figure 2. Arrhenius plot of the mean decay time of the signal at 785 nm, calculated using the parameters of the Kohlrausch, Williams and Watts (KWW) fit of the Transient Absorption Spectroscopy (TAS) data.The plots evidence a change in the transport mechanism around 200 K.The high temperature regime is described by an activation energy of [U ∼ 0.42] eV and is well matched by numerical simulations in sample Fe:LN 0.02 mol %.

Figure 4 .
Figure 4. (Left) Fraction of hops as a free polaron N FP−FP /N tot as a function of temperature and [Nb Li ] concentration for an undoped sample ([Fe] = 0). (Right) Total number of hops N FP−FP as a function of the temperature and [Fe] concentration for a stoichiometric sample ([Nb Li ] = 0).

Figure 5 .
Figure 5. (Left) N GP−GP /N tot as a function of temperature and antisite concentration for an undoped sample ([Fe] = 0). (Right) N GP−GP as a function of the temperature and [Fe] concentration for a congruent sample ([Nb Li ] = 19 × 10 25 m −3 ).Free polaron contribution has been disregarded for simplicity.Note the different temperature range in the two plots.

Table 1 .
Defect concentrations for the samples under study determined by optical absorption (see the text). (×10