# Optical second harmonic generation in Yttrium Aluminum Borate single crystals (theoretical simulation and experiment)

- Ali H Reshak
^{1}Email author, - S Auluck
^{2}, - A Majchrowski
^{3}and - IV Kityk
^{4}

**Received: **25 October 2007

**Accepted: **17 March 2008

**Published: **17 March 2008

## Abstract

Experimental measurements of the second order susceptibilities for the second harmonic generation are reported for YAl_{3}(BO_{3})_{4} (YAB) single crystals for the two principal tensor components xyz and yyy. First principle's calculation of the linear and nonlinear optical susceptibilities for Yttrium Aluminum Borate YAl_{3}(BO_{3})_{4} (YAB) crystal have been carried out within a framework of the full-potential linear augmented plane wave (FP-LAPW) method. Our calculations show a large anisotropy of the linear and nonlinear optical susceptibilities. The observed dependences of the second order susceptibilities for the static frequency limit and for the frequency may be a consequence of different contribution of electron-phonon interactions. The imaginary parts of the second order SHG susceptibility ${\chi}_{123}^{(2)}\left(\omega \right)$, ${\chi}_{112}^{(2)}\left(\omega \right)$, ${\chi}_{222}^{(2)}\left(\omega \right)$, and ${\chi}_{213}^{(2)}\left(\omega \right)$ are evaluated. We find that the 2*ω* inter-band and intra-band contributions to the real and imaginary parts of ${\chi}_{ijk}^{(2)}\left(\omega \right)$ show opposite signs. The calculated second order susceptibilities are in reasonable good agreement with the experimental measurements.

**PACS Codes:** 71.15. Mb; 71.15.-m

## I. Introduction

Yttrium Aluminium Borate YAl_{3}(BO_{3})_{4} (YAB) belongs to a family of double borates which crystallize in the trigonal structure of the mineral huntite CaMg_{3}(CO_{3})_{4} and belong to the space group R32 [1]. The general formula of these compounds is RX_{3}(BO_{3})_{4}, where R = Y^{3+}, Gd^{3+} or any other lanthanide, and X = Al^{3+}, Sc^{3+}, Ga^{3+}, Cr^{3+}, Fe^{3+} [2]. YAB is a non centro-symmetric crystal and as early as in 1974 it was reported to be a very effective second-harmonic generating material [3]. Furthermore, owing to its good chemical stability and possibility of substituting Y^{3+} ions with other lanthanide ions, namely Nd^{3+}, Yb^{3+}, Dy^{3+} and Er^{3+} [4] it is a promising material for laser applications. The nonlinear optical properties of this material in connection with lasing properties led to the construction of numerous systems generating red, green and blue light by self-frequency doubling effect [5]. YAB can be obtained as nano crystallite powders by simple technological approaches [6, 7]. They also possess relatively large two-photon absorption [6], which makes them promising third order optical materials. At the same time they are good matrices for different rare earth ions [8–11]. The existing data is restrained by consideration of the local crystalline fields and the influence of the rare earth ions [12–15]. Intrinsic defects also may play an important role in determining the optical susceptibilities [16, 17].

We feel that a reliable band structure will be of immense help in understanding the linear and nonlinear optical properties and show directions for technologists to obtain crystalline materials with desired optical properties. A reliable band structure calculation can help in determining the role of inter-band dipole matrix elements on the optical properties. It can give information on the dispersion of the bands in **k**-space and origin of the bands which are directly connected with the optical hyperpolarizabilities and susceptibilities [18]. In the present work we report first principle's calculation of the linear and nonlinear optical susceptibilities for YAB using the state-of-the-art full potential linear augmented plane wave method [19] which has proven to be one of the most accurate methods [20, 21] for the computation of the electronic structure of solids within density functional theory (DFT) [22]. One specific feature of the borate crystals is the coexistence of the strong covalent and ionic chemical bonds, which provide relatively flat **k**-dispersion of the bands [23]. Moreover there exists substantial anisotropy of the chemical bonds which substantially restrains the application of the pseudopotential method, even norm-conserving one [24]. The aim of this paper is to understand the origin of birefringence and the high *χ*^{2}(*ω*), using first principle's calculations.

