Bi3TeBO9: electronic structure, optical properties and photoinduced phenomena

Improvement of the nonlinear optical properties, mainly photoinduced enhancement of effective second harmonic generation was observed in Bi3TeBO9 (BTBO) crystals under irradiation of two external UV laser beams. The effect is caused by an interaction of two coherent photoinducing beams of nitrogen laser with wavelength 371 nm near the Bi3TeBO9 fundamental edge, resulting in formation of additional anisotropy. To explain the observed phenomenon and establish the general features of the electronic structure, the band structure calculations, as well as electronic density of states, elastic properties and principal optical functions were simulated by means of the density functional theory approach. The comparison of the nonlinear optical properties with the conventional α-BIBO crystal shows that BTBO is a promising material for the fabrication of UV laser optical triggers.


Introduction
The search for novel inorganic materials for nonlinear optics is an important direction of modern optoelectronics [1][2][3][4]. Representative applications of such materials are: second harmonic generation (SHG), optical parametric oscillation (OPO), optical parametric amplification (OPA) and self-frequency doubling (SFD). The materials should work in possibly broadest subsystem in connection with principal anionic borate clusters like [BO 3 O 5 ], etc., in order to achieve the highest secondorder nonlinear optical susceptibilities, which are determined by microscopic hyperpolarizabilities.
Recently, Bi 3 TeBO 9 (BTBO) crystal, having very large second-order nonlinear optical susceptibilities, was discovered [6,7]. The preliminary data indicate that its nonlinear optical coefficients are, respectively, higher than in the case of the exceptional a-BiB 3 O 6 (a-BIBO) compound [7].
Moreover, following the general considerations one can expect that the titled compound maybe a promising medium for the laser-operated nonlinear optical effects. For this reason, we will explore the changes of its nonlinear optical properties after illumination near the absorption edge. Additionally, the elastic properties of this material should play here crucial role. In this case, we will perform the bandstructure calculations.
In the present work, we report the synthesis of BTBO powders, determine crystal structure and finally show almost 10%-relative enhancement of the existing second order nonlinear optical susceptibilities under additional two-beam coherent UV laser light treatment with wavelength near the absorption edge of BTBO.

Experimental
Powder synthesis and XRD phase analysis BTBO powders were synthesized by means of conventional solid-state reaction route, similarly to the procedure described by Daub et al. [6]. The molar ratio of used precursor compounds: Bi 2 O 3 (SIAL, 99.95%), TeO 2 (MERCK, 99.99%) and B 2 O 3 (SIAL, 99.999%), was 3:2:1.1, respectively. The 5% excess of boron oxide was used to prevent the non-stoichiometry due to the presence of water and possible evaporation of B 2 O 3 at elevated temperatures. The precursor compounds were mixed and thoroughly ground in an agate mortar till the homogeneous mixture was obtained. It was initially preheated at 450°C for 12 h to evaporate the residuals of adsorbed water. In the following steps, the mixture was ground and calcined three times at 690°C for 12 h each time. After the final grounding, the pure, white powders of BTBO were obtained. Their phase uniformity was confirmed using the XRD analysis. The diffraction pattern of synthesized BTBO powder was collected using the BRUKER D8 Discover, equipped with standard CuK a radiator (k Ka1 = 1.54056 Å , k Ka2 = 1.54439 Å , Siemens KFL CU 2 K, 40 kV voltage and 40 mA current in operating mode). The radiator was equipped with Göbel FGM2 mirror. The Bragg-Brentano diffraction geometry was applied. The diffraction angle 2h B range was set to 10-70°with a step of 0.015°and acquisition time 2 s per step. The measurements were done at 25°C in the temperaturestabilized Anton Paar HTK 1200N chamber. For the data processing, the DIFFRAC.SUITE EVA application was used. The k Ka2 signal component was deconvoluted and removed, the background was subtracted and the data were smoothed using fast Fourier transform (FFT). The phase analysis was performed with support of Crystallography Open Database (COD). For crystallographic calculations, the MAUD ver. 2.32 program was used, allowing to perform the Rietveld refinement and calculate the precise unit cell parameters. The reference CIF file used for this analysis was received from the supplementary data given in Ref. [7]. The obtained BTBO powder diffraction pattern is shown in Fig. 1.

