Search tips
Search criteria 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2017; 7: 7360.
Published online 2017 August 4. doi:  10.1038/s41598-017-07164-1
PMCID: PMC5544691

Dielectric properties of semi-insulating Fe-doped InP in the terahertz spectral region


We report the values and the spectral dependence of the real and imaginary parts of the dielectric permittivity of semi-insulating Fe-doped InP crystalline wafers in the 2–700 cm−1 (0.06–21 THz) spectral region at room temperature. The data shows a number of absorption bands that are assigned to one- and two-phonon and impurity-related absorption processes. Unlike the previous studies of undoped or low-doped InP material, our data unveil the dielectric properties of InP that are not screened by strong free-carrier absorption and will be useful for designing a wide variety of InP-based electronic and photonic devices operating in the terahertz spectral range.


Semi-insulating (SI) iron-doped indium phosphide (InP:Fe) is widely used in electronic and photonic devices operating in the terahertz spectral range (THz range, 0.1~10 THz), including Schottky diode detectors1, high-electron mobility transistors2, photomixers3, and quantum cascade lasers (QCLs)46. Nominally updoped InP crystals always contain different unintentional impurities due to the growth processes with the concentrations up to 5·1015 cm−3 that results in shallow donor or acceptor energy levels within the energy gap. Iron doping provides acceptor levels in the mid-gap region of InP that compensate residual shallow donors and produce material with virtually no free carriers7. High-resistivity semiconductors are highly desired for devices operating in the THz spectral range as the optical losses at THz frequencies are typically dominated by the free-carrier absorption.

Given the importance of InP material to photonics and electronics, its optical properties have been extensively studied across the electromagnetic spectrum817. However, THz optical properties of SI InP have not been reported in the literature yet. The dielectric constants of nominally undoped and low-doped InP have been studied in the THz spectral range recently10, 11; however, they are strongly affected by a free carrier plasma that screens the intrinsic characteristics of SI InP material, especially at frequencies below 5 THz.

In this study, we present detailed investigation of the electrodynamic response of SI InP:Fe in the 2–700 cm−1 (0.06–21 THz) spectral region. The studied samples are obtained from commercial vendors and have resistivities exceeding 5 × 106 Ω·cm that corresponds to a free-carrier concentration below 109 cm−3. As expected for such low free carrier concentration, our optical data show no plasma effect in the entire spectral range of interest. Given the importance of SI InP to a variety of applications16, we believe that this report will be useful to a wide range of groups involved in the microwave and THz semiconductor devices research and development.

Experimental details


Two semi-insulating InP:Fe wafers with the nominal thicknesses of 350 ± 25 and 1000 ± 25 μm, obtained from two different vendors (AXT and Wafer Technology Ltd., respectively) were studied. Approximately 1 × 1 cm2 square sections were cut from the wafers and the thicknesses of the squares were measured at multiple locations to be 360 ± 2 μm and 991 ± 2 μm. The difference in the dielectric data obtained for the two samples was within the measurement error. Below we present and discuss the data obtained with the thicker sample.

Experimental setup and data processing

The measurements of the dielectric properties in the ν = 0.21–3.00 THz range (7–100 cm−1) were performed using a pulsed THz TeraView TPS-Spectra-3000 time-domain spectrometer. A spectrometer based on monochromatic and frequency-tunable continuous-wave (CW) backward-wave oscillators18 was used for the measurements in the 0.06–0.30 THz range (2–10 cm−1). Both spectrometers provide a possibility to derive the spectra of dielectric parameters - complex permittivity, dynamic conductivity, refractive index, etc., directly by measuring the values of the amplitude and the phase of the electric field of the radiation/electromagnetic wave passed through the plane-parallel sample. Additionally, we have performed high resolution (Δν < 0.3 cm−1) measurements of transmission and reflection coefficients of our samples in the 0.9–21 THz range (30–700 cm−1) using a standard vacuum Fourier-transform infrared spectrometer (FTIR) (Bruker Optics Vertex 80 v). For the data shown in the figures, we combined the data obtained with the backward-wave oscillators in the 2–10 cm−1 range with the data obtained by the TeraView system in the 10–100 cm−1 and the data obtained by the FTIR in the 100–700 cm−1 range. The results obtained by different measurement techniques in the overlapping spectral regions were the same within the experimental uncertainties.