In the Section 2 we present the computational and experimental details. Section 3 gives the results of the calculations and the measurements. The observed discrepancies are discussed following the band energy approach.

## II. Computational and experimental details

*α*=

*γ*= 90° and

*β*= 120°. YAB melts incongruently at 1280°C and decomposes into YBO

_{3}and AlBO

_{3}. Therefore, it cannot be crystallized from stoichiometric melts. The application of high temperature solution growth (HTSG) method allows lowering the temperature of YAB crystallization below the temperature of the peritectic transformation.

Self-consistent calculations of the electronic structure and optical properties based on the scalar relativistic full-potential linearized augmented plane wave method were carried out using the WIEN2K package [19]. This is a very accurate and efficient scheme to solve the Kohn-Sham equation of density functional theory (DFT) [22]. This is an implementation of the DFT with different possible approximations for the exchange-correlation (XC) potential. The XC is treated within the local density approximation (LDA) [25] and scalar relativistic equations are used to obtain self-consistency. The Kohn-Sham equations are solved using a basis of linear APW's. In the interstitial region the potential and the charge density are represented by Fourier series. In order to achieve energy convergence, the wave functions in the interstitial region were expanded in plane waves with a cut-off *K*_{max} = 9/*R*_{MT}, where *R*_{MT} denotes the smallest atomic sphere radius and *K*_{max} gives the magnitude of the largest *K* vector in the plane wave expansion. The *R*_{MT} are taken to be 2.14, 1.81, 1.28, and 1.28 atomic units (a.u.) for Y, Al, B and O respectively. The valence wave functions inside the spheres are expanded up to *l*_{max} = 10 while the charge density was Fourier expanded up to *G*_{max} = 14.

Self-consistency is obtained using 200 $\stackrel{\rightharpoonup}{k}$ points in the irreducible Brillouin zone (IBZ). We calculated the frequency dependent linear optical properties using 500 $\stackrel{\rightharpoonup}{k}$ points and nonlinear optical properties using 1400 $\stackrel{\rightharpoonup}{k}$ points in the IBZ. The self-consistent calculations are assumed to be converged when the total energy of the system is stable within 10^{-5} Ry.

## III. Results and Discussion

### A. First order optical susceptibilities and birefringence

The above expressions are written in atomic units with e^{2} = 1/m = 2 and *ħ* = 1. where *ω* is the photon energy and ${\text{P}}_{nn}^{X}\left(k\right)$ is the x component of the dipolar matrix elements between initial |*nk*⟩ and final |*n*'*k*⟩ states with their eigen-values E_{n} (*k*) and E_{n'}(*k*), respectively. *ω*_{nn'}(*k*) is the band energy difference *ω*_{nn'}(*k*) = E_{n}(*k*) - E_{n'}(*k*) and S_{k} is a constant energy surface S_{k =} {*k*; *ω*_{nn'}(*k*) = *ω*}.

*n*=

*n*

_{ e }-

*n*

_{0}, where

*n*

_{ e }is the index of refraction for an electric field oriented along the

**c**-axis and

*n*

_{0}is the index of refraction for an electric field perpendicular to the

**c**-axis. Figure 6, shows the birefringence Δ

*n*(

*ω*) for this single crystal. The birefringence is important only in the non-absorbing region, which is below the energy gap. In the absorption region, the absorption will make it difficult for these compounds to be used as nonlinear crystals in optic parametric oscillators or frequency doublers and triplers. We note that the spectral feature of Δ