Methods of calculation
Calculations of the electronic, optical and elastic properties of BTBO crystal were performed within the 10   framework of density functional theory (DFT) implemented in CASTEP (Cambridge Serial Total Energy Package) program-package of Biovia Materials Studio 8.0 [8]. The initial BTBO crystal structure was taken from Ref. [7]. We applied the plane-wave (PW) basis set. Eigen function and the eigenvalues were estimated from self-consistent solving of the oneelectron Kohn-Sham equation. Exchange and correlation effects were taken into account by using localdensity approximation (LDA) with Ceperley-Alder-Perdew-Zunger parameterization [9,10] and generalized gradient approximation (GGA) together with Perdew-Burke-Ernzerhof parameterization [11]. The cutoff energy of plane-wave basis was set to be equal to E cut = 830 eV. The Monkhorst-Pack k-points grid sampling was set at 2 9 2 9 2 points for the Brillouin zone used for electronic charge density calculations and geometry optimizations, respectively. Both the cell parameters and atomic positions were optimized (relaxed) using Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm [12], until the force on each ion was smaller than 0.01 eV/Å . The resulting structures were then used as input to the electronic structures and related properties calculations.

Structural properties
The crystal structure of BTBO crystal is presented in Fig. 3. BTBO belongs to hexagonal symmetry (space group P6 3 , no. 173) with Z = 2. The experimental and theoretically calculated principal band structure parameters, using DFT-LDA and DFT-GGA exchange-correlation functionals, together with NCPP, after the optimization of lattice geometry are presented in Table 1. Following Table 1, the calculated lattice constants are in good agreement with the experimental ones. Lattice constants calculated using the LDA functional are a bit smaller than the experimental ones, while in the case of GGA functional, these parameters are a bit higher. It can be explained by the renormalization of principal electronic  parameters within the LDA functional, which leads to underestimation of basic structural parameters.

Electronic properties
For calculations of electronic structure and other properties, the optimized crystal structure was used. Calculated band structure of BTBO crystal using LDA and GGA functionals for exchange-correlation effects description is shown in Fig. 4 In the spectral energy range within -15 to 10 eV, one can see two distinctive groups of peaks corresponding to principal electronic states of BTBO crystal. The first group corresponds to the energy states within the range of -10.6 to 0 eV. The 6s and 6p electrons of bismuth form bands in the range from -10.6 to -7 eV and from -6.5 to 0 eV. Also, a minor peak of 6s electrons is observed near 0 eV. Tellurium electrons form several peaks for 5s electrons, which are situated nearby -10 eV and at -4 eV, while the 5p electrons form peak at -8 and -1 eV, respectively. In the case of boron, the 2s peak at -6 eV is spectrally shifted toward relatively higher energies when  Figure 4 Band energy diagrams for BTBO crystal calculated using LDA and GGA exchange-correlation functionals. compared to the peak formed by tellurium. Wide band originating from -6.3 to -1 eV is formed by 2p electrons of boron. The electrons of oxygen atoms, originating from bands that are situated in the range from -10.3 to 0 eV, correspond to 2p subshell. As shown in Fig. 5, the top of the valence band originates from oxygen 2p states, bismuth 6s and tellurium 5p states. The bottom of the conduction band is formed by p states of oxygen and s states of tellurium. It can be assumed that the optical transitions in absorption spectra are likely formed by oxygen 2p to bismuth 7s or 7p states. Additional important information can be obtained from atomic charge bond population. In contrast to plane-wave basis, a linear combination of atomic orbitals (LCAO) basis set provides a natural way of specific quantities such as atomic charge, bond population and charge transfer. Population analysis in CASTEP is performed using a projection of the PW states onto a localized basis using a technique described by Sanchez-Portal et al. [13]. Population analysis of the resulting projected states is then performed using the Mulliken formalism [14]. This technique is widely used for the analysis of electronic structure calculations performed within a framework of LCAO basis sets. Calculated atomic populations, in units of the proton charge for BTBO crystal are presented in Table 3. One can see that the obtained charges calculated using LDA functional are close to the calculated using GGA.
The shortest bond lengths between constituent atoms in BTBO crystal, calculated using LDA and GGA functional together with overlap populations, are collected in Table 4. As can be seen, the calculated bond lengths and population values are in good agreement between one to another. All bonds of the titled crystal are of covalent type. It is easy to see that B-O bonds possess the highest degree of covalency, while Bi-O-the lowest one. Based on these results, we can conclude that BTBO crystal possesses covalent type of chemical bonding.