Due to low absorption, the spectra of the transmission and reflection coefficients of our samples contain pronounced maxima and minima that arise due to the interference of the radiation within plane-parallel slabs (a well-known Fabry-Perot effect19, 20). These spectra allow for the most precise determination of the dielectric parameters of the samples at the frequencies of the interference maxima where the interaction of the probing radiation with the material is the most effective. This approach was used in the present research: the dielectric parameters of the studied InP:Fe samples were determined by modeling the measured transmission coefficient maxima based on the Fresnel expressions that describe the optical properties of a plane-parallel layers. At frequencies where the samples were not transparent due to strong phonon or impurity absorption resonances, the dielectric parameters of the samples were determined by simultaneous processing both transmission and reflection spectra. To get the parameters of each resonance, we have modeled them with the Lorentzian lineshapes:


Here εʹ(ν) = n 2(ν)  κ 2(ν) and ε″(ν) = 2 are the real and imaginary parts of the complex dielectric permittivity ε *(ν) = εʹ(ν) + iε″(ν), n is the real and κ is the imaginary part of the complex refractive index n* = n + iκ, f j = Δε j ν j 2 is the oscillator strength of the j-th resonance, Δε j is its dielectric contribution, ν j represents the resonance frequency, and γ j is the damping factor.

Results and Discussion

The combined transmission and reflection spectra are shown in Fig. 1 (dots), together with the spectral modeling results (lines). Below 150–200 cm−1 the interferometric Fabry-Perot oscillations are clearly resolved with the maxima in the transmission spectrum corresponding to the minima in reflectivity. The interferometric effect is also seen at frequencies above 500 cm−1 where the absorption of the crystal is not too strong. A number of absorption bands are observed as pronounced minima in the transmission spectrum above 100 cm−1 and weaker absorption features are seen at lower frequencies. The origin of these absorption peaks is discussed below.

Figure 1
Room-temperature transmission (a) and reflection (b) spectra of semi-insulated Fe-doped InP wafer of 991 µm thickness. Black dots are the experimental data, red solid lines present theoretical modeling described in the main text.

Presented in Fig. 2 are the broad-band spectra of the real and imaginary dielectric permittivities and of the absorption coefficient of SI InP:Fe. The absorption coefficient was determined as α=ωcIm(ε+iε), where ω is the circular radiation frequency and c is the speed of light. The data reveals a complex structure with a set of absorption lines with the most intensive one located at 304 cm−1 (see Table 1). The strongest absorption line shown separately in the inset is a well-known transverse optical (TO) phonon21. This mode is responsible for the low transmission values between ≈150 cm−1 and ≈500 cm−1 in Fig. 1a and for the corresponding characteristic dispersion in the reflectivity spectrum in Fig. 1b. There might be weaker absorption bands in the frequency interval of the reststrahlen band, 300–350 cm−1 that are, however, not resolved in the spectra due to the dominant effect of the TO phonon absorption. Less intensive absorption bands on the left and right wings of the phonon resonance are assigned to multi-phonon (summation or differential22) processes (indicated in Fig. 2b with vertical arrows) or to the electron/hole transitions involving impurities.

Figure 2
Room temperature terahertz-infrared spectra of the real (a) and imaginary (left vertical axis in panel (b)) parts of the dielectric permittivity and of the absorption coefficient α (right vertical axis in panel (b)) of Fe-doped InP. Dots are the ...
Table 1
Parameters of the absorption lines (Eq. (1) in the text) identified in the dielectric data from the SI InP crystals.

