# Electronic structure, linear, nonlinear optical susceptibilities and birefringence of CuInX_{2}(X = S, Se, Te) chalcopyrite-structure compounds

- Ali Hussain Reshak
^{1}Email author and - S Auluck
^{2}

**Received: **21 December 2007

**Accepted: **19 May 2008

**Published: **19 May 2008

## Abstract

The electronic structure, linear and nonlinear optical properties have been calculated for CuInX_{2} (X = S, Se, Te) chalcopyrite-structure single crystals using the state-of-the-art full potential linear augmented plane wave (FP-LAPW) method. We present results for band structure, density of states, and imaginary part of the frequency-dependent linear and nonlinear optical susceptibilities. We find that these crystals are semiconductors with direct band gaps. We have calculated the birefringence of these crystals. The birefringence is negative for CuInS_{2} and CuInSe_{2} while it is positive for CuInTe_{2} in agreement with the experimental data. Calculations are reported for the frequency-dependent complex second-order non-linear optical susceptibilities ${\chi}_{123}^{(2)}\left(\omega \right)$. The intra-band and inter-band contributions to the second harmonic generation increase when we replace S by Se and decrease when we replace Se by Te. We find that smaller energy band gap compounds have larger values of ${\chi}_{123}^{(2)}\left(0\right)$ in agreement with the experimental data and previous theoretical calculations.

**PACS Codes:** 71.15.-m, 31.15.-p

## 1. Introduction

The ternary ${A}^{XI}{B}^{XIII}{C}_{2}^{XVI}$ semiconducting compounds which crystallize in chalcopyrite structure have drawn much attention in recent years [1–3, 31, 66]. They form a large group of semiconducting materials with diverse optical, electrical, and structural properties [4–11]. Ternary chalcopyrite appear to be promising candidates for solar-cells applications [12, 64, 65], 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 CO_{2} lasers [15–17].

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 ${A}^{XI}{B}^{XIII}{C}_{2}^{XVI}$ 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. CuInTe_{2} has positive birefringence whereas CuInS_{2} and CuInSe_{2} have negative birefringent [19, 20]. 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 [54, 55] 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 [21–23] and pressure-induced behaviors [24, 25]. Recently Eryigit et al. [53] have performed a first principles study of structural, dynamical, and dielectric properties of the chalcopyrite semiconductor CuInS_{2}. 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 [26–31] 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 CuInS_{2}, 0.01 to 0.416 eV for CuInSe_{2} and 0.18 to 0.424 eV for CuInTe_{2}. 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 CuInX_{2} 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

_{2}(X = S, Se, Te) semiconducting compounds crystallize in the chalcopyrite structure with tetragonal space group $I\overline{4}2d\left({D}_{2d}^{12}\right)$ having four formula units in each unit cell. CuInX

_{2}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 CuInS

_{2}, CuInSe

_{2}, CuInTe

_{2}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 [33, 18].

Lattice parameters, energy gaps, our calculated ${\epsilon}_{1}^{\perp}\left(0\right)$ and ${\epsilon}_{2}^{II}\left(0\right)$ and Δ*n*(0).

CuInS | CuInSe | CuInTe | |
---|---|---|---|

A (Å) | 5.52 | 5.78 | 6.179 |

C (Å) | 11.08 | 11.57 | 12.36 |

${E}_{g}^{\mathrm{exp}\phantom{\rule{0.5em}{0ex}}.}$ (eV) | 1.53 | 1.04 | 1.06 |

${E}_{g}^{Theory}$ (eV) | 0.01 | 0.01 | 0.18 |

${\epsilon}_{1}^{\perp}\left(0\right)$ | 11.3*, 9.71 | 14.8* | 13.5* |

${\epsilon}_{1}^{II}\left(0\right)$ | 11*, 9.68 | 14.6* | 13.8* |

Δ | -0.018 | -0.012 | 0.002* |

| 2.532 | 2.606 | 2.720 |

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 *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.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 *l*_{max} = 10 while the charge density was Fourier expanded up to *G*_{max} = 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 [67, 68]. The frequency dependent linear optical properties are calculated using 500 k-points and the nonlinear optical properties using 1500 $\stackrel{\rightharpoonup}{k}$ 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