Optical properties
The optical properties of solids are described by their dielectric function e = e 1 ? ie 2 , which real part e 1 corresponds to the refractive indices of the material, while its imaginary part e 2 describes the material absorption coefficients, which are directly related to the absorption spectra. The imaginary part of complex dielectric function e 2 , which corresponds to the absorption, can be obtained from the inter-band transition by numerical integration in k-space of elements of dipole matrix operator, between the filled VB states and empty CB levels (Eq. 1): where E-energy of probing optical beam; u-incident photon polarization vector; w c k and w v k -wave functions of the conduction and valence bands in kspace, respectively; X-unit cell volume; e-electron charge; e 0 -dielectric permittivity of vacuum; r-operator of electron position; and x-frequency of incident radiation.
Using Kramers-Kronig relation (Eq. 2) [15] from the imaginary part of dielectric function, its real part can be estimated as follows: Figure 6 depicts the dispersion of real e 1 and imaginary e 2 parts of dielectric functions calculated using GGA functional. The crystals of hexagonal symmetry have two principal directions that correspond to the direction along polar z-axis and direction perpendicular to z. In Fig. 6, we present the calculated dielectric function, which corresponds to photon polarization E || z and E ? z. The dielectric function shows insignificant anisotropy for BTBO crystal. The real part of dielectric function in the region from 7.5 to 18 eV has inversion of sign that can correspond to metallic reflection.

Elastic properties
In terms of crystal anisotropy studies, elastic properties of the crystal can be of special interest. The optimized crystal structure of investigated compound as input for the elastic constant simulations was used. The elastic constants are defined by means of the Taylor expansion of the total system energy E(V, d), with respect to a small lattice strain d of the primitive unit cell with volume V. The energy of a strained system is expressed as follows [16]: where E(V 0 , 0) is the energy of the unstrained system with equilibrium volume V 0 , s 0 is an element in the stress tensor and n i is a factor to take care of Voigt index.
For the hexagonal crystal, its elastic properties can be described by asymmetric 6 9 6 matrix of elastic coefficients. In case of hexagonal symmetry, we get the following relation between elastic constants: C 11 = C 22 ; C 13 = C 23 ; C 44 = C 55 ; C 66 = 1/2(C 11 -C 12 ); it is known that for practical application it is important to obtain a stable crystal structure. Knowing of the elastic constants and Young's modulus is important for technological aims, in the manufacturing process and crystal preparation. For crystal structure stability of hexagonal symmetry, the obtained elastic constants should satisfy the requirements of mechanical stability criteria [17], which are: C 33 [ 0; C 44 [ 0; C 12 [ 0; C 11 [ |C 12 |; (C 11 ? 2C 12 )C 33 [ 2C 13 2 ).
Following Table 5, the calculated elastic constants completely satisfy the mechanical stability criteria. Unfortunately, according to already published data concerning the BTBO crystal, the experimental results of its mechanical properties are not yet known. The values of C 11 and C 33 differ one from another, indicating the anisotropy of elastic properties. The universal anisotropy index for BTBO crystal is equal to A U = 0.412/0.254 from the calculation using GGA/ Re( ), Im( ) Figure 6 Calculated real Re(e) and imaginary Im(e) parts of the dielectric functions dispersion for the BTBO crystal using GGA functional. LDA functionals. Calculated elastic compliance matrix S ij for BTBO crystal is presented in Table 6. The Young's modulus components E x , E y and E z can be estimated from elastic compliance S ij using following relations: The axial compressibilities k a , k b and k c for hexagonal BTBO crystal are also calculated in the following way [18]: Both Young's modulus E and compressibility k, calculated using GGA/LDA functionals, are collected in Table 7.
Additional useful method to show mechanical anisotropy can be given by visualization of a threedimensional dependence of the bulk modulus and Young's modulus as a function of direction. For hexagonal crystal class, the directional dependencies of bulk modulus and Young's modulus are given by the following relations [19]: where B is the bulk modulus and E is the value of Young's modulus, in the direction determined by l 1 ; l 2 and l 3 -the direction cosines of [uvw] directions and S ij are the elastic compliance constants, which are inverse to the matrix of elastic constants C ij calculated for BTBO crystal. The surfaces of directional dependence of bulk modulus B and Young's modulus E for BTBO crystal, plotted in spherical coordinates, are depicted in Fig. 7. The surfaces of directional dependence of bulk and Young's modulus (Fig. 7) show that the BTBO crystal possesses low anisotropy. Both surfaces have the similar shape. In order to receive additional information and a better understanding of bulk modulus and Young's modulus relation along different directions, the projections of the moduli at ab and ac crystal plane are plotted (see Fig. 8).
It is easy to see that at ab plane both bulk modulus and Young's modulus are of a spherical shape, which means in this plane the BTBO crystal is almost isotropic. In the case of ac plane, some anisotropy of mechanical properties is shown.

