Unveiling the Defect Structure of Lithium Niobate with Nuclear Methods

X-ray and neutron diffraction studies succeeded in the 1960s to determine the principal structural properties of congruent lithium niobate. However, the nature of the intrinsic defects related to the non-stoichiometry of this material remained an object of controversial discussion. In addition, the incorporation mechanism for dopants in the crystal lattice, showing a solubility range from about 0.1 mol% for rare earths to 9 mol% for some elements (e.g., Ti and Mg), stayed unresolved. Various different models for the formation of these defect structures were developed and required experimental verification. In this paper, we review the outstanding role of nuclear physics based methods in the process of unveiling the kind of intrinsic defects formed in congruent lithium niobate and the rules governing the incorporation of dopants. Complementary results in the isostructural compound lithium tantalate are reviewed for the case of the ferroelectric-paraelectric phase transition. We focus especially on the use of ion beam analysis under channeling conditions for the direct determination of dopant lattice sites and intrinsic defects and on Perturbed Angular Correlation measurements probing the local environment of dopants in the host lattice yielding independent and complementary information.


Introduction
Present opto-electronic and future quantum-process based communication technologies rely strongly on optical materials with their inherent or artificially modified properties. Lithium niobate (LiNbO 3 , henceforth abbreviated LN) is one of the materials that has drawn most attention in this process. Its excellent linear and non-linear optical properties, facility to form waveguide using a variety of methods, ability to incorporate optically active dopants and the applications based on periodically poled crystals make it a favorite target for basic and applied research [1][2][3][4].
In order to achieve the modification of its properties with a certain goal in mind, the basic properties of LN have to be understood. Among them, the intrinsic defect structure of the congruent lithium-deficient crystals in comparison to stoichiometric ones, phenomena related with the phase ferroelectric-paraelectric phase transition, incorporation of dopants and their interaction with the lattice and lattice damage production and removal are topics for fundamental research fostering future application directed research.
Methods based on nuclear techniques played and will play an important role in materials research. Their application can clarify or-in combination with other methodsassist to clarify the topics mentioned in the preceding paragraph. While ion beam analysis methods are ideal for research on composition as well as lattice defect concentration and distribution, hyperfine interaction methods can yield complementary information on local environments in the lattice and dopant charge states. As will be demonstrated in the course of this review the combination of several of these techniques can unveil the nature of

Short Overview on Nuclear Methods for Materials Research
The nuclear techniques addressed in this review are characterized shortly indicating some peculiarities with regard to their application in research on LN. References to general introductory books and reviews for further in depth study are supplied for each of the methods.

Ion Beam Analysis
Ion beam analysis is a group of nuclear methods that are based on the interaction of energetic ion beams (typically few keV to tens of MeV) with materials [5][6][7][8]. Its main applications are: • Composition determination in bulk materials, surface layers and multilayer systems. Some of the methods attain a depth resolution down to sub-monolayer range; • Detection and quantitation of lattice damage in monocrystalline materials with high depth resolution up to a few μm from the surface; • Characterization of intrinsic and extrinsic defects in terms of lattice site location and defect type.
RBS relies on elastic backscattering of the incoming ions by the target atoms. Its origins can be traced back to the 1909 experiments performed by Geiger and Marsden [17] which led to the first modern model of the atom. The energy spectra of backscattered light ions (protons, He + and He ++ ) in the MeV range provide information on the composition of the material and its depth dependence up to several μm into the bulk. Since the

Short Overview on Nuclear Methods for Materials Research
The nuclear techniques addressed in this review are characterized shortly indicating some peculiarities with regard to their application in research on LN. References to general introductory books and reviews for further in depth study are supplied for each of the methods.

Ion Beam Analysis
Ion beam analysis is a group of nuclear methods that are based on the interaction of energetic ion beams (typically few keV to tens of MeV) with materials [5][6][7][8]. Its main applications are: • Composition determination in bulk materials, surface layers and multilayer systems. Some of the methods attain a depth resolution down to sub-monolayer range; • Detection and quantitation of lattice damage in monocrystalline materials with high depth resolution up to a few µm from the surface; • Characterization of intrinsic and extrinsic defects in terms of lattice site location and defect type.

Figure 2.
Backscattering spectrum (3.25 MeV He + ) for congruent LN with Rutherford Backscattering for niobium and Resonant Elastic Backscattering for oxygen using the 16 O(α,α) 16 O resonance (about 30 times higher than Rutherford scattering) at 3.045 MeV.
Another possibility for the detection and quantification of light elements in heavy matrices is to look for the light charged nuclei produced by ion beam induced reactions in the target material, within the so called NRA technique. Long measuring times are frequently necessary, since the reaction cross sections are typically much smaller than for backscattering. The analysis of the spectra is analogous to that of RBS and EBS. As an example, Figure 3 depicts a spectrum for LN irradiated with 1.75 MeV protons. Contributions of backscattering by the target and by the alpha particles produced in the 7 Li(p,α) 4 He reaction are clearly distinguishable due to the high Q-value (17.35 MeV) of the reaction employed.  (3.25 MeV He + ) for congruent LN with Rutherford Backscattering for niobium and Resonant Elastic Backscattering for oxygen using the 16 O(α,α) 16 O resonance (about 30 times higher than Rutherford scattering) at 3.045 MeV.
Another possibility for the detection and quantification of light elements in heavy matrices is to look for the light charged nuclei produced by ion beam induced reactions in the target material, within the so called NRA technique. Long measuring times are frequently necessary, since the reaction cross sections are typically much smaller than for backscattering. The analysis of the spectra is analogous to that of RBS and EBS. As an example, Figure 3 depicts a spectrum for LN irradiated with 1.75 MeV protons. Contributions of backscattering by the target and by the alpha particles produced in the 7 Li(p,α) 4 He reaction are clearly distinguishable due to the high Q-value (17.35 MeV) of the reaction employed.  The PIXE technique relies on the production of characteristic X-rays (mainly K-and L-lines) by interaction of the incoming ion (in most cases protons, for special purposes He + ) with inner shell electrons of the target atoms. It is most useful for the detection of medium and heavy elements in compound materials. The quantitative analysis of the elemental concentrations is significantly more complicated than in the previous cases since the technique integrates over the whole ion path in the material. This lack of depth The PIXE technique relies on the production of characteristic X-rays (mainly K-and L-lines) by interaction of the incoming ion (in most cases protons, for special purposes He + ) Crystals 2021, 11, 501 4 of 55 with inner shell electrons of the target atoms. It is most useful for the detection of medium and heavy elements in compound materials. The quantitative analysis of the elemental concentrations is significantly more complicated than in the previous cases since the technique integrates over the whole ion path in the material. This lack of depth resolution can be partially overcome by measuring the sample at different angles of incidence. Figure 4 shows an example of an X-ray energy spectrum for LN irradiated with 3.1 MeV He + ions. The PIXE technique relies on the production of characteristic X-rays (mainly K-and L-lines) by interaction of the incoming ion (in most cases protons, for special purposes He + ) with inner shell electrons of the target atoms. It is most useful for the detection of medium and heavy elements in compound materials. The quantitative analysis of the elemental concentrations is significantly more complicated than in the previous cases since the technique integrates over the whole ion path in the material. This lack of depth resolution can be partially overcome by measuring the sample at different angles of incidence. Figure 4 shows an example of an X-ray energy spectrum for LN irradiated with 3.1 MeV He + ions. Ions incident along major crystallographic directions (axes or planes) of monocrystalline materials show the channeling effect. Along a crystallographic axis or plane, the atoms of the crystal atoms appear perfectly aligned in strings or plates, respectively. These form, due to their repellent Coulomb potential, a "channel" in which the ions are guided and the ion fluence is concentrated in the center of the channel. Due to this guiding effect the probability for close encounters with lattice atoms are reduced diminishing the yield for ion scattering, nuclear reactions and X-ray emission.
In a typical channeling experiment the yield of backscattered ions, nuclear reactions products or characteristic X-rays is measured tilting the crystal stepwise from a major crystallographic direction to one of random ion beam incidence (in which the atomic arrangement appears disordered). These "scans" are essentially characterized by two parameters, depicted in Figure 5: Ions incident along major crystallographic directions (axes or planes) of monocrystalline materials show the channeling effect. Along a crystallographic axis or plane, the atoms of the crystal atoms appear perfectly aligned in strings or plates, respectively. These form, due to their repellent Coulomb potential, a "channel" in which the ions are guided and the ion fluence is concentrated in the center of the channel. Due to this guiding effect the probability for close encounters with lattice atoms are reduced diminishing the yield for ion scattering, nuclear reactions and X-ray emission.
In a typical channeling experiment the yield of backscattered ions, nuclear reactions products or characteristic X-rays is measured tilting the crystal stepwise from a major crystallographic direction to one of random ion beam incidence (in which the atomic arrangement appears disordered). These "scans" are essentially characterized by two parameters, depicted in Figure 5:

•
The minimum yield, χ min , i.e., the ratio of the reaction yield at perfect alignment with a crystallographic direction to that at random incidence. The lower χ min the better is the crystalline quality of the sample.

•
The half-width of the scan, ψ 1/2 , that is larger the less defective the material is.
Crystals 2021, 11, x FOR PEER REVIEW 5 of 57 • The minimum yield, χmin, i.e., the ratio of the reaction yield at perfect alignment with a crystallographic direction to that at random incidence. The lower χmin the better is the crystalline quality of the sample.

•
The half-width of the scan, ψ1/2, that is larger the less defective the material is. Possible applications of the channeling effect are [20]: • Determination of intrinsic lattice defect concentration (typically, genuine defects of the material or radiation induced defects) with high depth resolution. For the analysis, spectra are taken for the undamaged and the damaged sample aligned with a major crystallographic axis and for a random ion incidence angle. From their comparison a depth-resolved defect profile can be derived using appropriate • Determination of intrinsic lattice defect concentration (typically, genuine defects of the material or radiation induced defects) with high depth resolution. For the analysis, spectra are taken for the undamaged and the damaged sample aligned with a major crystallographic axis and for a random ion incidence angle. From their comparison a depth-resolved defect profile can be derived using appropriate calculation methods or programs, e.g., DICADA [21]; • Determination of intrinsic lattice defect type. Spectra are recorded in the same manner as in the item above but for different ion energies yielding a dependence of the minimum yield difference on E κ , where E is the ion energy. Determining the constant κ allows to distinguish e.g., between point defects, dislocations, inclusions, etc. • Lattice location of dopants in the crystal. While an angular scan for a dopant located on a regular lattice site of the crystal will follow that one of the host atoms, foreign atoms on interstitial sites will show scans with distinctly different shapes, depending on their relative position in the channel (e.g., a peak instead of a dip if located in the center). Measurements for different crystallographic directions allow to determine the lattice site by triangulation.
The channeling method was initially developed for the study of simple mono-or bi-elemental cubic crystals. Therefore, an extension to materials as complex as LN was not straightforward and it was only attempted in the 1980s. In contrast to simple materials, in most axial directions of LN the channels are highly asymmetric or split in several subchannels. In all cases they are dominated by the heaviest constituent (Nb). Figure 6 shows the projections of five important axial and a planar directions in LN indicating the regular and some interstitial lattice sites together with the respective calculated continuum potentials governing the channels.
The evaluation of the effect of a dopant on a non-Nb site-regular or interstitial-by simple reasoning can lead to serious errors. Therefore, the use of simulation software is mandatory for the lattice site determination of intrinsic and extrinsic defects in LN. Results presented and discussed in this review are based on two channeling simulation programs. FLUX [22,23] is a long-known code for RBS/C and NRA/C in cubic crystals that has been adapted to the trigonal LN lattice. CASSIS [24][25][26][27] is a multi-purpose code, i.e., it is applicable for all Bravais lattices and can also deal realistically with PIXE/C and high dopant concentrations. Figure 7 shows results of CASSIS for two-dimensional RBS/C calculations for the y-axis.
The SIMS method uses low energy heavy ions (mostly Ar + ) to remove and ionize surface atoms by sputtering and analyzes these secondary ions using a mass spectrometer. This allows determining the composition of materials with very good depth resolution. The sensitivity of the method-ranging from ppm to ppb depending on the material-is its main advantage. However, this method is destructive, in contrast to those mentioned above.
ERD uses typically He-ions with MeV energies to remove hydrogen from a sample surface. An absorber foil in front of the detection stops the backscattered He-ions and allows pass only the protons for analysis. In contrast, HIERD uses heavy ions with specific energy of several MeV/amu in order to sputter even heavy atoms from the surface. Various methods exist for the detection and analysis of the secondary ions emitted.
For better readability, technical details on the ion beam analysis such as ion species, charge state, current, detector geometry, installation used and provenience of samples are mostly omitted in this review. However, the information can always be recovered from the cited works. Figure 6. Channel projections with regular and some interstitial sites for important crystallographic directions and the corresponding calculated continuum potentials. The two-dimensional potentials for axial directions are plotted in logarithmic scale; the channeling occurs in the blue regions. The projection of the (0001)-plane and the corresponding one-dimensional continuum potential can be seen on the right side down.
The evaluation of the effect of a dopant on a non-Nb site-regular or interstitial-by simple reasoning can lead to serious errors. Therefore, the use of simulation software is mandatory for the lattice site determination of intrinsic and extrinsic defects in LN. Results presented and discussed in this review are based on two channeling simulation programs. FLUX [22,23] is a long-known code for RBS/C and NRA/C in cubic crystals that has been adapted to the trigonal LN lattice. CASSIS [24][25][26][27] is a multi-purpose code, i.e., it is applicable for all Bravais lattices and can also deal realistically with PIXE/C and high . Channel projections with regular and some interstitial sites for important crystallographic directions and the corresponding calculated continuum potentials. The two-dimensional potentials for axial directions are plotted in logarithmic scale; the channeling occurs in the blue regions. The projection of the (0001)-plane and the corresponding one-dimensional continuum potential can be seen on the right side down. dopant concentrations. Figure 7 shows results of CASSIS for two-dimensional RBS/C calculations for the y-axis. The SIMS method uses low energy heavy ions (mostly Ar + ) to remove and ionize surface atoms by sputtering and analyzes these secondary ions using a mass spectrometer. This allows determining the composition of materials with very good depth resolution. The sensitivity of the method-ranging from ppm to ppb depending on the material-is its main advantage. However, this method is destructive, in contrast to those mentioned above.
ERD uses typically He-ions with MeV energies to remove hydrogen from a sample surface. An absorber foil in front of the detection stops the backscattered He-ions and allows pass only the protons for analysis. In contrast, HIERD uses heavy ions with specific energy of several MeV/amu in order to sputter even heavy atoms from the surface. Various methods exist for the detection and analysis of the secondary ions emitted.
For better readability, technical details on the ion beam analysis such as ion species, charge state, current, detector geometry, installation used and provenience of samples are mostly omitted in this review. However, the information can always be recovered from the cited works. . RBS/C calculations with CASSIS. Upper row: Two-dimensional scans for the 0110 -axis, ψ 1 corresponds to a rotation around the c-axis and ψ 2 around the a-axis of LN. The calculated yields for Li, Nb and O are shown. Li and O exhibit a yield enhancement for the (0001)-or the (1120)-plane respectively. Lower row: One-dimensional scans through 0110 close to (1120), in which Li shows the typical asymmetric yield curve, and through the (0001)-plane, in which Li acts similar to a host atom located in the center of the planar channel.

