Pseudosymmetric Features and Nonlinear Optical Properties of Potassium Titanyl Phosphate Crystals

A number of publications containing structural data, characteristics of nonlinear optical properties of pure and doped crystals of potassium titanyl phosphate (KTP) family have been reviewed to analyze the structural and symmetry conditionality of nonlinear optical properties of these crystals. The pseudosymmetric features of KTP-type crystals with respect to inversion are investigated. Specifically, pseudoinversion distribution maps are calculated; pseudoinversion extrema and coordinates of pseudoinversion centres are found; and the distributions of pure and doped KTP-type structures and their individual atomic sublattices over the degree of pseudoinversion are analyzed. A correlation between the characteristics of nonlinear optical properties of a number of crystals belonging to the KTP family and the degree of pseudoinversion of their atomic structures is demonstrated.


Introduction
Study of the relationship of structural and symmetric features of crystals with their physical properties is an urgent problem of condensed-matter physics.Point symmetry determines the set of possible physical properties of crystals, primarily, in correspondence with the Neumann principle.The symmetric features of atomic structures of crystals can be characterized more completely taking into account the phenomenon of pseudosymmetry, which makes it possible to establish finer relationships of the structure-property type.Fedorov pseudosymmetry of crystals [1] is the phenomenon of invariance of a considerable part of the crystal atomic structure (part of electron density and (or) subsystem of atomic nuclei) with respect to some group of symmetry operations compatible with the lattice (with respect to some supergroup of the symmetry space group of crystal).
The pseudosymmetry of a specific structure can quantitatively be characterized by the degree of invariance (degree of pseudosymmetry) of its total electron density (r) with respect to some isometric operation [1,2]: Integration in Equation ( 1) is performed over the volume V of crystal unit cell.If is not a symmetry operation for the function (r), the degree of pseu- < 1; however, if (r) is symmetric with respect to the operation , reaches a maximally possible value: unity.
The second-order susceptibility of crystal determines the intensity of generation of the second optical harmonic and is a structure-and symmetry-sensitive property of crystal.For centrosymmetric crystals, the second-order susceptibility should be zero.One might suggest that reduction of symmetry will lead to some dependence of the second-order susceptibility of crystal on the degree of invariance of crystal structure with respect to inversion.leads to a dramatic decrease in the output second-harmonic intensity to zero.The intensity of second-harmonic generation (SHG) for KTi 1-x Zr x OPO 4 crystals reaches a maximum at x = 0.28, where it is more than doubled in comparison with the KTiOPO 4 sample [41].The SHG intensity increases with an increase in the zirconium content in RbTi 1-x Zr x OPO 4 crystals; at x = 0.034, it rises by 40% [50].The SHG intensity increases by approximately 35% in comparison with pure KTiOPO 4 samples after replacement of 6% titanium atoms with hafnium [47].RbTi 1-x Ta x OPO 4 and RbTi 1-x Nb x OPO 4 crystals were investigated in [52], as well as Yb-doped RbTi 1-x Ta x OPO 4 crystals and RbTi