*n*(

*ω*) shows strong oscillations around zero in the energy range up to 12.5 eV. Thereafter it drops to zero. We find that the calculated birefringence at zero energy is 0.025 in excellent agreement with our own measurement of 0.02. It is known that for the borates, the contribution of the electron-phonon interaction to the dielectric dispersion may be neglected for the SHG effects [29] contrary to the linear electro-optics Pockels effect. Comparing these dependences with the anisotropy for other borates [23, 24] one can conclude that the anisotropy caused by the chemical bonds is smaller than in the other borates [24].

### B. Second order susceptibilities

where *n* ≠ *m* ≠ *l*. Here *n* denotes the valence states, *m* the conduction states and *l* denotes all states (*l* ≠ *m*, *n*). There are two kinds of transitions which take place one of them *vcc'*, involving one valence band (*v*) and two conduction bands (*c* and *c'*), and the second transition *vv'c*, involving two valence bands (*v* and *v'*) and one conduction band (*c*). The symbols are defined as ${\Delta}_{nm}^{i}\left(\overrightarrow{k}\right)={\vartheta}_{nn}^{i}\left(\overrightarrow{k}\right)-{\vartheta}_{mm}^{i}\left(\overrightarrow{k}\right)$ with ${\overrightarrow{\vartheta}}_{nm}^{i}$ being the *i* component of the electron velocity given as ${\vartheta}_{nm}^{i}\left(\overrightarrow{k}\right)=i{\omega}_{nm}\left(\overrightarrow{k}\right){r}_{nm}^{i}\left(\overrightarrow{k}\right)$ and $\left\{{r}_{nm}^{i}\left(\overrightarrow{k}\right){r}_{ml}^{j}\left(\overrightarrow{k}\right)\right\}=\frac{1}{2}\left({r}_{nm}^{i}\left(\overrightarrow{k}\right){r}_{ml}^{j}\left(\overrightarrow{k}\right)+{r}_{nm}^{j}\left(\overrightarrow{k}\right){r}_{ml}^{i}\left(\overrightarrow{k}\right)\right)$. The position matrix elements between states *n* and *m*, ${r}_{nm}^{i}\left(\overrightarrow{k}\right)$, are calculated from the momentum matrix element ${P}_{nm}^{i}$ using the relation [35]: ${r}_{nm}^{i}\left(\overrightarrow{k}\right)=\frac{{P}_{nm}^{i}\left(\overrightarrow{k}\right)}{im{\omega}_{nm}\left(\overrightarrow{k}\right)}$, with the energy difference between the states *n* and *m* given by *ħ ω*_{
nm
}= *ħ*(*ω*_{
n
}- *ω*_{
m
}). *f*_{
nm
}= *f*_{
n
}- *f*_{
m
}is the difference of the Fermi distribution functions. *i*, *j* and *k* correspond to cartesian indices.

It has been demonstrated by Aspnes [36] that only one virtual-electron transitions (transitions between one valence band state and two conduction band states) give a significant contribution to the second-order tensor. Hence we ignore the virtual-hole contribution (transitions between two valence band states and one conduction band state) because it was found to be negative and more than an order of magnitude smaller than the virtual-electron contribution for these compounds. For simplicity we denote ${\chi}_{ijk}^{(2)}\left(-2\omega ;\omega ;\omega \right)$ by ${\chi}_{ijk}^{(2)}\left(\omega \right)$.

We have measured the second order susceptibilities of YAB single crystal using Nd-YAG laser at the fundamental wavelength 1064 nm. Since the investigated crystals belong to the point group R32 there are only five independent components of the SHG tensor, namely, the 123, 112, 222, 213 and 312 components (1, 2, and 3 refer to the x, y and z axes, respectively) [30]. These are ${\chi}_{123}^{(2)}\left(\omega \right)$, ${\chi}_{112}^{(2)}\left(\omega \right)$, ${\chi}_{222}^{(2)}\left(\omega \right)$, ${\chi}_{213}^{(2)}\left(\omega \right)$ and ${\chi}_{312}^{(2)}\left(\omega \right)$. Here ${\chi}_{ijk}^{(2)}\left(\omega \right)$ is the complex second-order nonlinear optical susceptibility tensor ${\chi}_{ijk}^{(2)}\left(-2\omega ;\omega ;\omega \right)$. The subscripts *i, j*, and *k* are Cartesian indices.