Hyperfine Interaction Methods
Hyperfine interaction methods such as Perturbed Angular Correlations (PAC), Mössbauer Effect (ME), Nuclear Magnetic Resonance (NMR) and Nuclear Quadrupole Resonance (NQR) are based on the interaction of a nucleus with an electric (or magnetic) field produced in a solid by the electrons and other nuclei near the nucleus. In several cases, techniques from the ion beam analysis and hyperfine interaction groups are used to study the same samples, obtaining complementary information and enlarging significantly the range of sensitivity for characterization of defects [28].
The PAC probe is an unstable nucleus decaying through a γ-γ cascade that has an intermediate state with suitable half-life and nuclear quadrupole moment (or dipole moment, in case one is interested in magnetic interactions). Conversion electrons from one of the transitions of the γ-γ cascade can also be used [29]. In the case of Hf-doped LN, the PAC probe is usually produced by the neutron capture reaction 180 Hf(n,γ) 181 Hf on natural Hf introduced in the samples during growth as HfO 2 . A typical activity of ≈15 µCi is obtained after reactor irradiation. Although thermal neutrons are the most important in the reactor irradiation, the LN crystal is also irradiated with fast neutrons, at a fluence 1-2 orders of magnitude smaller than the thermal flux [30] and these higher energy neutrons produce defects distinct from the ones produced by the thermal neutrons. Annealing in air at 700 • C is enough to restore the lattice completely [31]. The 111 In probe was implanted to a typical fluence of ≈10 13 cm −2 from which more than 90% were due to an unavoidable stable Cd contamination, resulting in a typical activity of 50 µCi. The crystals were annealed at temperatures in the 700-900 • C range under flowing wet oxygen to remove the defects introduced by ion implantation. PAC measurements were made using standard spectrometers with four BaF 2 detectors in a plane at 90 • intervals for γ-γ PAC or two BaF 2 detectors with two magnetic lenses in a plane at 90 • intervals for e --γ PAC [29].
Equation (1) is used to calculate the time differential anisotropy R(t) from the coincidence spectra N(θ,t), where θ is the angle between detectors and t is the time delay between events.
G 22 (t) is the perturbation function which describes the modulation of the angular correlation and A 22 is a parameter which depends on the nature of the radiations of the cascade [32]. For a cascade having an intermediate state with spin I = 5/2, such as the ones from the decay of 111 In and 181 Hf, three frequencies ω n = C n (η)ν Q are observed, corresponding to: The quadrupole interaction frequency (QIF), ν Q = eQV zz /h, and the asymmetry parameter η = (V xx − V yy )/V zz , which can be deduced from the frequency factors C n (η), contain information about the magnitude of the principal component V zz and the asymmetry of the EFG, respectively. The exponential factor allows for a Lorentzian distribution of EFGs around a mean value. Such a distribution can be caused, e.g., by different lattice defects in the vicinity of the probe atoms. The S 2n coefficients, which give the amplitudes for each frequency, can be calculated for polycrystalline as well as single crystalline samples. If fractions f i of the total number of radioactive probes are in different but unique sites experiencing different EFGs the observed spectra is described by a sum over several perturbation functions, f i G i 22 (t), with the sum of the f i normalized to unity. Figure 8 shows the PAC coincidence count rates N(θ,t) obtained with detectors at angles of θ = 90 • and θ = 180 • and the corresponding spectrum obtained in near-stoichiometric LN (doped with 4.6 wt % K 2 O) implanted with 111 In. The coincidence count rates N(θ,t) show the half-life of the intermediate state of the γ-γ cascade from the decay of 111 In (T 1/2 = 80 ns) with the perturbation resulting from the interaction of the quadrupole moment of the intermediate state of the cascade and the EFG produced by the surrounding electric charges. The anisotropy ratio R(t) (usually designated "PAC spectrum") calculated according to eq. 1 extracts the perturbation, exploiting the difference between minima and maxima in the experimental count rates at 90 • and 180 • .
Details on the other hyperfine interaction techniques can be found, e.g., in references [8,32]. In general one can extract information on the principal component V zz and on the asymmetry of the EFG, although the observable frequencies are dependent on the particular nuclear method and on the nuclear properties of the probe nucleus. One must emphasize that the knowledge of the EFG is not yet enough for a direct assignment of the probe to a given lattice site. Although much progress has been made in the calculation of EFG [33] this is still time-consuming and far from generalized. Nevertheless, one can compare EFG values obtained with different probes in LN [34]. of 111 In (T1/2 = 80 ns) with the perturbation resulting from the interaction of the quadrupole moment of the intermediate state of the cascade and the EFG produced by the surrounding electric charges. The anisotropy ratio R(t) (usually designated "PAC spectrum") calculated according to eq. 1 extracts the perturbation, exploiting the difference between minima and maxima in the experimental count rates at 90° and 180°. Details on the other hyperfine interaction techniques can be found, e.g., in references [8,32]. In general one can extract information on the principal component Vzz and on the asymmetry of the EFG, although the observable frequencies are dependent on the particular nuclear method and on the nuclear properties of the probe nucleus. One must emphasize that the knowledge of the EFG is not yet enough for a direct assignment of the probe to a given lattice site. Although much progress has been made in the calculation of EFG [33] this is still time-consuming and far from generalized. Nevertheless, one can compare EFG values obtained with different probes in LN [34].

Congruent Versus Stoichiometric Lithium Niobate
Already in the first growth studies it was noted that LN does not melt congruently with the first phase diagram published by Reisman and Holtzberg [35] in 1958. This nonstoichiometry implies and intrinsic defect structure and influences strongly many material properties which triggered intense research activities on the nature of the defects and the production of stoichiometric material.
Initial NMR measurements with 7 Li and 93 Nb in LN and lithium tantalate (LiTaO3, henceforth abbreviated LT) were reported in the late 1960s [36][37][38][39]. Measurements were performed with the then available magnetic fields below 1 T [40]. These measurements were characterized by wide resonance lines, whose width could be empirically related to the material's departure from stoichiometry. Measurements with NQR [41,42] and ME [43] have also shown broader than expected interaction lines in LN. Later, Peterson et al. [44] found a second 93 Nb NMR signal in congruent LN crystals, with a near zero QIF and

Congruent Versus Stoichiometric Lithium Niobate
Already in the first growth studies it was noted that LN does not melt congruently with the first phase diagram published by Reisman and Holtzberg [35] in 1958. This nonstoichiometry implies and intrinsic defect structure and influences strongly many material properties which triggered intense research activities on the nature of the defects and the production of stoichiometric material.
Initial NMR measurements with 7 Li and 93 Nb in LN and lithium tantalate (LiTaO 3 , henceforth abbreviated LT) were reported in the late 1960s [36][37][38][39]. Measurements were performed with the then available magnetic fields below 1 T [40]. These measurements were characterized by wide resonance lines, whose width could be empirically related to the material's departure from stoichiometry. Measurements with NQR [41,42] and ME [43] have also shown broader than expected interaction lines in LN. Later, Peterson et al. [44] found a second 93 Nb NMR signal in congruent LN crystals, with a near zero QIF and an intensity of approximately 6% of the total Nb signal that was assigned to Nb in Li sites (5.9%) thus providing a strong support to the Nb vacancy model, even if the authors assumed later that the interpretation of their NMR spectra was not so straightforward [45]. More recent measurements by Blümel et al. [40] using a higher magnetic field of ≈7 T did not confirm a second 93 Nb signal with such high intensity, thus favoring instead the Li vacancy model. A series of experimental studies complemented with simulations assuming O-, Li-or Nb-vacancies (but unfortunately omitting models based on ilmenite stacking faults) was published by Yatsenko et al. [46][47][48], favoring the formation of Li-vacancies in the vicinity of Nb located on regular Li sites.
Rebouta et al. [49] reported RBS/C and 181 Hf PAC measurements in LN doped with 1 mol% HfO 2 in the melt. The RBS/C measurements unambiguously determined that Hf replaced Li, but the obtained PAC spectra could only be fitted considering two QIF around 1200 MHz. Two frequencies of the same magnitude were also observed in 181 Hf PAC measurements done in polycrystalline LN samples [50]. These QIF were later detailed as ν Q1 = 1154(12) MHz and ν Q2 = 1213(12) MHz, with asymmetry parameters in the 0.2-0.3 range, from new measurements in single-crystals [34]. NMR studies [51,52] confirm an asymmetry parameter of 0.2 for congruent material. The observed frequency doublet was unexpected, as Hf was found only replacing Li and a single QIF should suffice. Significant progress was made when higher precision PAC measurements were performed using the 111 In probe, also replacing Li in congruent LN crystals [53]. Again, the new spectra could not be described by a single frequency even if considering an asymmetric EFG or a broad frequency distribution, both already indicative of crystal imperfections [54]. The frequency doublet now found was ν Q1 = 192(2) MHz and ν Q2 = 205(2) MHz, with asymmetry parameters in the 0.1-0.2 range. Since the frequencies were again very close, the simplest explanation is that the probes were in slightly different Li sites, one of them possibly linked to a defect located in the Li-Octahedron. Since the fraction of PAC probe atoms located at this additional site is significant (about one third) it should be energetically favorable for the inclusion of atoms causing a charge excess relative to Li. It is therefore unlikely that this defect is identical to a Li-position distorted by a nearby excess Nb atom, since the occupation of such a site by a dopant with a higher charge state than Li would increase further the local charge imbalance. A far more reasonable assumption is to identify this defect with a site, shifted form the regular lithium site that actually incorporates the excess Nb (or dopants) in congruent material. The broadening of quadrupole interaction lines seen by other techniques is also compatible with the existence of two nearby frequencies that cannot be resolved.
Growth procedures developed in the 1990s [55][56][57][58] made available high-quality crystals with stoichiometric composition, allowing a direct comparison of the two compositions.
The possible impacts of the Li-and the Nb-vacancy models on axial and planar channeling can be easily estimated qualitatively from the channel projections presented above. Eventually existing vacancies on regular Nb-sites should have only small effects on the Nb-channeling since even in the case of 5.9 mol% vacancies the steering effect is hardly diminished. On the other hand, anti-site Nb-atoms located on regular Li-sites should lower the minimum yields for Li-NRA measurements for congruent material in axial directions clearly since the channeling effect on the Li-rows is enhanced by the Nb. For a 1.2% fraction the average atomic number of the row increases from 3 (stoichiometric) to 3.5 (Nb-vacancy model). For planar channeling in the (0001)-direction Li-NRA-scans should show clear peaks in the aligned direction; but for the congruent case, it should be lower and narrower than for the stoichiometric.
Measurements with congruent and near-stoichiometric LN, grown with 4.6 wt % K 2 O in the melt ([Li]/[Nb] = 0.979 [59]), reported by Kling et al. [60,61] showed contrary effects. NRA/C using the 7 Li(p,α) 4 He nuclear reaction and RBS/C scans were performed simultaneously in order to study differences in Li-and Nb-sublattice between the two crystals. Several crystallographic axes and planes were investigated, as shown in Figure 9. The Li-NRA data corresponding to the (0001) plane, perpendicular to the c-axis, were particularly interesting. They showed a prominent flux peak that is considerably broadened in the congruent crystal in comparison to the near-stoichiometric one. The NRA-Li broadening suggests either a disordered Li sublattice or several well-defined defect sites for the Li ions. For the 0221 -axis, in which Li is located in the center of the channel, the minimum yield for Li-NRA/C drops drastically ruling out the possibility of simple Nb Li antisites, i.e., the Li-and the Nb-vacancy model [60] Besides the stronger reduction of the 7 Li(p,α) 4 He reaction yield at the surface the spectra for the stoichiometric extend also to higher depth indicating a reduced stopping power due to a diminished number of lattice defects in the channel compared to congruent material [61].
Further, a strong broadening of Nb-RBS/C in the 0110 -axis of the near-stoichiometric sample has been observed indicating an improvement of the Nb-sublattice structural ordering with regard to the congruent case [60]. These observations indicate that the excess Nb cannot be located in regular Li lattice sites, favoring the ilmenite stacking fault model but not ruling out the formation of small non-crystalline or Li-deficient (e.g., ilmenite-type LN or monoclinic lithium LiNb 3 O 8 ) regions in the material and on the other extreme, point defects with inversed cation occupation. particularly interesting. They showed a prominent flux peak that is considerably broadened in the congruent crystal in comparison to the near-stoichiometric one. The NRA-Li broadening suggests either a disordered Li sublattice or several well-defined defect sites for the Li ions. For the 〈022 1〉-axis, in which Li is located in the center of the channel, the minimum yield for Li-NRA/C drops drastically ruling out the possibility of simple NbLi antisites, i.e., the Li-and the Nb-vacancy model [60] Besides the stronger reduction of the 7 Li(p,α) 4 He reaction yield at the surface the spectra for the stoichiometric extend also to higher depth indicating a reduced stopping power due to a diminished number of lattice defects in the channel compared to congruent material [61]. The Li data are plotted as red circles and were obtained using NRA via the 7 Li(p,α) 4 He reaction and the Nb data are plotted as blue circles and were obtained using Rutherford Backscattering. The uncertainty bars for Nb-RBS are smaller than symbols. Adapted from Kling et al. [59,61].
Further, a strong broadening of Nb-RBS/C in the 〈011 0〉 -axis of the nearstoichiometric sample has been observed indicating an improvement of the Nb-sublattice structural ordering with regard to the congruent case [60]. These observations indicate that the excess Nb cannot be located in regular Li lattice sites, favoring the ilmenite stacking fault model but not ruling out the formation of small non-crystalline or Lideficient (e.g., ilmenite-type LN or monoclinic lithium LiNb3O8) regions in the material and on the other extreme, point defects with inversed cation occupation.
This point can be clarified by studying the Nb sublattice measuring the energy dependence of axial dechanneling induced by defects along several crystallographic directions. The difference in the minimum yields between near-stoichiometric and congruent, ΔXmin, at a fixed depth interval yields information about the type and size of defects in a non-perfect crystal [62]. The observed values for κ range between 0.69 to 0.78 Figure 9. RBS/C dips obtained for the (0001) plane in (a) congruent and (b) near-stoichiometric LiNbO 3 , as well as for the 0221 axis in (c) congruent and (d) near-stoichiometric LiNbO 3 . The Li data are plotted as red circles and were obtained using NRA via the 7 Li(p,α) 4 He reaction and the Nb data are plotted as blue circles and were obtained using Rutherford Backscattering. The uncertainty bars for Nb-RBS are smaller than symbols. Adapted from Kling et al. [59,61]. This point can be clarified by studying the Nb sublattice measuring the energy dependence of axial dechanneling induced by defects along several crystallographic directions. The difference in the minimum yields between near-stoichiometric and congruent, ∆X min , at a fixed depth interval yields information about the type and size of defects in a nonperfect crystal [62]. The observed values for κ range between 0.69 to 0.78 falls between the values for true point defects (κ = 1) and extended defect clusters (κ = 0.5) which indicates defects should have a size in the order of nanometers (as expected for local stacking fault sequences). However, a small admixture of true point defects cannot be excluded. In addition, the possibility that microcrystals of ilmenite-type or LiNb 3 O 8 (κ = 0) or dislocations (κ= −0.5) exist can be rejected.
A comparison was also made with PAC. Figure 10 shows PAC spectra obtained with the 111 In probe in the same near-stoichiometric crystal grown with 4.6 wt % K 2 O after implantation and annealing [63], compared with a spectrum previously obtained in a congruent crystal [53]. It is clear that the spectrum with the highest attenuation is the one obtained in the congruent crystal. The spectrum obtained in the near-stoichiometric crystal is described by a single QIF, ν Q = 191(2) MHz, i.e., only one of the two QIF previously found, and with a zero asymmetry parameter. It can also be seen that the pattern from the congruent crystal is not in phase with the one from the near-stoichiometric crystal, with the maxima in the congruent crystal occurring progressively a few nanoseconds after the corresponding ones for the stoichiometric crystal, due to the modulation imposed by the frequency doublet, which is also responsible for a higher attenuation. As the QIF have almost identical values, the PAC measurements suggest the existence of two non-equivalent Li sites in congruent crystals and only one Li site in stoichiometric material. No conclusions can be drawn on the Nb sublattice from these measurements, as 111 In replaces Li only [53]. crystal is described by a single QIF, νQ = 191(2) MHz, i.e., only one of the two QIF previously found, and with a zero asymmetry parameter. It can also be seen that the pattern from the congruent crystal is not in phase with the one from the nearstoichiometric crystal, with the maxima in the congruent crystal occurring progressively a few nanoseconds after the corresponding ones for the stoichiometric crystal, due to the modulation imposed by the frequency doublet, which is also responsible for a higher attenuation. As the QIF have almost identical values, the PAC measurements suggest the existence of two non-equivalent Li sites in congruent crystals and only one Li site in stoichiometric material. No conclusions can be drawn on the Nb sublattice from these measurements, as 111 In replaces Li only [53]. The joint analysis of the NRA and PAC data thus suggest the existence of two distinct Li sites in a congruent crystal, while only one of them remains in a stoichiometric crystal. Natural candidates for such defects are stacking faults. Several models have been proposed for LN in the literature. Both the Li and the Nb vacancy models include Nb on Li sites and thus a sequence of three Nb consecutive atoms along the c-axis would be created. Nassau et al. [64] proposed different stacking sequences in order to avoid such improbable arrangements and Smyth [65,66] proposed short stacking sequences with the ilmenite structure, which were later found energetically plausible [67]. The stacking sequence along the c-axis of ilmenite-like LN would be Nb, Li, V, Li, Nb, …, where V denotes the intrinsic vacancy. Since these inverse stacking sequences accommodate no Nb excess themselves they have to accompanied by Li vacancies. The "inverted" Li site in the ilmenite sequence is generally assumed to be one of this vacancies leading to the "complex The joint analysis of the NRA and PAC data thus suggest the existence of two distinct Li sites in a congruent crystal, while only one of them remains in a stoichiometric crystal. Natural candidates for such defects are stacking faults. Several models have been proposed for LN in the literature. Both the Li and the Nb vacancy models include Nb on Li sites and thus a sequence of three Nb consecutive atoms along the c-axis would be created. Nassau et al. [64] proposed different stacking sequences in order to avoid such improbable arrangements and Smyth [65,66] proposed short stacking sequences with the ilmenite structure, which were later found energetically plausible [67]. The stacking sequence along the c-axis of ilmenite-like LN would be Nb, Li, V, Li, Nb, . . . , where V denotes the intrinsic vacancy. Since these inverse stacking sequences accommodate no Nb excess themselves they have to accompanied by Li vacancies. The "inverted" Li site in the ilmenite sequence is generally assumed to be one of this vacancies leading to the "complex defect" model. This defect structure can be regarded as a Li-vacancy model in which the Nb is shifted from the regular Li-site and the inconvenient stacking of three consecutive Nb atoms along the c-axis is avoided. Further, it furnishes a straightforward explanation for the occurrence of the second electric field gradient in the Li-Octahedron observed in hyperfine interaction measurements. Figure 11 compares the regular LN stacking with other sequences, including short ilmenite stacking sequences and the "complex defect".
The channel projections of the lattice sites corresponding to the "inverted" Li and Nbsites in ilmenite stacking faults are plotted in Figure 12. The Nb on inverted locations shows up close but shifted from the regular Li sites and therefore should induce dechanneling near the Li-rows instead of an enhancement of the channeling effect for regular Li-site occupation. This effect is expected to be most visible in the case of the 0221 -axis and the (0001)-plane where these sites are located close to center of the channel.
defect" model. This defect structure can be regarded as a Li-vacancy model in which the Nb is shifted from the regular Li-site and the inconvenient stacking of three consecutive Nb atoms along the c-axis is avoided. Further, it furnishes a straightforward explanation for the occurrence of the second electric field gradient in the Li-octahedron observed in hyperfine interaction measurements. Figure 11 compares the regular LN stacking with other sequences, including short ilmenite stacking sequences and the "complex defect". The channel projections of the lattice sites corresponding to the "inverted" Li and Nb-sites in ilmenite stacking faults are plotted in Figure 12. The Nb on inverted locations shows up close but shifted from the regular Li sites and therefore should induce dechanneling near the Li-rows instead of an enhancement of the channeling effect for regular Li-site occupation. This effect is expected to be most visible in the case of the 〈022 1〉 -axis and the (0001)-plane where these sites are located close to center of the channel.  The channel projections of the lattice sites corresponding to the "inverted" Li and Nb-sites in ilmenite stacking faults are plotted in Figure 12. The Nb on inverted locations shows up close but shifted from the regular Li sites and therefore should induce dechanneling near the Li-rows instead of an enhancement of the channeling effect for regular Li-site occupation. This effect is expected to be most visible in the case of the 〈022 1〉 -axis and the (0001)-plane where these sites are located close to center of the channel. These results can be related to the lattice structure looking at the interaction potential in the (0001) plane seen by the incoming ions, shown in Figure 13, which was calculated using Molière's approximation to the Thomas-Fermi potential under the assumption of a static lattice [68]. Two pronounced minima are visible, showing the existence of two channels for the incoming ions, separated by a potential barrier associated to O ions. One of those minima, designated as Li-I, corresponds to the regular Li position, whereas the other one, named as Li-II, lies close to a regular Nb location. It could be considered as a Li that substitutes for an Nb and then shifts towards the nearest oxygen plane up to a position symmetrical to that of Li. The existence of those defective Li sites is a key result that provides a confirmation of theoretical analyses [67] proposing the occurrence of inverted Li-Nb sequences as they appear in the related ilmenite structure. channels for the incoming ions, separated by a potential barrier associated to O ions. One of those minima, designated as Li-I, corresponds to the regular Li position, whereas the other one, named as Li-II, lies close to a regular Nb location. It could be considered as a Li that substitutes for an Nb and then shifts towards the nearest oxygen plane up to a position symmetrical to that of Li. The existence of those defective Li sites is a key result that provides a confirmation of theoretical analyses [67] proposing the occurrence of inverted Li-Nb sequences as they appear in the related ilmenite structure. Taking as a basis this local ilmenite structure, the broadened flux peak in the NRA/channeling spectra has been simulated using the CASSIS code. The peak shape and particularly the broadening of the Li-NRA signal is reasonably well explained. The best fit is obtained assuming 2% of dispersed ilmenite stacking faults, Li-Nb-Nb-Li. On the contrary, other models exclusively involving Nb antisites and cation vacancies cannot account for this key experimental feature, as shown in Figure 14. Taking as a basis this local ilmenite structure, the broadened flux peak in the NRA/ channeling spectra has been simulated using the CASSIS code. The peak shape and particularly the broadening of the Li-NRA signal is reasonably well explained. The best fit is obtained assuming 2% of dispersed ilmenite stacking faults, Li-Nb-Nb-Li. On the contrary, other models exclusively involving Nb antisites and cation vacancies cannot account for this key experimental feature, as shown in Figure 14. The energy dependence of the axial dechanneling reported in [60] has also been simulated using CASSIS. Figure 15 shows a comparison of the experimentally observed energy dependence of the minimum yield Xmin in several axes for protons together with CASSIS simulations for (a) a defect free stoichiometric crystal; (b) congruent crystal containing 5.9 and 1.2% Nb antisites and (c) congruent cristal containing 2.0% ilmenite stacking inversions with and without 1.2% Nb antisites, from Kling at al. [26]. The best fit to the Nb data is obtained when the ilmenite stacking inversions with 1.2% Nb antisites, i.e., within the Li vacancy model.  (0001) Li-NRA angular scan assuming different defect models for congruent lithium niobate: (a) 5.9% niobium antisites, (b) 1.2% niobium antisites, and (c) 2% ilmenite type stacking faults. The Li data are plotted as red circles and were obtained using NRA via the 7 Li(p,α) 4 He reaction and the Nb data are plotted as blue circles and were obtained using Rutherford Backscattering. The uncertainty bars for Nb-RBS are smaller than symbols. Redrawn from Kling et al. [61].
The energy dependence of the axial dechanneling reported in [60] has also been simulated using CASSIS. Figure 15 shows a comparison of the experimentally observed energy dependence of the minimum yield X min in several axes for protons together with CASSIS simulations for (a) a defect free stoichiometric crystal; (b) congruent crystal containing 5.9 and 1.2% Nb antisites and (c) congruent cristal containing 2.0% ilmenite stacking inversions with and without 1.2% Nb antisites, from Kling at al. [26]. The best fit to the Nb data is obtained when the ilmenite stacking inversions with 1.2% Nb antisites, i.e., within the Li vacancy model. are smaller than symbols. Redrawn from Kling et al. [61].
The energy dependence of the axial dechanneling reported in [60] has also been simulated using CASSIS. Figure 15 shows a comparison of the experimentally observed energy dependence of the minimum yield Xmin in several axes for protons together with CASSIS simulations for (a) a defect free stoichiometric crystal; (b) congruent crystal containing 5.9 and 1.2% Nb antisites and (c) congruent cristal containing 2.0% ilmenite stacking inversions with and without 1.2% Nb antisites, from Kling at al. [26]. The best fit to the Nb data is obtained when the ilmenite stacking inversions with 1.2% Nb antisites, i.e., within the Li vacancy model.  The crystal structures of LiNbO 3 , perovskite and ilmenite are interrelated and some transitions between them are possible [69,70]. The ilmenite form of LN has been synthetized by ion-exchange from the ilmenite form of sodium niobate, NaNbO 3 [71][72][73]. This LN form transforms back to the regular structure above 900 K. The estimated enthalpy difference between the two LN forms [73] is in good agreement with the calculated difference of 0.1 eV in their lattice energies [67].
The channeling studies on near-stoichiometric LN clearly indicate the absence of any intrinsic disorder-including the occurrence of ilmenite stacking faults-suggesting that in this case the small energy difference is sufficient to suppress their existence. However, LN has a strong tendency to grow non-stoichiometric with the congruent melt showing an excess of Nb. The resulting need for an intrinsic defect structure may locally outbalance the extra energy needed for the formation of ilmenite-type stacking to accommodate the excess Nb and the Li vacancies leaving the "stoichiometric" remainder of the crystal virtually defect-free. This arrangement may result in a global energy minimum for the congruent material.
The current knowledge on the properties of ilmenite-like LN is still limited. Wiegel et al. [74] have reported the UV optical absorption edge in ilmenite-like LN moved towards longer wavelengths when compared to a congruent crystal. For comparison purposes, the UV absorption edge of LN moves towards shorter wavelengths when approaching stoichiometric composition [75]. Baran et al. [76] have reported a very strong Raman scattering peak at about 735 cm −1 for ilmenite-like LN. Kong et al. [77] have reported a weak peak at 738 cm −1 in congruent LN that is absent in near-stoichiometric crystals, in congruent crystals doped with MgO above 2 mol%, and in congruent crystals doped with 2.5 mol% TiO 2 and suggested that it comes from an ilmenite-like stacking defect in nonstoichiometric LN. Lengyel et al. [78] confirmed the initial decrease of the intensity of a band at around 740 cm −1 with MgO doping, but have seen this band reappear for heavily doped crystals.