Structure of KTP Crystals with Respect to Inversion
The complete characteristic of pseudosymmetry of any crystal under study with respect to inversion is a threedimensional distribution map of the degree of structural invariance (electron density) with respect to this operation (hereinafter, pseudoinversion), calculated for different positions of inversion points within their unit cell.These maps were obtained with scanning steps over the unit-cell axes a, b, and c chosen to be 0.05 of the cor responding unit-cell parameters.For 118 crystals (Figure 1) with known structure, we calculated three-dimensional pseudoinversion maps using Equation (1).The calculations were performed using the computer program and technique described in [54].Within this approach the electron density function is expanded in a Fourier series in structural amplitudes ( [1], see Formulas ( 5) and ( 6)). Figure 2 presents typical examples of cross sections of three-dimensional distribution maps of the degree of pseudoinversion for the structures of KTiOPO 4 , KSnOPO 4 , KTiOAsO 4 , Cs 0.625 K 0.375 TiOAsO 4 crystals (cuts by the plane z = 0.25).For the structures presented in Figure 2, the origin of coordinates is chosen on the two fold screw axis, and the coordinates of pseudoinversion peaks on the x and y axes are 0.25.We chose cuts by the plane z = 0.25 in Figure 2 because the z coordinate of the pseudoinversion peaks for the structures of the aforementioned crystals is also 0.25.
This situation is typical of most structures under study; in KTP crystals is accompanied by a phase transition to the centrosymmetric space group Pnan.Indeed, having    A value of 1100 was indicated in [3], with reference to [5], where a value of 2400 was reported.** A value of 140 was indicated in [25], with reference to [12], where a value of 40 was reported, and a value of 6 was indicated in [3], with reference to [5], where the corresponding value was found to be 60.*** Values of second-harmonic generation intensity for KTiOPO4 crystal were reported in [14] with reference to [6], where the corresponding values were given for quartz crystal.The values of [14] correspond to those of [6], when divided by I/I SiO2 value for KTiOPO4 crystal (also taken from [6]).added inversion to the set of symmetry operations of the space group Pna2 1 , which describes the structure of KTP crystals at room temperature, we obtain the group Pnan, where the inversion centre with respect to the twofold screw axis has coordinates (0.25, 0.25).Thus, in the polar phase of KTP structures, the pseudoinversion peaks are located specifically at inversion centres of these crystals in their high-symmetry nonpolar phase.With allowance for this circumstance, we will characterize the pseudosymmetry of the electron density in each crystal with a known structure by the maximum value of pseudoinversion  max in the three-dimensional map, and the point with the coordinates corresponding to the found  max values will be referred to as pseudoinversion centres.
Since the origin of coordinates is arbitrarily chosen in X-ray diffraction analysis, the coordinates of the pseudoinversion centres may differ from 0.25.In the Pna2 1 group [55], the origin of coordinates on the z axis can be chosen at any point, while in the directions of the x and y axes it may lie either on the twofold axis or at the intersection of mirror planes; therefore, the x and y coordinates of pseudoinversion centres can be either (0.25, 0.25) or (0, 0).
To refine the coordinates of pseudoinversion centres and  max values, we additionally calculated the distribution of the degree of pseudoinversion with a relative scanning step of 0.025 over the unit cell axes.Fixed refined x and y coordinates of pseudoinversion centres were used for repeated calculation of pseudoinversion distribution along the z axis with a relative scanning step of 0.001.3(b), which amounts to 31%, lies in the range of pseudoinversion values of 0.4 -0.5, which is followed by a sharp falloff.Therefore, the fraction of pseudo-centrosymmetric structures among doped KTP crystals is very small.The mean value <  max > for doped crystals is 0.490.Thus, doped KTP crystals are "less symmetric" with respect to inversion than pure compositions.
For the crystals listed in Table 2, along with the calculations of the pseudoinversion of their structures as a whole, pseudoinversion extrema for sublattices of individual types of atoms ( max (sublattices)) were also cal-culated.Pseudoinversion was calculated for the pure sublattices of all 118 crystals in Table 2; the distribution histogram of the corresponding extrema is shown in Figure 3(c).For 89 doped crystals in Table 2, the results of similar calculations for M-and M'-type sublattices containing doped atoms are presented as histograms in Figure 3(d).The histograms in Figure 3(c) indicate that the sublattices of Х, О, and F atoms are most pseudocentrosymmetric, sublattices of M-type atoms are least pseudo-centrosymmetric, and the pseudoinversion of the In the presence of impurities, the general view of the histogram for the M' sublattice (Figure 3(d)) barely differs from that in Figure 3(c); its characteristic maximum shifts to higher pseudoinversion values and the mean value <  max (M') > becomes 0.700.
The pseudoinversion histogram for the M-type sublettice changes more radically: the pronounced peak in the range of 0.3 -0.4 in Figure 3(c) disappears in Figure 3(d), and the distribution becomes more uniform in a wider pseudoinversion range; the fraction of crystals with ultimately acentric M sublattices increases.The mean pseudo inversion <  max (M) > becomes 0.384; i.e., it barely changes in comparison with <  max (M) > for M sublattices without impurities.Thus, the analysis of the pseudo inversion of individual sublattices suggests that the reductions of pseudoinversion of structures as a whole at a transition to doped KTP compositions, which is noted in Table 2 and Figure 3(b), is related to a great extent to the higher sensitivity of the pseudoinversion of M-type sublattice to the presence of doped atoms.Note that pseudosymmetry was previously studied [58] by the atomic displacement method [1] for 11 KTP-type structures.In particular, it was established that the potassium sublattice is less centrosymmetric in comparison with the TiO 6 -PO 4 subsystem, and its pseudosymmetry relative to inversion is more sensitive to introduction of impurities.

Comparison of the Nonlinear Optical Characteristics of KTP Crystals and the Pseudoinversion of Their Structures
A model was proposed in [8], according to which the second-order susceptibility of crystals is related to the symmetry of KTP-type structures and their pseudoinversion as follows: As can be seen in Table 1, the SHG data with respect to the reference sample (powder of pure SiO 2 crystal) differ by an order of magnitude in different studies for KTP [3,[5][6][7], and KTA [3,6,7] crystals.Based on this fact, we illustrated Equation ( 2) by selecting a group of com-   positions (from the aforementioned set of crystals) for which experimental SHG data were obtained either by the powder method [26], or directly with respect to a powder of pure KTiOPO 4 crystal, or the data can be recalculated with respect to it based on a specific publication.In addition, since most sources yield data on the ratio of second-harmonic intensities for the studied and reference samples (I/I KTP = I 2 /I 2 (KTP)), they were additionally recalculated into estimated values of the relative effective second-order susceptibility (it will be denoted as  2 / 2 (KTP)).In the first approximation, one can assume that where I 2 is the second harmonic intensity and I  is the primary radiation intensity.Therefore, the desired ratio  2 / 2 (KTP) was estimated to be T P coordinates (see [8] for the approximation model).data of Godfrey et al. [58], who established that the KSnOPO 4 structure can be partially described (in good approximation) by the Pnan group; exact description is obtained within the Pna2 1 group.

   
coordinates.The linear approximation of the dependence presented in Figure 4 within the model described in [8], is characterized by a correlation coefficient of 0.91, and the confidence interval boundaries are (0.76, 0.97) at a confidence probability of 0.95.

Conclusions
To date, despite the numerous publications on the structure and properties of KTP crystals, the question of the structural conditionality of the behavior of their nonlinear optical properties has not been completely clarified.
Equation ( 2) can be more pronounced within the concentration series of samples of the same qualitative composition.For example, the SHG intensity decreases with an increase in the tin fraction in the KTi 1-x Sn x OPO 4 series, and the calculation of pseudoinversion for a series of known structures of this composition indicates a monotonic increase in the latter (Figure 5).
In this paper, we reported the results of studying the pseudosymmetric features of known structures of KTP crystals with respect to inversion and tried to analyze the entire set of known nonlinear optical parameters of these crystals in view of the obtained pseudosymmetric characteristics.In particular, it was shown that doped structures of KTP crystals have on average a lower degree of pseudoinversion than "pure" compositions; in some cases this feature adequately correlates with the increase in the The boundary-composition crystal KSnOPO 4 has almost zero SHG intensity and the largest (in the KTi 1-x Sn x OPO 4 series) pseudoinversion: 0.886 (Table 1, no. 7; Table 2, no.4).This fact is in agreement with the relative intensity of the second optical harmonic.This correlation may manifest itself within the concentration series samples of the same qualitative composition.
We believe that, in order to establish the fundamental correlations between the structural and symmetric features of crystals (in particular, those belonging to the KTP family) and their nonlinear optical properties, for example, having the degree of pseudoinversion as a symmetric characteristic, it is necessary to primarily calculate this characteristic for the entire structure.This thesis is justified by the fundamental principles of symmetry in physical crystallography.The Neumann princeple, which sets a relationship between the symmetry of a medium (crystal) and the set of physical properties that are forbidden or allowed in this medium, deals with specifically the symmetry of the medium as a whole rather than the symmetry of its individual structural fragments within the unit cell.This approach was applied both in [8] and in this study.However, this does not depreciate the validity of the analysis of the characteristics of sublattices of individual types of atoms.Due to this analysis one can find sublattices with pseudosymmetric characteristics exhibiting a more significant sensitivity, for example, at a transition to doped compositions, and therefore, can determine to greater extent the behavior of the pseudosymmetric characteristics of crystal structures, as whole and physical properties of crystals.

Figure 1 .
Figure 1.Groups of KTP crystals under consideration.

Figure 3 .
Figure 3. Distribution histogram of the degree of pseudoinversion for (a) pure KTP crystals (29 structures), (b) doped KTP crystals (89 structures), (c) pure atomic sublattices of individual types for 118 KTP crystals fromTable 2, and (d) doped ublattices of individual types of atoms for 89 KTP crystals from Table 2. s

Figure 4 .
Figure 4. Correlation between the relative effective secondorder susceptibility 2/2(KTP) for a number of KTP crystals and their pseudoinversion in the

Table 2
The distribution of pure KTP crystals over pseudoinversion is fairly uniform.As is indicated in contains the maximum pseudoinversion values  max and coordinates of pseudoinversion centres z( max ) for a number of KTP structures.Figures 3(a) and 3(b) show the distribution histograms for the degree of pseudoinversion  max for pure and doped KTP crystals.

Table 2 ,
the mean value <  max > for them is 0.606.The situation for doped crystals is different: the pronounced maximum in the histogram in Figure

Table 1 ,
which are given in correspondence with the references to original studies.