Boron Inﬂuence on Defect Structure and Properties of Lithium Niobate Crystals

: Defect structure of nominally pure lithium niobate crystals grown from a boron doped charge have been studied by Raman and optical spectroscopy, laser conoscopy, and photoinduced light scattering. An inﬂuence of boron dopant on optical uniformity, photoelectrical ﬁelds values, and band gap have been also studied by these methods in LiNbO 3 crystals. Despite a high concentration of boron in the charge (up to 2 mol%), content in the crystal does not exceed 10 − 4 wt%. We have calculated that boron incorporates only into tetrahedral voids of crystal structure as a part of groups [BO 3 ] 3 − , which changes O–O bonds lengths in O 6 octahedra. At this oxygen–metal clusters MeO 6 (Me: Li, Nb) change their polarizability. The clusters determine optically nonlinear and ferroelectric properties of a crystal. Chemical interactions in the system Li 2 O–Nb 2 O 5 –B 2 O 3 have been considered. Boron, being an active element, structures lithium niobate melt, which signiﬁcantly inﬂuences defect structure and physical properties of a crystal grown from such a melt. At the same time, amount of defects Nb Li and concentration of OH groups in LiNbO 3 :B is close to that in stoichiometric crystals; photorefractive effect, optical, and compositional uniformity on the contrary is higher.


Introduction
Lithium niobate (LN, LiNbO 3 ) attracts attention due to its possible applications in integral and nonlinear optics, pure optics (generation of optical harmonics, lasing parametric generation, electro-optics, optical amplification, and conversion of optical radiation), acoustoelectronics (bandpass filters and SAW delay lines), quantum electronics, and solid state physics [1][2][3][4]. The equipment associated with modern optoelectronic and telecommunication technologies often includes LN crystals. Such applications claim LN crystals with high optical uniformity and optical damage resistance. Thus, a study of their defect structure and optical characteristics in dependence on obtaining conditions are highly relevant.
A LiNbO 3 crystal is a non-stoichiometric oxygen octahedral ferroelectric with high Curie temperature (1420 K), spontaneous polarization (P S = 5 × 10 −5 C/cm 2 ) and a wide homogeneity region on the phase diagram (44.5-50.5 mol % Li 2 O at 1460 K). LN should be considered as a solid solution LiNbO 3 :Nb [5,6]. Doping by a wide spectrum of metal elements is possible due to an octahedral coordination of metal ions in the LN structure. At this a significant however preserving symmetry distortion of MeO 6 (Me: Li + , Nb 5+ , dopant) octahedra can occur [1,2,7]. All these factors provide a possibility to control physical characteristics of the material. Being a phase of a variable composition, LN crystal has a highly developed defect structure. Optical damage (photorefractive effect) is determined by defects with localized electrons; such defects form photoelectric fields. LN is in general inapplicable for optical elements manufacturing. An increase in optical uniformity and optical damage resistance of SLN and NSLN crystals is achieved by two methods. The first one is HTTSSG (High Temperature Top Seeded Solution Growth) with addition~6 wt% of K 2 O flux (LN:K 2 O crystals) [29]. Recent papers [21][22][23][24][25] report the other know way: NSLN crystals with a high optical quality can be obtained from a congruent charge with a B 2 O 3 flux (LiNbO 3 :B crystals). Doping LN with boron allows us to combine approaching stoichiometric composition and Nb Li defects concentration decrease. As long as it is important to know technological details of obtaining of different types of highly useful LN crystals, it also is important to learn how different concentration of B exactly influence LN melt and thus structure and other properties.
In this work we bring together data from optical and atomic force microscopy, optical spectroscopy, Raman spectroscopy, laser conoscopy, photoinduced light scattering (PILS), IR-spectroscopy in the region of stretching vibrations of OH groups and computer simulation of the defect structure. Methods have been applied to nominally pure NSLN crystals grown by different technologies from a charge doped by boron (0.55-1.24 mol% B 2 O 3 in a charge). As long as grown crystals have approximately the same traces amount of boron in their structure, we consider that the designation in mol% means the B 2 O 3 content in the charge. We also present a review of chemical interactions taking place in systems Li 2 O-Nb 2 O 5 and Li 2 O-Nb 2 O 5 -B 2 O 3 . Such interactions determine fine features of the melt, crystallization, and crystals. Results for LiNbO 3 :B (0.55-1.24 mol%) crystals were compared with those of SLN and CLN crystals obtained due to a traditional technologies.
IR adsorption and Raman spectra are highly sensitive to changes in a crystal structure, thus, defectivity. Sites occupied by hydrogen in LN structure are the most sensitive ones to crystal field changes. The fundamental absorption edge position and features are sensitive towards the structural uniformity of the crystal [14,18,23,24,29,30]. These two methods also allow one to evaluate Li/Nb ratio, concentration of point defects Nb Li and V Li , type of complex defects that include OH-groups. PILS speckle-structure can be used to calculate photoelectric fields (photovoltaic E pv and diffuse E D ) using approach offered in [31].

Specific Features of Obtaining Nominally Pure LN Crystals from a Boron Doped Charge and the Method of Their Experimental Studies
Studied nominally pure CLN and LiNbO 3 :B (0.55-1.24 mol%) crystals were grown from a congruent melt [16]. Doping was carried out both by homogeneous doping of Nb 2 O 5 precursor [32] and direct solid-phase doping of LN charge [20,33]. Direct doping is a solid-phase synthesis of niobium pentoxide, lithium carbonate and boric acid with the subsequent production of a granular charge during annealing at pre-melting temperatures at 1240-1250 • C. This method was used to grow LiNbO 3 :B (0.55-0.83 mol%) crystals. Nevertheless, the method itself has several disadvantages: LN crystals have articulated secondary structure and a number of macro-defects [33,34] (Appendix A, Figures A1-A12).
Macro-and microdomain structure of doped LN crystals was studied by an image processing system «Thixomet». The system includes optical microscope Axio Observer D1m (Carl Zeiss, Oberkochen, Germany) connected through digital video camera Pix-eLink PL-B774U (Pixelink, Ottawa, Canada) with a computer equipped with the program ThixometPRO (Thixomet, Saint Petersburg, Russia). Studies were carried out in bright field and differential interference contrast. Fine crystal plates of Z-and X-orientation were pre-ground, polished and selectively etched at room temperature for 18 h in the HF:HNO 3 = 1:3 mixture. Nanostructures of doped crystals were studied by atomic force microscopes CMM-2000 (Zavod "Proton" MIET, Moscow, Russia) and Nano-R (Pacific Nano Technology, Santa Clara, CA, USA).
Earlier we have shown [19,20] that LN crystals with a high optical quality and a low photorefractive effect can be grown from a charge containing no more than~0.1 to 0.12 wt%. For example, LiNbO 3 :B crystals grown from melts with~0.12 to 0.25 wt% B 2 O 3 contained microstructure defects absent in other LN crystals types. The melt viscosity greatly increases at high boron content: a high complexing ability of boron leads to formation of an anti-seeding viscous film on the melt surface. This results in growth of LiNbO 3 :B crystals with numerous macro-and micro-defects irremovable by aftergrowth annealing. Thus, growing LN from boron-containing melts required the development of new approaches and solutions to adapt the commonly used technology. Changes affected technological parameters of growth, intrinsic equipment of a growth camera and synthesis of initial charge methods [19,20,25]. Papers [19,20] have revealed not only changes in the growth parameters, but also specific macrodefects, which were found only in LiNbO 3 :B crystals: "channels" and deviation of optical density (optical inhomogeneity of crystals). "Channels" are almost cylindrical long curved holes inside a LiNbO 3 :B boule with a diameter~1 to 200 µm. We have calculated the average density of such defects per 1 mm 2 of the investigated area,~0 to 7 mm −2 , Figure 1. The crystals under consideration are a LiNbO 3 :B series with an increase content of boron in the melt~0.12 to 0.25 wt%. The amount of "channels" decreases up to the total absence with a decrease in the boron concentration in the melt to~0.1 wt%, as paper [20] claims. Deviation of optical density due to deviations in the melt viscosity was discovered when examining polished LiNbO 3 :B plates by optical microscopy in the differential interference contrast (DIC) mode. DIC mode is actively used to study unpainted transparent objects. The resulting image is conventionally colored since it is the result of the polarized light beam interference on structural objects with different optical density. They look like colored spots of various sizes and are usually observed together with "channels", Figure 2. The latter is especially relevant for strongly doped melts. A colored with various pseudocolors structure without sharp boundaries is typical for all LiNbO 3 :B crystals that underwent only post-growth thermal annealing. Crystals turned to a single-domain state loose this type of a macrostructure, samples have uniform and monochrome image in a DIC mode, Figure 2e,f [19,20,25].
Crystals 2020, 10, x FOR PEER REVIEW 4 of 37 wt%. For example, LiNbO3:В crystals grown from melts with ~0.12 to 0.25 wt% B2O3 contained microstructure defects absent in other LN crystals types. The melt viscosity greatly increases at high boron content: a high complexing ability of boron leads to formation of an anti-seeding viscous film on the melt surface. This results in growth of LiNbO3:В crystals with numerous macro-and micro-defects irremovable by after-growth annealing. Thus, growing LN from boron-containing melts required the development of new approaches and solutions to adapt the commonly used technology. Changes affected technological parameters of growth, intrinsic equipment of a growth camera and synthesis of initial charge methods [19,20,25]. Papers [19,20] have revealed not only changes in the growth parameters, but also specific macrodefects, which were found only in LiNbO3:В crystals: "channels" and deviation of optical density (optical inhomogeneity of crystals). "Channels" are almost cylindrical long curved holes inside a LiNbO3:В boule with a diameter ~1 to 200 μm. We have calculated the average density of such defects per 1 mm 2 of the investigated area, ~0 to 7 mm -2 , Figure 1. The crystals under consideration are a LiNbO3:B series with an increase content of boron in the melt ~0.12 to 0.25 wt%. The amount of "channels" decreases up to the total absence with a decrease in the boron concentration in the melt to ~0.1 wt%, as paper [20] claims. Deviation of optical density due to deviations in the melt viscosity was discovered when examining polished LiNbO3:B plates by optical microscopy in the differential interference contrast (DIC) mode. DIC mode is actively used to study unpainted transparent objects. The resulting image is conventionally colored since it is the result of the polarized light beam interference on structural objects with different optical density. They look like colored spots of various sizes and are usually observed together with "channels", Figure 2. The latter is especially relevant for strongly doped melts. A colored with various pseudocolors structure without sharp boundaries is typical for all LiNbO3:B crystals that underwent only post-growth thermal annealing. Crystals turned to a single-domain state loose this type of a macrostructure, samples have uniform and monochrome image in a DIC mode, Figure 2e,f [19,20,25]. type in the LiNbO3:B (~6.20⋅10 -5 wt% in a crystal cone) crystal; (c), (d)-individual macrodefects of the "channel" type in the LiNbO3:B (~6.0⋅10 -5 wt% in a crystal cone) crystal; (e)-macrodefects of the "channel" type in the LiNbO3:B (0.1 wt% in the charge) crystal. We have noted that using growth parameters suitable for obtaining of pure LN, LN:Zn or Mg at obtaining of LiNbO3:B leads to cellular growth in the boule bottom, Figure  3a. Such conditions also make the crystallization front more sensitive towards smallest thermal conditions violations: the hardly notable change in water pressure in the cooling circuit of the growth unit leads to a change in the LiNbO3:B boule diameter, Figure 3b. These macro-and micro-defects decrease the useful volume of the crystal, suitable for manufacturing of optical devices units. We have noted that using growth parameters suitable for obtaining of pure LN, LN:Zn or Mg at obtaining of LiNbO 3 :B leads to cellular growth in the boule bottom, Figure 3a. Such conditions also make the crystallization front more sensitive towards smallest thermal conditions violations: the hardly notable change in water pressure in the cooling circuit of the growth unit leads to a change in the LiNbO 3 :B boule diameter, Figure 3b. These macroand micro-defects decrease the useful volume of the crystal, suitable for manufacturing of optical devices units.
Homogeneous doping essence is obtaining of Nb 2 O 5 :B precursor. Boric acid is added to a niobium strip product solution. The latter is obtained during extraction purification of niobium hydroxide [27,33]. Doped niobium pentoxide Nb 2 O 5 :B is used as a precursor during LN charge synthesis. A crystal LiNbO 3 :B (1.24 mol%) was grown from such a charge. Homogeneous doping allows us to grow crystals without growth stripes and other macro-and microdefects as those characteristic of direct doping [27]. Due to mass spectrometry, concentration of boron in crystals grown by both methods is close to traces (~10 −3 to 10 −5 wt%). Impurities concentrations were: Pb, Ni, Cr, Co, V, Ti, Fe, and Al less than 2 × 10 −4 , Ca, Si less than 1 × 10 −3 , F less than 1 × 10 −3 wt%. Homogeneous doping essence is obtaining of Nb2O5:B precursor. Boric acid is added to a niobium strip product solution. The latter is obtained during extraction purification of niobium hydroxide [27,33]. Doped niobium pentoxide Nb2O5:В is used as a precursor during LN charge synthesis. A crystal LiNbO3:В (1.24 mol%) was grown from such a charge. Homogeneous doping allows us to grow crystals without growth stripes and other macro-and microdefects as those characteristic of direct doping [27]. Due to mass spectrometry, concentration of boron in crystals grown by both methods is close to traces (~10 -3 to 10 -5 wt%). Impurities concentrations were: Pb, Ni, Cr, Co, V, Ti, Fe, and Al less than 2·10 -4 , Ca, Si less than 1·10 -3 , F less than 1·10 -3 wt%.