Ferroelectric to Paraelectric Phase Transition
The high-temperature phase of LN, above T C ≈ 1470 K for material with the congruent composition [79], is paraelectric and belongs to space group R3c. For LT, the transition temperature is significantly lower, T C ≈ 878 K [80]. In this phase, Nb (or Ta) a moves along the c-axis to the center of the octahedral oxygen cages, while for Li two models were proposed in the literature [81]. In the displacive model, Li moves into a triangularly coordinated site within the oxygen layer below the former Li site [82], while in the order-disorder model, the Li atoms are randomly distributed on both sides of the oxygen plane [83,84], as depicted in Figure 16. There is no agreement among the theoretical calculations on the phase transition. The calculations of Inbar et al. [85] suggest order-disorder features for the oxygen atoms, while the Li atoms are passive players. More recent molecular dynamics calculations of Phillpot at al. [86] and Sanna et al. [87] suggest that the phase transition presents displacive features for Nb/Ta and order-disorder features for Li. In contrast, the theoretical calculations of Toyoura et al. [88] suggest order-disorder features for both Li and Nb/Ta.
There are few measurements of the temperature dependence of the nuclear quadrupole interaction at the Nb/Ta site in LN and LT. Schempp et al. [41] reported 181 Ta NQR measurements in polycrystalline near-stoichiometric LN in a reduced temperature range from 21 to 515 K, where a continuous decrease of the EFG was observed. Moreover, the linewidths of the signal from near-stoichiometric samples were markedly narrower than the ones from Li-deficient material and became narrower with the increase in temperature. Figure 17 plots 181 Ta NQR and 181 Ta ME results in polycrystalline LT reported by Zhukov et al. [89] and by Löhnert et al. [90], respectively. The variation of the EFG derived from the NQR and ME measurements shows a continuous decrease with temperature in an extended temperature range which includes TC. This variation can be essentially reproduced using the point-charge model of De Wette et al. [91], taking into account thermal lattice expansion and atom displacement [90], using data from Abrahams et al. [83]. As thermal lattice expansion alone would make the EFG increase moderately with temperature, the main contribution to the EFG comes from the displacement of the atoms. More advanced ab initio cluster calculations of Shelyapina et al. [92] reproduced the essential features of the EFG temperature dependence up to TC. There is no agreement among the theoretical calculations on the phase transition. The calculations of Inbar et al. [85] suggest order-disorder features for the oxygen atoms, while the Li atoms are passive players. More recent molecular dynamics calculations of Phillpot at al. [86] and Sanna et al. [87] suggest that the phase transition presents displacive features for Nb/Ta and order-disorder features for Li. In contrast, the theoretical calculations of Toyoura et al. [88] suggest order-disorder features for both Li and Nb/Ta.
There are few measurements of the temperature dependence of the nuclear quadrupole interaction at the Nb/Ta site in LN and LT. Schempp et al. [41] reported 181 Ta NQR measurements in polycrystalline near-stoichiometric LN in a reduced temperature range from 21 to 515 K, where a continuous decrease of the EFG was observed. Moreover, the linewidths of the signal from near-stoichiometric samples were markedly narrower than the ones from Li-deficient material and became narrower with the increase in temperature. Figure 17 plots 181 Ta NQR and 181 Ta ME results in polycrystalline LT reported by Zhukov et al. [89] and by Löhnert et al. [90], respectively. The variation of the EFG derived from the NQR and ME measurements shows a continuous decrease with temperature in an extended temperature range which includes T C . This variation can be essentially reproduced using the point-charge model of De Wette et al. [91], taking into account thermal lattice expansion and atom displacement [90], using data from Abrahams et al. [83]. As thermal lattice expansion alone would make the EFG increase moderately with temperature, the main contribution to the EFG comes from the displacement of the atoms. More advanced ab initio cluster calculations of Shelyapina et al. [92] reproduced the essential features of the EFG temperature dependence up to T C .  181 Ta NQR and ME. Drawn using data from references [89,90]. The dashed line was drawn to guide the eye.
The temperature dependence of the nuclear quadrupole interaction at the Li site was measured using NMR, MS and PAC methods in LN and LT. 7 Li NMR measurements in single-crystalline LN were reported by Halstead [39] below TC, up to ≈950 K. A linear increase of the QIF was observed in this temperature range. The NMR signal had broad lines at low temperature, which became narrower with the increase in temperature, especially after ≈800 K. 7 Li NMR measurements in polycrystalline LT were reported by Slotfeldt-Ellingsen et al. [93,94], below and above TC. A linear increase of the QIF was observed until close to TC and then a decrease. PAC measurements in polycrystalline LN and LT doped with 181 Hf were reported by Catchen et al. [50,95]. They observed extensive line broadening at temperatures well below TC and significantly less line broadening above TC. The anisotropy coefficient was non-zero below TC and dropped to zero above TC. This result was not expected, because the axial symmetry at the Li site implies that η should vanish at temperatures both below and above TC. A second QIF was reported in the measurements above TC, with a value about an order of magnitude lower than the one of Hf probes in Li sites, which was tentatively attributed to Hf probes in Nb/Ta sites [95]. Moreover, in the case of a previous MS study in Fe doped LN, a second QIF was found at high temperature, but no more details were given [96]. Marques et al. [97] reported high precision PAC measurements in single-crystalline LT implanted with 111 In. While at RT two nearby QIF were found, νQ1 = 230(2) MHz and νQ2 = 242(2) MHz, with nonzero asymmetry parameters, above TC only one QIF and η = 0 was found. The observation of two nearby QIF is consistent with the previous measurements of Hauer et al. [53] in congruent LN and with the extensive line broadening reported in other studies well below TC. Ohkubo et al. [98] reported PAC measurements in polycrystalline LT with the 111m Cd and 117 Cd probes. Figure 18 shows the temperature dependence of the QIF for the measurements in LT.  181 Ta NQR and ME. Drawn using data from references [89,90]. The dashed line was drawn to guide the eye.
The temperature dependence of the nuclear quadrupole interaction at the Li site was measured using NMR, MS and PAC methods in LN and LT. 7 Li NMR measurements in single-crystalline LN were reported by Halstead [39] below T C , up to ≈950 K. A linear increase of the QIF was observed in this temperature range. The NMR signal had broad lines at low temperature, which became narrower with the increase in temperature, especially after ≈800 K. 7 Li NMR measurements in polycrystalline LT were reported by Slotfeldt-Ellingsen et al. [93,94], below and above T C . A linear increase of the QIF was observed until close to T C and then a decrease. PAC measurements in polycrystalline LN and LT doped with 181 Hf were reported by Catchen et al. [50,95]. They observed extensive line broadening at temperatures well below T C and significantly less line broadening above T C . The anisotropy coefficient was non-zero below T C and dropped to zero above T C . This result was not expected, because the axial symmetry at the Li site implies that η should vanish at temperatures both below and above T C . A second QIF was reported in the measurements above T C , with a value about an order of magnitude lower than the one of Hf probes in Li sites, which was tentatively attributed to Hf probes in Nb/Ta sites [95]. Moreover, in the case of a previous MS study in Fe doped LN, a second QIF was found at high temperature, but no more details were given [96]. Marques et al. [97] reported high precision PAC measurements in single-crystalline LT implanted with 111 In. While at RT two nearby QIF were found, ν Q1 = 230(2) MHz and ν Q2 = 242(2) MHz, with nonzero asymmetry parameters, above T C only one QIF and η = 0 was found. The observation of two nearby QIF is consistent with the previous measurements of Hauer et al. [53] in congruent LN and with the extensive line broadening reported in other studies well below T C . Ohkubo et al. [98] reported PAC measurements in polycrystalline LT with the 111m Cd and 117 Cd probes. Figure 18 shows the temperature dependence of the QIF for the measurements in LT.
The temperature dependence of the QIF doublet reported by Marques et al. [97] is shown in Figure 19. The temperature dependence of ν Q1 is close to the observed with 7 Li in LT [93], as expected since the 111 In probe atoms replace Li. The same is applicable for the 181 Hf, 117 Cd and 111m Cd probes which also replace Li. In contrast, ν Q2 that is associated to probes in Li-II sites, shows only a very weak temperature dependence. Close to and above T C , the derived ν Q1 and ν Q2 values are no longer distinguishable. This suggests that the differences between the Li-I and Li-II sites vanish in the high-temperature paraelectric phase. The temperature dependence of the QIF doublet reported by Marques et al. [97] is shown in Figure 19. The temperature dependence of νQ1 is close to the observed with 7 Li in LT [93], as expected since the 111 In probe atoms replace Li. The same is applicable for the 181 Hf, 117 Cd and 111m Cd probes which also replace Li. In contrast, νQ2 that is associated to probes in Li-II sites, shows only a very weak temperature dependence. Close to and above TC, the derived νQ1 and νQ2 values are no longer distinguishable. This suggests that the differences between the Li-I and Li-II sites vanish in the high-temperature paraelectric phase. Point-charge model EFG calculations reproduce the observed linear increase of νQ1 up to TC taking into account thermal lattice expansion only [97]. If the corresponding Li/O displacements of Abrahams et al. [83] are taken into account, then the EFG would be too   117 Cd, 111m Cd and 7 Li probes in lithium tantalate. Drawn using data from references [50,93,95,98].
The temperature dependence of the QIF doublet reported by Marques et al. [97] is shown in Figure 19. The temperature dependence of νQ1 is close to the observed with 7 Li in LT [93], as expected since the 111 In probe atoms replace Li. The same is applicable for the 181 Hf, 117 Cd and 111m Cd probes which also replace Li. In contrast, νQ2 that is associated to probes in Li-II sites, shows only a very weak temperature dependence. Close to and above TC, the derived νQ1 and νQ2 values are no longer distinguishable. This suggests that the differences between the Li-I and Li-II sites vanish in the high-temperature paraelectric phase. Point-charge model EFG calculations reproduce the observed linear increase of νQ1 up to TC taking into account thermal lattice expansion only [97]. If the corresponding Li/O displacements of Abrahams et al. [83] are taken into account, then the EFG would be too Point-charge model EFG calculations reproduce the observed linear increase of ν Q1 up to T C taking into account thermal lattice expansion only [97]. If the corresponding Li/O displacements of Abrahams et al. [83] are taken into account, then the EFG would be too large by a factor of two close to T C . The main contribution to the EFG at a Li site comes from the six surrounding O atoms, distributed in two planes, with Li-O distances of 2.08 Å and 2.29 Å at room temperature. At 940 K these distances are 2.05 Å and 2.51 Å [83], which results in a significantly increased EFG. To explain the observed decrease above T C it is necessary to decrease the difference between the two Li-O distances. Inbar and Cohen [85] proposed that the driving mechanism for the phase transition in LT is a displacement of the O atoms towards the Ta atoms. Since a displacement of the O only along the c-axis would result in too short Ta-O distances, they further proposed that the O atoms move also in the plane perpendicular to the c-axis, toward Li and lead to a coupled Li-O motion. EFG calculations using the above-mentioned point charge model show that such a Li-O coupled movement explains qualitatively the decrease of the EFG above T C [97].

