Dispersion of Phonon Surface Polaritons in ZnGeP2: Anisotropy and Temperature Impacts

Zinc germanium diphosphide (ZnGeP2) is an attractive and promising functional material for different devices of the nano- and optoelectronics. In this paper, dispersion of phonon surface polaritons (PSPs) in ZnGeP2 has been studied in the 200–500-cm−1 spectral range at 4 and 300 K. Dispersion of “real” and “virtual” PSPs were calculated for C-axis being normal and parallel to the surface. Anisotropy in ZnGeP2 leads to the different numbers of PSP dispersion branches for different orientations of the sample. The temperature-dependent phonon contributions in the dielectric permittivity shift dispersion of the surface polaritons in ZnGeP2 to the higher wavenumbers at 4 K. We have shown that experimental dispersion of PSP is in agreement with theory.


Introduction
Among other ternary A II B IV C V 2 compounds, ZnGeP 2 stands out for its properties, such as great values of the coefficients of the nonlinear susceptibility and birefringence, which attract interest of researches.
Mechanical stability, resistance to moisture, aggressive environments, and good thermal conductivity characterize ZnGeP 2 as an attractive and promising functional material for different devices of the nanoand optoelectronics. In [1], it has been discussed that diphosphides can be used in infrared converter systems. Experimental results of [2] prove that replacing CdTe by diphosphides can increase efficiency of the solar converters and simplify their mass deployment.
Therefore, the study of optical properties of ZnGeP 2 is crucial for industrial applications of this material.
In the previous paper [3], we have reported on the lattice vibration behaviors of ZnGeP 2 single crystals derived from the far-infrared (FIR) reflectance spectra in the temperature range of 4-300 K. From [3], one can notice that optical phonon behaviors in ZnGeP 2 at low temperatures are similar to those in other diphosphides [4,5]. These results show that surface polaritons in ZnGeP 2 can be excited in the number of spectral ranges. However, for ZnGeP 2 , the surface polaritons practically have not been studied yet. Results of [6] have shown that temperature has an impact on the dispersion of surface phonon polaritons in binary diphosphide Zn 3 P 2 .
Here, we present the results of our study of the dispersion of the surface polaritons in ZnGeP 2 at different temperatures.

Methods
In the experiments, a set of samples of single crystals of ZnGeP 2 cut into plates of size 2 × 4 × 0.5 mm was used. The samples of ZnGeP 2 were oriented to (001) and (100) crystallographic planes. The orientation of the crystallographic planes was controlled by X-ray diffraction.
As an input data for the calculations of the dispersion of phonon surface polariton (PSP), we used reflectance data of ZnGeP 2 from our study [3]. These spectra were measured at 300 K in the 200-500-cm −1 range, using a Bruker IFS 66v/s spectrometer with Hg lamp as the source of radiation with a resolution of 0.5 cm −1 , 256 scans per 20 s employing polarized radiation were collected in each experiment. The temperature of sample varied in the 4-300-K temperature range. The angle of incidence of radiation was less than 10°. Spectra were measured at two orientations of the electrical vector E of the IR radiation with respect to the crystal: E||c and E˩c.
The method of attenuated total reflection (ATR) was firstly used for the excitation of surface waves in metals [7]. Since general principles of excitation of PSP by ATR have been described in [8], this method remains the easiest and the most reliable for excitation of surface polaritons.
We obtained the ATR spectra of PSP for ZnGeP 2 in the 300-500-cm −1 frequency range using polarized radiation. To measure these data, we used FTIR spectrometer Bruker IFS 66 with accessory with grazing angle of the incidence equipped by Ge hemisphere as ATR element. The air gap between the investigated sample and the hemisphere varied from 1.5 to 2 μ.