*total Re*${\chi}_{ijk}^{(2)}\left(0\right)$ value compared to the other components (Table 1). A definite enhancement in the anisotropy on going from linear optical properties to the nonlinear optical properties is evident (Figures 7 and 8). It is well known that nonlinear optical susceptibilities are more sensitive to small changes in the band structure than the linear optical ones. Hence any anisotropy in the linear optical properties is enhanced more significantly in the nonlinear spectra.

Calculated total, intra-band and inter-band contributions of *Re* ${\chi}_{ijk}^{(2)}\left(0\right)$ in units of 10^{-7} esu, along with the measured ${\chi}_{ijk}^{(2)}\left(0\right)$ in units of pm/V.

Component | 123 | 112 | 222 | 213 |
---|---|---|---|---|

| -0.0011 | -0.009 | 0.01 | 0.0002 |

| -0.007 | 0.035 | -0.038 | 0.006 |

| 0.008 | -0.04 | 0.04 | -0.0085 |

| -0.5 | -0.7 | 0.8 | 0.5 |

| 0.12 | 1.03 |

*ω*inter-band and intra-band contributions to Im ${\chi}_{222}^{(2)}\left(\omega \right)$. We note the opposite signs of the two contributions throughout the frequency range. We have calculated the total complex susceptibility for ${\chi}_{123}^{(2)}\left(\omega \right)$,... and ${\chi}_{213}^{(2)}\left(\omega \right)$. The real part of the dominant component is shown in Figure 10. The zero-frequency limit of all components is listed in Table 1.

From above we can see the total second order susceptibility determining SHG is zero below half the band gap. The 2*ω* terms start contributing at energies ~1/2*E*_{
g
}and the *ω* terms for energy values above *E*_{
g
}. In the low energy regime (≤ 5 eV) the SHG optical spectra is dominated by the 2*ω* contributions. Beyond 5.7 eV (values of the fundamental energy gaps) the major contribution comes from the *ω* term.

One could expect that the structures in Im ${\chi}_{ijk}^{(2)}\left(\omega \right)$ could be understood from the features of *ε*_{2}(ω)_{.} Unlike the linear optical spectra, the features in the SHG susceptibility are very difficult to identify from the band structure because of the presence of 2*ω* and *ω* terms. But we use of the linear optical spectra to identify the different resonance leading to various features in the SHG spectra. The first spectral band in Im ${\chi}_{123}^{(2)}\left(\omega \right)$ between 0.0–5.0 eV is mainly originated from 2*ω* resonance and arises from the first structure in *ε*_{2}(ω). The second band between 5.0–7.0 eV is associated with interference between the *ω* resonance and 2*ω* resonance and is associated with high structure in *ε*_{2}(ω). The last structure from 7.0–8.0 eV is mainly due to *ω* resonance and is associated with the tail in *ε*_{2}(ω).

^{-7}esu, respectively. To evaluate the performed calculations we have done the measurements of the absolute value of ${\chi}_{123}^{(2)}\left(\omega \right)$ and ${\chi}_{222}^{(2)}\left(\omega \right)$ for the YAB single crystals for the Nd-YAG laser wavelength 1064 nm and we have revealed the corresponding values equal to about (0.042 and 0.061)×10

^{-7}esu, respectively confirming sufficiently good agreement. The calculated second order susceptibilities show substantially good agreement with the measured one.

## IV. Conclusion