Ion Implantation Induced Defects
Ion implantation matured to a standard technique for defect engineering and doping of all kinds of materials. In the case of LN the main applications that will be discussed in their respective section are the formation of optical waveguides (Section 5.3) as well as the introduction of optical active (Section 4.3) and insoluble dopants (Section 4.5). This section is dedicated solely to investigation of basic effects of ion implantation leading to the production of intrinsic defects and their removal by annealing.
The production of defects occurs in two ways: (i) by collisions of the incoming ion with a lattice atom that cause its displacement from the lattice site, eventually initiating a collision cascade (nuclear damage); (ii) by excitation of electrons in the material which, in the case of swift heavy ions, can cause the formation of amorphization tracks on the nanometer scale (electronic damage). Since the nuclear damage is mainly occurring at the end of the ion track, it is the dominating damage producer at low energy irradiation. On the other hand, the electronic damage dominates in swift ion irradiations on the major part of the trajectory. In both cases, increasing fluence causes the defects to accumulate until full amorphization of the material is achieved. Defects cause dechanneling similar to interstitials allowing to determine their concentration and depth distribution by ion channeling.

Nuclear Damage Dominated Cases
For irradiation, already the first systematic studies on ion implantation in LN by Jetschke et al. [99] and Götz and Karge [100,101] showed a clear anisotropy in the generation of lattice defects. For irradiation with 150 keV N + ions at room temperature amorphization at the surface was reached with fluences of 5 × 10 15 cm −2 and 2 × 10 16 cm −2 for x-and y-cut material. In contrast, z-cut material sustained a heavily damaged but not amorphized surface layer even at a fluence of 4 × 10 16 cm −2 . Three different stages in the damage evolution were distinguished: (i) In the pre-damage stage (fluence < 1 × 10 15 cm −2 ) mainly point defects were produced of which a significant fraction was attributed to Nb displaced from its regular lattice site to the free octahedron. This explains the lower dechanneling along the 0001 -axis in which these displaced Nb is still aligned with the guiding atomic rows, while for the 1120 -axis this lattice site appears in the center of the channel enhancing the dechanneling. (ii) The heavy damage stage (1 × 10 15 cm −2 < fluence < 5 × 10 15 cm −2 ) defect clusters form that cause 1120 -and 0110 -directions to be amorphous but still maintain some alignment with the 0001 -axis. (iii) In the final stage (fluence > 5 × 10 15 cm −2 ) saturation of defect formation is observed along the 1120 -and 0110 -axis. A recent study [102] revisited the issue employing 350 keV Ar + implantation into samples with x-and z-cut as well as a sample cut in a direction equally distanced from 1120 -and 0001 -axis by 45 • which allowed to study both axial directions in a single sample. Plotting the mean maximum damage concentration versus the ion fluence yielded perfect overlaps for 1120 -and 0001 channeling measurements for the x-z cut sample with those of the respective standard cuts. The same work presents an additional confirmation of displaced Nb occupying the intrinsic vacant by angular-resolved RBS/C for a 1 MeV I + implanted x-cut LN.
A study on ion irradiation at low temperature (15 K) was performed by Gischkat et al. [103] for various ions (50 keV H + , 160 keV Li + , 160 keV O + , 350 keV Ar + ) at fluences ranging from 5 × 10 11 cm −2 and 2 × 10 17 cm −2 . The analysis by RBS/C was done in situ at the same low temperature. Nevertheless, the analyzing He + -ion beam induced an annealing effect attributed to energy deposition by electronic processes. Taking into account this effect, the undisturbed damage concentrations were calculated indicating an enhanced defect formation with regard to room temperature.
Bianconi et al. [104] deduced, based on own experimental (C + , 0.5-5.0 MeV, O + , 0.44-8.0 MeV) and literature data, a universal curve for the Nb-defect fraction in the endof-track region as function of the energy deposited by nuclear stopping. The relation holds for all ions with Z ≤ 14 independent of the fluence.
While the Nb-related damage is easily detected by RBS/C the other two constituents of LN are more difficult to access. In a work by Schmidt et al. [105] the 7 Li(p,α) 4 He reaction was employed to study the Li-related damage in x-and z-cut material irradiated with 1 MeV I + (fluences between 2 × 10 13 cm −2 and 1 × 10 14 cm −2 ). The evolution of the Li-damage concentration with the fluence rate was found to accompany closely that for Nb in the z-cut material. Due to the low contribution of Li to the channeling effect in LN other crystallographic directions than along the 0001 -axis are difficult to investigate in terms of defect formation.
Götz and Karge [100] also found an enhancement of the etching rate in diluted hydrofluoric acid for room temperature implantation fluences above 5 × 10 14 N + cm −2 that increased with increasing fluence. Schrempel et al. [106] investigated the etching rates for z-cut LN implanted at 15 K and room temperature with 40 keV He + and 350 keV Ar + with fluences ranging from 5 × 10 12 to 5 × 10 16 cm −2 . The larger amount of damage created by low temperature irradiation as well as the higher nuclear damage produced by the irradiation by heavier ions increase the etching rate strongly.
Jentschke and Hehl [107] studied the influence of annealing procedures in air on the implantation induced by in P + (250 keV) and N + (150 keV) implanted LN. For the case of a partial damaged (5 × 10 14 N + cm −2 ) x-cut sample, a virtually complete recovery of the lattice was achieved after annealing at 300 • C for 30 min. On the other hand, x-and z-cut samples with RBS channeling spectra indicating full amorphization show only a moderate reduction of the dechanneling even for annealing temperatures up to 1100 • C. This is related to the loss of Li by Li 2 O outdiffusion from the amorphous layer during the annealing process. Random RBS spectra of these samples show a surface region with a higher yield of Nb than for the virgin crystal due to the formation of the Li-deficient LiNb 3 O 8 compound. The loss of Li during annealing at high temperature can be avoided by substituting air by flowing wet oxygen resulting in high-quality recrystallization of the material. This is demonstrated in Figure 20 for Cr-implanted LN after annealing in the two different atmospheres.

Electronic Damage Dominated Cases
For ions with specific energies above 0.1 MeV/amu, the electronic stopping power outweighs the nuclear along most of the path. For the formation of damage consisting of amorphous tracks the electronic stopping power of the ions has to exceed a material dependent threshold value.
The effect of electronic damage is seen most perfectly separated from nuclear damage for irradiation of materials with GeV ions. Canut et al. [108][109][110] performed a series of irradiation of LN using Sn-, Gd-and U-ions with energies in the GeV range and fluences ranging between 1 × 10 11 to 1 × 10 12 cm −2 . RBS/C investigations reveal even for the lowest fluences a significant damage in the observable depth range. The lattice disorder was found to increase with ion energy and fluence showing some saturation tendency for the latter at the maximum fluence. Plotting the damage fraction against the electronic stopping power for the different ions reveals a linear dependence.
Irradiations with lower ion energies, a few to tens of MeV, are of interest for optoelectronic applications and have therefore experienced an immense attention in the last two decades. In the following only a selective overview on the observed features can be given. Bentini et al. [111] investigated the implantation of 5 MeV O 3+ ions with fluences of 1 × 10 14 cm −2 and 6 × 10 14 cm −2 into x-cut LN. RBS/C spectra of the as-implanted samples show for all fluences a damage peak at the end-of-range (2.7 µm) of the implantation, while a second peak appears at the surface region for fluences above 4 × 10 14 cm −2 that has to be attributed to electronic damage. Samples annealed for 2 h at 235 • C in dry oxygen show a significant for the low fluence irradiation but a poor effect for higher fluences. The literature reports further similar RBS/C investigations on electronic damage for C [112], N [113], O [112][113][114][115], F [113,115,116], Si [117] and Mg [115].
Determination of the threshold electronic stopping power induced damage was the aim of several works. Rivera et al. [118] reported a value of about 6 keV/nm based on measurements for very low fluences (1 × 10 11 cm −2 to 1 × 10 12 cm −2 ) investigating a set of six different ion species. Ramos et al. [119] found a smaller value, 3.5 keV/nm, for irradiation with Cl-and Br-ions in the adjacent fluence range (1 × 10 12 cm −2 to 1 × 10 13 cm −2 ). An extensive study for the higher fluence range (1 × 10 13 cm −2 to 1 × 10 15 cm −2 ) using O-and Si-ions was recently presented by Wesch et al. [120]. A threshold value of 2.5 keV was extracted using the RBS/C data, in good agreement with the 2.2 keV reported in ref. [113] for 1 × 10 14 cm −2 O implantation. These data evidence the reduction of the threshold value with increasing fluence. A plot of the threshold electronic stopping power versus the logarithm of the ion fluence can be fitted with a straight line as demonstrated by Agulló-Lopez et al. [121].
The theoretical modelling of the processes leading to amorphization experimentally observed in the electronic stopping power regime is a very active field, steadily generating new approaches. A detailed discussion is therefore out of the scope of this paper. In brief, the most used description is based on thermal spike model that assumes that local melting of the crystal structure occurs above the threshold value [122]. However, this cannot explain the damage production observed even below the threshold. Other effects as the exciton based mechanisms have been invoked as an explanation for this effect [123,124]. Finally, the cumulative character of defect generation has been invoked to explain the damage enhancement and amorphization in sub-threshold irradiation [120,125].
The defect formation by Si + ion irradiation of LN along crystallographic directions has been investigated by Schrempel et al. [126,127]. For the case of low ion energies (<1 MeV for Si + ) it was observed that the nuclear damage maxima were shifted to higher depth for aligned irradiation with respect to those of random incidence. This effect is related to the reduced stopping power in axial direction due to the channeling effect that keeps the ion in regions with lower electron density and consequently less energy loss. Additionally, the aligned irradiation creates significantly less defects for the same fluence. For higher energies the electronic stopping passes the threshold. Surface region defect formation is observed and the shift of the damage maxima for channeling irradiation decreases again.

Neutron Irradiation Induced Defects
The interest in the behavior of LN-based acoustical devices close to the core of a nuclear fission reactor started in the 1970s [128][129][130][131]. Materials close to a reactor core are exposed to neutrons from the fission of 235 U and to gamma radiation [132]. While the neutrons directly released from fission have an average energy of ≈2 MeV (thus being considered fast neutrons), neutron moderation in most reactors extends the spectrum down to the meV energy range (thermal neutrons). Thus, most irradiations are done with neutrons with a wide energy range, unless specific filters are interposed [133]. Besides structural characterization performed with neutron scattering, LN crystals mostly with congruent composition were irradiated with fast and/or thermal neutrons with fluences up to 10 20 cm −2 [31,49,[134][135][136][137][138][139][140][141], although with scarce details on the actual neutron spectrum in many cases. An expansion along the c-axis and a shift towards longer wavelengths of the UV absorption edge were reported by Dowell et al. [128] for neutron fluences up to 4 × 10 17 cm −2 . Primak et al. [130] reported that irradiated LN becomes highly disordered and shows loss of piezoelectric response and optical birefringence at a fluence of 8 × 10 19 cm −2 . This process of destruction of LN was found to be significantly reduced when the LN sensors were kept at ≈800 K during the irradiation [131].
The basic physics of neutron damage production is well understood [142,143]. Fast neutrons with energy up to a few MeV produce displacement damage mostly through elastic collisions. In contrast, thermal neutrons produce displacement damage mostly through (n,γ) reactions, due to the recoil induced by the emission of the γ photons. In the case of LN one has to consider also the 6 Li(n,α) 3 H nuclear reaction, with a high cross-section for thermal neutrons [144] and whose products, α and 3 H, also produce displacement damage. Although in general fast neutrons are expected to create more damage than thermal neutrons, this is not always the case [145]. Figure 21 shows the number of collisions per neutron incident on a thin slab of LN (100 µm thick) as function of energy, covering a range from 0.1 meV to 20 MeV in 69 groups, calculated with the MCNPX Monte Carlo code [146]. Besides the total number of collisions, the contributions coming from collisions with Li, O and Nb are also shown. In the thermal neutron region, until about 1 eV, the number of collisions is essentially driven by reactions with Li (and within Li, with its less abundant isotope, 6 Li). After this energy, the number of collisions per incident neutron is essentially constant until about 1 MeV, when it starts to decrease.
In order to evaluate the damage caused by the above collisions, it is necessary to know the threshold energy for displacement of the atoms in LN. Seitz estimated the threshold energy for displacement of an ion or atom in a tightly bound solid, E d , to be 25 eV [147]. Calculations of E d values have been done using molecular dynamics for several technologically-relevant materials (see, e.g., references [148][149][150][151]), but are scarce for LN. The displacement energy values calculated by Petersen [152] are E d (Li) = 30 eV, E d (Nb) = 34 eV and E d (O) = 22 eV. The value of E d for oxygen has been estimated from irradiations with electrons as 53 eV [153,154], which is well within the normal range for oxides [155], but significantly higher than obtained by calculation. An estimate for the value of E d for Hf in Nb and Li sites can be obtained from the PAC measurements of Marques et al. [31]. The typical neutron fluence (all energies) for the preparation of the PAC measurements was less than 5 × 10 16 cm −2 , well below the fluence for which amorphization was observed [130]. Annealing in air at 573 K was enough to remove the coloration induced by the irradiation and annealing at 973 K was enough to remove all defects, as seen by RBS [31].  In order to evaluate the damage caused by the above collisions, it is necessary to know the threshold energy for displacement of the atoms in LN. Seitz estimated the threshold energy for displacement of an ion or atom in a tightly bound solid, Ed, to be 25 eV [147]. Calculations of Ed values have been done using molecular dynamics for several technologically-relevant materials (see, e.g., references [148][149][150][151]), but are scarce for LN. The displacement energy values calculated by Petersen [152] are Ed(Li) = 30 eV, Ed(Nb) = 34 eV and Ed(O) = 22 eV. The value of Ed for oxygen has been estimated from irradiations with electrons as 53 eV [153,154], which is well within the normal range for oxides [155], but significantly higher than obtained by calculation. An estimate for the value of Ed for Hf in Nb and Li sites can be obtained from the PAC measurements of Marques et al. [31]. The typical neutron fluence (all energies) for the preparation of the PAC measurements was less than 5 × 10 16 cm −2 , well below the fluence for which amorphization was observed [130]. Annealing in air at 573 K was enough to remove the coloration induced by the irradiation and annealing at 973 K was enough to remove all defects, as seen by RBS [31].
In a congruent LN crystal doped with 1 mol% HfO2 the spectra before and after thermal annealing show the same characteristic QIF around 1200 MHz, associated with Hf probes in Li sites. In contrast, in a near-stoichiometric crystal doped with 1 mol% HfO2 where Hf is distributed between Li (45%) and Nb (55%) sites, the spectrum after irradiation shows that the fraction of Hf probes in Nb sites is strongly reduced to 18%, with a full recovery only after annealing at 973 K. However, RBS data shows no change on the lattice site of Hf after irradiation [31]. Thus the change of lattice location after irradiation can only affect the radioactive 181 Hf atoms produced via the (n,γ) reaction. Since these atoms are much less than 1 ppm of the total amount of Hf in the crystal, no effect can be detected by RBS. It thus seems reasonable to assume that the displacement is due to the recoil induced by emission of the γ photons from the 180 Hf(n,γ) 181 Hf reaction, characterized by Alenius et al. [156]. As the maximum recoil that can be induced is 96 eV, In a congruent LN crystal doped with 1 mol% HfO 2 the spectra before and after thermal annealing show the same characteristic QIF around 1200 MHz, associated with Hf probes in Li sites. In contrast, in a near-stoichiometric crystal doped with 1 mol% HfO 2 where Hf is distributed between Li (45%) and Nb (55%) sites, the spectrum after irradiation shows that the fraction of Hf probes in Nb sites is strongly reduced to 18%, with a full recovery only after annealing at 973 K. However, RBS data shows no change on the lattice site of Hf after irradiation [31]. Thus the change of lattice location after irradiation can only affect the radioactive 181 Hf atoms produced via the (n,γ) reaction. Since these atoms are much less than 1 ppm of the total amount of Hf in the crystal, no effect can be detected by RBS. It thus seems reasonable to assume that the displacement is due to the recoil induced by emission of the γ photons from the 180 Hf(n,γ) 181 Hf reaction, characterized by Alenius et al. [156]. As the maximum recoil that can be induced is 96 eV, corresponding to the emission of a 5.6 MeV γ photon, only one displacement per recoil is expected.
The transfer of Hf form Nb to Li sites upon irradiation can be qualitatively explained by the existence of Li vacancies in the congruent and near-stoichiometric crystals. In this way, a recoiled 181 Hf from a Nb position finds easily a vacant Li site, leading to the observed change. In the case of a recoiled 181 Hf from a Li site, which finds another Li site, no change is seen. If the energy necessary to displace Hf from Li and Nb lattice sites is approximately the same and x represents the fraction of 181 Hf recoiled from Li and Nb sites, we can write:  [152].
Displacement of other dopants as a result of (n,γ) reactions is possible. Table 1 compiles data on the neutron capture cross section at 0.025 eV (1 b = 100 fm 2 ) [157], main activation reactions(s), respective cross section [158], maximum γ energy released [159] and maximum recoil energy. For the estimation of displacement damage, Dale et al. [160] have shown that is convenient to use a quantity named Non-Ionizing Energy-Loss (NIEL) to quantify the portion of energy lost by a particle that is used for displacement damage, with the remainder of the energy being dissipated as heat. It was shown that NIEL can be used to correlate displacement damage effects in materials induced by different types of radiation [160], even if some exceptions are known [161]. For a particle with energy E 0 incident on a material with a single type of atom with atomic weight A, it is defined as: where is N A is Avogadro's number, dσ/dT is the total differential cross section (elastic and inelastic) for energy transfer to an atom of the material, T min and T max are the minimal and maximal energies that can be transferred to a lattice atom, and Q(T) is a partition function which gives the fraction of T that is lost to NIEL. T min is usually taken as T min = 2E d , while T max = 4E 0 A/(A + A') 2 , where A' is the mass of the incident particle. The partition function Q(T) has a classical formulation given by Lindhard et al. [162] and fitted by Robinson et al. [143], as well as a more recent formulation by Akkerman et al. [163] with improved handling of low energy ions. Total NIEL for LN results from the application of Equation (4) to Li, Nb and O, taking stoichiometry in account, as: The product of the NIEL and the particle fluence gives the displacement damage energy deposition per unit mass of material. NIEL plays the same role to the displacement damage energy deposition as the stopping power to the total ionizing dose. Figure 22 shows the values of NIEL for electrons, neutrons, protons and alphas incident in LN, with energy up to 100 MeV, calculated with the tools available at the SR-NIEL platform [164], using the Akkerman partition function. Displacement energy values of 30 eV, 35 eV and 25 eV were used for Li, Nb and O, as the nearest to the values calculated by Petersen [152], due to the constraints of the calculations for neutrons. NIEL for incident neutrons has a weak dependence of the value of E d : a change from 20 to 50 eV in the E d of the 3 species leads to a neutron NIEL average decrease of 3% in the range 10 −4 -10 −2 MeV and less than 1% outside this range. Table 2 shows discrete values, for reference. The tabled values for neutrons were obtained by a simple linear interpolation between the two nearest values. It is worth noting that at 1 MeV protons are about 20 times more effective than neutrons in producing displacement damage in LN. energy up to 100 MeV, calculated with the tools available at the SR-NIEL platform [164], using the Akkerman partition function. Displacement energy values of 30 eV, 35 eV and 25 eV were used for Li, Nb and O, as the nearest to the values calculated by Petersen [152], due to the constraints of the calculations for neutrons. NIEL for incident neutrons has a weak dependence of the value of Ed: a change from 20 to 50 eV in the Ed of the 3 species leads to a neutron NIEL average decrease of 3% in the range 10 −4 -10 −2 MeV and less than 1% outside this range. Table 2 shows discrete values, for reference. The tabled values for neutrons were obtained by a simple linear interpolation between the two nearest values. It is worth noting that at 1 MeV protons are about 20 times more effective than neutrons in producing displacement damage in LN.