Nominally pure InP crystals contain different unintentional impurities due to the growth processes with the concentrations up to 5·1015 cm−3. The most common unintentional impurities in the SI InP are Si, S, Zn, C15, 17. In order to trap the free carriers produced by these impurities, InP crystals are doped with iron that creates acceptor levels in the mid-gap region of InP7. The activation energy of Fe impurity is 640 meV16 (frequencies above 5100 cm−1) and it is not expected to affect our measured spectra. The energies of the transitions related to the defects, i.e. indium and phosphorus vacancies and antisites are higher than 100 meV16 (frequencies above 800 cm−1), and are also not expected to show up in our spectra. However, the ionization energies of shallow donors (Si, S, Sn, Ge) are about 7.65 meV23, and of shallow acceptors (Zn, C, Si, Mn, Be, Mg) are about 25–40 meV13, 24. These transitions are likely to produce absorption features seen in our spectra.

The lowest-frequency absorption band at about 30 cm−1 in the spectra in Fig. 2 was previously observed in ref. 25. In accordance with the temperature dependence of the oscillator strength, this band was associated with the two-phonon absorption; however, the authors in ref. 23. did not assign the line to a specific phonon type and location. Koteles and Datars26 predicted an absorption line associated with the differential LO-TO phonon absorption that should be located either at the L or “hex” (see below) points of the Brillouin zone and that is expected to appear in the region 15–35 cm−1. The authors of ref. 26. also observed a number of two phonon summation and differential absorption lines which they attributed to phonons located at the X point, L point, and a point, located somewhere on the (111) hexagonal face of the Brillouin zone boundary, designated as “hex”.

The next absorption band at about 64 cm−1 could be a manifestation of differential phonon absorption LA-TA2 (hex)24 and/or an impurity transition. The energy of the residual donor impurities is reported to be about 7.65 meV (61.7 cm−1)23. The small difference between the observed location of the band and the literature data of the shallow donor impurity transition can be explained by the dependence of the energy gap of InP on the donor content and by the mutual influence of several impurities presented within its crystal lattice.

The two lines observed at 79 and 88 cm−1 could also be attributed to a differential phonon absorption. According to the data in ref. 26, the phonon positions are determined with an uncertainty of ±10 cm−1 which explains the difference between the peak positions obtained in this study and in the literature (see Table 1).

The line at 163 cm−1 is apparently some acceptor impurity contribution. It is unlikely that this line is a manifestation of an absorption due to a two-phonon or a compensated shallow donor process, since its dielectric contribution and oscillator strength are about an order of magnitude higher than those of the neighboring lines. According to the literature data24, Ge impurity has the ionization energy of 21 meV and could correspond to the line at about 163 cm−1.

At 235 cm−1 (29.1 meV), a relatively strong absorption band is registered. Due to its high intensity (oscillator strength), it is unlikely to be attributed to a multi-phonon absorption process. This line may be a manifestation of a shallow acceptor. Taking into account the position of the band, the most possible impurities could be Mn, Si, Be or Mg. We note that this line may also be produced by Zn contamination.

In Table 2 we present the values of the real and imaginary parts of the dielectric permittivity, conductivity, refractive index, extinction coefficient, loss tangent tan(δ) = ε″/ε′ and power absorption coefficient α at selected frequencies. The values are listed for various frequencies ν between 2 cm−1 and 100 cm−1. In Fig. 3 we show the same data for the real and imaginary permittivities, absorption coefficient and dynamical conductivity σ(ν) = νε″(ν)/2 (inset) in spectral form.

Table 2
Dielectric parameters of the measured SI InP:Fe crystals at selected frequencies.
Figure 3
Room temperature low-frequency spectra of the real εʹ (a) and imaginary ε″ (b) parts of the dielectric permittivity and a low-frequency spectrum of the absorption coefficient α (c) of Fe-doped InP. The dynamical ...

It is seen from Fig. 3 that, at THz frequencies below ≈15 cm−1, the measured absorption is larger than the absorption expected from impurities- or phonon-related processes (solid lines). The additional absorption at low THz frequencies should be related to the hopping conductivity of residual quasi-free charge carriers that is described by the Mott’s dispersion σ(ν) [proportional, variant] νs with s  127 and was previously observed in InP crystals (see, for example ref. 28). Dotted line in the inset of Fig. 3 shows that the low-frequency contribution to the conductivity is well reproduced by the Mott’s expression with s = 0.9.