We have performed experimental measurements of the second order susceptibilities for the second harmonic generation for the YAl_{3}(BO_{3})_{4} (YAB) single crystals for the two principal tensor components xyz and yyy. We have reported a first principle's calculation of the linear and non-linear optical susceptibilities using the FP-LAPW method within a framework of DFT. Our calculations show that YAB possesses a direct energy band gap of about 5.1 eV located at Γ point of the Brillouin zone. This is smaller than the experimental value of 5.7 eV. The calculated imaginary and real parts of the second order SHG susceptibility ${\chi}_{123}^{(2)}\left(\omega \right)={\chi}_{213}^{(2)}\left(\omega \right)$ and ${\chi}_{222}^{(2)}\left(\omega \right)={\chi}_{112}^{(2)}\left(\omega \right)$ were found to be in reasonable agreement with the measurements. We note that any anisotropy in the linear optical susceptibilities will significantly enhance the nonlinear optical susceptibilities. Our calculations show that the 2*ω* inter-band and intra-band contributions to the real and imaginary parts of ${\chi}_{ijk}^{(2)}\left(\omega \right)$ show opposite signs. This fact may be used in future for molecular engineering of the crystals in the desirable directions.

## Declarations

### Acknowledgements

The authors would like to thank the Institute of Physical Biology and Institute of System Biology and Ecology-Computer Center for providing the computational facilities. This work was supported from the institutional research concept of the Institute of Physical Biology, UFB (No. MSM6007665808), and the Institute of System Biology and Ecology, ASCR (No. AVOZ60870520).

## Authors’ Affiliations

## References