Extrinsic Point Defects
This section is dedicated to the determination of lattice sites occupied by foreign atoms in LN. The influence of dopant concentration, co-doping and stoichiometry of the material are addressed. The subsections are mainly organized according to the functionality of the dopant.

Dopants to Increase Photorefractive Damage Resistance
Although congruent LN exhibits excellent optical, electro-optical and non-linear optical properties, the relatively low optical damage resistance hampers its application in optical devices [165]. Besides the attempt to improve the optical damage resistance by growing stoichiometric LN, the introduction of divalent (Mg 2+ [166,167], Zn 2+ [168]), trivalent (Sc 3+ [169], In 3+ [170]) and tetravalent (Hf 4+ [171], Zr 4+ [172]) dopants has been proven to be successful. The observed minimum concentration for suppression of the optical damage (4.5 mol% for doping with MgO, 6 mol% for ZnO, 1.5 mol% for Sc 2 O 3 and In 2 O 3 and 4 mol% for HfO 2 and 6 mol% for ZrO 2 ) indicate that it is important to clarify the lattice location and interaction of these dopants with intrinsic defects of LN.
As the doping with Mg is by far the most technically most important way innumerous studies about the characteristics and application extensively studied and attempts were made. From the more fundamental viewpoint, the main interest is to understand the incorporation of this dopant in dependence on its concentration and to elucidate the underlying mechanisms.
Most of the theoretical approaches published so far depend largely on the model adopted to explain the non-stoichiometry of the congruent material. Zhang and Feng [173] proposed a series of overlapping concentration dependent mechanisms, based on the Nbvacancy model. At low concentration (<2.7 mol%) Mg substitutes Nb on antisites until they are consumed; for higher concentration up to the threshold Mg is incorporated at the Li site and accompanied by a trapped electron for charge compensation; above the threshold the incorporation occurs on both Li and Mg sites. These ideas were readopted by Liu et al. [174] and more recently by Abdi et al. [175]. Based on the Li-vacancy model, Iyi et al. [176] propose that Mg at first replaces Nb at Li sites up to 3 mol% MgO at which the removal of Nb Li is completed. In the range between 3 and 8 mol% MgO magnesium will only occupy regular Li sites while at higher concentration both cation sites are incorporating Mg. The studies by Donnerberg et al. [177,178] (further refined in [179,180]) are based on the occurrence of ilmenite-type stacking faults with Li vacancies but are formally treated as Nb antisites and Nb vacancies. In the following, we also adopt this procedure for simplicity. According to this model, Mg is replacing in the two concentration ranges < 1.5 mol% and 1.5 to 5.0 mol% simultaneously Nb on antisites and Li on regular sites employing different mechanisms. The preferential substitution of Nb Li by Mg in the second range results in their complete removal at the threshold level. For even higher Mg concentrations both Li and Nb sites should be occupied.
Since a direct determination of the lattice site of magnesium in dependence of its concentration has to be regarded as essential to decide which of the proposed models is the most adequate, an extended channeling study on a set of Mg-doped crystals was performed employing He + in order to improve the attainable minimum yield and minimize bremsstrahlung background [181]. As no adequate nuclear reaction for the study of magnesium is available, PIXE/C has been applied. Due to very low energy of the characteristic Mg-K X-rays a Si(Li)-detector with an extremely thin Be window was used to limit absorption. This set of crystals, containing one undoped and six doped (0.5, 1.0, 2.0, 4.0, 6.0 and 9.0 mol% MgO in melt), was grown at the Siemens Research Laboratory in Munich and thoroughly characterized regarding crystal composition, lattice constants, density, temperature, piezoelectricity and various optical properties [182][183][184][185][186][187].
Angular scans crossing through the 0001 -axis showed for all concentration dips with a width and minimum yield for the Mg-PIXE between those of Nb-RBS and Li-NRA. Computer simulations with CASSIS confirmed that this is consistent with the absence of non-axial lattice sites for Mg.
For MgO concentrations of 0.5 and 1.0 mol% the most striking observation is the strong peak in the Mg-PIXE yield for alignment with the 1120 -axis. This indicates that Mg is not occupying a regular cation lattice site but some interstitial site. The fact that the Mg-PIXE scans for the other axial directions show clear dips and a similar behavior to the Li-scans one can further exclude that Mg is located in the intrinsically vacant or in the Nb-octahedron. The computer simulation with the CASSIS code allows to identify the occupied position as the center of the Li-Octahedron which yields for Mg an oxygen environment similar to that of magnesium oxide. The epitaxial growth of LN on MgO substrates it is known that the (111)-plane of MgO and (001)-plane of LN have the same oxygen framework [188] which might explain this behavior.
It is tempting to identify the observed behavior with a substitution of Nb in the stacking but this would require a substitution of the following type: incorporating so much lithium that a near-stoichiometric crystal would be obtained. However, the measurements on the lithium fraction and Curie temperature for these samples exclude this scenario [182]. It is more likely that the Mg occupied sites associated to the ilmenite stacking faults (Li vacancies or Li sites in the vicinity) and the movement to the center of the Li-Octahedron can be interpreted as a relaxation mechanism. The Mg-scans along the 1120 -axis for concentrations above 2.0 mol% MgO in melt differ clearly from those at low concentration. For all axial directions clear minima are observed with scan curve that resemble those of Li. However, computer simulations demonstrate that Mg solely occupying regular Li sites would result in minimum yields that are significantly lower than the observed ones. The best fit to the experimental data was obtained assuming 2/3 of the Mg to occupy regular Li sites and 1/3 sites in the center of the Nb oxygen cage (occupation of regular Nb sites yielded even lower minima). Interestingly this observation holds for the whole concentration range from 2.0 to 9.0 mol% MgO, without any change at the threshold concentration for photorefractive damage suppression. Indeed, there exists a substitution reaction for divalent dopants that would be compatible with these findings: which on the other hand does not remove any antisites. A possible solution is that the reactions described by Equations (6) and (7) occur simultaneously until the threshold of 6 mol% Mg in the crystal is reached, causing a deviation in the Mg Li /Mg Nb ratio that is too small to be detectable by the PIXE/C measurements. A further study on this set of crystals focused on the influence of the Mg incorporation on the Li-NRA channeling minimum yield [189] providing further insight on the lattice modifications. The minimum yield for an atomic string in channeling directions depends strongly on repelling force of its continuum potential. The replacement of Li (Z = 3) by atoms with higher atomic number Nb (Z = 41) or Mg (Z = 12) increases the continuum potential of the Li-string and therefore reduces the minimum yield. Figure 23 shows the dependence of the Li-NRA minimum yield in various axial directions in dependence of the MgO-concentration. The main observation is that the minimum yield increases steeply for all but the 0001 -axis (in which all cations are aligned with each other) strongly for MgO melt concentrations up to 1 mol%. This indicates an increase in the disturbance of the lattice due to the interstitial location of magnesium rather than changes of the continuum potential due to removal of excess Nb. In the investigated axial directions ilmenite-like and regular Li site occupying Nb are located very close and do not produce distinguishable effects. For further increasing Mg concentrations the 1120 -and 0441axis no systematic alteration of the minimum yield is observed since the Li atomic strings are in close neighborhood to the-in terms of channeling-dominating Nb-rows. In addition, for both axial directions χ min does not decrease anymore to the value of the congruent case. On the other hand, Li is isolated in the 1120 -axis within the channel and the incorporation of the Mg reinforces the repelling force of the Li strings. This leads to a steady decrease of the Li-NRA minimum yields for this axial direction reaching at 9 mol% MgO a value below that of the congruent material.
In addition, one NMR study on 25 Mg (nuclear spin 5/2+) has been reported by Feng et al. [190] for LN doped with MgO with concentration ranging from 1.0 to 6.0 mol%. The NMR spectra for 25 Mg show clear differences for the various concentration but surprisingly no new line, that would indicate change in the lattice site, appears for concentrations above the threshold. Yatsenko et al. [191] used an indirect approach to obtain information on the incorporation of Mg using 7 Li-and 93 Nb-NMR in samples doped with 2 to 7 mol%. They conclude that the amount of vacancies is steadily increasing up to the threshold at which on the other hand no more Nb-antisites are present. Mg is supposed to occupy Li sites forming next-neighbor Mg Li -V Li complexes. yield is observed since the Li atomic strings are in close neighborhood to the-in terms of channeling-dominating Nb-rows. In addition, for both axial directions χmin does not decrease anymore to the value of the congruent case. On the other hand, Li is isolated in the 〈112 0〉 -axis within the channel and the incorporation of the Mg reinforces the repelling force of the Li strings. This leads to a steady decrease of the Li-NRA minimum yields for this axial direction reaching at 9 mol% MgO a value below that of the congruent material. In addition, one NMR study on 25 Mg (nuclear spin 5/2+) has been reported by Feng et al. [190] for LN doped with MgO with concentration ranging from 1.0 to 6.0 mol%. The NMR spectra for 25 Mg show clear differences for the various concentration but surprisingly no new line, that would indicate change in the lattice site, appears for concentrations above the threshold. Yatsenko et al. [191] used an indirect approach to obtain information on the incorporation of Mg using 7 Li-and 93 Nb-NMR in samples doped with 2 to 7 mol%. They conclude that the amount of vacancies is steadily increasing up to the threshold at which on the other hand no more Nb-antisites are present. Mg is supposed to occupy Li sites forming next-neighbor MgLi-VLi complexes.
Although less popular as dopant for photorefractive suppression in LN, Hf is highly interesting for studies on the basic incorporation mechanism of dopants since it can be accessed simultaneously by RBS/C and PAC. The advantage use of the combination of these two nuclear techniques in investigating has already been demonstrated in the case of intrinsic defects.
Hf-doped crystals LN were grown with 1.0 mol% [49,192] and 6.0 mol% [138] of HfO2 in congruent melt. For comparison, samples of near-stoichiometric material, to which 4.6 wt % K2O was added, were doped with 0.2 and 1.0 mol% [138,193] HfO2. Although less popular as dopant for photorefractive suppression in LN, Hf is highly interesting for studies on the basic incorporation mechanism of dopants since it can be accessed simultaneously by RBS/C and PAC. The advantage use of the combination of these two nuclear techniques in investigating has already been demonstrated in the case of intrinsic defects.
RBS enabled to determine the actual concentrations of the melt grown samples and derive the corresponding distribution coefficients (Table 3). Channeling measurements for various axial direction proved that in the 1.0 mol% congruent crystal Hf occupies exclusively the Li site [49,192]. For the samples doped with 6.0 mol% RBS measurements revealed that the Hf distribution in the samples is inhomogeneous (3.5-5.1 mol%) and the distribution coefficient lower than one [138]. Further, the higher minimum yield (10%) at 0221 -axis for Nb and Hf in comparison to the 1.0 mol% sample (2.5%) indicates an inferior crystal quality. Both findings indicate that the maximum solubility of Hf in LN has been exceeded in the sample. However, the quality of angular scans was still sufficient to obtain meaningful results for the incorporated Hf-fraction: the best fit was obtained assuming 55% of the Hf to occupy regular or slightly shifted (0.1 Å towards the octahedron center) Nb lattice sites and 45% to occupy regular Li sites. This ratio is not compatible with a charge compensation mechanism suggested in theoretical studies for which 25% of the Hf has to reside on Li sites and 75% on Nb sites. Besides the already known frequencies for Hf in Li sites, the Fourier analyses of the PAC spectra obtained in the crystal doped with 6 mol% HfO 2 and in the near-stoichiometric crystals has shown another peak at a frequency about 3.5 times smaller. The relative intensity agreed with the fraction for Hf in Nb sites determined by RBS and the EFG scaled with the 93 Nb NMR measurements previously reported by Peterson et al. [44]. It can thus be safely assigned to Hf probes occupying Nb sites at room temperature. A frequency doublet for probes in Nb sites, ν Q1 = 353 MHz and ν Q2 = 494 MHz, was necessary for a good fit of the PAC spectra [193]. However, the best fits for the 6 mol% congruent and 0.2 mol% near stoichiometric crystals, for which the fraction of probes in Nb sites was close to 50%, were obtained without a fraction of probes in Li II sites. This suggests that the defect associated Li II site has disappeared or is of too small importance to make an effect.
The influence of stoichiometry on the incorporation of Hf into LN was studied by RBS/C on samples grown from stoichiometric melt [138,193]. The sample with higher Hf concentration (1.0 mol%) shows clearly a mixture of Li-and Nb-lattice sites occupied by Hf. Computer simulations with FLUX yield a best fit for 66% Nb sites and 34% Li sites, i.e., roughly a 2:1 ratio. This deviates from the 3:1 ratio expected in the case of charge compensation during the incorporation of Hf into LN. From angular scans of 0.2 mol% HfO 2 sample, which can be regarded closer to a stoichiometric case, a fraction of more than 85% of Hf is located on the Li site. It can be concluded that Hf, despite being a tetravalent ion, shows a strong preference to incorporation on Li sites.
In order to demonstrate the influence of co-doping in LN on the incorporation mechanisms in LN, the case of Hf-Mg will be discussed in detail. A set of single crystals of congruent LN doped with 1 mol% of HfO 2 in melt and MgO concentrations ranging between 2.0 and 4.0 mol% in steps of 0.5 mol% as well as one with 6.0 mol% have been studied by RBS/C and PAC [194][195][196][197]. The evaluation of the first set of co-doped samples showed that the fraction of Hf located on Li sites decreases constantly with increasing MgO concentration from about 90% at 2.0 mol% MgO to about 20% at 4.0 mol% and changes to full occupation of Nb sites at 6.0 mol%. This reveals that the tendency of Mg to occupy Li sites is stronger than that of Hf but falls short on being a charge compensation effect. Unfortunately, only for the case of 6 mol% MgO data on the lattice site occupation by Mg are available from an additional study [198]. In this case the mixed occupation of regular Li and Nb-octahedron sites as in the single Mg-doped case has been confirmed indicating that Mg may keep its incorporation mechanisms despite co-doping.
PAC measurements were performed in the same crystals after irradiation with neutrons and annealing. The spectrum obtained in the 6 mol% co-doped crystal could only be fitted considering a doublet ν Q1 = 327(10) MHz and ν Q2 = 398(10) MHz for probes in Nb sites. This frequency doublet is slightly different from the one observed for Hf probes in Nb sites of near-stoichiometric crystals, although no shift was observed for probes in Li sites, ν Q1 = 1154(12) MHz and ν Q2 = 1213(12) MHz. Figure 24 shows the dependence of the fractions of Hf on the different Li and Nb sites on the MgO concentration. The fraction of Hf occupying the Nb antisite (Li-II site) is constant up to 2.5 mol% MgO in the crystal. In this concentration range, Mg is exclusively located at the center of the Li octahedron. It therefore appears that Mg substitutes in Li Li or vacant Li sites, decreasing the number of Li sites available for Hf. Further, full suppression of Hf on the Li-II occurs between 4 and 6 mol%, the region in which the suppression of the photorefractive damage is completed.
For indium only one PIXE/C study for co-doping with 6 mol% MgO in melt is available [199]. The results indicate that In mainly substitutes Li in its regular lattice site but a fraction occupying regular Nb sites cannot be ruled completely.
concentration. The fraction of Hf occupying the Nb antisite (Li-II site) is constant up to 2.5 mol% MgO in the crystal. In this concentration range, Mg is exclusively located at the center of the Li octahedron. It therefore appears that Mg substitutes in LiLi or vacant Li sites, decreasing the number of Li sites available for Hf. Further, full suppression of Hf on the Li-II occurs between 4 and 6 mol%, the region in which the suppression of the photorefractive damage is completed. For indium only one PIXE/C study for co-doping with 6 mol% MgO in melt is available [199]. The results indicate that In mainly substitutes Li in its regular lattice site but a fraction occupying regular Nb sites cannot be ruled completely.