In conclusion, we have performed room-temperature measurements of the dielectric properties of semi-insulating (SI) InP:Fe crystals in the 0.06–21 THz spectral range. Unlike the previous studies of undoped and low-doped InP material, our data unveil the dielectric properties of intrinsic InP that are not screened by the strong free-carrier absorption. A number of absorption resonances are discovered and their origin is analyzed. The values of the dielectric parameters of SI InP:Fe at frequencies between 2 and 700 cm−1 (0.06 and 21 THz) are presented. The data reported here is expected to be useful in designing and improving the performance of numerous microwave and terahertz semiconductor devices based on SI InP:Fe.


We acknowledge the assistance in the terahertz and infrared experiments of L. Kadyrov, M. Belyanchikov, Z. Bedran. The authors also express gratitude to Prof. Dan Botez and his group at the University of Wisconsin at Madison and Prof. Markus Amann and his group at the Technical University of Munich for providing the SI InP substrates for this study. The work was supported by the Russian Ministry of Education and Science (Program ‘5top100’). M.A.B. acknowledges financial support from the US National Science Foundation.

Author Contributions

Author Contributions

L.N.A. and E.S.Z. performed experiments, L.N.A. performed theoretical analysis and wrote the original manuscript with help of B.P.G., M.A.B. suggested the experiment and edited the manuscript, B.P.G. directed the project. All authors discussed results and reviewed the manuscript.


Competing Interests

The authors declare that they have no competing interests.


Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Contributor Information

L. N. Alyabyeva, ur.tpim@nl.avebaila.

M. A. Belkin, ude.saxetu.ece@niklebm.