- Mills AD: Inorg Chem. 1962, 1: 960-10.1021/ic50004a063.View ArticleGoogle Scholar
- Peterson GA, Keszler DA, Reynolds TA: Int J Inorganic Mat. 2000, 2: 101-10.1016/S1466-6049(00)00011-8.View ArticleGoogle Scholar
- Filimonov AA, Leonyuk NI, Meisner LB, Timchenko TI, Rez IS: Kristall Techn. 1974, 9: 63-10.1002/crat.19740090111.View ArticleGoogle Scholar
- Foldvari I, Beregi E, Baraldi A, Capelleti R, Ryba-Romanowski W, Dominiak-Dzik G, Munoz A, Sosa R: J Luminesc. 2003, 102–103: 395-10.1016/S0022-2313(02)00533-1.View ArticleGoogle Scholar
- Jaque D, Capmany J, Garcia Solé J: Appl Phys Lett. 1999, 75: 325-10.1063/1.124364. P. Decker, JM Dawes, Opt. Express, 12 (2004) 5922.View ArticleADSGoogle Scholar
- Majchrowski A, Kityk IV, Ebothe J: Phys Stat Solidi (B). 2004, 241: 3047-10.1002/pssb.200402071.View ArticleADSGoogle Scholar
- Li G, Cao O, Li Z, Huang Y: J Chinese Ceramic Society. 2007, 35: 381-Google Scholar
- Akhmetov SF, Akhmetova GL, Kovalenko VS, Leonuyk NI, Pashkova AV: Soviet-Phys Dokl. 1978, 23: 107-Google Scholar
- Iwai M, Mori Y, Sasaki T, Nakai S, Sarukura N, Liuand Z, Segawa Y: Jpn J Appl Phys. 1995, 34: 2338-10.1143/JJAP.34.2338.View ArticleADSGoogle Scholar
- Majchrowski A, Kityk IV: Ferroelectr Lett. 2002, 29: 31-10.1080/07315170210696.View ArticleGoogle Scholar
- Akhmetov SF, Akhmetova GL, Kovalenko VS, Leonyuk NI, Pashkova AV: Sov Phys-Dokl. 1978, 23: 107-Google Scholar
- Baraldi A, Capelletti R, Magnani N, Mazera M, Beregi E, Foldvari I: J Phys: Condensed Matter. 2005, 17: 6245-10.1088/0953-8984/17/39/010.ADSGoogle Scholar
- Zhang J-P, Chen G, Zhou H-B: Communic Theoreitlca OPhys. 2006, 45: 1121-ADSGoogle Scholar
- Brik MG: J Phys Chem Solids. 2006, 67: 738-10.1016/j.jpcs.2005.11.005.View ArticleADSGoogle Scholar
- Dekker P, Dawes J: Appl Phys B:Laser and Optics. 2006, 83: 267-10.1007/s00340-006-2168-3. OSA Trends in Optics and Photonic Series, 99, 267.View ArticleADSGoogle Scholar
- Zhao H, Li J, Wang J, H. Zhang, Wang X, Luiang X: J Rare Earths. 2006, 24: 152-Google Scholar
- Kityk IV, Makowska-Janusik M, Fontana MD, Aillerie M: Cryst Res Technol. 2001, 36: 577-10.1002/1521-4079(200107)36:6<577::AID-CRAT577>3.0.CO;2-3.View ArticleGoogle Scholar
- Kityk IV: Phys Solid State. 1991, 33: 1026-Google Scholar
- Blaha P, Schwarz K, Madsen GKH, Kvasnicka D, Luitz J: WIEN2k, An augmented plane wave plus local orbitals program for calculating crystal properties, Vienna University of Technology, Austria. 2001Google Scholar
- Gao S: Computer Physics Communications. 2003, 153: 190-10.1016/S0010-4655(03)00224-8.View ArticleADSGoogle Scholar
- Schwarz K: Journal of Solid State Chemistry. 2003, 176: 319-10.1016/S0022-4596(03)00213-5.View ArticleADSGoogle Scholar
- Hohenberg P, Kohn W: Phys Rev B. 1964, 136: 864-10.1103/PhysRev.136.B864.View ArticleADSMathSciNetGoogle Scholar
- Smok P, Kityk IV, Plucinski KJ, Berdowski J: Physical Review B. 2002, 65: 205103-10.1103/PhysRevB.65.205103.View ArticleADSGoogle Scholar
- Kityk IV, Mefleh A: Physica B: Condensed Matter. 1999, B262: 170-10.1016/S0921-4526(98)00654-1.View ArticleADSGoogle Scholar
- Kohn W, Sham LJ: Phys Rev A. 1965, 140: 1133-10.1103/PhysRev.140.A1133.View ArticleADSMathSciNetGoogle Scholar
- Levine BF: Phys Rev. 1973, B7: 2600-and references therein.View ArticleADSGoogle Scholar
- Hybertsen MS, Louie SG: Phys Rev B. 1986, 34: 5390-10.1103/PhysRevB.34.5390.View ArticleADSGoogle Scholar
- Tributsch H: Z Naturforsch A. 1977, 32A: 972-ADSGoogle Scholar
- Mefleh A, Benet S, Brunet S, Kaddouri H, Sahraoui B, Kityk IV, Makowska-Janusik M: Optical Materials. 1999, 13: 339-10.1016/S0925-3467(99)00078-6.View ArticleADSGoogle Scholar
- Boyd W: Nonlinear Optics. 1992, Academic Press, BostonGoogle Scholar
- Reshak AH, Auluck S: Phys Rev B. 2003, 68: 245113-10.1103/PhysRevB.68.245113.View ArticleADSGoogle Scholar
- Sharma S, Auluck S, Khan MA: Pramana Journal of Physics. 1999, 54: 431-10.1007/s12043-000-0135-9.View ArticleADSGoogle Scholar
- Sharma S, Dewhurst JK, Ambrosch-Draxl C: Phys Rev B. 2003, 67: 165332-10.1103/PhysRevB.67.165332.View ArticleADSGoogle Scholar
- Reshak AH: PhD thesis. 2005, Indian Institute of Technology-Rookee, IndiaGoogle Scholar
- Ambrosch-Draxl C, Sofo JO: Comput Phys Commun. 2006, 175: 1-10.1016/j.cpc.2006.03.005.View ArticleADSGoogle Scholar
- Aspnes DE: Phys Rev B. 1972, 6: 4648-10.1103/PhysRevB.6.4648.View ArticleADSGoogle Scholar
- Ebothe J, Kityk IV, Kisilewski J, Lukasiewicz T, Diduszko R: Appl Phys Lett. 2006, 89: 131106-10.1063/1.2357035.View ArticleADSGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.