Dopants to Enhance Photorefractive Damage
Despite the photorefractive effect being undesired for most optical applications it has an interesting application in the field of holographic recording of data [200]. Therefore also dopants that increase this effect have been the target of intense research. So far Fe has been identified as the single dopant with the most promising properties but other single dopants as Mn, Co, Cu, Ni, Cr, Mo, Ru, Rh, Ce and Tb [201][202][203] were investigated. Moreover, combinations of dopants, e.g., LiNbO3:Fe,Mn [204] or LiNbO3:Zr,Cu,Ce [205] were investigated intensively.
In contrast to the suppressing dopants the incorporation of photorefractive elements has gotten rather little attention by nuclear methods. Studies exist on the technically most important dopant, iron, on cobalt and on titanium. The latter enhances the photorefractive effect only very slightly but is of importance for optical waveguide applications.

Dopants to Enhance Photorefractive Damage
Despite the photorefractive effect being undesired for most optical applications it has an interesting application in the field of holographic recording of data [200]. Therefore also dopants that increase this effect have been the target of intense research. So far Fe has been identified as the single dopant with the most promising properties but other single dopants as Mn, Co, Cu, Ni, Cr, Mo, Ru, Rh, Ce and Tb [201][202][203] were investigated. Moreover, combinations of dopants, e.g., LiNbO 3 :Fe,Mn [204] or LiNbO 3 :Zr,Cu,Ce [205] were investigated intensively.
In contrast to the suppressing dopants the incorporation of photorefractive elements has gotten rather little attention by nuclear methods. Studies exist on the technically most important dopant, iron, on cobalt and on titanium. The latter enhances the photorefractive effect only very slightly but is of importance for optical waveguide applications.
It is generally presumed that the charge state of iron (Fe 2+ or Fe 3+ ) is of importance for the photorefractive efficiency. Various teams studied this topic applying Mössbauer Spectroscopy with 57 Co as radiation source. The first study by Vayashko et al. [206] used lithium deficient LN heavily Fe-doped (2.4 to 5.1 mol%) crystals and detected both charge states in the material. They also reported that the fraction of Fe 2+ increases with increasing Li-deficiency. Further studies with diffusion-doped [43] (0.5 to 1.0 mol% of Fe enriched to 86% in 57 Fe) and melt-doped [207,208] (0.22 mol% Fe 2 O 3 enriched to 90.7% in 57 Fe) by and confirmed that a major fraction (90%) of the iron is present as Fe 3+ in the virgin material. Annealing under reducing atmosphere (in Ar at 1000 • C for 140 h or under vacuum at 635 • C for 24 h), changes the charge state of the majority of the dopant to Fe 2+ . The measurements also indicate that the Fe dopant is located on an axial site [209].
Additional Mössbauer studies addressed the electron trapping in congruent, nearstoichiometric, Mg-doped and Fe-doped LN [210,211]. It was found that the electron trapping effect by Nb-antisite outweighs that attributed to Fe 3+ despite the larger cross section of the latter. For crystals doped with Mg up to threshold concentration the trapping is strongly reduced. Due to Li-loss during 57 Co diffusion at high temperature no quantitative results could be obtained for the near-stoichiometric material.
A first study on the lattice location of Fe was performed by Rebouta et al. [212] using samples of congruent LN doped in melt with 0.5 mol% Fe 2 O 3 . Due to the growth method, Fe was mainly incorporated as Fe 2+ (as was confirmed by optical absorption spectrometry) rather than in the form of Fe 3+ preferred in applications using the photorefractive effect However, previous studies suggest that the lattice location should not be affected by this [213]. PIXE/C with protons was used to determine the lattice location. Although the use of protons produces less pronounced dips and peaks in the angular scans the results definitely show that Fe is located in the Li site.
Zaltron et al. [214] presented a PIXE/C study on Fe-diffused samples. The material was produced by coating the LN substrate with a metallic Fe film and performing a diffusion treatment for 5 h at 1000 • C in oxygen atmosphere. Since the resulting samples contained Fe 3+ , an additional treatment under reducing atmosphere (96% Ar + 4% H 2 ) for different times at 500 • C was performed in order to obtain up to 100% Fe 2+ . The depth profiles of Fe before and after reduction were measured by SIMS revealing a distribution extending up µm into the bulk with a maximum concentration of about 0.9 mol% Fe at the surface. The distribution remained unaltered after reduction. Angular scans measured for the as-diffused (1% Fe 2+ ) and two reduced (20% and 100% Fe 2+ ) samples show fully overlapping Fe-PIXE curves. Despite the poor crystalline quality of the samples (Nb-RBS minimum yield of ≈45% in the 0221 -axis) the data indicate that Fe is located on the Li site. This confirms also the assumption that the lattice site does not depend iron being incorporated in the divalent or trivalent state.
The lattice site of Co in LN was studied by Szilagy et al. [215] in melt-doped (0.2 mol% CoO) congruent material. In this case, a He + ion beam was used to obtain PIXE/C scans for Co-and Nb-PIXE with pronounced dips for five major axial directions. The simultaneous recording of Nb-RBS axial scans allowed to demonstrate the high quality of the crystals grown (χ min ≈ 2%). Comparison of the experimental results with simulations using CASSIS demonstrated that Co is located exclusively on the regular Li site.
In contrast to its homologues Zr and Hf, Ti does not suppress the photorefractive effect. It is therefore of interest if its behavior is linked to some different incorporation mechanism. A first study was made by Buchal et al. [216] with a sample implanted with Ti at a fluence of 2.5 × 10 17 cm −2 and 360 keV energy and subjected to a heat treatment at 1000 • C. This caused an almost complete recovery of the crystal structure and a diffusion of Ti in the bulk material. The maximum concentration of Ti in the samples was reported to 12 mol%. The PIXE/C measurements for Ti-K characteristic X-ray showed an almost perfect overlap with that of Nb leading to the conclusion that Ti occupies the Nb lattice site. Kollewe et al. [217] studied two different cases: melt-doped material with 3 mol% Ti and Ti-indiffused waveguides with a peak concentration of about 9 mol%. PIXE/C scans were recorded for various axial directions. For measurements along the 0001 -axis both samples evidence the absence of off-axis lattice sites while for the other axial directions clear differences are visible. While for the melt-doped crystal an almost complete overlap of the Ti-PIXE scan with those of the Li-NRA is found, the indiffused sample shows a behavior that resembles mainly the Nb-PIXE signal. Finally, Hauer et al. [218] reported a combined PIXE/C and PAC study for LN doped in melt with 1 mol% of Ti. The angular PIXE/C scans for various axial directions confirm clearly that Ti occupies only Li sites at this concentration level and that no off-axis sites are occupied. The PAC measurements show a single QIF, ν Q = 15.4(2) MHz, which scales up very well with the results obtained with other probes occupying Li sites [34]. These findings led to the conclusion that Ti occupies solely Li-sites for lower concentrations while for higher concentrations a mixture of Li-and Nb-sites has to be assumed. Ti and Hf have obviously similar tendencies in the occupation of lattice sites although with opposite influence on the photorefractive index.

Optically Active Dopants
The excellent linear and non-linear properties turn LN into a highly desirable host material for optically active ions. Coherent stimulated emission from the rare earths Nd, Tm and Ho was observed in this material already in 1969 [219]. The interest on the laser applications of rare earths and Cr has led to innumerous studies on their optical properties in the last decades. Since these are strongly dependent on the dopant's crystallographic environment the determination of its lattice site is of outmost interest. Although some pertinent information can be obtained by spectroscopic methods the direct determination of the lattice location by ion beam methods is indispensable. The study of rare earths in LN Crystals 2021, 11, 501 32 of 55 by RBS/C is fostered by their larger mass in comparison to Nb as the heaviest component of the host material while chromium has to be investigated using PIXE/C.
For optical communications applications Er 3+ has attracted interest in a wide range of laser host materials due to its IR emission at about 1.5 µm coinciding with a minimum in the absorption in optical fibers. Erbium has been incorporated into LN by melt-doping, indiffusion and ion-implantation. RBS/C studies have been reported only for melt-doped LN. The first investigation used two crystals co-doped with 0.5 mol% Er 2 O 3 and 5.8 and 8.0 mol% MgO in melt (0.79 mol% Er, 4.8 and 7.5 mol% Mg in crystal), respectively [220,221]. Angular scans for axial directions show an exclusive occupation of the regular Li lattice site. For a second study, a series of singly-doped crystals with concentrations of 0.5, 1, 2, 3, 4 and 5 mol% in a congruent melt was grown and investigated [222]. The comparison of the experimental angular scans with simulations using the modified FLUX program confirmed that the main fraction Er is located in the Li-Octahedron but appears to be shifted form the regular site by 0.2 Å towards the closest oxygen layer. The Er concentration of these sites increases at first with increasing total Er concentration but reaches a saturation value at about 2.5 mol%. The high minimum yields for Er in axial directions point out that an additional fraction of Er occupies non-axial sites but cannot be identified with the known tetrahedral site in LN. The authors claim that these sites off the c-axis cannot be seen as a new crystalline structure but are rather amorphous precipitates or defect clusters.
Another important rare earth in optical communication application is praseodymium. Crystals doped with 0.2 and 0.5 mol% Pr 3+ in melt were studied by RBS/C [223][224][225]. The comparison of experimental scans and simulations with the FLUX code shows that as in the case of erbium Pr is located in the Li oxygen-cage but with a larger shift (0.45 Å) towards the closest oxygen layer. On the other hand, site selective spectroscopy performed on the same samples suggests that more than one site for Pr 3+ exists. This discrepancy can be easily understood taking into account the defect structure of the congruent material that causes (due to vacancies and stacking faults) different local surroundings for the Pr 3+ ion.
Also the lattice site of the probably most popular laser active ion, Nd, has been investigated using RBS/C. The history of the investigation of this impurity highlights the importance of computer simulations for the interpretation of lattice location measurements in complex non-cubic crystal structures like LN. The first studies of the Nd 3+ lattice sites for a 1 mol% Nd 3+ in melt (0.33 mol% in the crystal) where first interpreted based on assumptions adopted from the study of simpler crystal lattice and pointed out a mixed occupation of Li-and Nb-sites [192]. Another study even claimed the occupation of three different axial lattice sites within the Li-Octahedron [226,227]. A new interpretation of the experimental results using the FLUX code revealed that Nd 3+ is actually located in the Li-Octahedron displaced by 0.4 Å towards the closest oxygen layer [23].
In order to complete the overview on the lanthanides it has to be mentioned that three other rare earths were also found to occupy positions shifted from the Li site towards the closest oxygen plane: Eu (0.4 Å) [23,192,228], Ho (0.38 Å) [223,224] and Yb (0.3 Å) [223,224]. Plotting the five known displacement values against the ionic radius as suggested by Lorenzo et al. [223,224], shows an apparent linear dependence, depicted in Figure 25. The measured position of the main fraction of Er [222] (0.2 Å) has been added to the Lorenzo's data.
For Nd, Gd, Er, Tm, and Lu (1 mol% in melt) RBS/C studies on samples co-doped with Mg (6 mol%) have been performed [199,221]. For Nd, Gd and Tm the angular scans show only a partial occupation of the Li octahedron. In contrast to the singly-doped crystals, remarkably large fractions of randomly incorporated ions (almost complete for Nd, 70% for Gd, 20% for Tm) are observed. Rare earths with low or no random fractions are Er and Lu. Interestingly, also the distribution coefficient of Nd in the MgO co-doped crystal, k = 0.03, is an order of magnitude lower than in the singly-doped material [192]. This indicates that the incorporation of lanthanides with large ionic radius into the LN lattice is strongly hampered by Mg. The lattice location study for Mg on the LiNbO 3 :Lu,Mg sample by Kling et al. [198] demonstrates that, in contrast to Hf-co-doping, the distribution of Mg in the lattice is identical to the singly 6 mol% MgO-doped crystal. In the competition for the occupation of sites in the Li-Octahedron the divalent Mg may be favored with regard to the trivalent lanthanides since it requires less charge compensation by formation of additional Li vacancies. mixed occupation of Li-and Nb-sites [192]. Another study even claimed the occupation of three different axial lattice sites within the Li-octahedron [226,227]. A new interpretation of the experimental results using the FLUX code revealed that Nd 3+ is actually located in the Li-octahedron displaced by 0.4 Å towards the closest oxygen layer [23].
In order to complete the overview on the lanthanides it has to be mentioned that three other rare earths were also found to occupy positions shifted from the Li site towards the closest oxygen plane: Eu (0.4 Å) [23,192,228], Ho (0.38 Å) [223,224] and Yb (0.3 Å) [223,224]. Plotting the five known displacement values against the ionic radius as suggested by Lorenzo et al. [223,224], shows an apparent linear dependence, depicted in Figure 25. The measured position of the main fraction of Er [222] (0.2 Å) has been added to the Lorenzo's data. For Nd, Gd, Er, Tm, and Lu (1 mol% in melt) RBS/C studies on samples co-doped with Mg (6 mol%) have been performed [199,221]. For Nd, Gd and Tm the angular scans show only a partial occupation of the Li octahedron. In contrast to the singly-doped crystals, remarkably large fractions of randomly incorporated ions (almost complete for Nd, 70% for Gd, 20% for Tm) are observed. Rare earths with low or no random fractions are Er and Lu. Interestingly, also the distribution coefficient of Nd in the MgO co-doped crystal, k = 0.03, is an order of magnitude lower than in the singly-doped material [192]. This indicates that the incorporation of lanthanides with large ionic radius into the LN lattice is strongly hampered by Mg. The lattice location study for Mg on the LiNbO3:Lu,Mg sample by Kling et al. [198] demonstrates that, in contrast to Hf-co-doping, the distribution of Mg in the lattice is identical to the singly 6 mol% MgO-doped crystal. In the competition for the occupation of sites in the Li-octahedron the divalent Mg may be favored with regard to the trivalent lanthanides since it requires less charge compensation by formation of additional Li vacancies. Trivalent chromium has been already used in the very first laser as optically active ion. The potential application as a compact tunable laser initiated an intense research on the properties of Cr-doped LN. The introduction of chromium in LN with chromium can be achieved by melt doping [229,230], ion implantation [231] and diffusion [232]. Since chromium as a transition metal shows the tendency to increase the photorefractive effect the co-doping with Mg at high concentration is desired. However, the co-doping alters Trivalent chromium has been already used in the very first laser as optically active ion. The potential application as a compact tunable laser initiated an intense research on the properties of Cr-doped LN. The introduction of chromium in LN with chromium can be achieved by melt doping [229,230], ion implantation [231] and diffusion [232]. Since chromium as a transition metal shows the tendency to increase the photorefractive effect the co-doping with Mg at high concentration is desired. However, the co-doping alters the absorption (Cr-doped LN is green, Cr-Mg-co-doped pink) and luminescence properties. Optical [233], ESR [234] and EXAFS [235] studies suggest that these effects are related to changes in the lattice site occupied by the Cr-only Li-sites for singly-doped material, Liand Nb-sites for Mg-co-doped.
The lattice location of Cr was studied by PIXE/C in a congruent sample doped with 0.1 mol% Cr and one with 6 mol% Mg and 0.1 mol% Cr in melt [236]. The very low Cr dopant concentration turned to out to be extremely challenging for the experiment and a copper monolayer had to be deposited on the samples to suppress low energy X-ray emission due to charging effects by the ion beam. Nevertheless, it was possible to determine the lattice location by comparing the angular scans obtained with computer simulations using CASSIS. Surprisingly, it was found that already the singly doped crystal contains about 40% of the Cr on niobium sites and 60% on Li sites. This implies that chromium, in contrast to all other trivalent dopants studied so far, shows a tendency for charge compensation Cr 3+ Li -Cr 3+ Nb replacing Li 1+ -Nb 5+ ). The results align well with those of the ESR study by Jaque et al. [234] suggesting that in singly-doped material mainly Cr 3+ Li -Cr 3+ Nb occur with some Cr 3+ on Li-sites. For the case of the co-doped material the fraction of Cr located on Li sites is reduced to about 20% while 80% occupy the Nb octahedron although shifted about 0.1 Å towards the center of the octahedron. Waveguide lasers require only a thin near-surface with the optically active ion, instead of bulk-doping. Ion implantation is an advantageous tool to dope a material with a controllable concentration and depth distribution. In order to verify the applicability of this method y-cut LN samples were implanted with 100 keV Cr to fluences of 1 × 10 16 and 5 × 10 16 cm −2 , respectively [115]. This resulted in a peak concentration of 2.2 mol% Cr for the latter fluence with a projected range of about 55 nm. In order to remove the implantation induced lattice damage annealing was performed under wet oxygen atmosphere for various temperatures and times. The success of the annealing process was verified by RBS/C showing that for 1000 • C and 4 h the almost completely amorphized surface layer was reconstructed and the RBS/C spectra hardly distinguishable from those for virgin material. In order to verify the incorporation of chromium in the crystal lattice angular scans were taken with PIXE/C. An almost perfect overlap with the results for the melt-doped case demonstrate the successful incorporation of Cr into the LN crystal.