Results and Discussion
In the tetragonal ZnGeP 2 , B 2 (can be observed employing the radiation polarized parallel to the c-axis) and E (the radiation polarized normal to the c-axis) modes are infrared active. According to the data from [3], in the 200-500-cm −1 range, one can observe two B 2 and four E modes. Spectral analysis of the reflectance data has been performed employing factorized model for the dielectric permittivity in the following form: where ε ∞ is high-frequency dielectric permittivity and ν jl and γ jl (ν jt and γ jt ) are the frequency and damping coefficient of the jth longitudinal (transversal) mode. Dispersion of the real part of the dielectric permittivity of ZnGeP 2 obtained using model (1) at 4 and 300 K is presented in Fig. 1.
Results of spectral analysis of the reflectance data of ZnGeP 2 at 4 and 300 K are shown in Table 1. Observed temperature shift of the phonon modes is in agreement with theoretical calculations from [9] and has been also observed in the reflectance spectra of Al 2 O 3 [10] and Raman spectra of Si, Ge, and α-Sn [11].
Before discussing the obtained results, we make a short review of influence of the anisotropy on the dispersion relation of PSP. The coordinate axes x and y lay in the sample's plane, and z is normal to the surface. Let us denote the component of the dielectric tensor which is normal (parallel) to the optical axis of Fig. 1 Real part of the dielectric function of ZnGeP 2 single crystal. Ε||c: 1 4 K, 2 300 K; Ε˩c: 3 4 K, 4 300 K the crystal as ε ⊥ (ε || ). Neglecting the damping, we have analyzed this dielectric function applying conditions described in [8] and found spectral ranges of existence of the "real" (type 1) and "virtual" (type 2) PSP in ZnGeP 2 at 4 and 300 K, which are shown in Table 2.
In the case of uniaxial media, according to the terminology of [6], "ordinary" PSP occurs when the optical axis C is parallel to the surface (C||y). The dispersion relation in this case (ε x = ε z = ε ⊥ ) is as follows: This dispersion relation involves only ε ⊥ , and it is of the same form as for PSP dispersion in isotropic media.
"Extraordinary" PSP occur in the two other configurations. The optical axis is perpendicular to the surface (C||z). The dispersion relation can be written down as follows: where ε z = ε || , ε x = ε ⊥ . The optical axis is parallel to the surface of the sample and is parallel to the direction of propagating of PSP (C||x) with ε x = ε || , ε z = ε ⊥ .
Experimental evaluation of the dispersion of surface polaritons has been performed applying the method of ATR, which was used for excitation of the surface waves in ZnGeP 2 for C||x and C||y orientation of the sample at several angles of the incidence. Figure 2 presents experimental ATR spectra recorded at 30°. In the presence of the SP damping and dissipation of the electromagnetic wave energy, spectra show certain number of minima, whose wavenumbers correspond to the PSP. In Fig. 3, we show three experimental and five calculated ATR spectra (C||y) as the socalled ATR surface R(ν,α), which is a threedimensional presentation of the system transmission that depends on the radiation frequency ν and angle α. This ATR surface R(ν,α) = I(ν,α)/I 0 (ν,α) has three "canyons" connected with a "pass." I(ν, α) is the intensity of radiation passing through the "ATR    element-gap-sample" system; I 0 (ν,α) is the intensity of the incident radiation onto the ATR unit. The depth of the "canyon" depends on the following system parameters: gap d between the ATR element and sample, radiation frequency ν, complex permittivity ε(ν,k) of the sample, and the refractive indexes of the ATR unit and gap. The surface polariton dispersion curves ν s (k) correspond to the "canyons", i.e., to the set of ATR spectra minima. (Here, ν s is the surface polariton frequency and k is the SP wave vector). Thus, we can distinguish three branches of PSP dispersion in ZnGeP 2 within the measured spectral range. From the minima of ATR spectra, the PSP dispersion has been evaluated: where ν is the frequency of the ATR spectrum minimum, c is the speed of light in a vacuum, and n is the refractive index of material of the ATR element.
We have calculated dispersion of PSP in ZnGeP 2 for three possible configurations applying Eqs. 2-4 in the spectral ranges shown in Table 2 at 4 and 300 K. We present the obtained PSP dispersion branches in Figs. 4, 5, and 6. To emphasize the influence of the temperature, we plot 4 K PSP dispersion curves in blue and 300 K in orange.
For configuration C||x, we obtained one branch of the dispersion of the "real" PSP and two branches of the dispersion of the "virtual" PSP. The "real" PSP branch goes to the infinity wave vector within the specified spectral range. The "virtual" dispersion branches terminate at the cross points of the polariton dispersion with the dispersions of the longitudinal optical phonons: 1. 412 cm −1 , 4 K; 409 cm −1 , 300 K; 2. 362 cm −1 , 4 K; 361 cm −1 , 300 K.
As already mentioned, for C||y, the dispersion relation is similar to the isotropic case, since only ε ⊥ is involved   in Eq. (2). Therefore, the quantity of dispersion branches is equal to the quantity of the Reststrahlen bands for E c. We obtained three dispersion branches of "real" PSP for C||y orientation.
Due to the temperature shift of the phonon frequencies of ZnGeP 2 to the higher wave numbers [3], almost all branches of the 4 K dispersion lay above the corresponding branches at 300 K. The temperature-dependent phonon contributions in the dielectric permittivity, which are described by frequency, damping, and strength in frames of Lorentz model [3], affect dispersion of the surface polaritons in ZnGeP 2 .

Conclusions
Thus, in this paper, dispersion of PSP in ZnGeP 2 has been studied in the 200-500-cm −1 spectral range at 4 and 300 K. Dispersion of "real" and "virtual" PSP were calculated for C-axis being normal and parallel to the surface. Anisotropy in ZnGeP 2 leads to the different numbers of PSP dispersion branches for different orientations of the sample. The temperature-dependent phonon contributions in the dielectric permittivity shift dispersion of the surface polaritons in ZnGeP 2 to the higher wave numbers at 4 K. Experimental dispersion of PSP is in agreement with theory.