The Growing and Characterization Methods of the LN Crystals
An SLN crystal was grown from a melt with 58.6 mol% Li2O. Growing of stoichiometric LN crystals is a complex technical task due to a crystal's wide homogeneity area on a phase diagram [4,7,35]. A great difference between the melt and crystal compositions require a decrease in growth velocity V several orders of magnitude. This suppresses concentration over-cooling that can lead to changes in the crystal composition at different stages of the growth. Thus, even at a very small growth velocity (much less than V = 0.1 mm/h), one can obtain a crystal 10 mm in diameter and 4 mm long from 200 g of melt. Its optical uniformity Δn (where n is a refractive index) is compared with that of CLN crystals of the same diameter grown at a velocity V = 3-5 mm/h [36]. Traditional technology provides very little possibility to obtain SLN single crystals of a big enough diameter to create optical elements suitable in optical instrumentation, quantum electronics, etc.
Nominally pure CLN was grown using original ICT KSC RAS charge that allows us to obtain water white crystals [37]. All crystals were grown by Czochralski on Crystal-2 (Zavod Kristall Ltd, Saint Petersburg, Russian Federation) with automatic control of a crystal diameter.
Crystals were turned to a single-domain state by a high temperature electro-diffuse annealing: constant current was applied to crystals during their cooling in a temperature interval 1240-880 °С with a speed 20 grad/h. The grade of single domain state was determined by analysis of the frequency dependence of electrical impedance and by determining the value of the static piezomodulus (d333st).
Samples for optical experiments were cut from grown boules in the shape of parallelepipeds so that edges coincided with crystallographic axes X, Y, and Z (the polar axis) and in the shape of Z-oriented plane-parallel plates ~1 mm thick. Parallelepiped samples had sizes ~7 × 6 × 5 mm 3 . Faces of all studied samples were thoroughly polished.
Raman spectra were excited with a 514.5 nm line of laser 2018-RM (SpectraPhysics, Milpitas, CA, USA) and detected by a spectrograph T64000 (Horiba Jobin Yvon, Palaiseau,

The Growing and Characterization Methods of the LN Crystals
An SLN crystal was grown from a melt with 58.6 mol% Li 2 O. Growing of stoichiometric LN crystals is a complex technical task due to a crystal's wide homogeneity area on a phase diagram [4,7,35]. A great difference between the melt and crystal compositions require a decrease in growth velocity V several orders of magnitude. This suppresses concentration over-cooling that can lead to changes in the crystal composition at different stages of the growth. Thus, even at a very small growth velocity (much less than V = 0.1 mm/h), one can obtain a crystal 10 mm in diameter and 4 mm long from 200 g of melt. Its optical uniformity ∆n (where n is a refractive index) is compared with that of CLN crystals of the same diameter grown at a velocity V = 3-5 mm/h [36]. Traditional technology provides very little possibility to obtain SLN single crystals of a big enough diameter to create optical elements suitable in optical instrumentation, quantum electronics, etc.
Nominally pure CLN was grown using original ICT KSC RAS charge that allows us to obtain water white crystals [37]. All crystals were grown by Czochralski on Crystal-2 (Zavod Kristall Ltd, Saint Petersburg, Russian Federation) with automatic control of a crystal diameter.
Crystals were turned to a single-domain state by a high temperature electro-diffuse annealing: constant current was applied to crystals during their cooling in a temperature interval 1240-880 • C with a speed 20 grad/h. The grade of single domain state was determined by analysis of the frequency dependence of electrical impedance and by determining the value of the static piezomodulus (d 333st ).
Samples for optical experiments were cut from grown boules in the shape of parallelepipeds so that edges coincided with crystallographic axes X, Y, and Z (the polar axis) and in the shape of Z-oriented plane-parallel plates~1 mm thick. Parallelepiped samples had sizes~7 × 6 × 5 mm 3 . Faces of all studied samples were thoroughly polished.
Raman spectra were excited with a 514.5 nm line of laser 2018-RM (SpectraPhysics, Milpitas, CA, USA) and detected by a spectrograph T64000 (Horiba Jobin Yvon, Palaiseau, France) equipped with a confocal microscope. Low excitation power (not higher than 3 mW under the microscope) decreased an influence of photorefractive effect on Raman spectrum. Spectra resolution was 1 cm −1 . Spectra were treated using Horiba LabSpec 5.0 and Origin 8.1 programs. Instrumental error in determination of frequency (ν), width (S) and intensity (I) of Raman bands was ±1.0, ±3.0 cm −1 and 5%, respectively.
PILS and laser conoscopy methods are described in detail in papers [11,13,38]. Conoscopic patterns and PILS patterns were excited by laser Nd:YAG (MLL-100, Changchun New Industries Optoelectronics Tech. Co. Ltd, Changchun, China), wavelength λ 0 = 532 nm Crystals 2021, 11, 458 7 of 37 and power p = 1 and 90 mW. The studied sample was located on a movable two-coordinate plate, which allowed us to obtain several conoscopic patterns from different areas of a crystal. Conoscopic pattern formed on a translucent screen and was registered by a digital photo camera.
Photoelectric fields studies were carried out using PILS patterns excited by argon laser 2018-RM (SpectraPhysics, Mountain View, CA, USA) lines 476.5 nm (p = 216 mW); 488.0 nm (p = 98 mW); 514.5 nm (p = 282 mW); 530.9 nm (p = 160 mW). Laser beam diameter was 1.8 mm. Scattered radiation fell on a translucent screen and was registered by a digital photo camera. PILS in both cases was registered in ee-type geometry: laser radiation that fells on the crystal and scattered radiation have the same polarization, at the same time E vector is oriented parallel to the crystals polar axis Z.
PILS indicatrix parameters of studied crystals provided values of strengths of photovoltaic and diffusion electric fields. The error of our experiment calculations was 1.5-2%. The calculation of photovoltaic and diffusion field strengths was carried out in Mathcad 15.0 using approach from [31]. Refractive indexes of ordinary and extraordinary rays were determined from empirical equations from [39]. PILS indicatrix angle was calculated due to a formula θ = arctg(a/b), where a is the indicatrix size in the positive direction of the polar axis, b is the distance between the crystal and the screen.
Optical absorption spectra of LN crystals were determined using spectrophotometer Varian Cary 2300 (Varian Inc., Palo Alto, CA, USA), SF-256 UVI (Granat, Saint Petersburg, Russian Federation) and monochromator MDR-41 (OKB Spectr, Saint Petersburg, Russian Federation). A deuterium lamp was used as a radiation source. Band gap was determined using crystal transmission spectrum. An invert spectrum was obtained from the latter-an absorption spectrum. The obtained absorption spectrum in the decreasing linear part of the graph was approximated by a straight line until it crossed the abscissa axis. The point of intersection of this straight line and the abscissa axis is the wavelength corresponding to the absorption edge of the crystal. A bandgap was determined due to a formula E = hc/λ, where λ is a wavelength of absorption edge, h-Planck's constant, c-light speed in vacuum. The band gap determination error was ±1.0 nm.
Differential thermal analysis (DTA) was carried out on an installation with a thermal block made of ruby single crystal, which provides a gradient-free zone, and high sensitivity and resolution [40]. Error in thermal effects temperature determination was ±0.5 • C. The detailed description of the method is available from [40].
Computer simulation of B 3+ ion localization and electrostatic interaction of point charges in LiNbO 3 :B crystals was carried out by a calculation of sum energy of Coulomb interaction of point charges (U, eV) in an oxygen-octahedral cluster of a LN structure (Li + , Nb 5+ , O 2− ) with B 3+ element.
The chemical analysis was carried out by atomic emission spectrometry (spectrometer ICPS-9000 by Shimadzu), with an error~1%.

Study of LiNbO 3 :B Crystals Defective Structure by Raman Spectroscopy
Raman spectra of LiNbO 3 :B (0.55-1.24 mol%) crystals were previously studied in papers [21,22,41,42]. The papers focus on an influence of a disorder in cation sublattice on A 1 (TO) symmetry type phonons in Y(ZZ)Ȳ scattering geometry and on photorefractive effect manifestation in Raman spectra. Figure 4 demonstrates Raman spectra of SLN, CLN and LiNbO 3 :B (0.55-1.24 mol%) in Y(ZX)Ȳ and Y(ZZ)Ȳ scattering geometries. These geometries correspond to fundamental vibration of the lattice of E(TO) and A 1 (TO) symmetry types. Spectra fragments with the most articulated changes at a crystal composition variation are given on Figure 5. Table 1 demonstrates main parameters (ν-frequency, S-band width, I-band intensity) of bands appearing in studied crystals. Figure 6 demon-strates changes in half-widths and intensities of 576 and 630 cm −1 Raman bands and in PILS speckle-structure indicatrix opening angle θ. SLN Raman spectrum contains all nine bands corresponding to E(TO) symmetry type phonons permitted by selection rules in Y(ZX)Ȳ scattering geometry, Figure 4, Table 1. Low-intense bands 179 and 611 cm −1 are clearly observed in SLN Raman spectrum. However, they dissolve by disordering effects in spectra of CLN and LiNbO 3 :B (0.55-1.24 mol%). Frequencies of all bands stay the same within error, Figures 4 and 5, Table 1. This indicates that secondary structure changes (they appear at changes in Li/Nb ratio and B 3+ cation concentration) affect very little quasi-elastic lattice constants in studied LN crystals. At the same time, Figures 4-6 and Table 1 reveal that width and intensities of LiNbO 3 :B Raman bands change significantly in all Raman spectra areas: two-particle states of acoustic phonons (100-150 cm −1 ); vibration of cations located in oxygen octahedra BO 6 (B: Nb, Li, dopant) (200-300 cm −1 ); vibrations of oxygen octahedra atoms (500-900 cm −1 ). Table 1  А1(ТО) symmetry type phonons in Y(ZZ)Ȳ scattering geometry and on photorefractive effect manifestation in Raman spectra. Figure 4 demonstrates Raman spectra of SLN, CLN and LiNbO3:B (0.55-1.24 mol%) in Y(ZX)Ȳ and Y(ZZ)Ȳ scattering geometries. These geometries correspond to fundamental vibration of the lattice of Е(ТО) and А1(ТО) symmetry types. Spectra fragments with the most articulated changes at a crystal composition variation are given on Figure 5. Table 1 demonstrates main parameters (ν-frequency, Sband width, I-band intensity) of bands appearing in studied crystals. Figure 6 demonstrates changes in half-widths and intensities of 576 and 630 cm -1 Raman bands and in PILS speckle-structure indicatrix opening angle θ. SLN Raman spectrum contains all nine bands corresponding to Е(ТО) symmetry type phonons permitted by selection rules in Y(ZX)Ȳ scattering geometry, Figure 4, Table 1. Low-intense bands 179 and 611 cm -1 are clearly observed in SLN Raman spectrum. However, they dissolve by disordering effects in spectra of CLN and LiNbO3:B (0.55-1.24 mol%). Frequencies of all bands stay the same within error, Figures 4 and 5, Table 1. This indicates that secondary structure changes (they appear at changes in Li/Nb ratio and В 3+ cation concentration) affect very little quasielastic lattice constants in studied LN crystals. At the same time, Figures 4-6 and Table 1 reveal that width and intensities of LiNbO3:B Raman bands change significantly in all Raman spectra areas: two-particle states of acoustic phonons (100-150 cm -1 ); vibration of cations located in oxygen octahedra ВО6 (В: Nb, Li, dopant) (200-300 cm -1 ); vibrations of oxygen octahedra atoms (500-900 cm -1 ). Table 1 shows that LiNbO3:В Raman bands in the area 150-300 cm -1 are narrower than those of CLN; 152 and 240 cm -1 widths of LiNbO3:В Raman bands coincide with those of SLN. The 578 cm -1 band corresponds to doubly degenerate Е(TO) vibrations of oxygen atoms in О6 octahedra. It is much wider in spectra of LiNbO3:B (0.55-1.24 mol%) than in SLN and CLN, Table 1, Figure 6. Obtained data show that even very a faint change in boron concentration in the charge and melt (0.55-1.24 mol% В2О3) leads to a significant ordering of cation sublattice along the polar axis at both doping methods. At the same time, oxygen octahedra О6 are distorted. New bands in the area of Raman spectra corresponding to oxygen octahedra vibrations would mean a change in their geometry. However, the area contains no new bands.     S  152  7  152  12  152  7  152  9  152  9  152  10  179  -----------240  9  240  11  41  9  241  11  240  10  240  11  268  10  268  14  270  13  271  12  270  13  270  13  324  10  324  13  325  12  324  14  324  14  323  16  371  17  371  23  371  24  370  24  370  26  371  26  434  10  434  14  432  9  432  10  432  11  432  14  578  16  576  15  575  32  576  33  576  33  575  24  611  -----------630 1  20  626  25  626  41  628  42  626  46  626  28  876 1  20  876  30  875  25  877  25  875  26  875  31  1 Note. 630 and 876 cm -1 bands correspond to А1(ТО) and А1(LO) symmetry types phonons, respectively. The phonons are usually inactive in Y(ZX)Ȳ scattering geometry, in our case they appear due to the photorefractive effect. We do not show widths of low-intense 180 and 610 cm -1 bands due to a great error in their determination.   1 Note. 630 and 876 cm −1 bands correspond to A 1 (TO) and A 1 (LO) symmetry types phonons, respectively. The phonons are usually inactive in Y(ZX)Ȳ scattering geometry, in our case they appear due to the photorefractive effect. We do not show widths of low-intense 180 and 610 cm −1 bands due to a great error in their determination.
Due to selection rules [2,43], fundamental vibrations of E(TO) symmetry type should manifest in Y(ZX)Ȳ scattering geometry, A 1 (TO) symmetry type-in Y(ZZ)Ȳ. A 1 (TO) symmetry type vibrations are forbidden in Y(ZX)Ȳ scattering geometry [2]. However, due to photorefractive effect, they manifest in this scattering geometry. The bands corresponding to A 1 (TO) symmetry type vibrations have intensities proportional to photorefractive effect value. Papers [2,44] have shown that the 630 cm −1 band is the most suitable for photorefractive effect evaluation. The band corresponds to A 1 (TO) vibrations of oxygen atoms in O 6 octahedra. Due to selection rules [2,43], fundamental vibrations of E(TO) symmetry type should manifest in Y(ZX)Ȳ scattering geometry, А1(ТО) symmetry type-in Y(ZZ)Ȳ. А1(ТО) symmetry type vibrations are forbidden in Y(ZX)Ȳ scattering geometry [2]. However, due to photorefractive effect, they manifest in this scattering geometry. The bands corresponding to А1(ТО) symmetry type vibrations have intensities proportional to photorefractive effect value. Papers [2,44] have shown that the 630 cm -1 band is the most suitable for photorefractive effect evaluation. The band corresponds to А1(ТО) vibrations of oxygen atoms in О6 octahedra.
The Е(ТО) 580 cm -1 band does not change with changes in photorefractive effect value, at least in crystals with a low photorefractive effect. Thus, photorefractive effect can be evaluated due to the formula Irel = (I630/I580)⋅100% [11,13,44,45]. Structure disorder contributes the most to LN Raman bands widening during dopant concentration increase at a constant temperature [2,43]. Structure disorder and photorefractive effect contribute the most to LN Raman bands intensity changes [2,43]. Figure 6b shows that the least relative intensity of the 630 cm -1 band is observed in CLN and SLN crystals spectra, the greatestin LiNbO3:В (0.55-0.83 mol%) crystals. In LiNbO3:В spectra the 630 cm -1 band's intensity first increases, than decreases with an increase in В 3+ concentration, Figure 6b. In addition, concentration behavior of 630 cm -1 band intensity and widths of 630 cm -1 (А1(ТО)) and 576 cm -1 (Е(ТО)) bands in the Y(ZX)Ȳ Raman spectra of LiNbO3:В well coincide with the behavior of PILS opening angle θ, Figure 6. Thus, oxygen octahedra distortions that occur due to a change in cation sublattice order and photorefractive effect depend on boron concentration in the charge in studied LiNbO3:В crystals.
Obtained data allow us to conclude that boron structures the melt, which actually means that a congruent LN crystal is growth from a boron-doped charge. Such a crystal The E(TO) 580 cm −1 band does not change with changes in photorefractive effect value, at least in crystals with a low photorefractive effect. Thus, photorefractive effect can be evaluated due to the formula I rel = (I 630 /I 580 ) × 100% [11,13,44,45]. Structure disorder contributes the most to LN Raman bands widening during dopant concentration increase at a constant temperature [2,43]. Structure disorder and photorefractive effect contribute the most to LN Raman bands intensity changes [2,43]. Figure 6b shows that the least relative intensity of the 630 cm −1 band is observed in CLN and SLN crystals spectra, the greatest-in LiNbO 3 :B (0.55-0.83 mol%) crystals. In LiNbO 3 :B spectra the 630 cm −1 band's intensity first increases, than decreases with an increase in B 3+ concentration, Figure 6b. In addition, concentration behavior of 630 cm −1 band intensity and widths of 630 cm −1 (A 1 (TO)) and 576 cm −1 (E(TO)) bands in the Y(ZX)Ȳ Raman spectra of LiNbO 3 :B well coincide with the behavior of PILS opening angle θ, Figure 6. Thus, oxygen octahedra distortions that occur due to a change in cation sublattice order and photorefractive effect depend on boron concentration in the charge in studied LiNbO 3 :B crystals.
Obtained data allow us to conclude that boron structures the melt, which actually means that a congruent LN crystal is growth from a boron-doped charge. Such a crystal has noticeable differences in fine structural features and physical characteristics from a nominally pure CLN crystal, grown from a congruent charge. Obtained results can have the following explanation.
We believe that boron structures the melt so that crystals grown from it have a decreased amount of Nb Li . Raman spectra confirm this conclusion. Raman spectra of both nominally pure and doped CLN crystals contain a low-intense band 120 cm −1 (A 1 (TO)), Figure 5. The band corresponds to two-particle states of acoustic phonons with a total wave vector equal to zero [2]. Due to [46], it cannot correspond to a pseudoscalar mode A 2 forbidden in a point group C 3v by the selection rules. This band is split in two, 105 and 118 cm −1 bands, in a spectrum of a CLN crystal, Figure 5. This happens due to a refinement of the selection rules in the wave vector of two-particle states of A 1 (TO) acoustic phonons [2,47,48]. It is a well-known fact [2,47,48], that spectrum of a perfect SLN crystal is free of the 120 cm −1 band, Figure 5. In addition, SLN crystals are characterized by the perfect order in a cation sublattice, it is free of Nb Li defects. Intensity of this band first decreases, than increases in LiNbO 3 :B with an increase in a dopant concentration, Figure 5. The band is split in two components 112 and 123 cm −1 in a spectrum of a crystal LiNbO 3 :B (1.24 mol%), Figure 5. This indicates a high structure perfection of the crystal. A decrease in the 120 cm −1 band intensity also indicates a decrease in Nb Li defects amount at its doping by B 3+ ions. Cation sublattice orders (Raman bands in the area 150-300 cm −1 narrow, Table 1) simultaneously with a decrease in resonant anharmonic interaction between the lowest-frequency fundamental A 1 (TO) vibrations 254, 274 cm −1 and two-particle A 1 (TO) acoustic excitations. The latter is discussed in detail in [49]. A degree of mixing between A 1 (TO) single-phonon and multi-phonon states depends on the value of this interaction, which also influences the spectrum in the 120 cm −1 area. Thus, LN crystals doping can effectively govern interaction between single-phonon and two-phonon states. Paper [50] has revealed that 120 cm −1 band intensity is connected with an acoustic Q-factor of LN crystal. The higher Q-factor is, the smaller the band intensity is, due to a decrease in the amount of Nb Li defects. Thus, our data prove a higher structure perfection of nominally pure LN crystals grown from boron doped charge compared to that of a CLN crystal. The crystal LiNbO 3 :B (1.24 mol%) is the most structurally perfect one.

Study of LiNbO 3 :B Crystal Structure by IR Spectroscopy in the Area of OH-Groups Vibration
LN single crystals obtained at air atmosphere always contain OH groups incorporated into the structure during growth [14]. Localization of hydrogen atoms in LN structure greatly influences charge distribution and polarizability of oxygen octahedra clusters MeO 6 (Me: Li, Nb, dopant). Hydrogen atoms bond with oxygen of LN structure by hydrogen bonds. OH groups increase LN crystal conductivity and decrease optical damage and coercive field value [10,29]. Doping of LN crystal leads to a total restructure of the whole hydrogen bonds system.
Stretching and deformation vibrations of OH groups are supposed to manifest in Raman and IR absorption spectra of LN crystals in the areas 1600-1800 and 3450-3550 cm −1 , respectively. The amount of bands in these areas and their basic parameters depend on localization of hydrogen atoms in the crystal structure. The localization in its turn depends on stoichiometry, doping particularities and growth details of a single crystal [14]. A perfect strictly stoichiometric LN crystal has no site for a hydrogen atom [51]. IR spectrum of a real NSLN crystal contains only one narrow (S = 3 cm −1 ) 3466 cm −1 band [26,52]. A slight shift of stoichiometry leads to widening and splitting of this band in IR absorption spectrum into 3466 and 3480 cm −1 bands. This is caused by a formation of complex defects Nb Li -OH, V Li , etc., in a structure of a non-stoichiometric LN crystal [51]. CLN crystal IR spectrum contains three components of the band: and intense 3466 cm −1 band, an average intense 3481 cm -1 band and a low-intense 3489 cm −1 band [53]. Papers [29,52] have demonstrated that as LN structure approaches stoichiometric composition, width of the 3466 cm −1 band and intensity of 3481 and 3489 cm −1 bands decrease. Note that SLN crystal has a ratio Li/Nb = 1, point defects, such as the same cations in the neighbor sites are absent. Thus, studies of parameters of IR absorption bands in the OH-groups vibrations area contains important information about the nature of defects and their localization in the structure; the structural rearrangements that occur in LN crystals with a change in composition; and the peculiarities of a single crystal growth. Such information is also important at refinements of technologies of growing optically and structurally perfect LN single crystals. Figure 7 demonstrates IR absorption spectra of studied LN crystals in the area 3420-3550 cm −1 . This area contains stretching vibrations of OH groups. Table 2 demon-strates parameters of spectral bands. It is obvious that spectra of all studied crystals are different. An SLN crystal spectrum contains three bands with frequencies 3465, 3480, and 3488 cm −1 , Figure 7. The bands in an SLN spectrum are much narrower than in other crystals, Table 2. It is believed that a band in the frequency range 3465-3466 cm −1 is characteristic of a highly perfects SLN crystals with only one hydrogen atom position, [29,52]. Figure 7 clearly shows that a CLN spectrum is different from the one of an SLN crystal. IR spectra of a CLN is a blurred wide absorption band consisting of several components with the same polarization with frequencies 3467 (ν 1 ), 3483 (ν 2 ), and 3486 (ν 3 ) cm −1 and a weak band with a frequency~3490 cm −1 (ν 4 ). A 3480-3485 cm −1 band is considered characteristic of a CLN crystal and is attributed to stretching vibrations of a V Li -OH complex [54]. V Li defects are absent from SLN crystal structure. This means that IR spectrum should not contain bands with frequencies 3480-3485 cm −1 , which is confirmed experimentally [29,52]. IR spectra of LiNbO 3 :B (0.55-0.83 mol%) is similar to those of CLN: an OH-groups absorption band is split in three components of the same polarization~3470, 3483, and~3486 cm −1 . Thus, hydrogen bonds in CLN and LiNbO 3 :B crystals structure should be similar. Three components of CLN crystal are considered to be connected with stretching vibrations of OH groups located near Nb Li 4+ -V Li − defects [53,55]. However, intensities of all three components are higher in LiNbO 3 :B than in CLN spectrum, Table 2. At the same time widths of two bands are narrower than in CLN spectrum, widths of the band 3485 cm −1 is not narrower, Table 2. This indicates ordering of OH groups distribution in LiNbO 3 :B crystals structure. Figure 7 demonstrates IR absorption spectra of studied LN crystals in the area 3550 cm -1 . This area contains stretching vibrations of OH groups. Table 2 demon parameters of spectral bands. It is obvious that spectra of all studied crystals are di An SLN crystal spectrum contains three bands with frequencies 3465, 3480, and 34 1 , Figure 7. The bands in an SLN spectrum are much narrower than in other crystals 2. It is believed that a band in the frequency range 3465-3466 cm -1 is characteris highly perfects SLN crystals with only one hydrogen atom position, [29,52]. Fi clearly shows that a CLN spectrum is different from the one of an SLN crystal. IR of a CLN is a blurred wide absorption band consisting of several components w same polarization with frequencies 3467 (ν1), 3483 (ν2), and 3486 (ν3) cm -1 and a wea with a frequency ~3490 cm -1 (ν4). A 3480-3485 cm -1 band is considered characteris CLN crystal and is attributed to stretching vibrations of a VLi-ОН complex [54]. VLi are absent from SLN crystal structure. This means that IR spectrum should not c bands with frequencies 3480-3485 cm -1 , which is confirmed experimentally [29 spectra of LiNbO3:B (0.55-0.83 mol%) is similar to those of CLN: an OH-groups abso band is split in three components of the same polarization ~3470, ~3483, and ~348 Thus, hydrogen bonds in CLN and LiNbO3:B crystals structure should be similar components of CLN crystal are considered to be connected with stretching vibrat OH groups located near NbLi 4+ -VLi − defects [53,55]. However, intensities of all thre ponents are higher in LiNbO3:B than in CLN spectrum, Table 2. At the same time of two bands are narrower than in CLN spectrum, widths of the band 3485 cmnarrower, Table 2. This indicates ordering of OH groups distribution in LiNbO3:B c structure.   A computer simulation has been carried out in a paper [56]. The paper considers possible sites of hydrogen atom in SLN crystals. The paper has shown that all hypothetical hydrogen sites on O-O bonds located near V Li defects are unstable and should be excluded from further consideration. Moreover, the hydrogen atom site on C axis has not been confirmed. Analysis of possible hydrogen atoms trajectories in LN crystal structure shows that it can be associated in SLN structure only with one oxygen ion in the upper oxygen triangle of NbO 6 octahedron. Dipole moment orientation of OH-group is slightly different from a direction of a short O-O bond. Thus, creation of a hydrogen bond leads not only to a cardinal change in the wave functions of the outer electron orbitals of the oxygen ion and the parameters of its electronic polarizability, but also to a strong distortion of the entire NbO 6 octahedron [56]. Such changes in a crystal structure should manifest in 850-900 cm −1 area of Raman spectra. Bonds in this area correspond to stretching bridge vibrations B-O-B (B: Nb or dopant) of oxygen in NbO 6 octahedra along the polar axis [2]. Table 1 shows that spectra of LiNbO 3 :B and nominally pure SLN, CLN crystals contain only one band in the area 850-900 cm −1 . The width of this band depends on the Li/Nb value. This band width is minimal in SLN and maximal in CLN crystal spectrum, Table 1 Widening of bands in the 850-900 cm −1 area is determined by the fact that Nb cations are located in different sites in octahedra, their concentration is quite high and they interact weakly between each other [2]. Changes in LN crystals Raman spectra in the area of stretching bridge vibrations correlate with changes observed in IR spectra, Figure 7. This result confirms conclusions of a theoretical work [56] and experimental data [29,52].
Papers [57,58] have clearly indicated that the ratio I rel = I 3480 /I 3465 between intensities of 3480 and 3465 cm −1 IR bands almost linearly decreases with an increase in R = Li/Nb ratio. Note that I rel = 0 if the 3480 cm −1 band is absent from the spectrum. Evaluation of I rel allows one to determine R with higher precision (0.01 mol% Li 2 O) than other methods [57,58]. For example, determination of Li 2 O concentration due to Curie temperature has an error 0.1 mol% Li 2 O. A ratio I rel = I 3480 /I 3465 = 0.78 for an SLN crystal studied in our research. This means that R = Li/Nb = 0.98. This result is different from data obtained by Raman spectroscopy: at R = 1 a low-intense 120 cm −1 band should be absent from SLN spectrum, Figure 5. Changes can be explained by the following: 120 cm −1 Raman band is sensitive only towards the R = Li/Nb ratio (which automatically include all point and complex defects); 3466 cm −1 IR band is sensitive only towards Nb Li -OH complex defects. This is why R determined due to IR spectra is smaller, than that determined due to Raman spectra.
A Li/Nb ratio can also be calculated with a high precision (up to 0.01 mol%) from a fundamental absorption edge [59]. However, this kind of experiment should be carried out at fine plates (0.5-1.0 mm) of LN crystals. Moreover, studied samples must not have additional absorption bands in the fundamental absorption edge, as it is in crystals doped by multiply charged cations Fe, Cu, etc. Empirical formulae for Li/Nb evaluation in nominally pure LN crystals were suggested in [59]. The formulae apply data on fundamental absorption edge: where λ 20 -fundamental optical absorption edge corresponding to a wavelength at an absorption coefficient α = 20 cm −1 .
Knowing Li/Nb value, one can calculate concentration of point defects in LN crystal. Due to a Li-vacancies compensation model [2], concentration of Nb Li 4+ and V Li − point defects in mol% in a crystal lattice of nominally pure LN crystals can be calculated due to a formula [59]: Paper [59] offers a model for determination of Li/Nb ratio in crystals that have absorption bands directly in the region of the fundamental absorption edge. Ratio of the integrated intensity of the band with a maximum at 3466 cm −1 (A int, I I ) to the total integrated intensity of the IR absorption spectrum (A int, I ) (A int, I I /A int, I = α 1 (ν)dv/ α(ν)dv) was determined for each crystal's spectrum. In the formula α(ν)-absorption coefficient in dependence on wavelength in cm −1 .
A Li/Nb ratio can be determined in crystals doped by photochrome or photorefractive cations absorbing in visible or UV spectral area due to the obtained A int, I I /A int, I value and Li/Nb ratio [59]. This method can obviously be applied for other cases; for example, for LN crystal samples with a width >1 mm, when fundamental optical absorption edge λ 20 is impossible to determine. In this paper we have used samples with a width Z = 4 mm. Our calculations concerning Li/Nb and Nb Li 4+ and V Li − point defects concentration obtained due to the abovementioned formulae are demonstrated in Table 2.
A Klauer [60] method can be applied to calculate concentration of OH-groups: where C(OH)-concentration of OH defects in a sample, (cm −3 ), A int, I -a total integral intensity of IR spectrum in the OH stretching vibrations area (A int, I = α(ν)dv, α(ν)-absorption coefficient in dependence on wavelength in cm −1 ), a OH -intensity if ion absorption = (9.125 ± 1.369) × 10 −18 cm [60,61]. Calculation results are demonstrated in Table 2. Concentration of OH-groups is minimal in SLN and maximal in LiNbO 3 :B (0.55 and 0.83 mol%) crystal, A decrease in the amount of bonded OH groups apparently leads to an increase in free protons in LN crystals, which can influence electrical conductivity of the crystals and decrease photorefractive effect.

Study of LiNbO 3 :B Crystals Optical Uniformity by PILS and Laser Conoscopy
We have studied optical uniformity of crystals by laser conoscopy at laser radiation (λ 0 = 532 nm) power p = 1 and 90 mW, Figure 8. When conoscopic patterns of LN crystals at low power differ from pattern of a perfect crystal, it means that changes are caused by defects, mechanical stresses, or compositional heterogeneity. Conoscopic pattern distortions at high-power radiation are caused by laser-induced static and fluctuating defects. Laser conoscopy and PILS method are unable to provide direct information about crystals intrinsic structure and defects type. However, they detect information on the optical homogeneity and photorefractive properties of crystals.

Study of LiNbO3:B Crystals Optical Uniformity by PILS and Laser Conoscopy
We have studied optical uniformity of crystals by laser conoscopy at laser radiation (λ0 = 532 nm) power p = 1 and 90 mW, Figure 8. When conoscopic patterns of LN crystals at low power differ from pattern of a perfect crystal, it means that changes are caused by defects, mechanical stresses, or compositional heterogeneity. Conoscopic pattern distortions at high-power radiation are caused by laser-induced static and fluctuating defects. Laser conoscopy and PILS method are unable to provide direct information about crystals intrinsic structure and defects type. However, they detect information on the optical homogeneity and photorefractive properties of crystals.   Figure 8. A "Maltese cross" is deformed in the center, fragments are shifted from its center; angle between branches is not 90°; isochromes are slightly oval, however their geometry is not distorted. Appearance of anomalous biaxiality can be caused by the existence of local micro-regions and clusters with birefringence in the crystal. In addition, observed distortion in LiNbO3:В (0.83 mol%) crystal conoscopic patterns can be associated with an increase in the effect of photorefraction, Figure 8.   Figure 8. A "Maltese cross" is deformed in the center, fragments are shifted from its center; angle between branches is not 90 • ; isochromes are slightly oval, however their geometry is not distorted. Appearance of anomalous biaxiality can be caused by the existence of local micro-regions and clusters with birefringence in the crystal. In addition, observed distortion in LiNbO 3 :B (0.83 mol%) crystal conoscopic patterns can be associated with an increase in the effect of photorefraction, Figure 8.
PILS is a direct consequence of the photorefractive effect. It arises in a ferroelectric crystal on defects induced by exciting radiation with fluctuating refractive index, dielectric constant, and other physical parameters [62,63]. Laser radiation scattered by such defects interferes with pumping, forming a complex pattern of minima and maxima of the scattered light intensity (speckle structure), Figure 8. It causes strong laser beam destruction in the crystal; it is an interfering factor for the radiation generation and conversion [4]. The shape and features of the PILS indicatrix speckle structure depend on the crystal structure, the state of its defectiveness, as well as on the radiation polarization and the experiment geometry [62,64]. The time and opening angle of the PILS indicatrix determine the response speed of the electro-optical modulators and gates. PILS parameters are determined by the depth of traps, the mobility of electrons responsible for the magnitude of the photorefractive effect, as well as the interaction of laser radiation with a defective crystal [62]. Figure 8 demonstrates PILS patterns of studied crystals. The time of full PILS indicatrix speckle structure opening for the above crystals is 60 s. According to PILS data, a crystal LiNbO 3 :B (1.24 mol%) possesses the highest optical damage resistance. PILS speckle structure indicatrix of this crystal does not open even at a high excitation radiation power 160 mW, despite other crystals, Figure 8. Conoscopic and PILS patterns show that structure uniformity of LiNbO 3 :B is similar to that of a CLN crystal and is much higher than that in a SLN crystal, Figure 8. Note that photorefractive effect is much higher in SLN than in CLN and LiNbO 3 :B (0.55-1.24 mol%) crystals. An SLN crystal PILS opening angle θ is 56 • ; it does not exceed 22 • in LiNbO 3 :B (0.55-1.24 mol%) crystals, Figure 6, Table 3. It is a known fact that SLN crystals have non-uniform refractive index along the growth axis, while in CLN crystals it is stable. Photorefractive effect manifests itself much stronger in SLN crystals, than in CLN crystals [2,65]. A change in the structure and properties of a melt, as well as crystal doping, adds levels in the band gap and thereby changes the photorefractive properties and electrical conductivity of the crystal. The crystal LiNbO 3 :B (0.83 mol%, direct doping) is photorefractive, has a pronounced PILS and is characterized by a uniform refractive index along the growth axis, as evidenced by a clear conoscopic pattern, Figure 8. At the same time, the crystal LiNbO 3 :B (1.24 mol%, homogeneous doping) has a reduced effect of photorefraction; a more blurred conoscopic pattern indicates the presence of fluctuations in the refractive index along the polar axis. A similar, but much more diffuse conoscopic pattern is observed for SLN crystal, Figure 8. At the same time, a discontinuous structure of the laser beam is observed when it propagates along the polar axis in SLN, (Figure 21 in [66]). This indicates a much higher inhomogeneity of the refractive index in SLN in comparison with the LiNbO 3 :B (1.24 mol%, homogeneous doping) crystal. For other crystals under study, a discontinuous structure of the laser beam has not been observed yet.
A decrease in the photorefractive effect can be due to a change in the band structure and electro-optical properties of the crystal due to the structuring of the melt with boron. A low photorefractive effect value in the LiNbO 3 :B (1.24 mol%, homogeneous doping) crystal can be explained by a change in the electro-optical properties of the crystal. Due to works [67,68], a change in electro-optical coefficients in nominally pure LN crystals at a change in Li/Nb ratio leads to a change in ionic contribution to the magnitude of the electro-optical effect. We assume that boron influence leads to such a structuring of a melt that crystals grown from the melt have a more strict crystal structure compared to CLN crystals, in which MeO 6 octahedra have less ability to deform.  change in Li/Nb ratio leads to a change in ionic contribution to the magnitude of the electro-optical effect. We assume that boron influence leads to such a structuring of a melt that crystals grown from the melt have a more strict crystal structure compared to CLN crystals, in which МеO6 octahedra have less ability to deform. Figure 9 demonstrates optical absorption spectra of SLN, CLN, and LiNbO3:В (0.55, 0.83 mol%) crystals. The absorption edge of LiNbO3:B crystals is shifted to longer wavelengths compared to the absorption edge of SLN, CLN crystals. At this, absorption spectra of LiNbO3:B crystals has a steeper rise in comparison with SLN. This indicates a higher optical uniformity of crystals grown from a boron-doped charge. It is believed that NbLi defects bring the main contribution to the change in the position of the fundamental absorption edge in nominally pure CLN crystals [2,59]. Thus, we can draw conclusions about the mechanism of non-photorefractive impurities incorporation into the crystal lattice. Coincidence of the concentration dependences of NbLi defects and the Urbach parameter on the Li/Nb ratio in the crystal has been spotted in a paper [59]. This agrees with the conclusions of [69] that the Urbach absorption occurs in the LN crystal as a result of the transition of electrons from the filled states "2p" of oxygen to the empty states "4d" of niobium. A decrease in the concentration of NbLi defects as the Li/Nb ratio approaches unity leads to a fundamental absorption edge shift toward longer wavelengths as a result of a decrease in the density of local states associated with NbLi defects near the bottom of the conduction band. The transition of electrons from the valence band to these states is possible in LN under the light. Figure 9 demonstrates that CLN absorption edge is maximally shifted to short-wavelength region. Table 2 shows that for CLN Li/Nb = 0.942, С(VLi) = 3.90 mol% and С(NbLi) = 0.98 mol%. It is a known fact that CLN crystals have Li deficiency near 6 mol% (Li/Nb = 0.946) [2]. According to the compensation of Li-vacancies model, the crystal lattice of a CLN crystal contains ~1 mol% NbLi 4+ point defects and ~4 -4+ Coincidence of the concentration dependences of Nb Li defects and the Urbach parameter on the Li/Nb ratio in the crystal has been spotted in a paper [59]. This agrees with the conclusions of [69] that the Urbach absorption occurs in the LN crystal as a result of the transition of electrons from the filled states "2p" of oxygen to the empty states "4d" of niobium. A decrease in the concentration of Nb Li defects as the Li/Nb ratio approaches unity leads to a fundamental absorption edge shift toward longer wavelengths as a result of a decrease in the density of local states associated with Nb Li defects near the bottom of the conduction band. The transition of electrons from the valence band to these states is possible in LN under the light. Figure 9 demonstrates that CLN absorption edge is maximally shifted to short-wavelength region. Table 2 shows that for CLN Li/Nb = 0.942, C(V Li ) = 3.90 mol% and C(Nb Li ) = 0.98 mol%. It is a known fact that CLN crystals have Li deficiency near 6 mol% (Li/Nb = 0.946) [2]. According to the compensation of Li-vacancies model, the crystal lattice of a CLN crystal contains~1 mol% Nb Li 4+ point defects and 4 mol% V Li − point defects [54,70]. At the same time, Nb Li 4+ is completely absent from a perfect SLN crystal structure. LiNbO 3 :B (0.55 and 0.83 mol%) crystals have ordering of structural units of the cation sublattice and the content of Nb Li and V Li point defects between SLN and CLN [24]. Approximation to stoichiometry in LiNbO 3 :B crystals is possible by binding excess niobium in a congruent melt due to the complexing ability of boron compounds [24]. Obviously, the position of the fundamental absorption edge is determined not only by the concentration of Nb Li and V Li point defects. The number of shallow electron traps depends significantly on the Li/Nb ratio [2,7,12,71]. Under laser radiation a spatial charge separation occurs and an internal electric field arises, leading to a photoinduced change in refractive indices in LN crystal as a result of photoexcitation processes (photovoltaic and diffusion current) [7,12,63,71]. Photovoltaic E pv and diffusion E D fields presented in Table 3 were calculated due to formulae [31]:

Investigation of Photovoltaic Fields and Bandgap in LiNbO3:В Crystals by PILS and Optical Spectroscopy
where E pv -photovoltaic field, E D -diffusion field, λ-wavelength, θ in s -scattered radiation angle, Γ −c and Γ +c -amplification factors (indices "−" and "+" indicate the direction of the scattered radiation against and along the direction of the crystal polar axis, respectively), n e and n o -refractive indices of the extraordinary and ordinary rays, respectively, r 33 and r 51 -LN electro-optical coefficients.
An amplification factor Γ θ in s can be calculated in dependence on PILS opening angle due to a formula [31]: where I s -scattered radiation intensity, I Ω so -primary scattering intensity (incident beam) , l e f f -The effective interaction interval, which is calculated depending on the scattering angle, according to the following formulas [31]: where d-crystal thickness, w p -laser beam diameter. Figures 10 and 11 demonstrate angular distribution of the scattered radiation intensity at different wavelengths and dependences E D (a) and E pv (b) on the wavelengths in studied LN crystals. The band gap for the studied crystals is in the range 3.37-3.78 eV, Table 3. The least band gap was observed in LiNbO 3 :B and SLN, Table 3, the greatest-in CLN (3.78 eV). Due to [72], a pure CLN crystal band gap is 3.72 eV, which is close to the value characteristic of wide-gap semiconductors.
Moreover, LN crystal grown at air always contain 10 16 -10 18 cm -3 protons bonded with oxygen by a hydrogen bond. Hydrogen atoms form such complex defects as VLi-OH, NbLi-OH, etc., [7,14,15]. OH-groups play important role in formation of a secondary defect   A predominant photorefraction mechanism in a LN crystal is photovoltaic, which means that the value of the photovoltaic field (Epv) is much greater than that of the diffusion field (ED) [71]. Table 3, Figures 10 and 11 show: an increase in photovoltaic field in a LN crystal is accompanied with an increase in PILS scattering angle θ and an asymmetry of its indicatrix. The photovoltaic field Epv determines the magnitude of the induced birefringence (the effect of photorefraction). SLN has the least Еpv at 476.5, 488.0, 514.5, and 530.9 nm excitation. Thus, a diffusion charge transport mechanism is manifested in SLN stronger than in other crystals. This picture is characteristic of crystals with a large number of shallow electron traps.
A dependence of Еpv has a maximum at excitation by 514.5 nm for all studied crystals, but LiNbO3:В (0.55 mol%), Figure 11. At the same time, Еpv values of LiNbO3:В (0.55 and 0.83 mol%) crystals are close: 5458 and 5554 V/cm, respectively. Curious that ED of a LiNbO3:В (0.83 mol%) crystal is close to that of a CLN. This may indicate an almost equal number of shallow electron traps in these crystals. ED value of a LiNbO3:В (0.55 mol%) crystal is larger than that of CLN and LiNbO3:В (0.83 mol%) crystals; at the same time, its PILS indicatrix opening angle is minimal, Table 3. A predominant photorefraction mechanism in a LN crystal is photovoltaic, which means that the value of the photovoltaic field (E pv ) is much greater than that of the diffusion field (E D ) [71]. Table 3, Figures 10 and 11 show: an increase in photovoltaic field in a LN crystal is accompanied with an increase in PILS scattering angle θ and an asymmetry of its indicatrix. The photovoltaic field E pv determines the magnitude of the induced birefringence (the effect of photorefraction). SLN has the least E pv at 476.5, 488.0, 514.5, and 530.9 nm excitation. Thus, a diffusion charge transport mechanism is manifested in SLN stronger than in other crystals. This picture is characteristic of crystals with a large number of shallow electron traps.
A dependence of E pv has a maximum at excitation by 514.5 nm for all studied crystals, but LiNbO 3 :B (0.55 mol%), Figure 11. At the same time, E pv values of LiNbO 3 :B (0.55 and 0.83 mol%) crystals are close: 5458 and 5554 V/cm, respectively. Curious that E D of a LiNbO 3 :B (0.83 mol%) crystal is close to that of a CLN. This may indicate an almost equal number of shallow electron traps in these crystals. E D value of a LiNbO 3 :B (0.55 mol%) crystal is larger than that of CLN and LiNbO 3 :B (0.83 mol%) crystals; at the same time, its PILS indicatrix opening angle is minimal, Table 3.

Chemical Interactions in Systems Li 2 O-Nb 2 O 5 , Li 2 O-B 2 O 3 -Nb 2 O 5 ; Particularities of LN Crystallization from Melts with Non-Metal Dopants
Solid-phase chemical reactions with the formation of various borates occur in the process of high-temperature synthesis of the charge in the Li 2 O-B 2 O 3 -Nb 2 O 5 system. To evaluate the melt-crystal system, it is advisable to use the distribution coefficients K D and the ∆C parameter. K D is the ratio of the impurity concentration in the crystal at the initial moment of growth and the impurity concentration in the melt, ∆C characterizes the compositional uniformity of the crystal; it is defined as the difference between the dopant concentration at the cone and the end of the crystal. A crucial particularity of boron is its ability to change LiNbO 3 :B melt and crystal structure without incorporation to the crystal. Its effective distribution coefficient (K D ) is only~3 × 10 −4 [19,20,25]. Such small values of the boron distribution coefficient in the Li 2 O-Nb 2 O 5 -B 2 O 3 system are due to the fact [73] that boron should not incorporate into the LiNbO 3 :B crystal, since the LN phase has no region of boron solubility in the solid state, Figure 12. The small amount of boron that does incorporate into the LN structure can be either be mechanically captured or incorporate into the tetrahedral voids of the crystal structure. Note that boron B 3+ has a small ion radius. Quasi-double diagrams are polythermal sections of the general quasi-ternary diagram of the state of systems Li 2 O-Nb 2 O 5 -B 2 O 3 : LiNbO 3 -LiBO 2 , and LiNbO 3 -Li 3 B 2 O 4.5 . These cross sections have a simple form, since they contain only one two-phase eutectic transformation each and have no solubility regions, Figure 12, [73]. This type of phase diagram is generally favorable for the quality of the growing LiNbO 3 :B crystal since crystallization proceeds from the two-phase region and the only crystallizing phase is LN. If we disregard the radical change in the structure of the melt upon doping with boron, crystal grown from such a melt should chemically and structurally correspond to a CLN grown from an undoped melt. Wherein studied LiNbO 3 :B crystals have Nb Li defects amount close to that of an SLN crystal. Moreover, the LiNbO 3 :B crystal has a more ordered structure than CLN [25]. Note that Nb Li are the deepest electron traps responsible for the photorefractive effect. The reason for this is apparently an electronic structure of boron: it is a strong complexing element due to 1 electron at a p-orbital. The electron provides high ionization energies and electronegativity values at a small ionic radius, which leads to a significant change in the structure of the boron-containing melt.

Chemical Interactions in Systems Li2O-Nb2O5, Li2O-B2O3-Nb2O5; Particularities of LN Crystallization from Melts with Non-Metal Dopants
Solid-phase chemical reactions with the formation of various borates occur in the process of high-temperature synthesis of the charge in the Li2O-B2O3-Nb2O5 system. To evaluate the melt-crystal system, it is advisable to use the distribution coefficients КD and the ∆С parameter. КD is the ratio of the impurity concentration in the crystal at the initial moment of growth and the impurity concentration in the melt, ∆С characterizes the compositional uniformity of the crystal; it is defined as the difference between the dopant concentration at the cone and the end of the crystal. A crucial particularity of boron is its ability to change LiNbO3:В melt and crystal structure without incorporation to the crystal. Its effective distribution coefficient (КD) is only ~3·10 −4 [19,20,25]. Such small values of the boron distribution coefficient in the Li2O-Nb2O5-B2O3 system are due to the fact [73] that boron should not incorporate into the LiNbO3:В crystal, since the LN phase has no region of boron solubility in the solid state, Figure 12. The small amount of boron that does incorporate into the LN structure can be either be mechanically captured or incorporate into the tetrahedral voids of the crystal structure. Note that boron В 3+ has a small ion radius. Quasi-double diagrams are polythermal sections of the general quasi-ternary diagram of the state of systems Li2O-Nb2O5-B2O3: LiNbO3-LiBO2, and LiNbO3-Li3B2O4.5. These cross sections have a simple form, since they contain only one two-phase eutectic transformation each and have no solubility regions, Figure 12, [73]. This type of phase diagram is generally favorable for the quality of the growing LiNbO3:В crystal since crystallization proceeds from the two-phase region and the only crystallizing phase is LN. If we disregard the radical change in the structure of the melt upon doping with boron, crystal grown from such a melt should chemically and structurally correspond to a CLN grown from an undoped melt. Wherein studied LiNbO3:В crystals have NbLi defects amount close to that of an SLN crystal. Moreover, the LiNbO3:В crystal has a more ordered structure than CLN [25]. Note that NbLi are the deepest electron traps responsible for the photorefractive effect. The reason for this is apparently an electronic structure of boron: it is a strong complexing element due to 1 electron at a p-orbital. The electron provides high ionization energies and electronegativity values at a small ionic radius, which leads to a significant change in the structure of the boron-containing melt.  Table 4 demonstrates ∆С and КD for exact B, Mg, Zn, and Ce dopants concentrations in the melt. The necessity to use КD and ∆С is substantiated in detail in the works [27,74]. The distribution coefficient КD for doped (B, Zn, Mg, and Ce) LN crystals is clearly a function of the electronic structure of the dopant; the ∆С case is not so obvious, Table 4. The ∆С parameter value is quite small for all studied dopants except Zn. In fact, the value is comparable to the error of the method for determining the dopant concentration, Table 4.  Table 4 demonstrates ∆C and K D for exact B, Mg, Zn, and Ce dopants concentrations in the melt. The necessity to use K D and ∆C is substantiated in detail in the works [27,74]. The distribution coefficient K D for doped (B, Zn, Mg, and Ce) LN crystals is clearly a function of the electronic structure of the dopant; the ∆C case is not so obvious, Table 4. The ∆C parameter value is quite small for all studied dopants except Zn. In fact, the value is comparable to the error of the method for determining the dopant concentration, Table 4.  Table 4 data reveal high concentration uniformity in dopant distribution along the polar axis at a K D coefficient strongly different from the unity. This is especially characteristic of LiNbO 3 :B crystals, in which K D << 1 and ∆C~10 −6 wt%. The explanation of this fact can be the following. The initial LN charge composition ceases to be congruent with the addition of dopant, Figure 12. This, in particular, means that the spectrum of variations of ionic complexes in the melt in terms of structure and components greatly increases. The melt is trapped in a limited zone near the crystallization front, which has a constant temperature. Consequently, those ionic complexes for which this temperature is solidus (T c1 ) will crystallize, naturally, taking into account some supercooling. As part of the volume of the melt is consumed, the concentration of the dopant (in our case, boron) in the remaining melt increases and the ratio of the concentrations of various ionic complexes changes. Accordingly, the fraction of ionic complexes for which the temperature T c1 is the solidus temperature will also change, it will become noticeably less. This will continue until the system reaches a certain critical state, at which the concentration of complexes with solidus temperature T c1 is insufficient for the growth of a crystal with a constant dopant concentration. This circumstance limits the fraction of the melt that can be crystallized to obtain a compositionally uniform LiNbO 3 :B crystal. Upon reaching the described critical state, the further behavior of the system can be different: from a significant change in the dopant concentration distribution along the growth axis to defects such as cellular growth and crystallization of a different composition phase, Figure 3 [19,27,74]. Many factors affect the conditions for reaching such a critical state, and hence the possible dimensions of a LiNbO 3 :B crystal of high optical quality with a constant dopant concentration distribution in the crystal boule volume. This is the composition of the melt, the thermodynamics of the initial components, the graphical expression of which is the phase diagram, the structure of the melt consisting of ionic complexes with different thermodynamic and kinetic characteristics, and even the technical capabilities (sensitivity and reaction time constant) of the control and monitoring system of the growth process.
An important conclusion follows from the above reasoning for the technology of doped lithium niobate crystals: LN crystals with a uniform distribution of dopant can be grown by Czochralski at K D that is noticeably different from unity during crystallization of only a well-defined part of the melt. Consequently, the length and diameter of such a LN crystal boule is limited. In different systems, the limiting sizes of such structurally and compositionally uniform crystals will differ markedly. Moreover, for different systems, there may be different physical and chemical reasons for this, expressed in the difference in the thermodynamic parameters of the systems, that is, in the form of phase diagrams. LN has no regions of homogeneity with either boron or its compounds in the boron-containing system Li 2 O-Nb 2 O 5 -B 2 O 3 . As the crystal grows, only the pure LN phase crystallizes, Figure 12. However, at the same time, there will be a significant increase in the boron content in the melt and, as a consequence, a decrease in the crystallization temperature and a rather radical increase in the melt viscosity, which limits convective flows, Figure 12. C 0 , Crystals 2021, 11, 458 22 of 37 C 1 and C 2 -the composition of the melt or the concentration of boron in the melt, T 1 and T 2 are the crystallization temperature corresponding to the composition of the C 1 and C 2 . All this, most likely, limits the maximum possible boron concentration in the melt to 0.12 wt%, since an increase in the boron concentration above this critical concentration in the melt leads to cellular growth and other irreparable defects (for example, defects in the form of "channels") in LiNbO 3 :B crystals, Figures 1-3, [19,25]. Our studies have shown that the fraction of the melt, the crystallization of which leads to the production of an optically uniform LiNbO 3 :B crystal, does not exceed 17-18% [19,25].
The above reasoning is valid for evaluating the technological growth parameters when doping lithium niobate with any dopants. For example, the phase diagram of the Li 2 O-Nb 2 O 5 -MgO system is less complex than that of the system Li 2 O-Nb 2 O 5 -ZnO [75,76]. The Zn diagram has a greater number of phases and they are less stable; the concentration regions of existence of different phases and mixtures of phases are less extended [75]. These physical and chemical particularities of a Zn system explain a greater number of concentration thresholds at LN doping with Zn rather than Mg. The particularities also explain LiNbO 3 :Zn crystals tend to cracking and separation of impurity phases during growth [27]. Phase stresses in LiNbO 3 :Zn crystals are also a consequence of Zn system particularities. The stresses are the reason why these crystals claim special modes of aftergrowth treatment and turning them to a single domain state [27]. Thus, a volume of a structurally and chemically uniform boule and the fraction of the crystallized melt for LiNbO 3 :Zn will be less than for LiNbO 3 :Mg crystals. The fraction of crystallized melt in the case of Li 2 O-Nb 2 O 5 -ZnO system is ≤~20%, Li 2 O-Nb 2 O 5 -MgO system ≤~30% [16,27,74]. The fraction in the case of Li 2 O-Nb 2 O 5 -CeO 2 system strongly depends on the initial dopant concentration in the melt; the crystallized melt fraction is significantly lower than in systems with magnesium and zinc. These facts are caused by significantly lower distribution coefficients K D , Table 4. This leads to a fast increase in the dopant concentration in the melt and achieving of the critical state at which the complexes concentration with solidus temperature T c1 is insufficient for crystal growth with a constant dopant concentration. This is why only ≤~12% fraction of the melt should crystallize to turn into a structurally and chemically uniform LiNbO 3 :Ce crystal [76].
Results obtained in [77] based on the analysis of the dopant electronic structure are important for predicting the technological conditions of growth and the quality control of doped LN crystals. At this, p-elements (non-metals, in our case, boron), which have a higher chemical activity due to the greater number of valence electrons than that of metals, make the melt more homogeneous at the level of ionic complexes. It is exactly this fact that allows us to grow LN crystals from a congruent melt with a defective structure close to one in SLN. Metals (s and d elements, Zn and Mg) have a similar effect on the LN melt and crystals properties. Nevertheless, Zn-doped crystals have an increased ununiformity of dopant incorporation to the crystals, tend to cracking due to phase stresses; they also are characterized by additional claims towards growth conditions and after-growth electrothermal treatment of boules. The f-element metals, due to their electronic structure, form the structure of the melt in such a way that only a relatively small fraction of the melt can be crystallized to obtain a compositionally uniform crystals, for example, LiNbO 3 :Ce.
LN is characterized by a high melting temperature (~1526 K) and a high chemical activity of the melt, which requires platinum crucibles. This is the reason for the extremely small number of works devoted to the study of the structure of even nominally pure LN melts. In recent years, the need has arisen to obtain optically highly perfect LN single crystals (nominally pure and doped). The existence in the melt of strongly bound groups of atoms or ions (clusters) of a certain structure was considered as the main concept in the literature [78]. In [79][80][81], high-temperature Raman spectroscopy was effectively used to study the crystallization from a melt. In the pre-crystallization temperature of a LN melt, the Raman spectra revealed a noticeable change in the structure of the melt and a mismatch between the structure of anionic motifs in the melt and crystal. The mismatch effect near the melting point creates serious obstacles to the nucleation of equilibrium structures [79][80][81]. Metastable phases have an advantage in crystallization under these conditions. The phases crystal lattice is very far from the structure of the short-range order of the melt. The mismatch in the anionic structure of the melt and crystal can also affect the growth of crystal faces or the formation of structural defects [79,80].
LN structure is a framework of oxygen octahedra O 6 , articulated by vertices and faces, with the densest hexagonal packing. Inside the octahedra there are intrinsic (Li + and Nb 5+ ) or doping cations. The bonds in the octahedra occupied by niobium are predominantly covalent [1]. At the same time, the lithium ion is bound to oxygen atoms only by electrostatic interaction. Due to the predominance of the covalent bond type, the niobium cation tends to form in the melt anionic motifs consisting of NbO 4 tetrahedra. That is, when LN melts, anionic structure rearranges significantly: the coordination of niobium atoms changes from octahedral to tetrahedral, which is confirmed by Raman spectroscopy data [79][80][81]. Strong covalent bond implies the ability to preserve the structure of oxyanions in the molten state. Thus, the melt can contain not only isolated tetrahedral groups, but also complexes with a stable structure [79]. Raman studies of Li 2 O-Nb 2 O 5 system melts have revealed not only bands corresponding to vibrations of terminal NbO 3− and middle NbO 2− groups (815-870 cm −1 ) but also bands corresponding to symmetric stretching bridge vibrations of Nb-O-Nb bonds (670-690 cm −1 ). Bridge bands are located at the junction of the tetrahedra [79,80]. Probably, octahedra NbO 6 exist in the melt in the pre-crystallization temperature range along with NbO 4 tetrahedra. Earlier paper [82] studied Raman spectra of Nb 2 O 5 , intended for charge synthesis and LN crystal growth. The study has revealed five bands 814, 845, 902, 965, and 995 cm −1 in the region of stretching bridge vibrations of oxygen atoms along the polar axis (800-1100 cm −1 ). This indicates a variety of island structures consisting of fragments of chains of octahedra and tetrahedra in niobium pentoxide.
Boron is a highly chemically active dopant. Its fundamental properties are change in crystallization temperature, viscosity, and surface tension of melts [83]. We [19] have qualitatively revealed an increase in melt viscosity, in the melting temperature by~10 K, and the Curie temperature by~50 K compared to nominally pure CLN crystal during growth of studied LiNbO 3 :B (0.55-1.24 mol%) crystals. The crystal chemistry of boron oxide compounds is extremely diverse: it is determined by the possibility of double hybridization of the boron atom, sp 2 -and sp 3 -, with the formation of [BO 3 ] 3− triangles and [BO 4 ] 5− tetrahedra, respectively [84]. Boron neutral atom radius is 0.88 Å; B 3+ -0.15 Å for B(III), 0.25 Å for B(IV). In the ground state, boron atoms have the 2s 2 2p 1 configuration with one unpaired electron and two vacant p-orbitals. Due to the presence of a free orbital in a small boron atom, boron is one of the strongest acceptors of unshared electron pairs. The use of this orbital in the donor-acceptor interaction allows the coordination number of the boron atom to increase to four. The atoms and atomic groups surrounding the boron atom are located at three corners of the tetrahedron, the fourth corner of which remains free, and a positive electric field of significant intensity is created in this corner. This circumstance explains the pronounced ability of boron compounds to provide stable molecular complexes. Boron triangles and tetrahedra can be present in isolation or polymerize among themselves through a common oxygen atom [84]. The tendency towards formation of bulky polyanions explains the high viscosity of melts [85]. It is well known that alkali metal oxides are modifiers of the boron-oxygen network of the melt and transfer the boron atom to a four-coordinated state, i.e., [BO 4 ] tetrahedra are formed [85]. When boron-oxygen triangles and tetrahedra are combined, the absolute values of specific negative charges for complex anions decrease monotonically in the series of ortho-, pyro-, meta-, and polyborates. This ensures sufficient stability of numerous condensed compounds with monovalent as well as large divalent elements. In [80], temperature changes in lithium metaborate melts were investigated from the Raman spectra. The work has shown the presence of a chain boron-oxygen anion in the melt. All these crystal-chemical laws predetermine the reaction potential of melts containing borates.
The triple system Li 2 O-B 2 O 3 -Nb 2 O 5 has also been studied in [73,86]. The authors analyzed over 30 different molar ratios of the components. Paper [73] has revealed in the system six different lithium borates (LiB 3 O 5 , Li 2 B 4 O 7 , LiBO 2 , Li 6 B 4 O 9 , Li 4 B 2 O 5 and Li 3 BO 3 ), one niobium borate (Nb 3 BO 9 ) and three different lithium niobates (LiNb 3 O 8 , LiNbO 3 and Li 3 NbO 4 ). According to this work, lithium borates exhibit remarkable solubility in LN at temperatures below 1100 • C, and thus are suitable for crystal growth as a flux. In addition, the concentration range of LN crystallization, determined by the spontaneous nucleation method, is rather wide in the Li 2 O-B 2 O 3 -Nb 2 O 5 ternary system.
There is a known HTTSSG method of growing SLN single crystals with a reduced photorefractive effect from a congruent melt in the presence of K 2 O flux [36,87,88]. The absence of potassium in the crystal structure becomes evident from a comparison of the ionic radii of Li + and K + (0.68 Å and 1.38 Å, respectively). Isomorphic substitution of related elements of the alkaline group (lithium with potassium) in the cationic sublattice of the crystal seems unlikely with such a significant difference in ionic radii. Thus, potassium, like boron, does not incorporate in the crystal structure [89][90][91].
We assume that boron-containing polyanions, forming stable covalent bonds in the melt with niobium-containing polyanions, thereby bind the excess of niobium and increase the Li/Nb ratio in the melt. As a result, the grown crystal approaches a stoichiometric composition in terms of the cation sublattice ordering degree, in the same way as when using a K 2 O flux. Such a crystal is characterized by a reduced content of Nb Li defects and lithium vacancies (V Li ) in comparison with a CLN crystal, Table 2. This is confirmed by an increase in Curie temperature (T C ) by~47 K in comparison with a CLN crystal [19]. For a comparable in magnitude change in the T C of a LN crystal doped with metallic impurities (Mg, Zn, etc.), significantly higher impurity concentrations (~2 to 3 wt%) are required.
Paper [92] has shown that melts of inorganic polymers that form chain and ring structures with the inclusion of ions of a solute (B 2 O 3 , Na 2 B 4 O 7 , Li 6 B 4 O 9 , etc.), which tend to form glasses, are good solvents. For example, in [93], the formation of a hightemperature poorly soluble borate Al 5 BO 9 is reported; a fair amount of Al 2 O 3 is removed from the melt. It can be assumed that metal cations in the melt that incorporate into the crystal structure as uncontrolled impurities [87] will also be removed from the melt, as a result of which the LN crystal will be more perfect.
Due to Raman data (Table 1) oxygen octahedra are less deformed in the structure of a LiNbO 3 :B (homogeneous doping) than in LiNbO 3 :B (direct doping) crystal. This obviously happens due to an increase in the ordering of Li + and Nb 5+ cations and vacancies along the polar axis. We have concluded this because widths of 432 and 875 cm −1 bands in a LiNbO 3 :B (direct doping) crystal increase in the Raman spectrum. The bands correspond to deformation and stretching vibrations of oxygen atoms in the bridge Me-O-Me (Me: Li + , Nb 5+ ). It is important to stress that 875 cm −1 band parameters are sensitive towards the magnitude of the dipole moment and, thus, the spontaneous polarization of the crystal [2]. LiNbO 3 :Mg crystals were studied in a paper [32]: the crystals were obtained by direct and homogeneous doping. Raman spectroscopy, optical study of macro-and microstructure have revealed higher ordering of homogeneously doped crystals compared to directly doped crystals. Magnesium incorporates directly to niobium pentoxide structure at a homogeneous doping. It has inter-polyhedra coordination with a uniform distribution in the melt volume. Clusters formed in a melt with different doping methods have different structures and sizes; therefore, melts, other things being equal, should crystallize in different ways [50,79,94]. Crystallization occurs through the attachment of the clusters to a growing crystal.
[BO 4 ] 5− tetrahedra can be both of correct and distorted form. Boron forms three covalent bonds with oxygen atoms, and the fourth bond in the tetrahedron is formed by the donor-acceptor mechanism. For this reason, the lengths of B-O bonds in the tetrahedron vary from 1.462 to 1.512 Å, the spread of angles <O-B-O>-104-115 • [84]. During homogeneous LN crystal doping, the boron-containing reagent is introduced directly into the niobium stripping solution. Only one type of borates (Nb 3 BO 9 ) forms initially. Thus, B-O and Nb-O are less different in the length than at direct doping, when seven different borates form [73]. Obviously, due to this, tetrahedra, and then the octahedra formed in the precrystallization region, will be less distorted during homogeneous doping. Therefore, our Raman data confirmed that oxygen octahedra are less deformed in the structure of the LiNbO 3 :B (homogeneous doping) than in LiNbO 3 :B (direct doping) crystals, Table 1.

Computer Simulation of the B 3+ Ion Localization in the LN Crystal Structure
Calculations of the total energy of the Coulomb interaction of point charges (U, eV) were carried out for the oxygen-octahedral lithium niobate cluster (Li + , Nb 5+ , O 2− ) with the B 3+ ion considered in the sp 2 -hybrid state as part of the [BO 3 ] 3plane triangles. Calculations were performed using the formula: where q 1 and q 2 -charges in fractions of an electron, r 12 -distance between the centers of interacting charges [Å], k-constant expressed by the formula (eV·Å): where e-charge of an electron, ε 0 -dielectric constant.
The system (cluster) consists of two Li + cations, two Nb 5+ cations, one B 3+ cation, and twenty O 2− oxygen anions; it is not electrically neutral. We consider only a cluster consisting of six O 6 oxygen octahedra, "torn out" from a large electrically neutral system, in order to study how the energy of the B 3+ interaction with the surrounding fragment of the LN crystal structure depends on the B 3+ position in tetrahedral voids. Taking a larger fragment of the structure is impractical, since the electrostatic interaction of point charges strongly decreases with distance. Li + and Nb 5+ cations sites correspond in our calculations to the structure of the ferroelectric LN phase: lithium is shifted to the lower oxygen plane, niobium-to the upper oxygen plane [2]. The calculations were based on the structural data of the CLN crystal [95]. The work simulates two processes. In the first case, the energy of the Coulomb interaction of the B 3+ ion with a fragment of the CLN crystal structure (6 O 6 octahedra) is calculated. In this case, we took coordinates of the cell basis with constant lattice parameters (a and c = 5.1489 and 13.8631 Å, respectively), though changing with increasing temperature [2]. In the second case, we took the same the LN structure fragment, but with different lattice parameters a and c. For this, the cell parameters were normalized to a temperature of 297 K [96]. Due to similar results obtained in the first and second cases, only experimental data with constant parameters a and c are discussed below.
A fragment of the structure (six oxygen octahedra) taken for modeling is shown on Figure 13. We considered seven possible locations of the boron cation: at the centers of the tetrahedral faces of the first and second octahedral layers bordering the corresponding octahedra (lithium, niobium, and vacant), and also in the plane of the oxygen triplet separating the octahedral layers. In this case, the calculations did not take into account the replacement of lithium with niobium. The calculation results are shown on Figure 14. According to Figure 14, the maximum energy value for both cases corresponds to the presence of a boron cation in the [BO 3 ] 3group bordering on niobium octahedra. For other possible arrangements of boron, the sum of the energy of the Coulomb interaction is much lower, which can be considered as a theoretically possible location of B 3+ cations in the crystal structure.  We did not take into account the influence of boron in the tetrahedral face on point charges outside the system under consideration, since this influence significantly decreases with distance. It should be noted that the B 3+ cation, embedding in the tetrahedral voids of the crystal structure during the growth process introduces an excess positive  We did not take into account the influence of boron in the tetrahedral face on point charges outside the system under consideration, since this influence significantly decreases with distance. It should be noted that the B 3+ cation, embedding in the tetrahedral voids of the crystal structure during the growth process introduces an excess positive We did not take into account the influence of boron in the tetrahedral face on point charges outside the system under consideration, since this influence significantly decreases with distance. It should be noted that the B 3+ cation, embedding in the tetrahedral voids of the crystal structure during the growth process introduces an excess positive charge into the system. In this case, to achieve the minimum energy of the system, it will be advantageous for boron to occupy positions, first of all, in tetrahedral voids bordering lithium and vacant octahedra, or in the oxygen plane separating the oxygen-octahedral layers, Figure 14. The formation of Nb Li point defects when boron is in the tetrahedral face is unlikely, since in this case it leads to an increase in the total energy due to the high localization of positive charges.
Thus, according to calculations, the presence of the B 3+ ion in the crystal structure can prevent the formation of a deep electron trap Nb Li , at least within the considered system. The results indicate that, in contrast to the metallic K 2 O flux, when the potassium element does not incorporate into the crystal structure, the nonmetallic B 2 O 3 flux has a complex effect on the LN single crystals structure. It structures the melt, changes the microscopic crystallization mechanisms, thereby, aligns and brings K Li and K Nb closer to unity (the distribution coefficients of lithium and niobium in the LiNbO 3 :B crystal). It reduces the formation of Nb Li point defects by at least the number of B 3+ cations incorporated into the tetrahedral voids of the structure, corresponding to the boron concentration in the LiNbO 3 :B crystal (~4 × 10 −4 mol%). As shown above, in the Li 2 O-Nb 2 O 5 -B 2 O 3 system, there are no regions of homogeneity of boron and its compounds with lithium niobate, and CLN should be the only crystallizing phase. However, due to the binding of an excess of niobium in the melt with boron, the fraction of lithium in the structure of the growing crystal increases, which increases the Li/Nb ratio and brings the structure of LiNbO 3 :B crystals closer to that of SLN. The increase in the lithium content is additionally confirmed by the calculation of the Li 2 O concentration in the crystals from the increase in the Curie temperature for CLN and LiNbO 3 :B (0.83 mol%) crystals-1145 and 1189 • C, respectively [97].
Thus, due to Raman spectroscopy, structural units of the cation sublattice are more ordered in LiNbO 3 :B crystals, Table 1. LN crystals doped with certain concentrations of metal cations also have an increased cation sublattice order [2,98]. Metal dopants do incorporate into the crystal structure. In the case of boron, such ordering is impossibleboron incorporates into crystal tetrahedral voids and due to our data (see paragraph 3.5.) changes O-O bonds length. This changes the O 6 oxygen octahedra size and the distribution of Nb 5+ and Li + cations in octahedra. As a result, a more energetically favorable distribution of cations over octahedra is established in the crystal, the arrangement of cations along the polar axis orders. Indeed, the anionic sublattice of the LiNbO 3 :B crystal turns out to be less distorted than of CLN crystal, Table 1. At this, the polarizability of the oxygen-octahedral MeO 6 clusters, which determines the nonlinear optical properties of the crystal, changes. This confidently manifests in the LiNbO 3 :B crystals Raman spectrum. 576 and 630 cm −1 bands correspond to doubly degenerate E(TO) and totally symmetric A 1 (TO) vibrations of oxygen atoms in O 6 oxygen octahedra, respectively. They substantially widen~2 times in LiNbO 3 :B spectrum in comparison with SLN and CLN, Table 1. Changes in O-O bonds length caused by trace amounts of boron localizing in tetrahedral voids can be achieved only by a much higher concentration of metal elements Zn 2+ , Mg 2+ [2,25,26,98]. Bands corresponding to vibrations of oxygen atoms in oxygen octahedra widen smoothly in LiNbO 3 :Zn and LiNbO 3 :Mg Raman spectra and sharply-in LiNbO 3 :B spectrum. This can be explained by the following: tetrahedral voids of the LN crystal structure act as a kind of "buffer", they compensate for the deformation of the oxygen framework of the structure [99]. They compensate for various influences on the anionic sublattice, including the distortion of oxygen octahedra due to the presence of metal cations Zn, Mg, etc., in them. In the case of LiNbO 3 :B crystals, some of the tetrahedra are already filled with B 3+ cations, which negatively affects their "buffer" ability, which explains the sharp increase in the 576 and 630 cm −1 Raman bands widths.  ). A decrease in the 880 cm −1 band width in LiNbO 3 :B spectrum may be associated with an increase in the cation sublattice order of these crystals. Aligning k 0 Li and Nb distribution coefficients during crystal growth explains the approach of the Li/Nb ratio to 1 in LiNbO 3 :B; boron in tetrahedral voids limits the Nb Li structural defect formation.
Thus, the ordering of the cation sublattice of the crystal manifests itself in the Raman spectrum not only in the range 200-300 cm −1 , corresponding to the vibrations of metal cations in oxygen octahedra along the polar axis, but is also indirectly traced by the parameters of the 880 cm −1 band. Due to this data, even traces amounts of B 3+ cations incorporating to LN oxygen tetrahedra noticeably distort the anionic framework. This leads to a noticeable asymmetry of oxygen octahedra (compared with that in CLN crystals) and a change in their polarizability. The deformation and polarizability of oxygen-octahedral MeO 6 clusters, in turn, determine the LN crystal electro-optical properties [1,2]. At this, the cationic sublattice of LiNbO 3 :B crystals, on the contrary, turns out to be more ordered in comparison with the that of the CLN crystal, which is also confirmed by Raman data, Table 1.

Conclusions
Optical microscopy, absorption spectroscopy, Raman spectroscopy, laser conoscopy, optical spectroscopy, optical and atomic force microscopy, photoinduced light scattering, and infrared spectroscopy in the field of stretching vibrations of OH groups and computer simulation were applied to study nominally pure CLN and SLN crystals and a series of LiNbO 3 :B crystals grown from a charge of various genesis, doped with a non-metallic element boron. A quantitative assessment of the strengths of the photovoltaic and diffusion fields has been carried out according to the PILS characteristics for all studied crystals. It was shown that the value of the diffusion field responsible for the diffusion mechanism of charge transfer in the LiNbO 3 :B crystal has an intermediate value between the values CLN and SLN; it depends on the boron concentration in the charge. The PILS indicatrix opening angle also depends on the boron concentration in the charge. The angle in LiNbO 3 :B crystal is smaller than that in SLN crystal. A LiNbO 3 :B band gap corresponds to the value for a SLN crystal. At the same time, LiNbO 3 :B crystals have greater optical uniformity compared to SLN crystals. Positions of OH groups is more ordered in LiNbO 3 :B crystal structure than in CLN crystal. Thus, growing nominally pure LN crystals from a melt structured with non-metal boron makes it possible to control the features of the secondary structure, optical uniformity, the magnitude of photoelectric fields and the width of the band gap in the crystal.
Macro-and micro-structure was studied in LiNbO 3 :B (0.55-1.24 mol%) crystal grown from a different charge. The charge contained different concentrations of boron and was prepared by solid phase synthesis-granulation and by homogeneous doping of Nb 2 O 5 :B precursor. Boron is a chemically active element; it is shown to have a strong effect on the melt structure, making it more homogeneous and viscous and changing the size and structure of clusters in the melt. This leads to a significant change in the secondary structure (defect sublattice) and properties of LiNbO 3 :B crystals, in comparison with nominally pure LN crystals. The melt temperature is slightly increased (1264 • C) compared to nominally pure CLN (1257 • C). Boron concentration should not exceed~0.12 wt% in the melt for the successful growth of optically and structurally uniform LiNbO 3 :B crystals. At the same time, the boron concentration in the grown LiNbO 3 :B crystal is~10 −3 to 10 −5 wt%, which corresponds to trace amounts of impurities. Calculations have revealed that boron can incorporate in trace amounts into the faces of oxygen tetrahedra of the LN crystal structure, bordering on the lithium or vacant oxygen octahedra, or into the oxygen plane separating the oxygen-octahedral layers. Trace amounts of boron in the LiNbO 3 :B structure, apparently obeys a certain mechanism of incorporation into oxygen tetrahedra, like the threshold mechanisms for doping metals. Boron introduces an additional positive charge into the system, thereby preventing the formation of Nb Li point defects. On the other hand, incorporating into oxygen tetrahedra, boron noticeably distorts the anionic framework of the crystal structure, changing the O-O bond lengths and changing the polarizability of the oxygen-octahedral MeO 6 clusters, which determines the nonlinear optical properties of the crystal. Moreover, when the O 6 octahedra are distorted, the structural units of the cation sublattice are simultaneously ordered along the polar axis.
Thus, a new approach was formulated to the preparation of nominally pure LN single crystals. Composition and structure of such crystals approach SLN. Flux in trace amounts has a multi-stage and complex effect on the structural and optical properties of LN single crystals. The results obtained make it possible to expand the concept of "doping". Doping usually means the introduction of a noticeable amount of dopants directly into the crystal structure in order to change the physical properties of the material. Despite the high content of boron in the charge (up to 2.0 mol%), its concentration in the crystal is at the level of trace amounts, i.e., orders of magnitude is lower than the concentration of doping metals. Moreover, the B 3+ ion occupies the tetrahedral voids of the structure. However, boron oxide compounds, as strong complexing agents and solvents, have a noticeable effect on the structure and physical characteristics of the melt, and, consequently, on the structure of the grown crystal. Boron structures the melt in a certain way and decreases the amount of Nb Li defects. Probably, it also reduces the content of uncontrolled impurities. Thus, boron reduces the effect of photorefraction in a LN single crystal, increases its structural and optical uniformity. This approach allows one to obtain nominally pure LN crystals by directed structuring of the melt. The crystals have the same properties as doped crystals, but the ordering of the structural units of the cation sublattice and Nb Li defects approach the SLN crystal. At the same time, LiNbO 3 :B crystal has a significantly lower photorefractive effect than SLN crystals. In this case, the degree of the O 6 octahedra distortion and the photorefractive properties of LiNbO 3 :B crystals can be varied using different methods of the charge doping.  Data Availability Statement: Raw data of this paper will be available from corresponding author, N.V., on a reasonable request. appearance of micro-( Figure A11) and nanodefects ( Figure A12). The crack on Figure A1 developed in accordance with the symmetry of the crystal at angles between the individual branches of~120 degrees. Probably, the second phase on Figure A7a can be in a quasicrystalline state since the spatial regions of its separation contain 5th-order symmetry axes.                     Crystals 2020, 10, x FOR PEER REVIEW 2 of 38 characterized by high values of photoelectric fields and photorefraction effect. The latter can be varied in a very wide range [1,2,8].
Optical damage resistance can be increased in congruent LN (CLN, R = Li/Nb = 0.946) crystals by their doping with non-photorefractive (Me: Zn, Mg, In, etc.) cations [2]. Unlike multiply charged photorefractive cations, they do not change their charge state in the crystal (they are not electron donors) under the action of optical radiation. The influence of such dopants on crystal properties is caused by their ability to change the amount of point defects and linked molecular complexes in the crystal cation sublattice. The molecular complexes in question can be caused by OH groups in the crystal structure [1,2,[9][10][11]. Point defects NbLi are Nb 5+ cations in the Li + sites of a perfect stoichiometric (SLN, R = 1) LN composition. They, along with transition metal impurities (for example, Fe), are deep electron traps and influence photorefractive effect the most [1,2]. Moreover, a LN structure contains a lot of shallow electron traps besides NbLi that influence photorefractive effect [12]. The complexity of LN doping task increases provided that significant concentrations of metal dopants inevitably lead to a disorder in optical and structural uniformity of a single crystal [1,2,[9][10][11]13].
Moreover, LN crystal grown at air always contain 10 16 -10 18 cm -3 protons bonded with oxygen by a hydrogen bond. Hydrogen atoms form such complex defects as VLi-OH, NbLi-OH, etc., [7,14,15]. OH-groups play important role in formation of a secondary defect