Pentavalent and Hexavalent Dopants
Dopants in the pentavalent state are of interest since they can replace Nb without need of charge compensation. However, none of the elements in the pnictogen group incorporates into the LN lattice. An alternative are the homologues of Nb, V and Ta. The latter is known to be able to form mixed crystals covering the full stoichiometry range [237].
Since Ta has about double the atomic mass of Nb it was easy to study this impurity by RBS/C. The results of congruent crystal doped with 1 mol% of Ta showed that the angular scans for Nb and Ta fully overlap for all investigated axial and planar direction evidencing that Ta fully replaces Nb on its lattice site [228].
Similar to the case of pnictogens, also no doping of LN with chalcogenides could be achieved. Nevertheless, three elements that readily form hexavalent oxides-Mo [238], W [239][240][241] and U [242] have been successfully incorporated into LN. The lattice location of latter two was subject to studies by RBS/C.
Tungsten doping was performed in three different ways: ion implantation, doping in melt and diffusion. In the first case a fluence of 1 × 10 14 cm −2 of 186 W was implanted into x-cut congruent LN resulting in maximum W concentration of 1.25 mol% peaking at about 50 nm depth [239]. In order to remove the implantation induced damage the sample was subjected to a two-step annealing process (600 • C for 1 h, 800 • C for 4 h) in a wet flowing oxygen atmosphere. The angular scans of the recrystallized material showed a complete overlap of the Nb-and W signals for all crystallographic directions proving the complete incorporation of tungsten on the niobium site. In a second study, LN crystals were grown from melt doped with 1 and 2.5 mol% WO 3 [240]. The concentration in the grown samples were determined by RBS to be 0.2 and 0.5 mol%, respectively, fixing the distribution coefficient of tungsten to 0.2. The RBS/C measurements confirmed the results of the ion implantation work on W replacing solely Nb. In the third approach a y-cut LN crystal was placed on a bed of metallic tungsten and kept at 1000 • C for 2.5 h under vacuum which yielded an extremely W-rich (about 35 mol%) layer at the surface [241]. A subsequent annealing treatment at 1000 • C for 2 h reduced the W concentration to 0.75 mol% and improved the crystal quality sufficiently to perform RBS/C scans. As in the previous two cases, the exclusive occupation of Nb sites by W was observed.
For the investigation of the lattice site of uranium a crystal grown from a congruent melt doped with uranylacetate-dihydrate (0.5 at.% U equivalent to roughly 2.5 mol%) was used, in which uranium should be incorporated in the hexavalent state. The samples obtained contained about 0.18 mol% U (distribution coefficient ≈ 0.07) as determined by RBS. Optical absorption spectroscopy confirmed that uranium had been incorporated as U 6+ [242]. As in the case of W the angular scans for U-RBS and W-RBS overlapped fully confirming also for uranium the occupation of the Nb lattice site [240]. Calvo et al. [242] observed that a heat treatment under vacuum, i.e., reducing conditions for LN, turns the U-doped crystals from orange to red. Optical absorption studies revealed the existence of trivalent (and eventually tetravalent) uranium after this treatment. In order to check if this alteration of the valence state has any influence on the lattice site of uranium RBS/C with a sample reduced for 16 h at 1000 • C under high vacuum were performed. The angular scans, shown in Figure 26 for various crystallographic directions, show that part of the U has been pushed out of the Nb site but is still on an axial position [243]. Employing simulations using CASSIS suggest that about 20% of the U is now on Li sites while the remainder still occupies the Nb site. These findings corroborate the assumption that the cation lattice site occupied by ad dopant is mainly governed by its valence state. scans, shown in Figure 26 for various crystallographic directions, show that part of the U has been pushed out of the Nb site but is still on an axial position [243]. Employing simulations using CASSIS suggest that about 20% of the U is now on Li sites while the remainder still occupies the Nb site. These findings corroborate the assumption that the cation lattice site occupied by ad dopant is mainly governed by its valence state.

Noble Metals
The introduction of noble metals into glass in order to change its color is a method applied since Roman times [244]. The coloration is due the formation of noble metal nanoparticles for which the existence of surface plasmon resonances causes strong optical absorption peaks. The associated non-linear optical properties raised the interest on formation of such nanoparticles in modern optical materials. LN can incorporate noble metals as Cu and Au only at very low concentrations in melt-doped crystals causing the Figure 26. Angular RBS/C scans for the 0221 -axis for Nb and U in the as-grown and reduced sample. In the reduced case the narrowing of the uranium-dip near the axial direction is due to the partial occupancy of Li-sites by U.

Noble Metals
The introduction of noble metals into glass in order to change its color is a method applied since Roman times [244]. The coloration is due the formation of noble metal nanoparticles for which the existence of surface plasmon resonances causes strong optical absorption peaks. The associated non-linear optical properties raised the interest on formation of such nanoparticles in modern optical materials. LN can incorporate noble metals as Cu and Au only at very low concentrations in melt-doped crystals causing the absence of surface plasmon related effects. Since the first successful formation of Cunanoparticles and the observation of their typical optical absorption behavior by Saito et al. [245] in 1987, this method has been applied for the production of Cu [246], Ag [247], Au [248] and mixed [249] nanoparticles. The influence of ion energy (20 keV-3 MeV), fluence (4 × 10 13 -2 × 10 17 cm −2 ), implantation temperature (room temperature-700 K), post-implantation treatment (annealing at various temperatures and atmospheres) on the optical properties has been studied in detail [178][179][180][181]. RBS and RBS/C was employed in several studies to investigate the distribution of the implanted species, the implantationinduced damage and its removal by annealing as well as some structural characteristics of the implanted elements.
The best investigated case is that of Au. Due its high mass it can be easily detected by RBS and the implantation-related peak in the spectrum is well separated from niobium. Shi et al. [250] demonstrated that an implantation at 1 MeV with a fluence of 4 × 10 13 Au + cm −2 only partially damages the implanted layer while for 5 × 10 14 Au + cm −2 already full amorphization is achieved. Taking RBS/C spectra they were able to calculate damage depth profiles with nanometric resolution. A comparison with the popular TRIM code [251] for damage profile computation revealed that the experimentally determined damage reaches far deeper into the crystal than computed slowing down recovery of the crystal lattice during annealing. The study by Takahira et al. [252] was the first one to investigate the recrystallization of an Au-implanted (550 keV, 1 × 10 16 to 1 × 10 17 cm −2 ) LN crystals in detail. They showed that annealing at 600 • C for 1 h in air efficiently removes a large part of the damage for the sample implanted with 3.7 × 10 16 Au + cm −2 -corresponding to a peak concentration of 7 at.%. It is noteworthy that the Au implantation profile is not affected by the heat treatment.
Dev et al. [253] performed RBS/C measurements for the 0001 -axis of LN implanted with 1 × 10 15 cm −2 Au (1.43 MeV) and annealed at 450 • C in flowing oxygen atmosphere for 40 min. The curves for the Nb-and Au-RBS overlap mostly but show a distinct kink at about 0.5 • from the axial direction, ruling out a formation of nanoclusters due to a low Au concentration (0.5 mol%). This indicates that a major part of the Au is in axial lattice site but in contrast to the claims of the authors, not necessarily the Nb-site. Due to its high atomic mass the Au-RBS scan would also overlap for Li-, free octahedron or other aligned sites. The kink-related fraction is definitely related to a fraction on a non-axial site, possibly the tetrahedral site (although this would have to be verified by additional studies).
For a higher fluence (1 × 10 16 Au + cm −2 , 800 keV, implanted at room temperature) Kling et al. [254] investigated the recrystallization of the sample after subsequent annealing at 600 • C (2 h) and 800 • C (1 h) in flowing wet oxygen. RBS/C measurements showed that even at 800 • C only a partial recovery of the crystal lattice was achieved. Further, it was found that the yield for Au-RBS is reduced for the 0110 -axial direction. Since optical absorption demonstrated the formation of Au-nanoparticles in the sample an incorporation into the LN lattice was ruled out leaving only the hypothesis that the Au-cluster are aligned the with the LN lattice. Complementary studies by TEM of this sample revealed that gold did not form spherical nanoparticles as deduced from various optical investigations but thin flakes parallel to the surface with lateral dimension of up to 200 nm. An electron diffraction measurement confirmed that these flakes are highly coherent with the LN lattice.
In another study by the same authors, the interaction between Au-and He-implants was in the focus of the interest [255]. All samples were implanted with 20 keV He + to a fluence of 1 × 10 16 cm −2 . Au + -implantations with 50 keV (5 × 10 15 cm −2 ) were used to obtain two well-separated profiles while others using 800 keV Au + (1 × 10 16 cm −2 ) achieved an overlap with the profile of He. In RBS/C spectra for the samples with separate profiles the two damage-related peaks for He and Au appear separately even after annealing at 270 • C. Annealing at 600 • C causes the Au to migrate towards the bulk and broadens the profile while at 800 • C a movement of the gold towards the surface is observed. This mobility of Au contrast clearly with the static Au profiles observed for single implantation at the same temperature discussed above. Most likely the diffusion of Au at 600 • C is facilitated by the presence of He in the material but not showing a real gettering effect as in Si [256]. As in the previous case, diffusion towards the surface and an alignment of the Au-nanoparticles was found for the samples with overlapping profiles. TEM observations revealed again the existence of Au-flakes but with drastically reduced lateral sizes (<100 nm) compared to singly-implanted LN opening a path for the control of nanoparticle properties by co-implantation of noble gases.
For the case of Ag-implantation, Milz et al. [257] studied the influence of implantation and annealing temperatures on the damage production and removal. X-cut LN was implanted with 380 keV Ag + to a high fluence (1 × 10 17 cm −2 ) at liquid nitrogen (−196 • C) temperature, room temperature and 700 K (423 • C). RBS/C spectra show that in the first two cases the implanted surface layer was completely amorphized while in the hightemperature case the layer was heavily damaged and thinner retaining some crystallinity. Rapid thermal annealing at temperatures ranging from 300-900 • C was employed for recrystallization. In the sample implanted at room temperature an increasing reduction of the amorphized layer with increasing temperature was observed but only at 900 • C a partial recovery of the crystal lattice at the surface was observed.
The interaction of Ag and O implants was the topic of a study by Williams et al. [258]. Z-cut samples were implanted with Ag and O at energies of 160 and 35 keV and fluences of 4 × 10 16 cm −2 and 8 × 10 16 cm −2 , respectively at RT and 500 • C. Both implantation sequences were used and additionally a singly Ag + implanted samples (1.5 MeV, 2 × 10 16 cm −2 and 1.7 × 10 17 cm −2 , same temperatures) were produced. The implanted material was annealed for 1 h at 500 • C in air. From RBS/C spectra recorded for the Ag-implanted sample it was concluded that in case of high-temperature implantation is located closer to the surface than in that for room temperature. On the other hand, during the annealing the Ag moves towards the bulk in the high-temperature sample while it approaches the surface in the RT-sample. Co-implanted material shows some difference between the two implantation sequences when the high-temperature procedure is applied. For Ag-then O-implantation the Ag-profile stays practically unaltered with most of the Ag remaining at the surface.
For the opposite sequence, a significant Ag-diffusion towards the bulk was observed that seems to stabilize the Ag-nanoparticles formed.
The incorporation of Pt into LN was achieved by 50 keV Pt + ion implantation at room temperature at a fluence of 5 × 10 15 cm −2 [259]. A maximum concentration of about 1.5 mol% of Pt was attained. The samples were annealed under flowing wet oxygen atmosphere at 600 • C and 800 • C for 1-6 h in steps of one hour or at 1000 • C for 60, 90 and 120 min. Further, an incremental annealing was done at 300 • C, 400 • C, 450 • C, 500 • C, 600 • C and 700 • C for 30 min followed by one at 1000 • C for up to 5 h. Studying the Pt-RBS yield in the surface region (0-10 nm) showed no effect for 600 • C, a slight diffusion towards the bulk after 6 h at 800 • C and a very rapid one at 1000 • C leaving the surface region virtually without Pt. RBS/C spectra showed no alignment effect for Pt except for the sample that underwent incremental annealing. A minimum yield of 53% is observed for alignment with 0110 -axis while for all other crystallographic directions random behavior is found. It can be concluded that Pt forms precipitates that are preferentially aligned with the 0110 -axis of LN as in the case of Au.
Iridium implanted into LN (130 keV, 5 × 10 16 cm −2 ) with a high peak concentration has been studied by Kling et al. [260] The samples were annealed under flowing wet oxygen at 500 • C, 800 • C and 1000 • C for 30 min. In contrast to Pt, significant diffusion towards the bulk was found by RBS already at 500 • C. After heat treatment at 1000 • C a practically uniform concentration of 0.015 mol%, i.e., about two orders of magnitude lower than in the as-implanted material, is reached. RBS/C scans show an alignment effect with the LN for all investigated axial directions with a width comparable to that of Nb-RBS. However, the high minimum yields (40-80%) do not allow attributing positively an occupation of the Nb-site by Ir.

Optical Waveguides
Photonic applications using optical materials like LN receive a constantly increasing interest. A key part in any photonic device is an optical waveguide. It consists of a thin surface layer (in the order of the wavelength of the light guided) with a refractive index larger than the underlying substrate. Total reflection at the interfaces with the bulk and the air confines the light propagation to the waveguide. Many different approaches for the fabrication of waveguides in LN and other optical materials have been developed in last decades. For a comprehensive review see Bazzan and Sada [261].

Indiffusion
The formation of LN waveguides by in-diffusion of Ti from a thin metallic layer deposited on the surface was reported for the first time by Kaminow and Carruthers [262] in 1973 and became with time the most common method. It combines a good guiding behavior with the preservation of the electro-optic properties. The downside is that Ti, in contrast to its homologues Zr and Hf, it does not suppress the photo-refractive effect.
For applications, it is of important to determine the depth profile of the in-diffused Ti that controls the refractive index profile of the optical waveguide. The method of choice is in this case the SIMS. The first Ti-profile measurement was performed in 1979 by Burns et al. [263]. The samples were prepared from x-cut LN coated with 17.5 nm of Ti and diffusing at 900-1100 • C for 6h in Ar atmosphere. SIMS concentration profiles for Ti and Li ions for a sample treated at 1000 • C show pronounced peaks with half-width of 200-300 nm near the sample surface. A clear correlation between decreasing Ti-and increasing Li-concentration with depth is observed. These results suggest that the diffusion treatment time was too short for a sufficient redistribution of the Ti. Arizmendi et al. [264,265] investigated in detail this first phase of the diffusion process combining ion beam analysis, X-ray diffraction and electron microscopy techniques for samples indiffused in Ar and O 2 atmospheres between 500 and 950 • C. RBS/C measurements for samples annealed at 500 • C showed the formation of an amorphous TiO 2 -layer on the substrate and a lack of oxygen in the LN surface. For the technically more important case of 950 • C (30 min, dry oxygen) a Ti x Nb 1-x O 2 compound with x = 0.65 in the surface region is inferred from the RBS results in accordance with x-ray diffraction. Angular scans for several axial directions revealed that this new compound is aligned along the c-axis of the LN substrate but not with other axial directions. The authors argue that the compound formed in this step constitutes the actual source for the further in-diffusion of Ti.
Based on these observations for technical applications the diffusion treatment at high temperature and for extended time periods (up to 16 h) has been adopted. Several groups reported SIMS studies on the Ti-diffusion process. Bremer et al. [266] reported a Gaussian profile peaking at the surface for Ti after indiffusion for 16 h under wet Ar atmosphere. A vast investigation of the main determining factors that govern Ti-diffusion (Ti-film thickness, diffusion temperature) has been presented by Caccavale et al. [267,268]. With the help of SIMS measurements on samples annealed at 900-1000 • C they could confirm that only for the highest temperature a complete indiffusion of Ti takes place while for lower temperatures an enhanced Ti-concentration near the surface is found. Further, they were able to detect a significant anisotropy in the diffusion process leading to a diffusion coefficient twice as high for z-cut material than for x-cut. Caccavale et al. [269] have given an impressive example of the lateral resolution of the SIMS technique by studying the lateral diffusion of Ti in waveguides. It was possible to determine the lateral Ti-profiles of two waveguides separated only by a few µm. For near-stoichiometric LN produced by the vapor transport equilibration Zhang et al. [270] observed a diffusion profile that can rather be described by an error-function than by a Gaussian as in the congruent material.
To overcome the photorefractive effect in Ti-diffused waveguides diffusion into Mgdoped material or Mg co-diffusion can be applied. The first SIMS study by Bremer et al. [271] used already Ti-diffused waveguides as described above [266]. MgO was deposited as Mg-source on the waveguide and indiffused at 950 • C in order to avoid further Ti-diffusion and taking into account the more rapid diffusion of Mg. A single Mgdiffusion treatment leads to a Gaussian distribution of the Mg peaking at the surface with some reduction of the Li content at the waveguide surface. Repetition of the Mg-indiffusion causes a drastic (20%) loss of Li in the surface region accompanied by high Mg-peak in this region. The authors interpreted this as hint to the substitution of Li by Mg during doping. Caccavale et al. [272] concluded that the Ti-profile suffers only minor additional diffusion by Mg co-diffusion but that on the other hand Mg-diffusion is clearly enhanced in comparison to an undoped crystal [273].
A comparative study on the Ti diffusion in undoped and Mg-melt-doped congruent and near-stoichiometric LN was presented by Kumar et al. [274]. Samples were indiffused at 1080 • C for 23 h in dry air atmosphere. The values of the Ti-diffusion constant determined using SIMS profiles show that Ti-diffusion is slowed down in MgO-doped material for near-stoichiometric (1 mol% MgO) and congruent material (5 mol% MgO) with respect to their undoped counterparts.
An important alternative is the indiffusion of Zn from vapor phase into LN, reported for the first time by Herreros and Lifante [275]. Nevado et al. [276,277] produced samples indiffused at 550 • C and 700 • C with a diffusion times of 2 h and 8 h, respectively. The samples were subsequently annealed for 4 h at the 700 • C or 800 • C, respectively, to form the actual waveguide. As-diffused and annealed samples were investigated by SIMS and RBS. In both cases the as-diffused samples showed Zn-rich surface layers followed by a zone with decreasing Zn concentration. After the annealing Zn diffused from the Zn-rich superficial layer several µm deep into the bulk forming the actual waveguide. A further quantitative analysis was performed using HIERD [278]. This high depth-resolution quantitative technique allowed to confirm that in the as-diffused case the top layer consists of LN and ZnNb 2 O 6 phases, i.e., a partial substitution of Li by Zn occurs. On the other hand, in the annealed sample a homogenous mixture of ZnO and LN was found on top and a Zn-diffusion profile for the remainder of the waveguide.
For completeness, it should be mentioned, that also other elements, Sc [279] and Zr [280], which are candidates for waveguide formation were studied with regard to their diffusion properties using SIMS. Further, several studies focused on the co-diffusion of two elements: Zn/Ni [281], Zr/Ti [282,283], Sc/Ti [283] as well as for waveguide laser applications, the important Er/Ti [283,284].

