1. Introduction
The ternary AXIBXIIIC2XVI
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabdgeabnaaCaaaleqabaGaemiwaGLaemysaKeaaOGaemOqai0aaWbaaSqabeaacqWGybawcqWGjbqscqWGjbqscqWGjbqsaaGccqWGdbWqdaqhaaWcbaGaeGOmaidabaGaemiwaGLaemOvayLaemysaKeaaaaa@3A12@ semiconducting compounds which crystallize in chalcopyrite structure have drawn much attention in recent years 1233166. They form a large group of semiconducting materials with diverse optical, electrical, and structural properties 4567891011. Ternary chalcopyrite appear to be promising candidates for solar-cells applications 126465, light-emitting diodes 13, nonlinear optics 14, and optical frequency conversion applications in all solid state based tunable laser systems. They have potentially significant advantages over dye lasers because of their easier operation and the potential for more compact devices. Tunable frequency conversion in mid infrared (IR) is based on optical parametric oscillators (OPO) using pump lasers in near IR 15. On the other hand frequency doubling devices also allow one to expand the range of powerful lasers to far IR such as CO2 lasers 151617.
In last decade, first-principle calculations have been successfully used to obtain different properties of materials. The structural parameters and dynamical properties of crystals determine a wide range of microscopic and macroscopic behavior: diffraction, sound velocity, elastic constants, Raman and infrared absorption, inelastic neutron scattering, specific heat, etc. The AXIBXIIIC2XVI
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabdgeabnaaCaaaleqabaGaemiwaGLaemysaKeaaOGaemOqai0aaWbaaSqabeaacqWGybawcqWGjbqscqWGjbqscqWGjbqsaaGccqWGdbWqdaqhaaWcbaGaeGOmaidabaGaemiwaGLaemOvayLaemysaKeaaaaa@3A12@ chalcopyrites concerned in this paper use copper as the group XI element, indium as the group XIII element, and sulphur, selenium and tellurium as the group XVI element. Since these compounds display large birefringence 18, they are potentially interesting as nonlinear optical materials, as well as semiconductors. CuInTe2 has positive birefringence whereas CuInS2 and CuInSe2 have negative birefringent 1920. So far, however, the trends of the coupling coefficients in these materials are not well understood. Rashkeev et al. 17 used the linear muffin-tin orbitals (LMTO) method within muffin-tin approximation (MTA) to predict enhancement of χ(2) (ω) by substitution of S by Se and Te. It is well known that MTA works reasonably well in highly coordinated systems, such as face-centered cubic (FCC) metals 56. For covalently bonded solids or layered structures, such as chalcopyrites, the MTA is a poor approximation and leads to discrepancies with experiments 52. The more general treatment of the potential, such as provided by a full potential (FP) method has none of the drawbacks of the atomic sphere approximation (ASA) and MTA based methods. In full potential methods the potential and charge density are expanded into lattice harmonics inside each atomic sphere and as a Fourier series in the interstitial region. In the present work we use the full potential linear augmented plane wave (FP-LAPW) method which has proven to be one of the most accurate methods 5455 for computation of the electronic structure of solids within density functional theory (DFT). Hence the effect of full potential on the linear and nonlinear optical properties can be ascertained.
The dynamical properties of chalcopyrites are well known and have been described by many workers. There is much research interests in their hardness 212223 and pressure-induced behaviors 2425. Recently Eryigit et al. 53 have performed a first principles study of structural, dynamical, and dielectric properties of the chalcopyrite semiconductor CuInS2. Yet up to now there is no comprehensive work that concerns the electronic structure, frequency dependent dielectric function, birefringence and frequency dependent nonlinear optical (NLO) properties of these compounds, although their potential NLO applications has been emphasized. A comprehensive understanding of the origin of the optical nonlinearity and high χ(2) (ω) of these materials are very interesting subjects.
There exist a number of calculations of electronic band structure and optical properties using different methods starting from Miller's empirical rule to the current first principles 262728293031 methods. We are not aware of any full-potential calculations for these compounds. Most the existing ab initio calculations are based on MTA. We therefore think it worthwhile to perform ab initio calculations using a full potential method. The calculated energy gaps vary from 0.01 to 0.812 eV for CuInS2, 0.01 to 0.416 eV for CuInSe2 and 0.18 to 0.424 eV for CuInTe2. Thus there is a large variation in the energy gaps, suggesting that the energy band gap depends on the method of band structure calculation. Also some of the calculated energy gaps are equal to the measured energy gap which is not expected from calculations based on the local density approximation. There is a dearth of theoretical calculations for the frequency dependent dielectric function ε(ω). Also there is no experimental data for the NLO. There is only one report 31 in which they have analyzed the second order of NLO from the chemical bond viewpoint. We think it is worthwhile to perform these calculations. Our calculations will highlight the effect of replacing S by Se and Se by Te on the electronic and optical properties in the CuInX2 compounds. Our motivation in this paper is to understand the origin of the high χ(2) (ω) and the degree of birefringence in these materials as well as to study the trends with moving from S to Se to Te.
2. Details of Calculations
The ternary CuInX2 (X = S, Se, Te) semiconducting compounds crystallize in the chalcopyrite structure with tetragonal space group I4¯2d(D2d12)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabdMeajjqbisda0yaaraGaeGOmaiJaemizaq2aaeWaaeaacqWGebardaqhaaWcbaGaeGOmaiJaemizaqgabaGaeGymaeJaeGOmaidaaaGccaGLOaGaayzkaaaaaa@3659@ having four formula units in each unit cell. CuInX2 is a ternary analog of diamond structure and essentially a superlattice (or superstructure) of zinc blende. Like the atoms in diamond and zinc blende structures, each constituent atom in these ternary compounds, XI, XIII, and XVI, is tetrahedrally coordinated to four neighboring atoms. The Cu atom is located at (0, 0, 0), (0, 0.5, 0.25), In atom at (0, 0, 0.5), (0, 0.5, 0.75), and X atoms at (x, 0.4, 0.125), (-x, 0.75, 0.125), (0.75, x, 0.875), (0.4,- x, 0.875), where x is equal to 0.20, 0.22, 0.225 for CuInS2, CuInSe2, CuInTe2 respectively. In the present work we used the experimental lattice parameters 32 as listed in Table 1. These ternary compounds have a small energy gap and this proved to be of great interest in the nonlinear optical properties 3318.
<p>Table 1</p>
Lattice parameters, energy gaps, our calculated ε1⊥(0)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIXaqmaeaacqGHLkIxaaGcdaqadaqaaiabicdaWaGaayjkaiaawMcaaaaa@31ED@ and ε2II(0)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaeHbbjxAHXgaiqaacaWFjbGaa8xsaaaakmaabmaabaGaeGimaadacaGLOaGaayzkaaaaaa@3B85@ and Δn(0).
CuInS2
CuInSe2
CuInTe2
A (Å)
5.52a
5.78a
6.179a
C (Å)
11.08a
11.57a
12.36a
Egexp.
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabdweafnaaDaaaleaacqWGNbWzaeaacyGGLbqzcqGG4baEcqGGWbaCcqGGUaGlaaaaaa@32A6@ (eV)
1.53b, r, 1.55e, 1.48i, 1.54d
1.04b, 0.96g, 1.01h, 0.98i, r, 0.88d
1.06b, j, 0.96f
EgTheory
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabdweafnaaDaaaleaacqWGNbWzaeaacqWGubavcqWGObaAcqWGLbqzcqWGVbWBcqWGYbGCcqWG5bqEaaaaaa@35BA@ (eV)
0.01c, -0.14r, 0.812*
0.01c, -0.2r, 0.416*
0.18c, 0.424*
ε
1
⊥
(
0
)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIXaqmaeaacqGHLkIxaaGcdaqadaqaaiabicdaWaGaayjkaiaawMcaaaaa@31ED@
11.3*, 9.71p, 8.6q
14.8*
13.5*
ε
1
I
I
(
0
)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIXaqmaeaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaeHbbjxAHXgaiqaacaWFjbGaa8xsaaaakmaabmaabaGaeGimaadacaGLOaGaayzkaaaaaa@3B83@
11*, 9.68p, 8.4q
14.6*
13.8*
Δn(0)
-0.018k, -0.019*
-0.012m, -0.013*
0.002*
N
0
2.532n, o, 2.71*
2.606n, 2.83*
2.720n, 2.91*
aRef.32 bRef.36 cRef.17 dRef.37 eRef.38 fRef. 11 gRef.39 hRef.40 iRef.41 jRef.42 kRef.19 mRef.20 nRef.31 oRef.57 pRef.53 qRef.62 rRef.63
*This work
In our calculations we use the state-of-the-art full potential linear augmented plane wave (FPLAPW) method in a scalar relativistic version as embodied in the WIEN2k code 34. This is an implementation of the DFT with different possible approximations for the exchange-correlation (XC) potential. Exchange and correlation are treated within the local-density approximation (LDA), and scalar relativistic equations are used to obtain self-consistency. Kohn-Sham equations are solved using a basis of linear APW's. The LDA is known to understand energy gaps by up to 40%. Attempts to improve on this by using LDA+U, SIC (self interaction methods) and GW methods. These methods are very time consuming and are a challenge by itself. As the calculation of the nonlinear optical properties are very time consuming. Hence we have taken a simple approach of using LDA and simply correcting the energy by a rigid shift of the valence/conduction bands.
In order to achieve energy eigenvalues convergence, the wave functions in the interstitial region were expanded in plane waves with a cut-off Kmax = 9/RMT, where RMT denotes the smallest atomic sphere radius and Kmax gives the magnitude of the largest K vector in the plane wave expansion. The RMTare taken to be 2.00 atomic units (a.u.) for Cu, In, S, Se and Te respectively. The valence wave functions in side the spheres are expanded up to lmax = 10 while the charge density was Fourier expanded up to Gmax = 14. Self-consistent calculations are considered to be converged when the total energy of the system is stable within 10-4 Ry. The integrals over the Brillouin zone are performed using 250 k-points in the irreducible Brillouin zone (IBZ). The BZ integrations are carried out using the tetrahedron method 6768. The frequency dependent linear optical properties are calculated using 500 k-points and the nonlinear optical properties using 1500 k⇀
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiqbdUgaRzaaoaaaaa@2C6B@ points in the IBZ. Both the plane wave cutoff and the number of k-points were varied to ensure total energy convergence.
3. Results and Discussion
3.1. Band structure and density of states
The band structure, total density of states (TDOS) along with the X-p/s, and Cu-d/p/s partial DOS for CuInX2 (X = S, Se, Te) compounds are shown in Figure 1. In all cases, the valence band maximum (VBM) and the conduction band minimum (CBM) are located at Γ resulting in a direct energy gap of 0.812, 0.416, and 0.424 eV for CuInS2, CuInSe2 and CuInTe2, respectively. CuInSe2 and CuInTe2 have a reduction of the gap in comparison to CuInS2. The reduction in the energy gap can be attributed due to the fact that the conduction bands shift towards Fermi energy (EF) when we move from S to Se to Te. In the conduction bands shifting Cu-s states has a small effect whereas shifting In-s states has a strong effect in opening the energy band gap, while leaving the valence bands unchanged. The overall reduction in the energy band gap is consistent with an overall weakening of the bonds, and, therefore, with a smaller bonding antibonding splitting. A comparison of the experimental and theoretical band gaps is given in Table 1. The calculated energy gaps are smaller than the experimental gaps as expected from an LDA calculation 17. We note that the energy gap decreases when S is replaced by Se/Te in agreement with the experimental data. This trend in the energy band gaps is not present in previous LMTO calculations 17, suggesting that the MTA and ASA are poor approximations and lead to discrepancies on comparison with the experimental data 52. The band structure and hence the DOS can be divided into six groups. From the PDOS we are able to identify the angular momentum character of the various structures. The lowest group has mainly In-d states. The second group between -11.0 to -14.0 eV has significant contributions from X-s states. The third group -6.0 to -7.0 eV is mainly In-X bond. The groups from -5.5 eV up to Fermi energy (EF) are due to Cu-d states with some contribution from X-p states. The electronic structure of the upper valence band is dominated by Cu-d and X-p states. We note that most of Cu-d character is concentrated in the upper valence band. The Cu-s states are pushed from the valence bands into conduction bands. The last group from 0.5 eV and above has contributions from X-p, Cu-spd, and In-sp states. The trends in the band structures (as we move from S to Se to Te) can be summarized as follows: (1). The second group in CuInSe2 is shifted towards lower energies by around 0.5 eV in comparison with CuInS2, while in CuInTe2 it is shifted towards higher energies by around 1 eV, with reduces the bandwidth of both CuInSe2 and CuInTe2 with respect to CuInS2. (2). The bandwidth of third group is increased with moving from S to Se to Te, this group is shifted towards higher energies by around 0.2 eV. (3). The bandwidth of conduction band reduces slightly by around 0.5 eV on going from S to Se to Te causing to increase the gap between the third and fourth groups. From the PDOS, we note a strong hybridization between Cu-d and X-p states around -4 eV. Following Yamasaki et al. 43 we can define degree of hybridization by the ratio of Cu-d states and X-p states within the muffin tin sphere. Based on this we can say that the hybridization between Cu-d and X-p states becomes weak when we move from S to Se to Te. Also we note that Cu-s strongly hybridized with Cu-p, and In-s with In-p states. Cu-s, Cu-p, In-s, and In-p show strong hybridization with In-d states.
<p>Figure 1</p>
Band structure and total density of states (states/eV-unit cell), along with Cu-p/s, X-p, and Cu-d partial densities of states
Band structure and total density of states (states/eV-unit cell), along with Cu-p/s, X-p, and Cu-d partial densities of states.
![]()
3.2. Linear optical response and birefringence
Since the investigated compounds have tetragonal symmetry, we need to calculate two dielectric tensor components, corresponding to electric field E→
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiqbdweafzaalaaaaa@2C1C@ perpendicular and parallel to c-axis, to completely characterize the linear optical properties. These are ε2⊥(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaacqGHLkIxaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@32CE@ and ε2II(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaeHbbjxAHXgaiqaacaWFjbGaa8xsaaaakmaabmaabaGaeqyYdChacaGLOaGaayzkaaaaaa@3C64@, the imaginary parts of frequency dependent dielectric function. We have performed calculations of the frequency dependent dielectric function for these compounds using the expressions 4445
ε
2
I
I
(
ω
)
=
12
m
ω
2
∫
B
Z
∑
|
P
n
n
′
Z
(
k
)
|
2
d
S
k
∇
ω
n
n
′
(
k
)
ε
2
⊥
(
ω
)
=
6
m
ω
2
∫
B
Z
∑
[
|
P
n
n
′
X
(
k
)
|
2
+
|
P
n
n
′
Y
(
k
)
|
2
]
d
S
k
∇
ω
n
n
′
(
k
)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=MipeYlH8Hipec8Eeeu0xXdbba9frFj0xb9Lqpepeea0xd9q8qiYRWxGi6xij=hbbc9s8aq0=yqpe0xbbG8A8frFve9Fve9Fj0dmeaabaqaciGacaGaaeqabaqabeGadaaakeaafaqaaeGabaaabaGaeqyTdu2aa0baaSqaaiabikdaYaqaamXvP5wqSXMqHnxAJn0BKvguHDwzZbqegeKCPfgBaGabaiaa=LeacaWFjbaaaOWaaeWaaeaacqaHjpWDaiaawIcacaGLPaaacqGH9aqpjuaGdaWcaaqaaiabigdaXiabikdaYaqaaiabd2gaTjabeM8a3naaCaaabeqaaiabikdaYaaaaaGcdaWdrbqaamaaqaeajuaGbaWaaSaaaeaadaabdaqaaiabdcfaqnaaDaaabaGaemOBa4MafmOBa4MbauaaaeaacqWGAbGwaaWaaeWaaeaacqWGRbWAaiaawIcacaGLPaaaaiaawEa7caGLiWoadaahaaqabeaacqaIYaGmaaGaemizaqMaem4uam1aaSbaaeaacqWGRbWAaeqaaaqaaiabgEGirlabeM8a3naaBaaabaGaemOBa4MafmOBa4MbauaaaeqaamaabmaabaGaem4AaSgacaGLOaGaayzkaaaaaaWcbeqab0GaeyyeIuoaaSqaaiabdkeacjabdQfaAbqab0Gaey4kIipaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaacqGHLkIxaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaiabg2da9KqbaoaalaaabaGaeGOnaydabaGaemyBa0MaeqyYdC3aaWbaaeqabaGaeGOmaidaaaaakmaapefabaWaaabqaKqbagaadaWcaaqaamaadmaabaWaaqWaaeaacqWGqbaudaqhaaqaaiabd6gaUjqbd6gaUzaafaaabaGaemiwaGfaamaabmaabaGaem4AaSgacaGLOaGaayzkaaaacaGLhWUaayjcSdWaaWbaaeqabaGaeGOmaidaaiabgUcaRmaaemaabaGaemiuaa1aa0baaeaacqWGUbGBcuWGUbGBgaqbaaqaaiabdMfazbaadaqadaqaaiabdUgaRbGaayjkaiaawMcaaaGaay5bSlaawIa7amaaCaaabeqaaiabikdaYaaaaiaawUfacaGLDbaacqWGKbazcqWGtbWudaWgaaqaaiabdUgaRbqabaaabaGaey4bIeTaeqyYdC3aaSbaaeaacqWGUbGBcuWGUbGBgaqbaaqabaWaaeWaaeaacqWGRbWAaiaawIcacaGLPaaaaaaaleqabeqdcqGHris5aaWcbaGaemOqaiKaemOwaOfabeqdcqGHRiI8aaaaaaa@A3A2@
The above expressions are written in atomic units with e2 = 1/m = 2 and ħ = 1. where ω is the photon energy and PnnX(k)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabbcfaqnaaDaaaleaacqWGUbGBcqWGUbGBaeaacqWGybawaaGcdaqadaqaaiabdUgaRbGaayjkaiaawMcaaaaa@3340@ is the x component of the dipolar matrix elements between initial |nk⟩ and final |n'k⟩ states with their eigenvalues En(k) and En' (k), respectively. ωnn' (k) is the energy difference ωnn' (k) = En (k) - En'(k) and Sk is a constant energy surface Sk = {k; ωnn' (k) = ω}
Figure 2 shows the calculated imaginary part of the anisotropic frequency dependent dielectric function ε2⊥(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaacqGHLkIxaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@32CE@ and ε2II(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaeHbbjxAHXgaiqaacaWFjbGaa8xsaaaakmaabmaabaGaeqyYdChacaGLOaGaayzkaaaaaa@3C64@. A considerable anisotropy is found between ε2⊥(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaacqGHLkIxaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@32CE@ and ε2II(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaeHbbjxAHXgaiqaacaWFjbGaa8xsaaaakmaabmaabaGaeqyYdChacaGLOaGaayzkaaaaaa@3C64@. Broadening is taken to be 0.02 eV. It is well known that LDA calculations underestimate the energy gaps. A very simple way to overcome this drawback is to use the scissors correction, which merely makes the calculated energy gap equal to the experimental gap. Our calculated optical properties are scissors corrected 1435 by 0.716, 0.622, and 0.635 eV for CuInS2, CuInSe2, and CuInTe2 respectively. These values are the difference between the calculated and measured energy gap (Table 1). We note that the magnitude of ε2(ω) is increased when move from S to Se to Te. Our calculated ε2⊥(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaacqGHLkIxaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@32CE@ and ε2II(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaeHbbjxAHXgaiqaacaWFjbGaa8xsaaaakmaabmaabaGaeqyYdChacaGLOaGaayzkaaaaaa@3C64@ show reasonable agreement with the experimental data 20.
<p>Figure 2</p>
Calculated ε2⊥(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaacqGHLkIxaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@32CE@ and ε2II(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaeHbbjxAHXgaiqaacaWFjbGaa8xsaaaakmaabmaabaGaeqyYdChacaGLOaGaayzkaaaaaa@3C64@ (dark curve) along with the experimental data20 (light curve) for CuInS2 and CuInSe2, and ε2⊥(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaacqGHLkIxaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@32CE@ (dark curve) and ε2II(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaeHbbjxAHXgaiqaacaWFjbGaa8xsaaaakmaabmaabaGaeqyYdChacaGLOaGaayzkaaaaaa@3C64@ (light curve) for CuInTe2
Calculated ε2⊥(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaacqGHLkIxaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@32CE@ and ε2II(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaeHbbjxAHXgaiqaacaWFjbGaa8xsaaaakmaabmaabaGaeqyYdChacaGLOaGaayzkaaaaaa@3C64@ (dark curve) along with the experimental data20 (light curve) for CuInS2 and CuInSe2, and ε2⊥(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaacqGHLkIxaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@32CE@ (dark curve) and ε2II(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaeHbbjxAHXgaiqaacaWFjbGaa8xsaaaakmaabmaabaGaeqyYdChacaGLOaGaayzkaaaaaa@3C64@ (light curve) for CuInTe2.
![]()
It would be worthwhile to attempt to identify the transitions that are responsible for the structures in ε2⊥(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaacqGHLkIxaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@32CE@ and ε2II(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaeHbbjxAHXgaiqaacaWFjbGaa8xsaaaakmaabmaabaGaeqyYdChacaGLOaGaayzkaaaaaa@3C64@ using our calculated band structure. Generally the peaks in the optical response are caused by the electric-dipole transitions between the valence and conduction bands. These peaks in the linear optical spectra can be identified from the band structure. Figure 3, show the transitions which are responsible for the structures in ε2⊥(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaacqGHLkIxaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@32CE@ and ε2II(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaeHbbjxAHXgaiqaacaWFjbGaa8xsaaaakmaabmaabaGaeqyYdChacaGLOaGaayzkaaaaaa@3C64@. In order to identify these structures we need to look at the optical matrix elements. We mark the transitions, giving the major structure in ε2⊥(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaacqGHLkIxaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@32CE@ and ε2II(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaeHbbjxAHXgaiqaacaWFjbGaa8xsaaaakmaabmaabaGaeqyYdChacaGLOaGaayzkaaaaaa@3C64@ in the band structure diagram. These transitions are labeled according to the optical transitions in Figure 2. For simplicity we labeled the transitions in Figure 3, as A, B, and C. The transitions (A) are responsible for the structures of ε2⊥(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaacqGHLkIxaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@32CE@ and ε2II(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaeHbbjxAHXgaiqaacaWFjbGaa8xsaaaakmaabmaabaGaeqyYdChacaGLOaGaayzkaaaaaa@3C64@ in the energy range between 0.0 eV and 2 eV, the transitions (B) are responsible for the structures of ε2⊥(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaacqGHLkIxaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@32CE@ and ε2II(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaeHbbjxAHXgaiqaacaWFjbGaa8xsaaaakmaabmaabaGaeqyYdChacaGLOaGaayzkaaaaaa@3C64@ in the energy range 2.0 eV and 4.0 eV, and the transitions (C) are responsible for the structures of ε2⊥(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaacqGHLkIxaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@32CE@ and ε2II(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaeHbbjxAHXgaiqaacaWFjbGaa8xsaaaakmaabmaabaGaeqyYdChacaGLOaGaayzkaaaaaa@3C64@ between 4.0 eV and 6.0 eV.
<p>Figure 3</p>
The optical transitions depicted on a generic band structure
The optical transitions depicted on a generic band structure.
![]()
The calculated ε1⊥(0)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIXaqmaeaacqGHLkIxaaGcdaqadaqaaiabicdaWaGaayjkaiaawMcaaaaa@31ED@ and ε2II(0)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaeHbbjxAHXgaiqaacaWFjbGaa8xsaaaakmaabmaabaGaeGimaadacaGLOaGaayzkaaaaaa@3B85@ is listed in Table I. We note that a smaller energy gap yields a larger ε1⊥(0)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIXaqmaeaacqGHLkIxaaGcdaqadaqaaiabicdaWaGaayjkaiaawMcaaaaa@31ED@ and ε2II(0)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabew7aLnaaDaaaleaacqaIYaGmaeaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaeHbbjxAHXgaiqaacaWFjbGaa8xsaaaakmaabmaabaGaeGimaadacaGLOaGaayzkaaaaaa@3B85@ value. This could be explained on the basis of the Penn model 46, where ε(0) is related to Eg by the equation ε(0) ≈ 1 + (ħωP/Eg)2. Hence smaller Eg yields a larger ε(0).
These compounds show considerable anisotropy which is lead to an important quantity in SHG and OPO. The important quantity is the birefringence. The birefringence is important to fulfill phase-matching condition. The birefringence can be calculated from the linear response functions from which the anisotropy of the index of refraction is obtained. Using the expression 47
n
(
ω
)
=
(
1
/
2
)
[
ε
1
(
ω
)
2
+
ε
2
(
ω
)
2
+
ε
1
(
ω
)
]
1
/
2
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=MipeYlH8Hipec8Eeeu0xXdbba9frFj0xb9Lqpepeea0xd9q8qiYRWxGi6xij=hbbc9s8aq0=yqpe0xbbG8A8frFve9Fve9Fj0dmeaabaqaciGacaGaaeqabaqabeGadaaakeaacqWGUbGBdaqadaqaaiabeM8a3bGaayjkaiaawMcaaiabg2da9maabmaabaGaeGymaeJaei4la8YaaOaaaeaacqaIYaGmaSqabaaakiaawIcacaGLPaaadaWadaqaamaakaaabaGaeqyTdu2aaSbaaSqaaiabigdaXaqabaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaamaaCaaaleqabaGaeGOmaidaaOGaey4kaSIaeqyTdu2aaSbaaSqaaiabikdaYaqabaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaamaaCaaaleqabaGaeGOmaidaaaqabaGccqGHRaWkcqaH1oqzdaWgaaWcbaGaeGymaedabeaakmaabmaabaGaeqyYdChacaGLOaGaayzkaaaacaGLBbGaayzxaaWaaWbaaSqabeaacqaIXaqmcqGGVaWlcqaIYaGmaaaaaa@5005@
we can determine the value of the extraordinary and ordinary refraction indices. The birefringence is the difference between the extraordinary and ordinary refraction indices, Δn = ne - n0, where ne is the index of refraction for an electric field oriented along c-axis and n0 is the index of refraction for an electric field perpendicular to c-axis.
Figure 4 shows the birefringence Δn(ω) for CuInX2 compounds. The birefringence is important only in the non-absorbing region, which is below the gap. In general we note that the shape of Δn(ω) for the three compounds is rather similar. This attributed due to the fact that these compounds have similar band structures. This curve shows strong oscillations around zero in the energy range up to 12.5 eV. Thereafter it drops to zero. We find that the birefringence is negative for CuInS2 and CuInSe2 and positive for CuInTe2 in agreement with the experimental data 1920 (Table 1). We have compared our calculated refractive index with that obtained by D. Xue et al. 31 in Table 1. Good agreement is found. This agreement will be in the favor of our calculated χijk(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqWGPbqAcqWGQbGAcqWGRbWAaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@36F6@. Previous work of D. Xue et al. 58, discussed the correlation between the nonlinear tensor dij and refractive index (n0). They compared the theoretical predictions of the nonlinear optical tensor coefficient d11 of K4Gd2(CO3)3Fe4 with the experimental data and found reasonable agreement with the experimental data of Mercier et al. 59. This means that there is a correlation between the values of the refractive indices (n0) and nonlinear tensor coefficient dij. Based on these results we can say that if our calculated refractive indices (n0) are in good agreement with that obtained by D. Xue et al. 31, then we expect our nonlinear optical susceptibilities χ312(2)(0)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqaIZaWmcqaIXaqmcqaIYaGmaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabicdaWaGaayjkaiaawMcaaaaa@34D6@ to be in good agreement with the nonlinear tensor coefficient d36 which is obtained by D. Xue et al. 31.
<p>Figure 4</p>
Calculated Δn(ω) for CuInX2 compounds
Calculated Δn(ω) for CuInX2 compounds.
![]()
3.3. Non-linear optical response
The complex second-order nonlinear optical susceptibility tensor χabc(2)(−2ω;ω;ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqWGHbqycqWGIbGycqWGJbWyaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabgkHiTiabikdaYiabeM8a3jabcUda7iabeM8a3jabcUda7iabeM8a3bGaayjkaiaawMcaaaaa@3E3B@ can be generally written as the sum of three physically different contributions in the form 48:
χ
int
e
r
a
b
c
(
−
2
ω
;
ω
;
ω
)
=
e
3
ℏ
2
Ω
∑
c
v
n
,
k
r
v
c
a
{
r
c
n
b
r
n
v
c
}
(
ω
n
v
−
ω
c
n
)
[
2
f
v
c
ω
c
v
−
2
ω
+
f
n
c
ω
n
c
−
ω
+
f
v
n
ω
v
n
−
ω
]
χ
int
r
a
a
b
c
(
−
2
ω
;
ω
;
ω
)
=
i
e
3
2
ℏ
2
Ω
∑
c
v
,
k
f
v
c
[
2
ω
c
v
(
ω
c
v
−
2
ω
)
r
v
c
a
(
r
v
c
;
c
b
+
r
c
v
;
b
c
)
+
1
ω
c
v
(
ω
c
v
−
ω
)
(
r
v
c
;
c
a
r
c
v
b
+
r
v
c
;
b
a
r
c
v
c
)
+
1
ω
c
v
2
(
1
ω
c
v
−
ω
−
4
ω
c
v
−
2
ω
)
r
v
c
a
(
r
c
v
b
Δ
c
v
c
+
r
c
v
c
Δ
c
v
b
)
−
1
2
ω
c
v
(
ω
c
v
−
ω
)
(
r
v
c
;
a
b
r
c
v
c
+
r
v
c
;
a
c
r
c
v
b
)
]
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=MipeYlH8Hipec8Eeeu0xXdbba9frFj0xb9Lqpepeea0xd9q8qiYRWxGi6xij=hbbc9s8aq0=yqpe0xbbG8A8frFve9Fve9Fj0dmeaabaqaciGacaGaaeqabaqabeGadaaakeaafaqaaeWabaaabaGaeq4Xdm2aa0baaSqaaiGbcMgaPjabc6gaUjabcsha0jabdwgaLjabdkhaYbqaaiabdggaHjabdkgaIjabdogaJbaakmaabmaabaGaeyOeI0IaeGOmaiJaeqyYdCNaei4oaSJaeqyYdCNaei4oaSJaeqyYdChacaGLOaGaayzkaaGaeyypa0tcfa4aaSaaaeaacqWGLbqzdaahaaqabeaacqaIZaWmaaaabaGaeS4dHG2aaWbaaeqabaGaeGOmaidaaiabfM6axbaakmaaqafabaqcfa4aaSaaaeaacqWGYbGCdaqhaaqaaiabdAha2jabdogaJbqaaiabdggaHbaadaGadaqaaiabdkhaYnaaDaaabaGaem4yamMaemOBa4gabaGaemOyaigaaiabdkhaYnaaDaaabaGaemOBa4MaemODayhabaGaem4yamgaaaGaay5Eaiaaw2haaaqaamaabmaabaGaeqyYdC3aaSbaaeaacqWGUbGBcqWG2bGDaeqaaiabgkHiTiabeM8a3naaBaaabaGaem4yamMaemOBa4gabeaaaiaawIcacaGLPaaaaaGcdaWadaqaaKqbaoaalaaabaGaeGOmaiJaemOzay2aaSbaaeaacqWG2bGDcqWGJbWyaeqaaaqaaiabeM8a3naaBaaabaGaem4yamMaemODayhabeaacqGHsislcqaIYaGmcqaHjpWDaaGccqGHRaWkjuaGdaWcaaqaaiabdAgaMnaaBaaabaGaemOBa4Maem4yamgabeaaaeaacqaHjpWDdaWgaaqaaiabd6gaUjabdogaJbqabaGaeyOeI0IaeqyYdChaaOGaey4kaSscfa4aaSaaaeaacqWGMbGzdaWgaaqaaiabdAha2jabd6gaUbqabaaabaGaeqyYdC3aaSbaaeaacqWG2bGDcqWGUbGBaeqaaiabgkHiTiabeM8a3baaaOGaay5waiaaw2faaaWcbaGaem4yamMaemODayNaemOBa4MaeiilaWIaem4AaSgabeqdcqGHris5aaGcbaGaeq4Xdm2aa0baaSqaaiGbcMgaPjabc6gaUjabcsha0jabdkhaYjabdggaHbqaaiabdggaHjabdkgaIjabdogaJbaakmaabmaabaGaeyOeI0IaeGOmaiJaeqyYdCNaei4oaSJaeqyYdCNaei4oaSJaeqyYdChacaGLOaGaayzkaaGaeyypa0tcfa4aaSaaaeaacqWGPbqAcqWGLbqzdaahaaqabeaacqaIZaWmaaaabaGaeGOmaiJaeS4dHG2aaWbaaeqabaGaeGOmaidaaiabfM6axbaakmaaqafabaGaemOzay2aaSbaaSqaaiabdAha2jabdogaJbqabaaabaGaem4yamMaemODayNaeiilaWIaem4AaSgabeqdcqGHris5aOWaamqaaKqbagaadaWcaaqaaiabikdaYaqaaiabeM8a3naaBaaabaGaem4yamMaemODayhabeaadaqadaqaaiabeM8a3naaBaaabaGaem4yamMaemODayhabeaacqGHsislcqaIYaGmcqaHjpWDaiaawIcacaGLPaaaaaaakiaawUfaaiabdkhaYnaaDaaaleaacqWG2bGDcqWGJbWyaeaacqWGHbqyaaGcdaqadaqaaiabdkhaYnaaDaaaleaacqWG2bGDcqWGJbWycqGG7aWocqWGJbWyaeaacqWGIbGyaaGccqGHRaWkcqWGYbGCdaqhaaWcbaGaem4yamMaemODayNaei4oaSJaemOyaigabaGaem4yamgaaaGccaGLOaGaayzkaaGaey4kaSscfa4aaSaaaeaacqaIXaqmaeaacqaHjpWDdaWgaaqaaiabdogaJjabdAha2bqabaWaaeWaaeaacqaHjpWDdaWgaaqaaiabdogaJjabdAha2bqabaGaeyOeI0IaeqyYdChacaGLOaGaayzkaaaaaOWaaeWaaeaacqWGYbGCdaqhaaWcbaGaemODayNaem4yamMaei4oaSJaem4yamgabaGaemyyaegaaOGaemOCai3aa0baaSqaaiabdogaJjabdAha2bqaaiabdkgaIbaakiabgUcaRiabdkhaYnaaDaaaleaacqWG2bGDcqWGJbWycqGG7aWocqWGIbGyaeaacqWGHbqyaaGccqWGYbGCdaqhaaWcbaGaem4yamMaemODayhabaGaem4yamgaaaGccaGLOaGaayzkaaGaey4kaScabaWaamGaaeaajuaGdaWcaaqaaiabigdaXaqaaiabeM8a3naaDaaabaGaem4yamMaemODayhabaGaeGOmaidaaaaakmaabmaabaqcfa4aaSaaaeaacqaIXaqmaeaacqaHjpWDdaWgaaqaaiabdogaJjabdAha2bqabaGaeyOeI0IaeqyYdChaaOGaeyOeI0scfa4aaSaaaeaacqaI0aanaeaacqaHjpWDdaWgaaqaaiabdogaJjabdAha2bqabaGaeyOeI0IaeGOmaiJaeqyYdChaaaGccaGLOaGaayzkaaGaemOCai3aa0baaSqaaiabdAha2jabdogaJbqaaiabdggaHbaakmaabmaabaGaemOCai3aa0baaSqaaiabdogaJjabdAha2bqaaiabdkgaIbaakiabfs5aenaaDaaaleaacqWGJbWycqWG2bGDaeaacqWGJbWyaaGccqGHRaWkcqWGYbGCdaqhaaWcbaGaem4yamMaemODayhabaGaem4yamgaaOGaeuiLdq0aa0baaSqaaiabdogaJjabdAha2bqaaiabdkgaIbaaaOGaayjkaiaawMcaaiabgkHiTKqbaoaalaaabaGaeGymaedabaGaeGOmaiJaeqyYdC3aaSbaaeaacqWGJbWycqWG2bGDaeqaamaabmaabaGaeqyYdC3aaSbaaeaacqWGJbWycqWG2bGDaeqaaiabgkHiTiabeM8a3bGaayjkaiaawMcaaaaakmaabmaabaGaemOCai3aa0baaSqaaiabdAha2jabdogaJjabcUda7iabdggaHbqaaiabdkgaIbaakiabdkhaYnaaDaaaleaacqWGJbWycqWG2bGDaeaacqWGJbWyaaGccqGHRaWkcqWGYbGCdaqhaaWcbaGaemODayNaem4yamMaei4oaSJaemyyaegabaGaem4yamgaaOGaemOCai3aa0baaSqaaiabdogaJjabdAha2bqaaiabdkgaIbaaaOGaayjkaiaawMcaaaGaayzxaaaaaaaa@8D6F@
In which {rnvbrvcc}=(1/2)(rnvbrvcc+rnvcrvcb)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaamaacmaabaGaemOCai3aa0baaSqaaiabd6gaUjabdAha2bqaaiabdkgaIbaakiabdkhaYnaaDaaaleaacqWG2bGDcqWGJbWyaeaacqWGJbWyaaaakiaawUhacaGL9baacqGH9aqpdaqadaqaamaalyaabaGaeGymaedabaGaeGOmaidaaaGaayjkaiaawMcaamaabmaabaGaemOCai3aa0baaSqaaiabd6gaUjabdAha2bqaaiabdkgaIbaakiabdkhaYnaaDaaaleaacqWG2bGDcqWGJbWyaeaacqWGJbWyaaGccqGHRaWkcqWGYbGCdaqhaaWcbaGaemOBa4MaemODayhabaGaem4yamgaaOGaemOCai3aa0baaSqaaiabdAha2jabdogaJbqaaiabdkgaIbaaaOGaayjkaiaawMcaaaaa@56A0@ is a symmetrized combination of the dipole matrix elements rcna=δcnPcna/imωcn
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabdkhaYnaaDaaaleaacqWGJbWycqWGUbGBaeaacqWGHbqyaaGccqGH9aqpdaWcgaqaaiabes7aKnaaBaaaleaacqWGJbWycqWGUbGBaeqaaOGaemiuaa1aa0baaSqaaiabdogaJjabd6gaUbqaaiabdggaHbaaaOqaaiabdMgaPjabd2gaTjabeM8a3naaBaaaleaacqWGJbWycqWGUbGBaeqaaaaaaaa@430F@, which are in turn obtained from the momentum matrix elements Pcna
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabdcfaqnaaDaaaleaacqWGJbWycqWGUbGBaeaacqWGHbqyaaaaaa@304C@. Superscripts (a, b, c) refer to the Cartesian coordinates. Here v stands for a valence band state and c for the conduction band state and n the intermediate band, which is either in the valence or in the conduction band. It has been demonstrated by Aspnes 49 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. Here we have ignored the virtual-hole contribution (transitions between two valence band states and one conduction band state) because it was shown to be negative and more than an order of magnitude smaller than the virtual-electron contribution for these compounds. For simplicity we call χabc(2)(−2ω;ω;ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqWGHbqycqWGIbGycqWGJbWyaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabgkHiTiabikdaYiabeM8a3jabcUda7iabeM8a3jabcUda7iabeM8a3bGaayjkaiaawMcaaaaa@3E3B@ as χabc(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqWGHbqycqWGIbGycqWGJbWyaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@36C6@. Since the CuInX2 compounds are belong to the point group 4¯2m
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiqbisda0yaaraGaeGOmaiJaemyBa0gaaa@2E5A@ there are only two independent components of the SHG tensor, namely, 123 and 312 components (1, 2, and 3 refer to x, y and z axes, respectively) 50. These are χ123(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqaIXaqmcqaIYaGmcqaIZaWmaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@35B5@ and χ312(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqaIZaWmcqaIXaqmcqaIYaGmaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@35B5@. In the static limit, these two components are equal according to the Kleninman "permutation" symmetry, which dictates additional relations between tonsorial components beyond the purely crystallographic symmetry.
The second-order nonlinear optical susceptibility is very sensitive to the scissors correction. The scissors correction has a profound effect on magnitude and sign of χijk(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqWGPbqAcqWGQbGAcqWGRbWAaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@36F6@. That is attributed to the fact that local density approximation calculations underestimate the energy gaps. Moreover, it is clear from the formulae of χijk(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqWGPbqAcqWGQbGAcqWGRbWAaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@36F6@ that, χijk(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqWGPbqAcqWGQbGAcqWGRbWAaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@36F6@ depends on the energy gap. Hence when we used the scissors correction we found a considerable effect on χijk(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqWGPbqAcqWGQbGAcqWGRbWAaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@36F6@. It is well known that nonlinear optical properties are more sensitive to small changes in the band structure than the linear optical properties. Hence any anisotropy in the linear optical properties is enhanced in the nonlinear spectra. This is attributed to the fact that the second harmonic response χijk(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqWGPbqAcqWGQbGAcqWGRbWAaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@36F6@ involves 2ω resonance in addition to the usual ω resonance. Both ω and 2ω resonances can be further separated into inter-band and intra-band contributions.
The calculated imaginary part of SHG susceptibility χ123(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqaIXaqmcqaIYaGmcqaIZaWmaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@35B5@ is shown in Fig. 5. We show the 2ω inter-band and intra-band contributions for CuInX2 compounds. We note the opposite signs of the two contributions throughout the frequency range. Both these contributions and the total Im χ123(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqaIXaqmcqaIYaGmcqaIZaWmaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@35B5@ increase on moving from S to Se and decreases on moving from Se to Te. As we note from Figure 1, CuInSe2 and CuInTe2 has a reduction of the energy gap relative to CuInS2. The bandwidth of main groups in band structure and DOS of CuInSe2 and CuInTe2 increases over the value in CuInS2 resulting in a reduction of the gap. Im χ123(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqaIXaqmcqaIYaGmcqaIZaWmaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@35B5@ inCuInS2 is smaller than in CuInSe2 and CuInTe2. From this we conclude that the region around the conduction band minimum does not make a significant contribution to χ123(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqaIXaqmcqaIYaGmcqaIZaWmaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@35B5@. Also the VBM and CBM in CuInSe2 and CuInTe2 seem to be more parallel than in CuInS2. This gives an additional amplifying factor when one takes the integrals over the Brillouin zone, making the oscillator strength of the entire peak larger. This clearly shows that trends in χ123(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqaIXaqmcqaIYaGmcqaIZaWmaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@35B5@ cannot be based only on the minimum band gap. This indicates that the average band gap plays a more significant role.
<p>Figure 5</p>
Calculated Im χ123(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqaIXaqmcqaIYaGmcqaIZaWmaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@35B5@ along with the (2ω) intra-band and inter-band contributions for CuInX2 compounds
Calculated Im χ123(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqaIXaqmcqaIYaGmcqaIZaWmaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@35B5@ along with the (2ω) intra-band and inter-band contributions for CuInX2 compounds.
![]()
In Fig. 6 we present the inter-band and intra-band contributions to ω and 2ω resonances for CuInTe2. We find that ω resonance is smaller than 2ω resonance. As can be seen the total SHG susceptibility is zero below half the band gap. The 2ω terms start contributing at energies ~1/2 Eg and the ω terms for energy values above Eg. In the low energy regime (≤ eV) the SHG optical spectra is dominated by 2ω contributions. Beyond 3 eV the major contribution comes from ω terms.
<p>Figure 6</p>
Calculated total Im χ123(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqaIXaqmcqaIYaGmcqaIZaWmaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@35B5@, along with the (2ω)/(1ω) intra-band and inter-band contributions for CuInTe2 compound
Calculated total Im χ123(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqaIXaqmcqaIYaGmcqaIZaWmaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@35B5@, along with the (2ω)/(1ω) intra-band and inter-band contributions for CuInTe2 compound.
![]()
The structures in Im χ123(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqaIXaqmcqaIYaGmcqaIZaWmaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@35B5@ can be understood from the structures in ε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 can make use of the linear optical spectra to identify the different resonance leading to various features in the SHG spectra. The first structure in Im χ123(2)(ω)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqaIXaqmcqaIYaGmcqaIZaWmaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabeM8a3bGaayjkaiaawMcaaaaa@35B5@ between 1.0 – 3.0 eV for CuInS2, 0.5 – 2.0 eV for CuInSe2, and 1.0 – 2.5 eV for CuInTe2 is associated with interference between a ω resonance and 2ω resonance and arises from the first hump in ε2(ω). The second structure between 3.0 – 4.5 eV for CuInS2, 2.0 – 3.0 eV for CuInSe2, and 2.5 – 4.5 eV for CuInTe2 is due mainly to ω resonance and associated with high peak in ε2(ω). The last structure from 4.5 – 5.5 for CuInS2, 3.0 – 5.5 eV for CuInSe2, and 4.5 – 5.5 eV for CuInTe2 is manly due to ω resonance and associated with the tail in ε2(ω).
In Table 2, we present the values of χ123(2)(0)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqaIXaqmcqaIYaGmcqaIZaWmaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabicdaWaGaayjkaiaawMcaaaaa@34D6@. These values clearly increase on going from S to Se to Te in agreement with experiment and theoretical calculations. We note that there is a large variation in the theoretical and experimental values of χ123(2)(0)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqaIXaqmcqaIYaGmcqaIZaWmaeaacqGGOaakcqaIYaGmcqGGPaqkaaGcdaqadaqaaiabicdaWaGaayjkaiaawMcaaaaa@34D6@, suggesting that these values depend on the method of calculations/measurements. Also some of the calculated vales are equal to the measured one. However our calculated χ123(2)(0)
MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=wiFfYlOipiY=Hhbbf9v8qqaqFr0xc9vqpe0di9q8qqpG0dHiVcFbIOFHK8Feei0lXdar=Jb9qqFfeaYRXxe9vr0=vr0=LqpWqaaeaabiGaciaacaqabeaabeqacmaaaOqaaiabeE8aJnaaDaaaleaacqaIXaqmcq