Effective SHG measurements
Following the performed band-structure calculations and the elastic features, one can expect that the titled crystal may be operated by an external laser light of wavelength near the fundamental absorption edge. For this reason, we used two coherent nitrogen 371 nm 9 ns laser pulses with frequency repetition about 10 Hz in order to study the possible changes of the effective second-order nonlinear optical susceptibilities. As a probing beam, an Er:glass 20 ns laser (at k = 1540 nm) was used. The detection was performed using interferometric filter at 770 nm and the photomultiplier connected to the oscilloscope. In order to monitor the possible photoluminescence near 371 nm laser treatment, the additional control of the spectral lines and their time dependencies were performed in the range of 330-905 nm by spectrometer. There was any wider spectral noticeable photoluminescence in this spectral range. Additionally, the time shift between the excitation and the induced nonlinear signals did not exceed few nanoseconds, while for the fluorescence it should be of an order higher. The use of two coherent beams interacting at different angles should stimulate the occurrence of additional birefringence [20], which also changes the nonlinear optical susceptibilities. Moreover, during the photoinduced changes in the common borate crystals, the contributions of anharmonic phonons begin to play an additional role [21,22]. They may be responsible for the photoinduced second-order nonlinear optical susceptibility changes versus the cationic content, as shown in Fig. 9.
It is crucial that, contrary to the well-studied a-BIBO crystal, for the titled BTBO there exists a huge opportunity to operate the output SHG using external photoinducing UV coherent nanosecond laser pulses. The arrow indicates an excitation of the prevailing number of anharmonic phonons by the nitrogen pulsed laser beams. Because they are anharmonic, their symmetry is the same as for the second-order nonlinear optical susceptibility. As a consequence, it is manifested in Fig. 9 and the optical functions are changed significantly. The effect disappears immediately after the switching off the photoinducing beam treatment.

Conclusions
An opportunity to influence the second-order nonlinear optical susceptibility of BTBO powder using two coherent beams treatment of 371 nm laser near the absorption edge was established. The performed band-structure calculations and the related elastic properties show that BTBO possesses relatively improved elastic features with respect to other borates. This may favor laser operation with the output second-order nonlinear optical susceptibilities due to changes in microscopical hyperpolarizability. This phenomenon may be a consequence of both the elastooptical photoinduced anisotropy and possible contribution of the phonon anharmonicities. The comparison with the conventional a-BIBO crystal shows that BTBO is a promising material for the fabrication of UV laser optical triggers. Figure 9 Dependence second-order nonlinear optical susceptibility of BTBO (red plot) at fundamental wavelength 1540 nm versus the power density of the 371-nm pulsed 9 ns laser. The latter is incident under the optimal angle between photoinducing beams equal to about 20 degrees. The arrow corresponds to the critical photoinducing power density at which it gets maximal changes of the second-order susceptibility. The measurement was taken with comparison to a-BIBO (black plot).