Proton-Exchange
A simple alternative to metal-diffusion based waveguide production is the protonexchange (PE) technique reported in 1981 by Jackel and Rice [285], trailed by the production of the first optical waveguide in the following year [286]. Proton exchange is achieved by placing the LN crystal in a liquid source (most commonly melted benzoic acid) that acts as hydrogen source at elevated temperature. For a review see Cabrera et al. [287]. Unfortunately, the as-exchanged layer suffers from lattice stress and damage, as found in early RBS/C studies [288,289], which degrades the optical properties. Therefore, a subsequent annealing procedure at about 300 • C in air is employed in order to obtain a high-quality waveguide. Figure 27 shows the effect for LN proton-exchanged using benzoic acid. Studies of the as-exchanged layer thickness are performed by analyzing the damage profile in axially aligned directions by RBS. Changing process temperatures and times as well as the cut of the substrate allowed to determine their influence on the exchange-depth and the hydrogen diffusion coefficient. Results were reported for benzoic acid in the low (110-220 °C) [290] and high temperature range (200-300 °C) [291,292]. In accordance with the expectations, the exchange depth was found to be proportional to the square root of the exchange time while the diffusion constant follows an Arrhenius law [290] Analogue findings were reported for samples treated with sulfuric [293] and pyrophosphoric acid [294]. The comparison of channeling spectra taken before and after annealing shows a strong but incomplete recovery of the LN crystal structure [291,292,294,295].
Nb-disorder profiles extracted from the RBS/C spectra show differences between xand z-cut LN [292]. The as-exchanged z-cut material shows a low damage (up to 25% of the Nb-atoms displaced) increasing from surface to bulk with an abrupt end. On the other hand, for x-cut material the damage was found to be constant at a high level (80%). After annealing the damage profiles extended for both cases deeper into the bulk and interestingly, showed a smaller degree of recovery for the z-cut substrate. Olivares et al. [291] measured RBS/C spectra and angular scans for Nb-RBS along several axes in asexchanged waveguides demonstrating high-temperature treatment leads to a superior crystal quality. Moreover, virtually no differences between samples prepared under slightly different conditions (addition of lithium benzoate, annealing atmosphere) were observed. The authors also performed ERD angular scans for the lattice location of hydrogen at several planar directions [291,296]-due to the need of large ion beam incidence angles (θ > 70°) for this technique no appropriate axial directions were available. Studies of the as-exchanged layer thickness are performed by analyzing the damage profile in axially aligned directions by RBS. Changing process temperatures and times as well as the cut of the substrate allowed to determine their influence on the exchange-depth and the hydrogen diffusion coefficient. Results were reported for benzoic acid in the low (110-220 • C) [290] and high temperature range (200-300 • C) [291,292]. In accordance with the expectations, the exchange depth was found to be proportional to the square root of the exchange time while the diffusion constant follows an Arrhenius law [290] Analogue findings were reported for samples treated with sulfuric [293] and pyrophosphoric acid [294]. The comparison of channeling spectra taken before and after annealing shows a strong but incomplete recovery of the LN crystal structure [291,292,294,295].
Nb-disorder profiles extracted from the RBS/C spectra show differences between x-and z-cut LN [292]. The as-exchanged z-cut material shows a low damage (up to 25% of the Nb-atoms displaced) increasing from surface to bulk with an abrupt end. On the other hand, for x-cut material the damage was found to be constant at a high level (80%). After annealing the damage profiles extended for both cases deeper into the bulk and interestingly, showed a smaller degree of recovery for the z-cut substrate. Olivares et al. [291] measured RBS/C spectra and angular scans for Nb-RBS along several axes in as-exchanged waveguides demonstrating high-temperature treatment leads to a superior crystal quality. Moreover, virtually no differences between samples prepared under slightly different conditions (addition of lithium benzoate, annealing atmosphere) were observed. The authors also performed ERD angular scans for the lattice location of hydrogen at several planar directions [291,296]-due to the need of large ion beam incidence angles (θ > 70 • ) for this technique no appropriate axial directions were available. A comparison of the experimental results with computer simulations based on a theoretical model proposed by Kóvacs et al. [297] (shift of H in c-direction by 0.1 Å, in y-direction ≈ 1Å) revealed a good agreement. Although somewhat ambiguous, these results support clearly support the formation of O-H suggested by theory.
As in the case of Ti-diffusion the actual elemental profiles within the waveguide are of outmost interest for establishing a correspondence with refractive index profile. Suitable nuclear depth profiling techniques employed in proton-exchanged LN for hydrogen are ERD, NRA and SIMS, for lithium NRA and SIMS and for oxygen SIMS.
Canali et al. [289] measured the H and Li depth profiles of as-exchanged layers using nuclear reactions. In the case of hydrogen profiling was done with the 15 N(p,αγ) 12 C reaction that allows due to its resonant character a very high depth resolution; Li-profiles were taking using the non-resonant 7 Li(p,α) 4 He reaction that allows only moderate depth resolution. The hydrogen profiles overlap with the damage depth observed in the RBS/C. The Li-NRA spectra reveal a 70% reduction of lithium in the treated layer corroborating that indeed an exchange reaction between Li and H takes place. Additional SIMS studies by the authors confirmed these results. Hsu et al. [295] verified the depletion in as-exchanged layers and demonstrated that after annealing a partial recovery of the Li concentrations occurs although a thin surface layer stays poor in lithium. Resonant NRA depth profiling by Ito and Kawamoto [298] revealed a reduction of the H-concentration in the surface after annealing and further diffusion into the bulk out of the depth region accessible by this method.
Rottschalk et al. [299] employed ERD for the H depth profiling of PE waveguides (benzoic acid, 180 • C, 5 h or 250 • C, 1 min) in as exchanged state and after annealing at various temperatures. Despite the extremely short treatment time, it was clearly visible that higher temperature yielded a higher hydrogen surface concentration. As expected, the hydrogen diffuses more rapidly with increasing annealing temperature reducing its concentration at the surface. The work of Hagner and Bachmann [300] investigated the H-concentration in the as-exchanged waveguide in dependence of the amount of lithium benzoate admixture to the benzoic acid used in the process. The ERD spectra showed a reduction of the hydrogen incorporation into the layer with otherwise identical conditions. Hydrogen profiles in PE-LN layers (benzoic acid, 185 • C, 70 min) annealed at 400 • C for 6-180 min were obtained by Casey et al. [301] and Zavada et al. [302] using SIMS. The measured H concentration profiles had to be fit using two Gaussian diffusion profiles which suggest the existence of two different types of hydrogen diffusion with distinct diffusion coefficients. The authors suggest that the fast component is related to interstitial diffusion while the slower represents substitution of lithium by hydrogen. De Souza et al. [303,304] confirmed and quantified the fractions of these two different species of hydrogen with NMR measurements for several exchange conditions. In a second study [305] the authors succeeded to separate the signal origination from the "static" hydrogen and compared the results with calculations for four defect models. The best agreement was found for H occupying a position within an oxygen plane along a side of an oxygen triangle equidistant from two oxygen atoms.
Yamamoto and Taniuchi [306] as well as Cheng et al. [307] reported SIMS-based hydrogen profiles for the alternative pyrophosphoric acid exchange method and show that the profiles obtained are very similar to those of samples using the common benzoic acid.

Ion Implantation
Among the waveguide fabrication techniques ion implantation is the one that can be applied to largest variety of materials. For comprehensive reviews see Chen et al. [308] and Peña-Rodriguez et al. [309]. As discussed in Section 3.3 defects are produced by electronic and nuclear interactions. In many cases, these defects cause an alteration of the refractive index forming an optical waveguide. The first successful waveguide formation using ion implantation was reported by Destefanis et al. [310]. For waveguides based on nuclear damage induced defects the quality of the waveguides is improved by an annealing procedure to remove electronic damage in the light guiding layer but leaving the buried damage layer intact. In the meantime, waveguide formation based on electronic damage production gained importance with many ions reported to be useful: B, C, N, O, F, Si, P, Ti, Cu, Ni, and Ag.
Barfoot et al. [311] studied the damage built-up during implantation with 1 MeV He + (2.0 × 10 15 and 2.4 × 10 16 cm −2 ) and annealing induced recovery by in-situ RBS/C. Varying the implantation temperature between 40 and 140 • C resulted in a decrease of the Nb-lattice disorder with increasing temperature. In all cases a saturation of the damage was observed at about 1 × 10 16 He + cm −2 coinciding with the value found for optimal waveguide performance. Further, the best annealing temperature was determined to be about 180 • C. Williams et al. [312] investigated similar waveguides (implanted with 1.7 MeV He + , 1.5 × 10 15 and 2.0 × 10 16 cm −2 , RT; annealed at 175 • C and 400 • C) by NRA using the 7 Li(p,α) 4 He reaction. A comparison indicated among the NRA/C spectra indicated a lithium movement towards the nuclear damage region during annealing. The authors suggest that this effect is responsible for the increase in the extraordinary index in the surface region.
Low-energy He + -implantation (50 keV, 2.0 × 10 15 and 2.4 × 10 16 cm −2 ) performed at temperatures of −80 • C, −30 • C, 25 • C and 80 • C has been studied by Al-Chalabi et al. [313] employing RBS/C and HIERD. The damage in the sample implanted at −80 • C appeared strongly enhanced while no difference between those implanted at −30 • C and 80 • C was detected. Li-and H-profiles obtained from the HIERD data indicate a loss of lithium in the implanted layer that is compensated by the incorporation of hydrogen. The lithium concentration reduction was found to increase with increasing fluence and implantation temperature.
Kling et al. [314] studied the annealing temperature dependence of the lattice damage in y-cut LN using RBS/C. Samples were implanted at room temperature with 20 keV He + to a fluence of 1 × 10 16 cm −2 or with 25 keV Li + (1 × 10 16 cm −2 )-that yield the same damage profile-in order to check for the chemical influence of the implant. Annealing was performed under wet oxygen atmosphere for 10 min at five different temperatures (150-350 • C). At first, an increase of the damage during the heat treatment was observed with increasing temperature reaching a maximum at 270 • C with an apparent activation energy of 0.14 eV, as depicted in Figure 28. For Li-implanted samples an almost complete recovery of the lattice has been observed at this temperature indicating that He-implantation induces fundamentally different damage structures. A further increase of temperature induced a decrease in the number of defects with a much larger activation energy of 1.3 eV. This two-phase behavior resembles closely that observed for He-implantation in silicon [315]. An additional TEM study on the same samples annealed at 270 • C revealed the existence of platelet voids in the depth region at which the He ions were stopped. This region forms naturally a layer of low refractive index.
The annealing behavior of LN implanted with 800 keV He + ions at very low temperature (−173 • C) to fluences of 5 × 10 15 and 1 × 10 16 cm −2 was the topic of a study of Gischkat et al. [316]. RBS/C measurements were performed at implantation temperature, at room temperature and after heat treatments at 250 • C and 300 • C in air. An apparently higher defect concentration observed in x-cut samples with respect to z-cut material was observed. As discussed in the Section 3.3 this anisotropy is due to the dislocation of Nb to the free octahedron site during implantation leading to different defect geometries in terms of channeling for the two directions. The dechanneling effect due to defects was observed to decrease with heat treatment for the low-fluence case and to increase for the high-fluence case. In the latter case, the enhancement was attributed to the formation of dislocation loops. recovery of the lattice has been observed at this temperature indicating that Heimplantation induces fundamentally different damage structures. A further increase of temperature induced a decrease in the number of defects with a much larger activation energy of 1.3 eV. This two-phase behavior resembles closely that observed for Heimplantation in silicon [315]. An additional TEM study on the same samples annealed at 270 °C revealed the existence of platelet voids in the depth region at which the He ions were stopped. This region forms naturally a layer of low refractive index.  Besides helium, the implantation of nitrogen has already been explored in 1980s with view to its applications in optical waveguides [99,100]. Damage profiles for samples implanted with 200 keV N + with fluences between 5 × 10 14 and 2 × 10 16 cm −2 were determined by RBS/C and correlated with the refractive index of the material. The formation of optically isotropic waveguide by amorphization due to electronic damage formation of a surface layer by implantation of high-energy (5 MeV and 7.5 MeV) Si + ions has been studied by Olivares et al. [117,317]. From RBS/C spectra for samples implanted with fluences ranging from 5 × 10 13 to 1 × 10 15 cm −2 the authors determined threshold values for the amorphization and the associated electronic stopping power. Similar studies have also been performed for H [318], O [111,[319][320][321] and F [116].
Optical waveguides can also be formed by high-fluence implantation of Ti with subsequent annealing/diffusion treatment. Bremer et al. [322] implanted y-cut LN with 200 keV Ti to fluences of 2.5 × 10 17 cm −2 and 4 × 10 17 cm −2 at liquid nitrogen temperature, room temperature and 350 • C. The samples underwent subsequent heat treatment at 1000 • C under wet oxygen atmosphere for 2 to 8 h. RBS/C studies demonstrated a good crystal lattice recovery was achieved in all cases for long annealing times but that a higher implantation temperature favors the recrystallization process. Optical waveguide properties were observed for all samples.
A technique related to implantation is ion beam mixing. A LN sample plated with a 40 nm metallic Ti layer has been irradiated with 3 MeV Ti + ions at a fluence of 1 × 10 16 cm −2 by Bremer et al. [266]. After removal of excess Ti from the surface the sample was annealed at 1000 • C under wet oxygen atmosphere to achieve recrystallization. The SIMS profile measurement showed a Ti profile with a steady slope reaching deep into the bulk material which differs strongly from the Gaussian profiles found in indiffused waveguides.

Combination of Techniques
The combination of the proton-exchange and the implantation technique produces interesting results. For samples subjected first to a He + implantation (200 keV, 3 × 10 16 cm −2 ) and then undergoing a proton exchange treatment the hydrogen indiffusion is strongly enhanced [323]. It appears that the damage introduced by ion implantation process facilitates the incorporation of the hydrogen. For the inverse treatment sequence, it was observed that a proton-exchanged sample subsequently implanted with 6 MeV O + at a fluence of 6 × 10 14 cm −2 is fully amorphized while PE-only and implantation-only samples show only minor damage [324]. Obviously, the proton-exchange process leads to a more fragile lattice increasing the susceptibility to lattice disorder.

Interaction with Optical Active Dopants
For optical waveguide lasers the absence of interferences between the waveguide formation process and the optically active dopants is mandatory. Therefore lattice site location and crystal structure studies for several rare earths in bulk and processed waveguide materials have been performed. Herreros et al. [325,326] compared the lattice sites for LiNbO 3 :Tm,Mg and LiNbO 3 :Er,Mg in as-grown and ion-implanted (2.6 and 2.0 MeV He + , respectively, 2 × 10 16 cm −2 ; annealing: 250 • C, 30 min, air) samples. Angular scans for an axial direction demonstrated that the lattice sites of Tm and Er stay unaltered after waveguide fabrication. The same authors also studied the structural differences in waveguides fabricated by proton-exchange (benzoic acid, 300 • C, 20 h) and by He + implantation (2 MeV, 2 × 10 16 cm −2 ) in LiNbO 3 :Nd,Mg using RBS/C. Again, the results show that ion-implanted waveguides essentially maintain their properties while proton-exchanged present a larger deterioration [327]. Cajzl et al. [328] reported on the influence of the three waveguide fabrication methods (Ti-diffusion, proton-exchange and ion implantation) on the lattice site of erbium. Angular RBS scans confirm that none of the methods has any measurable influence on the position of the Er. These results indicate that rare-earth ions are successfully incorporated of into LN waveguides.

Conclusions
The relevance of nuclear techniques in the investigation of defect structures in modern materials has been demonstrated for the case of LN, presently one of the most important optoelectronic materials. It has been shown that the techniques originally developed for a simple crystal system can be successfully transferred to complex cases.
In the case of the intrinsic defect structure of congruent LN the combination of Ion Beam Analysis and Perturbed Angular Correlation studies delivered evidence for an additional lattice site in the Li-Octahedron. The features of this site observed by nuclear techniques support the interpretation that the excess Nb is incorporated through ilmenitetype stacking faults dispersed in the crystal lattice.
Channeling measurements are found to be the method of choice for the characterization of ion implantation induced damage and its recovery laying the basis for the theoretical interpretation of damage formation in oxides.
The determination of lattice sites for a large number of elements, including the role of co-doping furnishes pertinent information for future engineering of LN by introduction of foreign atoms. It was demonstrated that the main factor determining the occupation of a Lior a Nb-octahedron by a dopant is the valence state. Divalent, trivalent and tetravalent ions occupy-with the single exception of Cr, the Li-Octahedron, pentavalent and hexavalent the Nb.
Depth profiling and defect characterization in optical waveguides by ion beam methods is essential for the development of opto-electronic devices.