_{2}(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 CuInS

_{2}, CuInSe

_{2}and CuInTe

_{2}, respectively. CuInSe

_{2}and CuInTe

_{2}have a reduction of the gap in comparison to CuInS

_{2}. The reduction in the energy gap can be attributed due to the fact that the conduction bands shift towards Fermi energy (E

_{F}) 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 (E

_{F}) 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 CuInSe

_{2}is shifted towards lower energies by around 0.5 eV in comparison with CuInS

_{2}, while in CuInTe

_{2}it is shifted towards higher energies by around 1 eV, with reduces the bandwidth of both CuInSe

_{2}and CuInTe

_{2}with respect to CuInS

_{2}. (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.

### 3.2. Linear optical response 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 eigenvalues E_{n}(*k*) and E_{
n'
}(*k*), respectively. *ω*_{
nn'
}(*k*) is the energy difference *ω*_{
nn'
}(*k*) = E_{n} (*k*) - E_{
n'
}(*k*) and S_{
k
}is a constant energy surface S_{k} = {*k*; *ω*_{
nn'
}(*k*) = *ω*}

_{2}, CuInSe

_{2}, and CuInTe

_{2}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 ${\epsilon}_{2}^{\perp}\left(\omega \right)$ and ${\epsilon}_{2}^{II}\left(\omega \right)$ show reasonable agreement with the experimental data [20].

The calculated ${\epsilon}_{1}^{\perp}\left(0\right)$ and ${\epsilon}_{2}^{II}\left(0\right)$ is listed in Table I. We note that a smaller energy gap yields a larger ${\epsilon}_{1}^{\perp}\left(0\right)$ and ${\epsilon}_{2}^{II}\left(0\right)$ value. This could be explained on the basis of the Penn model [46], where *ε*(0) is related to *E*_{
g
}by the equation *ε*(0) ≈ 1 + (*ħω*_{
P
}/*E*_{
g
})^{2}. Hence smaller *E*_{
g
}yields a larger *ε*(0).

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* = *n*_{
e
}- *n*_{
0
}, where *n*_{
e
}is the index of refraction for an electric field oriented along **c**-axis and *n*_{0} is the index of refraction for an electric field perpendicular to **c**-axis.

*n*(

*ω*) for CuInX

_{2}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 CuInS

_{2}and CuInSe

_{2}and positive for CuInTe

_{2}in agreement with the experimental data [19, 20] (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 ${\chi}_{ijk}^{(2)}\left(\omega \right)$. Previous work of D. Xue et al. [58], discussed the correlation between the nonlinear tensor

*dij*and refractive index (

*n*

_{0}). They compared the theoretical predictions of the nonlinear optical tensor coefficient

*d*

_{ 11 }of K

_{4}Gd

_{2}(CO

_{3})

_{3}Fe

_{4}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 (

*n*

_{0}) and nonlinear tensor coefficient

*dij*. Based on these results we can say that if our calculated refractive indices (

*n*

_{0}) are in good agreement with that obtained by D. Xue et al. [31], then we expect our nonlinear optical susceptibilities ${\chi}_{312}^{(2)}\left(0\right)$ to be in good agreement with the nonlinear tensor coefficient

*d*

_{ 36 }which is obtained by D. Xue et al. [31].

### 3.3. Non-linear optical response

In which $\left\{{r}_{nv}^{b}{r}_{vc}^{c}\right\}=\left(1/2\right)\left({r}_{nv}^{b}{r}_{vc}^{c}+{r}_{nv}^{c}{r}_{vc}^{b}\right)$ is a symmetrized combination of the dipole matrix elements ${r}_{cn}^{a}={\delta}_{cn}{P}_{cn}^{a}/im{\omega}_{cn}$, which are in turn obtained from the momentum matrix elements ${P}_{cn}^{a}$. 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 ${\chi}_{abc}^{(2)}\left(-2\omega ;\omega ;\omega \right)$ as ${\chi}_{abc}^{(2)}\left(\omega \right)$. Since the CuInX_{2} compounds are belong to the point group $\overline{4}2m$ 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 ${\chi}_{123}^{(2)}\left(\omega \right)$ and ${\chi}_{312}^{(2)}\left(\omega \right)$. 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 ${\chi}_{ijk}^{(2)}\left(\omega \right)$. That is attributed to the fact that local density approximation calculations underestimate the energy gaps. Moreover, it is clear from the formulae of ${\chi}_{ijk}^{(2)}\left(\omega \right)$ that, ${\chi}_{ijk}^{(2)}\left(\omega \right)$ depends on the energy gap. Hence when we used the scissors correction we found a considerable effect on ${\chi}_{ijk}^{(2)}\left(\omega \right)$. 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 ${\chi}_{ijk}^{(2)}\left(\omega \right)$ involves 2*ω* resonance in addition to the usual *ω* resonance. Both *ω* and 2*ω* resonances can be further separated into inter-band and intra-band contributions.

*ω*inter-band and intra-band contributions for CuInX

_{2}compounds. We note the opposite signs of the two contributions throughout the frequency range. Both these contributions and the total Im ${\chi}_{123}^{(2)}\left(\omega \right)$ increase on moving from S to Se and decreases on moving from Se to Te. As we note from Figure 1, CuInSe

_{2}and CuInTe

_{2}has a reduction of the energy gap relative to CuInS

_{2}. The bandwidth of main groups in band structure and DOS of CuInSe

_{2}and CuInTe

_{2}increases over the value in CuInS

_{2}resulting in a reduction of the gap. Im ${\chi}_{123}^{(2)}\left(\omega \right)$ inCuInS

_{2}is smaller than in CuInSe

_{2}and CuInTe

_{2}. From this we conclude that the region around the conduction band minimum does not make a significant contribution to ${\chi}_{123}^{(2)}\left(\omega \right)$. Also the VBM and CBM in CuInSe

_{2}and CuInTe

_{2}seem to be more parallel than in CuInS

_{2}. 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 ${\chi}_{123}^{(2)}\left(\omega \right)$ cannot be based only on the minimum band gap. This indicates that the average band gap plays a more significant role.

*ω*and 2

*ω*resonances for CuInTe

_{2}. 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

*E*

_{ g }and the

*ω*terms for energy values above

*E*

_{ g }. In the low energy regime (≤ eV) the SHG optical spectra is dominated by 2

*ω*contributions. Beyond 3 eV the major contribution comes from

*ω*terms.

The structures in Im ${\chi}_{123}^{(2)}\left(\omega \right)$ 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 ${\chi}_{123}^{(2)}\left(\omega \right)$ between 1.0 – 3.0 eV for CuInS_{2}, 0.5 – 2.0 eV for CuInSe_{2}, and 1.0 – 2.5 eV for CuInTe_{2} 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 CuInS_{2}, 2.0 – 3.0 eV for CuInSe_{2}, and 2.5 – 4.5 eV for CuInTe_{2} is due mainly to *ω* resonance and associated with high peak in *ε*_{2}(*ω*). The last structure from 4.5 – 5.5 for CuInS_{2}, 3.0 – 5.5 eV for CuInSe_{2}, and 4.5 – 5.5 eV for CuInTe_{2} is manly due to *ω* resonance and associated with the tail in *ε*_{2}(*ω*).

_{2}with the available experimental data (Table 2) we found that our result appears to be very close to the experimental data [36] whereas the other theoretical value [17] appears to be more than twice the experimental data. This is attributed to our use of the FP. Thus from Table 2, we can conclude that the smaller band gap compounds gives higher values of ${\chi}_{123}^{(2)}\left(0\right)$ in agreement with the experiment [36, 51] and theory [17]. The lack of experimental data prevents any conclusive comparison with experiment over a large energy range. It is well known that [60, 61]${\chi}_{ijk}^{(2)}\left(\omega \right)$ = 2

*d*

_{ ij }, based on this expression we compared our calculated ${\chi}_{312}^{(2)}\left(0\right)$ with the values of

*d*

_{ 36 }which was obtained by D. Xue et al. [31] at the wavelength 10.6

*μ*m. Good agreement is found especially both calculations show the negative sign and same trends for both CuAlSe

_{2}and CuAlSe

_{2}, and positive sign for the CuAlTe

_{2}(Table 2).

Experimental, calculated total and intra-inter- band of the zero frequency of the real part of the ${\chi}_{123}^{(2)}\left(\omega \right)$. *χ*^{123}(0) is expressed in units of 1 × 10^{-8} esu.

CuInS | CuInSe | CuInTe | |
---|---|---|---|

| 14 | ||

| 31.7 | 72.5 | 126.0 |

| -51.8 | -115.0 | -187.2 |

| 83.5 | 187.5 | 313.2 |

| -45.3* | -28.2* | 8.5* |

d | -22.48 | -14.51 | 3.85* |

| -44.96** | -29.02** | 7.7** |

## 4. Conclusion

We have performed calculations of band structure, DOS, frequency dependent linear, birefringence and nonlinear optical response for CuInX_{2} (X = S, Se, Te) compounds, using FP-LAPW method. Our results for band structure and DOS show that these compounds have similar structures with direct energy band gap. These energy band gap changes when S is replaced by Se/Te. This is attributed to the fact that the bandwidth of conduction bands increases on going from S to Se and decreases on going from Se to Te in agreement with the experimental data. Our calculated energy band gaps show better agreement with the experiment than that obtained by using LMTO method [17]. This is attributed to the use of a FP method. As a result of this the LMTO method is unable to predict correctly the optical transitions, even when the scissor correction is applied. All structures in the imaginary part of the frequency dependent dielectric function *ε*_{2}(*ω*) are shifted towards lower energies when S is replaced by Se and Te. We have calculated *ε*_{2}(*ω*) and find a considerable anisotropy between ${\epsilon}_{2}^{\perp}\left(\omega \right)$ and ${\epsilon}_{2}^{II}\left(\omega \right)$. We find that the values of *ε*_{1}(0) increases with decreasing energy gap. This could be explained on the basis of the Penn model. We have calculated the birefringence of these compounds and find it is negative for CuInS_{2} and CuInSe_{2} while it is positive for CuInTe_{2} in agreement with the experimental data. Our calculations of SHG susceptibility show that the intra-band and inter-band contributions are significantly increased when S is replaced by Se/Te. Our calculations show that smaller band gap materials have higher ${\chi}_{123}^{(2)}\left(0\right)$ values. Since ${\chi}_{123}^{(2)}\left(0\right)$ is roughly inversely proportional to the band gap, one might think that this would lead to a possible route to further enhancement of Im ${\chi}_{123}^{(2)}\left(\omega \right)$. The enhancement of SHG is considerable when S is replaced by Se/Te. We have compared our calculated linear and nonlinear optical properties with the previous calculations of D. Xue et al. [31]. Good agreement is found.

## Declarations

### Acknowledgements

The author would like to thank the Institute Computer Center for providing the computational facilities. Also to thank Prof. Claudia Ambrosch-Draxl, 000 Montanuniversität Leoben, 8700 Leoben, Franz-Josef-StraBe 18, Austria, and Dr. Sangeeta Sharma, Institut für Theoretische Physik, Freie Universität Berlin, Arnimallee 14, D-14195 Berlin, Germany. 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

- Deb SK, Zunger A: Ternary and Multinary compounds. 1987, Pittsburgh, Pennsylvania: Materials Research SocietyGoogle Scholar
- Jaffe GE, Zunger A: Phys Rev B. 1983, 28: 5822-1984, 29: 1882View ArticleADSGoogle Scholar
- Rincon C, Bellabarba C: Phys Rev B. 1986, 33: 7160-View ArticleADSGoogle Scholar
- Parthe E: Crystal Chemistry of Tetrahedral Structures. 1964, New York: Gordon and BreachGoogle Scholar
- Goryunova NA: The Chemistry of Diamond-like semiconductors. 1965, New York: Chapman and HallGoogle Scholar
- Shay JL, Wernick JH: Ternury Chalcopyrite Semiconductors. Growth, Electronic properties and Applications. 1974, Oxford: PergamonGoogle Scholar
- Kaufmann U, Schneider J: Festkorperprobleme XIV. Edited by: Treusch J. 1974, Braunschweig: Vieweg, 229-View ArticleGoogle Scholar
- Wagner J: Electroluminescence. Edited by: Pankove JO. 1977, Berlin: Springer, 171-View ArticleGoogle Scholar
- Mackinnon A: Festkorperprobleme XXI. Edited by: Treusch J. 1981, Dortmund: Vieweg, 149-View ArticleGoogle Scholar
- Miller A, Mackinnon A, Weaire D: solid State Physics. Edited by: Ehrenreich H, Seitz F, Turubull D. 1981, New York: Academic, 36:Google Scholar
- Pamplin BR, Kiyosawa T, Mastumoto K: Prog Cryst Growth Charact. 1979, 1: 331-View ArticleGoogle Scholar
- Kazmerski LL: Nuovo Cimento, D. 1983, 2: 2013-View ArticleADSGoogle Scholar
- Shay JL, Schiavone LM, Buehier E, Wernick JH: J Appl Phys. 1972, 43: 2805-Wagner S, Shay JL, Tell B, Kasper HM: Appl Phys Lett 1973, 22: 351View ArticleADSGoogle Scholar
- Levine BF: Phys Rev B. 1973, 7: 2600-and references thereinView ArticleADSGoogle Scholar
- Hopkins FK: Laser Focus World. 1995, 31: 87-Google Scholar
- Catella GC, Burlage D: MRS Bull. 1998, 23: 28-View ArticleGoogle Scholar
- Rashkeev SN, Lambrecht WRL: Phys Rev B. 2001, 63: 165212-View ArticleADSGoogle Scholar
- Chemla DS, Kupecek PJ, Robertson DS, Smith RC: Opt Commun. 1971, 3: 29-View ArticleADSGoogle Scholar
- Boyd GD, Kasper H, McFee JH: IEEE J Quantum Electronics. 1971, QE-7: 563-View ArticleADSGoogle Scholar
- Alonso MI, Wakita K, Pascual J, Garriga M, Yamamoto N: Phys Rev B. 2001, 63: 075203-View ArticleADSGoogle Scholar
- Garbato L, Rucci A: Philos Mag. 1977, 35: 1685-View ArticleADSGoogle Scholar
- Neumann H: Cryst Res Technol. 1989, 24: 815-View ArticleGoogle Scholar
- Rincon C, Villareal I, Galindo H: J Appl Phys. 1999, 86: 2355-View ArticleADSGoogle Scholar
- Choi I, Yu PY: Phys Status Solidi B. 1999, 211: 143-View ArticleADSGoogle Scholar
- Roa L, Chervin JC, Itie JP, Polian A, Gauthier M, Chevy A: Phys Status Solidi B. 1999, 211: 455-Status Solidi B 1999, 213:,349View ArticleADSGoogle Scholar
- Miller RC: Appl Phys Lett. 1964, 5: 17-View ArticleADSGoogle Scholar
- Flytzanis C: Nonlinear Optics. Edited by: Rabin H, Tang CL. 1975, New York: Academic Press, 1: 9-Google Scholar
- Chemla DS, Jerphagnon J: Handbook on Semiconductors. Edited by: Moss TS. 1980, Amsterdam: North-Holland, 2: 545-Google Scholar
- Chen C, Ye N, Lin J, Jiang J, Zeng W, Wu B: Adv Mater. 1999, 11: 1071-View ArticleGoogle Scholar
- Xue D, Zhang S: Physica B. 1999, 262: 78-View ArticleADSGoogle Scholar
- Xue D, Betzler K, Hesse H: Phys Rev B. 2000, 62: 13546-View ArticleADSGoogle Scholar
- Pearson WB: Handbook of Lattice Spacing and Structures of Metals and Alloys. 1967, Oxford: Pergamon Press, 2:Google Scholar
- Boyd GD, Kasper HM, McFee JH: IEEE J Quantum Electron. 1971, 7: 563-View ArticleADSGoogle Scholar
- Blaha P, Schwarz K, Madsen GKH, Kvasnicka D, Luitz J: WIEN2k, An augmented plane wave plus local orbitals program for calculating crystal properties. 2001, Vienna University of Technology, AustriaGoogle Scholar
- Nastos F, Olejnik B, Schwarz K, Sipe JE: Phys Rev B. 2005, 72: 045223-View ArticleADSGoogle Scholar
- Shay JL, Wernick JH: Ternary Chalcopyrite Semiconductors: Growth, Electronic properties, and applications. 1975, Oxford: Pergamon PressGoogle Scholar
- Poplavnoi AS, Polygalov YI: Izv Akad Nauk SSSR, Neorg Mater. 1971, 7: 1706-1711. [Inorg Mater (USSR) 1971, 7: 1527–1531]Google Scholar
- Tell B, Shay JL, Kasper HM: Phys Rev B. 1971, 4: 2463-View ArticleADSGoogle Scholar
- Shay L, Tell B, Kasper HM, Schiavone LM: Phys Rev B. 1973, 7: 4485-View ArticleADSGoogle Scholar
- Shay L, Kasper HM: Phys Rev Lett. 1972, 29: 1162-View ArticleADSGoogle Scholar
- Gonzalez J, Rincon C: J Apply Phys. 1989, 65: 2031-View ArticleADSGoogle Scholar
- Twaites MJ, Tamlinson RD, Hampshire MJ: Solid State Commun. 1977, 23: 905-View ArticleADSGoogle Scholar
- Yamasaki T, Suzuki N, Motizuki K: J Phys C. 1987, 20: 395-View ArticleADSGoogle Scholar
- Reshak AH, Auluck S: Phys Rev B. 2003, 68: 245113-View ArticleADSGoogle Scholar
- Sharma S, Auluck S, Khan MA: Pramana Journal of Physics. 1999, 54: 431-View ArticleADSGoogle Scholar
- Penn DR: Phys Rev. 1962, 128: 2093-View ArticleADSMATHGoogle Scholar
- Wooten F: Optical properties of solids. 1972, New York and London; Academic pressGoogle Scholar
- Rashkeev SN, Lambrecht WRL, Segall B: Phys Rev B. 1998, 57: 3905-View ArticleADSGoogle Scholar
- Aspnes DE: Phys Rev B. 1972, 6: 4648-View ArticleADSGoogle Scholar
- Boyd W: Nonlinear optics. 1992, Boston: Academic PressGoogle Scholar
- Kurtz SK, Jerphagon J, Choy MM: In numerical data and functional relationships in science and technology. Landolt-Bornstein, new series. group III. Edited by: Hellwege KH, Hellwege AM. 1979, Berlin: Springer, 11: 671-Google Scholar
- Reshak AH, Auluck S: Physica B. 2005, 358: 158-View ArticleADSGoogle Scholar
- Eryigit R, Parlak C, Eryigit R: Eur Phys J B. 2003, 33: 251-View ArticleADSGoogle Scholar
- Shiwu Gao: Computer Physics Communications. 2003, 153: 190-View ArticleGoogle Scholar
- Karlheinz Schwarz: Journal of Solid State Chemistry. 2003, 176: 319-View ArticleGoogle Scholar
- Reshak AH: Eur Phys J B. 2005, 47: 503-View ArticleADSGoogle Scholar
- Boyd GD, Kasper H, McFee JH: IEEE J Quantum Electron. 1971, QE-7: 563-View ArticleADSGoogle Scholar
- Xue D, Betzler K, Hesse H, Lammers D: J Appl Phys. 2000, 87: 2849-View ArticleADSGoogle Scholar
- Mercier N, Leblanc M, Durand J: Eur J Solid State Inorg Chem. 1997, 34: 241-Google Scholar
- [http://www.rp-photonics.com/nonlinear_polarization.html]
- Boyd RY: Principles of Nonlinear Optics. 1982, NY: Acdemic Press, 420-Google Scholar
- Marquez R, Rincon C: Phys Status Solidi (b). 1995, 191: 115-View ArticleADSGoogle Scholar
- Jaffe JE, Zunger A: Phys Rev B. 1983, 28: 5822-View ArticleADSGoogle Scholar
- Shao LX, Chang KH, Chung TH, Tseng BH, Hwang HL: J Physics and Chemistry. 2003, 64: 1505-ADSGoogle Scholar
- Ramanathan K, Hasoon FS, Smith S, Young DL, Contreras MA, Johnson PK, Pudov AO, Sites JR: Journal of Physics and Chemistry of Solids. 2003, 64: 1495-View ArticleADSGoogle Scholar
- Lambrecht WRL, Rashkeev SN: Journal of Physics and Chemistry of Solids. 2003, 64: 1615-View ArticleADSGoogle Scholar
- Jepsen O, Andersen OK: Solid State Commun. 1971, 9: 1763-Lehmann G and Taut M, Phys Status Solidi B 1972, 54: 496View ArticleADSGoogle Scholar
- Wilson JA, Yoffe AD: Adv Phys. 1969, 18: 193-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.