1. Han SP, et al. InGaAs Schottky barrier diode array detector for a real-time compact terahertz line scanner. Opt. Express. 2013;21:25874–25882. doi: 10.1364/OE.21.025874. [PubMed] [Cross Ref]
2. Ostinelli O, Alt AR, Lövblom R, Bolognesi CR. Growth and characterization of iron-doped semi-insulating InP buffer layers for Al-free GaInP/GaInAs high electron mobility transistors. J.Appl.Phys. 2010;108:114502. doi: 10.1063/1.3516490. [Cross Ref]
3. Chimot N, et al. Photomixing at 1.55 μm in ion-irradiated In0.53Ga0.47As on InP. Opt. Express. 2006;14:1856–1861. doi: 10.1364/OE.14.001856. [PubMed] [Cross Ref]
4. Ajili L, Scalari G, Hoyler N, Giovannini M, Faist J. InGaAs–AlInAs/InP terahertz quantum cascade laser. Appl. Phys. Lett. 2005;87:141107. doi: 10.1063/1.2081122. [Cross Ref]
5. Vijayraghavan K, et al. Terahertz sources based on Čerenkov difference-frequency generation in quantum cascade lasers. Appl. Phys. Lett. 2012;100:251104. doi: 10.1063/1.4729042. [Cross Ref]
6. Belkin MA, Capasso F. New frontiers in quantum cascade lasers: high performance room temperature terahertz sources. Phys. Scr. 2015;90:118002. doi: 10.1088/0031-8949/90/11/118002. [Cross Ref]
7. Clerjaud B. Transition-metal impurities in III-V compounds. J. Phys. C. 1985;18:3615–3661. doi: 10.1088/0022-3719/18/19/005. [Cross Ref]
8. Joyce HJ, et al. Ultralow surface recombination velocity in InP nanowires probed by terahertz spectroscopy. Nano lett. 2012;12:5325–5330. doi: 10.1021/nl3026828. [PubMed] [Cross Ref]
9. Joyce HJ, et al. Electronic properties of GaAs, InAs and InP nanowires studied by terahertz spectroscopy. Nanotechnology. 2013;24:214006. doi: 10.1088/0957-4484/24/21/214006. [PubMed] [Cross Ref]
10. Zhang C, Jin B, Chen J, Wu P, Tonouchi M. Noncontact evaluation of nondoped InP wafers by terahertz time-domain spectroscopy. J. Opt. Soc. Am. B. 2009;26:A1–A5. doi: 10.1364/JOSAB.26.0000A1. [Cross Ref]
11. Xu X, Zhang L, Zhang C. Ultra-broadband terahertz spectroscopy of InP wafer using coherent heterodyne time-domain spectrometer. Proc. SPIE. 2015;9625:96250B-1–96250B-7. doi: 10.1117/12.2214351. [Cross Ref]
12. Hahn S, Dornich K, Hahm T, Kohler A, Niklas JR. Contact free defect investigation of wafer annealed Fe-doped SI-InP. Proc. Materials Science in Semiconductor. 2006;9:355–358. doi: 10.1016/j.mssp.2006.01.077. [Cross Ref]
13. Messner C, Kostner H, Höpfel RA, Unterrainer K. Time-resolved THz spectroscopy of proton-bombarded InP. J. Opt. Soc. Am. B. 2001;18:1369–1371. doi: 10.1364/JOSAB.18.001369. [Cross Ref]
14. Němec H, et al. Ultrafast carrier dynamics in Br+-bombarded InP studied by time-resolved terahertz spectroscopy. Phys. Rev. B. 2008;78:235206. doi: 10.1103/PhysRevB.78.235206. [Cross Ref]
15. Mayer KM, et al. Characterization of processing‐induced defects in high‐purity InP by photoluminescence. J. Appl. Phys. 1992;72:1080–1085. doi: 10.1063/1.351783. [Cross Ref]
16. Kaminski P, Kozlowski R, Strzelecka S, Piersa M. High-resolution photoinduced transient spectroscopy as a new tool for quality assessment of semi-insulating InP. J. Phys.- Condens. Mat. 2004;16:S225–S233. doi: 10.1088/0953-8984/16/2/027. [Cross Ref]
17. Palik ED. Handbook of Optical Constants of Solids. San Diego: Academic press.; 1998.
18. Gorshunov B, et al. Terahertz BWO-spectroscopy. Int. J. Infrared Milli. 2005;26:1217–1240. doi: 10.1007/s10762-005-7600-y. [Cross Ref]
19. Born M, Wolf E. Principles of Optics. Oxford: Pergamon; 1980.
20. Dressel, M. & Gruner, G. Electrodynamics of Solids (Cambridge University Press, Cambridge, 2002).
21. Lockwood DJ, Guolin Y, Rowell NL. Optical phonon frequencies and damping in AlAs, GaP, GaAs, InP and InSb studied by oblique incidence infrared spectroscopy. Solid State Commun. 2005;136:404–409. doi: 10.1016/j.ssc.2005.08.030. [Cross Ref]
22. Johnson EJ, Willardson R, Beer AC. Semiconductors and Semimetals. New York: Academic Press; 1967. Optical Properties of III–V Compounds in.
23. Kubota E, Katsui A, Sugii K. Photoluminescence study of residual shallow acceptors in liquid‐encapsulated Czochralski‐grown InP. J. Appl. Phys. 1986;59:3841–3846. doi: 10.1063/1.336725. [Cross Ref]
24. The Ioffe Institute electronic archive. New Semiconductor Materials. Characteristics and Properties (2015).
25. Gorshunov BP, Kalinushkin VP, Kozlov GV, Sirotinskij OI. Dielectric properties of compensated indium phosphide in submillimeter range. Fiz. Tverd. Tela. 1990;32:1530–1532.
26. Koteles ES, Datars WR. Two-phonon absorption in InP and GaP. Solid State Commun. 1976;19:221–225. doi: 10.1016/0038-1098(76)90854-1. [Cross Ref]
27. Mott NF, Davis EA. Electronic Processes in Noncrystalline Solids. London: Oxford University; 1971.
28. Mareš JJ, Štofik JK, ŠSmí V, Zema J. Ac properties of semi-insulating InP. Philos. Mag. B. 1991;64:233–243. doi: 10.1080/13642819108207616. [Cross Ref]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group