The optical properties of a single-phase corundum-structured In2O3 epilayer grown by a mist chemical vapor deposition method have been studied. Raman scattering measurements on a c face and on a lateral face reveal all of the seven Raman-active modes of the corundum structure, with good adherence to the Raman selection rules. Three out of the four infrared-active modes are observed in the spectroscopic ellipsometry measurements. The phonon frequencies obtained from Raman and ellipsometry measurements are in excellent agreement with density functional perturbation theory calculations. No trace of bixbyite phase was detected in the spectra. In the UV region, the imaginary part of the dielectric function shows two distinct onsets of strong absorption associated with direct band-to-band transitions at 3.38 and 3.86 eV.

Over the last decade, In2O3 and other semiconducting oxides such as Ga2O3 have been attracting a growing interest as technologically relevant wide gap materials.1,2 When doped with tin, In2O3 is a transparent conductive oxide exhibiting high electron mobilities with applications in transparent contacts for solar cells,3 thin film transistors,4,5 and light-emitting diodes and lasers.6 While the stable bixbyite structure of In2O3 has been extensively investigated,1,7–10 studies of a metastable rhombohedral polymorph α-In2O3 are scarce and preliminary.11 

Bandgap engineering in the α-In2O3α-Ga2O3α-Al2O3 system opens the possibility of tuning the bandgap from 3.8 to 8.8 eV.12 A corundum-structured Ga2O3 on the c-plane of a sapphire substrate was obtained by a mist chemical vapor deposition (MCVD).13 Later, growth of α-In2O3 by MCVD on sapphire substrates with Fe2O3 buffer layers was reported,14 and more recently, α-In2O3 directly grown on sapphire substrate by MCVD has been achieved.15,16 The vibrational and dielectric properties of high-quality single-crystal α-Ga2O3 films grown by MCVD were thoroughly investigated,17–21 providing insight into the band structure in the form of effective masses.22 A similarly good knowledge of the optical properties of the corundum-structured In2O3 phase is required to advance toward novel devices based on this III-oxide system.

Whereas the Raman-active phonons of bixbyite In2O3 have been extensively investigated and their symmetry assignment has been performed by polarized Raman scattering,7 most Raman spectra of α-In2O3 reported in the literature have been obtained from nanocrystalline powders with limited crystallinity,23,24 which precluded the symmetry analysis of the Raman modes. An unpolarized Raman spectrum of an α-In2O3 thin film grown by metal organic chemical vapor deposition (MOCVD) showing a residual bixbyite component was reported,11 but the observed modes were assigned in reference to the modes of the rhombohedric BiFeO3 crystals, which belongs to a different space group (R3c) than α-In2O3 (R3¯c). On the other hand, the very nature of the bandgap of In2O3 has been controversial owing to its strong optical absorption in the ultraviolet at 3.8 eV and weaker absorption below 3 eV, which early interpretations ascribed to an indirect gap. More recent studies have attributed the gradual absorption onset to dipole-forbidden or low-dipole-intensity direct transitions in bixbyite In2O3 and to low density of states and nonparabolicity of the conduction band in the case of the corundum-structured In2O3.25,26 Some dispersion of the gap values reported in the literature exists owing to the characteristics of the samples and the methods of analysis.26,27

In this Letter, we present a study of the optical properties of a high-quality, single-phase α-In2O3 epilayer by means of Raman scattering and spectroscopic ellipsometry measurements. An epilayer of thickness around 800 nm was grown on a c-plane sapphire substrate by MCVD. Mist particles were formed by a 2.4 MHz ultrasonic transducer and carried by a flow of O2 gas to the reaction chamber heated at 550 °C, where the (0001) sapphire substrate was placed on a sample holder. The details of the growth are described in Ref. 16. X-ray characterization of the sample revealed a high crystalline quality, showing 2θ/θ scans with a distinct (0006) α-In2O3 peak and no trace of bixbyite reflections. Hall effect experiments yielded a free electron concentration of 7.6×1018 cm−3 and a mobility of 151cm2/Vs. All characterization experiments were carried out at room temperature.

Raman scattering measurements were carried out using a Jobin-Yvon T64000 spectrometer equipped with a LN2-cooled charge coupled detector (CCD). Confocal micro-Raman measurements with a pinhole aperture of 100 μm at the back focal plane were performed to minimize the signal from the sapphire substrate. The Raman spectra were excited with the 532 nm line of a Nd:YAG laser. The light was collected in backscattering configuration through an 80× objective with a numerical aperture of 0.75, and the polarized Raman spectra were recorded both on the c-face and on a lateral fracture face with a spectral bandwidth of 2.2 cm−1. The frequency and symmetry of the α-In2O3 modes were calculated using density functional perturbation theory (DFPT) as implemented in the abinit package.28 The calculations were made in the local density approximation (LDA) with Perdew–Wang exchange correlation, a 6 × 6 × 6 k-point grid, and a plane wave basis set with an energy cutoff of 70 Ha. The frequencies of the phonon modes and their symmetry are listed in Table I.

TABLE I.

Frequencies of the optically active modes of α-In2O3 as determined from Raman and spectroscopic ellipsometry measurements, compared with DFPT calculations.

Raman activeInfrared active
ω (cm–1)ω (cm–1)
SymmetryExperimentalTheoreticalSymmetryExperimentalTheoretical
A1g 163.1 161.0 Eu  170.5 
Eg 177.9 174.3 A1u  265.2 
Eg 219.2 219.3 Eu 299.5 295.3 
Eg 271.5 269.3 Eu 432.6 429.9 
Eg 384.7 383.6 Eu 490.8 483.4 
A1g 500.6 498.4 A1u  488.8 
Eg 591.3 585.2    
Raman activeInfrared active
ω (cm–1)ω (cm–1)
SymmetryExperimentalTheoreticalSymmetryExperimentalTheoretical
A1g 163.1 161.0 Eu  170.5 
Eg 177.9 174.3 A1u  265.2 
Eg 219.2 219.3 Eu 299.5 295.3 
Eg 271.5 269.3 Eu 432.6 429.9 
Eg 384.7 383.6 Eu 490.8 483.4 
A1g 500.6 498.4 A1u  488.8 
Eg 591.3 585.2    

Spectroscopic ellipsometry was performed using two different instruments spanning together a spectral range from 250 cm–1 (30 meV) up to 6.5 eV. The spectral resolution of the Fourier-transform infrared ellipsometer was set to 4 cm–1 while for the visible and ultraviolet (UV) spectral ranges, an instrument equipped with an autoretarder was used. The spectral ranges of both ellipsometers overlap. Experiments were performed at three different angles of incidence Φ of 50°, 60°, and 70°.

Figure 1(a) shows the Raman spectra of the α-In2O3 epilayer in different backscattering configurations, where z is along the c axis and x is taken to be along the (101¯0) direction.17 All spectra are dominated by the strong A1g mode at 163.1 cm–1, which mainly involves In atoms motion along the c axis. The intensity ratio of this low-frequency A1g mode to all other modes is much higher than in α-Ga2O3,17 suggesting a significantly higher In–O bond polarizability. The A1g modes exhibit stronger polarizability than the Eg modes, typically more than one order of magnitude higher, resulting in the Raman spectra being dominated by the A1g modes. This is in good agreement with the Raman intensities derived from the calculated polarizabilities, which are shown in Fig. 1(b) for the z(xx)z¯ configuration. Although the polarizability of the high-frequency A1g mode, which predominantly involves motion of the O atoms perpendicularly to the c axis, is notably higher than that of the low-frequency A1g mode, its Raman intensity is lower because the Raman scattering cross section is inversely proportional to the mode frequency,29 and the occupation factor is also lower. Furthermore, the low-frequency A1g Raman peak exhibits a very narrow linewidth, which reflects the long lifetime of this mode as a consequence of the lack of anharmonic decay channels into lower energy phonons. In contrast, the linewidth of the high-frequency A1g mode is 17 cm−1 and the spectral power of the mode is distributed over a wider spectral range.

FIG. 1.

(a) Raman spectra of the corundum-structured α-In2O3 for different backscattering configurations. The intensities of the spectra have been normalized to the dominant A1g mode. The dotted line indicates the phonon overtone distribution (POD) according to the calculated phonon density of states. The residual Raman peaks from the sapphire substrate are marked by asterisks. (b) Integrated intensities of the z(xx)z¯ Raman peaks relative to the low-frequency A1g intensity (red bars) compared with the Raman intensities derived from the calculated Raman polarizabilities (green bars).

FIG. 1.

(a) Raman spectra of the corundum-structured α-In2O3 for different backscattering configurations. The intensities of the spectra have been normalized to the dominant A1g mode. The dotted line indicates the phonon overtone distribution (POD) according to the calculated phonon density of states. The residual Raman peaks from the sapphire substrate are marked by asterisks. (b) Integrated intensities of the z(xx)z¯ Raman peaks relative to the low-frequency A1g intensity (red bars) compared with the Raman intensities derived from the calculated Raman polarizabilities (green bars).

Close modal

It should be noted that the diagonal component of the polarizability tensor for the Eg (219.2 cm–1) mode is exceedingly small, which makes this mode barely visible in parallel polarization spectra. In contrast, the off diagonal component is comparable to those of the other modes and, consequently, the Eg (219.2 cm–1) mode can be distinctly observed in the x(yz)x¯ spectrum. The weak signal from the forbidden A1g modes detected in the crossed polarization spectra may arise from leakage from other polarizations in slightly misaligned samples, particularly in the case of measurements on the lateral face, where fine control of the film orientation is less precise.

The parallel polarization spectra display weak broad bands that are especially noticeable in the 150–300 cm–1 range. These bands are attributed to second-order scattering. Overtone scattering can be estimated by considering the phonon density of states with a doubled frequency axis and weighted by the statistical factor (1+nBE)2, where nBE is the Bose–Einstein distribution function.30 Using a smoothed phonon density of states derived from our DFPT calculations, we obtain the phonon overtone distribution (POD) shown as a dotted line in Fig. 1. The calculated profile exhibits a prominent band in the low frequency range that is in reasonable agreement with the background signal of the Raman spectra.

The phonon dispersion and density of states obtained from our DFPT calculations are shown in Fig. 2. As can be seen from the plots of the projected PDOS, phonons below 250 cm–1 correspond predominantly to In motion, whereas above 250 cm–1, the vibrations mostly involve O motion, without any significant phonon gap between both kinds of modes. There is no phonon gap between acoustic and optical modes either, as low-frequency optical modes dominated by In motion cross the acoustic phonon branches and form nondispersive bands. This may contribute to enhance phonon–phonon scattering and limit the thermal conductivity.31 These nondispersive bands also play a major role in the second order scattering observed in the spectra [see Fig. 1(a)].

FIG. 2.

Phonon dispersion of α-In2O3 along high symmetry lines. Raman (IR) active modes at zone center are labeled according to their respective symmetry in green (red) letters. The phonon density of states and the projected density of states on In (red shaded area) and O atoms (blue shaded area) are also displayed.

FIG. 2.

Phonon dispersion of α-In2O3 along high symmetry lines. Raman (IR) active modes at zone center are labeled according to their respective symmetry in green (red) letters. The phonon density of states and the projected density of states on In (red shaded area) and O atoms (blue shaded area) are also displayed.

Close modal

Spectroscopic ellipsometry yields the ellipsometric angles Ψ and Δ, corresponding to amplitude ratio and phase difference of the different polarization orientations of the reflected light. The experimental results were analyzed using a multi-layer model shaped to the sample structure to determine finally the complex dielectric function (DF). In this case, the layer model consists of α-Al2O3 to describe the substrate and on top, the α-In2O3 layer of interest. In the UV range, an additional effective medium layer that accounts for surface roughness using Bruggeman's32 formalism was added. As demonstrated previously for cubic and hexagonal InN33,34 as well as for bixbyite In2O3,8 this low refractive index layer also accounts for the electron accumulation at the surface of α-In2O3.26 From UV-ellipsometry, the thickness of the α-In2O3 film and the effective medium layer were determined to be 796 and 9.1 nm, respectively.

In the infrared, on top of a sapphire substrate model,19 the DF of the α-In2O3 layer was numerically fitted point-by-point (pbp) until best agreement with the experimental ellipsometric angles was achieved. This pbp DF was line shape fitted by an analytical model to determine optical parameters such as phonon wavenumbers in the infrared or transition energies in the UV. Due to the c-plane orientation of the sample, only the ordinary DF εord is accessible by spectroscopic ellipsometry, i.e., with the electric field vector perpendicular to the c-axis of the crystal.

The characteristic phonon energies are extracted from the pbp DF by assuming an analytic form of the DF containing the dielectric limit ε and the IR allowed Eu phonons of the corundum-structured crystal structure, using Lorentzian broadened phonon oscillators with the phonon frequency ωTO and the broadening parameter γTO. To account for the free-electron contribution, a Drude term is added to the DF:35 

(1)

It should be noted that out of the four infrared-active TO phonons (Eu) expected in the ordinary DF, only three are visible within our spectral range. According to our DFPT calculations, the fourth mode is expected at lower frequency (see Table I) and it should exhibit an extremely low oscillator strength.

Figure 3 depicts the real (blue) and imaginary (red) parts of the pbp DF together with the best fit of Eq. (1). The phonon frequencies of the Eu phonons detected are in good agreement with the theoretical values (see Table I). Furthermore, the observed intensities also agree qualitatively very well with the calculated oscillator strengths of the respective modes. Similar to α-Ga2O3 and α-Al2O3, the lowest and the highest wavenumber phonon have a very low intensity.22,44 The A1u phonons are only allowed in the extraordinary DF, and thus, they are not accessible experimentally in our sample. The fitted plasma frequency ωP=1413 cm−1 and the plasma broadening γP=406 cm−1 are averaged over the high-electron density accumulation surface layer and the lower carrier concentration bulk α-In2O3. Therefore, this result is not directly comparable to the Hall-effect experiment in such a case. The dielectric limit is found to be ε=3.9±0.1.

FIG. 3.

Point-by-point fitted complex dielectric function and line shape model in the spectral range of infrared-active optical phonons. Artifacts resulting from noise in the point-by-point fit are marked by an asterisk.

FIG. 3.

Point-by-point fitted complex dielectric function and line shape model in the spectral range of infrared-active optical phonons. Artifacts resulting from noise in the point-by-point fit are marked by an asterisk.

Close modal

For the UV spectral range, the agreement of a model DF to ellipsometric angles Ψ and Δ is shown in Fig. 4. The difference between the experimental data and the model in the range from 0.5 eV to ca. 3.5 eV is only of minor importance, since this is below the absorption edge and the optical properties are, therefore, not significantly affected. The imperfect match of the Fabry–Pérot oscillations leads to the layer thickness value becoming less precise. A possible explanation for the difference could be layer thickness inhomogeneities. The resulting real (blue) and imaginary parts (red) of the DF in the UV are shown in Fig. 5 (continuous curves). The used model function consists of a substrate model37 and a modified version of Elliott's model38 for the In2O3 layer, with the imaginary part

(2)

including two direct allowed interband transitions, very similar to the case of the ordinary DF in α-Ga2O3,36 which is also shown in Fig. 5 as a dashed gray curve for comparison. The real part is obtained by Kramers–Kronig transformation. α-In2O3 is a uniaxial material, so especially the region around the onset of strong absorption has to be considered carefully.39 Advanced ellipsometry data analysis frequently makes use of an analytical model to account for the extraordinary dielectric function on c-plane samples.40,41 In the current work, the line shape of the extraordinary dielectric function of α-In2O3 was chosen similar to the case of α-Ga2O3, but with the absorption onset lowered to the region of the ordinary DF of α-In2O3. While the agreement of model and data was improved, the overall result of the ordinary dielectric function only changed marginally. Thus, the same procedure was not necessary for the IR dielectric function. In the case of α-In2O3, the Elliott model yields two distinct absorption steps at characteristic energies of 3.38 and 3.86 eV, respectively. A number of theoretical and experimental values for the absorption onset reported in the literature lie in the region around 3–3.3 eV.25–27,42 Nevertheless, more recent experiments by Nishinaka and Yoshimoto43 found a value of 3.75 eV. Our result of two absorption steps provides a refinement of the current state of knowledge and fits very well to theoretical findings,25 where also two steps can be seen.

FIG. 4.

Ellipsometric angles Ψ and Δ and the corresponding fit in the ultraviolet range for three angles of incidence.

FIG. 4.

Ellipsometric angles Ψ and Δ and the corresponding fit in the ultraviolet range for three angles of incidence.

Close modal
FIG. 5.

Dielectric function of α-In2O3 in the spectral region around the onset of strong interband absorption. The imaginary part of the dielectric function is additionally decomposed into the contributions of different transitions (dotted lines). For comparison, the ε2 of the ordinary DF of Ga2O3 from Ref. 36 is reproduced as well (dashed gray curve).

FIG. 5.

Dielectric function of α-In2O3 in the spectral region around the onset of strong interband absorption. The imaginary part of the dielectric function is additionally decomposed into the contributions of different transitions (dotted lines). For comparison, the ε2 of the ordinary DF of Ga2O3 from Ref. 36 is reproduced as well (dashed gray curve).

Close modal

In summary, we have presented a comprehensive study of the optical properties of a single-phase α-In2O3 thin film grown by the MCVD method. The Raman spectra attest the high quality of the crystal and they do not show any trace of the stable bixbyite phase. All Raman active modes of the rhombohedric phase are observed and their symmetries identified by means of polarized Raman measurements. The infrared-active modes accessible from the spectroscopic ellipsometry measurements are determined from the pbp-fitted complex DF. The frequencies of both the Raman and infrared-active modes are in excellent agreement with DFPT calculations in the LDA approximation. In the UV region, the imaginary part of the DF exhibits two distinct onsets of absorption at 3.38 and 3.86 eV.

See the supplementary material for the ellipsometric angles Ψ and Δ measured in the IR spectral range at three different angles and the corresponding pbp fits. A comparison of pbp fits based on isotropic and anisotropic models is performed. The extraordinary DF used in the anisotropic model employed in the UV spectral range is plotted in the spectral region around the onset of strong interband absorption.

This work has been partially supported by the Generalitat Valenciana under Grant No. Prometeu/2021/066 and was funded in part by the Leibniz Science Campus GraFOx II. The authors would like to thank Mr. H. Yokoo of Kogakuin University for sample preparation.

The authors have no conflicts to disclose.

Ramon Cuscó: Formal analysis (equal); Investigation (equal); Writing – original draft (equal). Tomohiro Yamaguchi: Resources (equal); Writing – review and editing (equal). Elias Kluth: Formal analysis (equal); Investigation (equal); Writing – original draft (equal). Rüdiger Goldhahn: Supervision (equal); Writing – review and editing (equal). Martin Feneberg: Formal analysis (equal); Investigation (equal); Supervision (equal); Writing – original draft (equal).

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1.
O.
Bierwagen
,
Semicond. Sci. Technol.
30
,
024001
(
2015
).
2.
J. A.
Spencer
,
A. L.
Mock
,
A. G.
Jacobs
,
M.
Schubert
,
Y.
Zhang
, and
M. J.
Tadjer
,
Appl. Phys. Rev.
9
,
011315
(
2022
).
3.
J.-H.
Kim
,
H.-J.
Seok
,
H.-J.
Seo
,
T.-Y.
Seong
,
J. H.
Heo
,
S.-H.
Lim
,
K.-J.
Ahn
, and
H.-K.
Kim
,
Nanoscale
10
,
20587
(
2018
).
4.
Z.
Yuan
,
X.
Zhu
,
X.
Wang
,
X.
Cai
,
B.
Zhang
,
D.
Qiu
, and
H.
Wu
,
Thin Solid Films
519
,
3254
(
2011
).
5.
J. S.
Lee
,
Y.-J.
Kwack
, and
W.-S.
Choi
,
ACS Appl. Mater. Interfaces
5
,
11578
(
2013
).
6.
M. T.
Hardy
,
C. O.
Holder
,
D. F.
Feezell
,
S.
Nakamura
,
J. S.
Speck
,
D. A.
Cohen
, and
S. P.
DenBaars
,
Appl. Phys. Lett.
103
,
081103
(
2013
).
7.
C.
Kranert
,
R.
Schmidt-Grund
, and
M.
Grundmann
,
Phys. Status Solidi RRL
8
,
554
(
2014
).
8.
M.
Feneberg
,
J.
Nixdorf
,
C.
Lidig
,
R.
Goldhahn
,
Z.
Galazka
,
O.
Bierwagen
, and
J. S.
Speck
,
Phys. Rev. B
93
,
045203
(
2016
).
9.
D.
Kim
and
S.
Lee
,
Sci. Rep.
10
,
4957
(
2020
).
10.
K. O.
Egbo
,
A. E.
Adesina
,
C. V.
Ezeh
,
C. P.
Liu
, and
K. M.
Yu
,
Phys. Rev. Mater.
5
,
094603
(
2021
).
11.
C. Y.
Wang
,
Y.
Dai
,
J.
Pezoldt
,
B.
Lu
,
T.
Kups
,
V.
Cimalla
, and
O.
Ambacher
,
Cryst. Growth Des.
8
,
1257
(
2008
).
12.
S.
Fujita
and
K.
Kaneko
,
J. Cryst. Growth
401
,
588
(
2014
).
13.
D.
Shinohara
and
S.
Fujita
,
Jpn. J. Appl. Phys.
47
,
7311
(
2008
).
14.
N.
Suzuki
,
K.
Kaneko
, and
S.
Fujita
,
J. Cryst. Growth
364
,
30
(
2013
).
15.
T.
Yamaguchi
,
S.
Takahashi
,
T.
Kiguchi
,
A.
Sekiguchi
,
K.
Kaneko
,
S.
Fujita
,
H.
Nagai
,
M.
Sato
,
T.
Onuma
, and
T.
Honda
,
Appl. Phys. Express
13
,
075504
(
2020
).
16.
A.
Taguchi
,
S.
Takahashi
,
A.
Sekiguchi
,
K.
Kaneko
,
S.
Fujita
,
T.
Onuma
,
T.
Honda
, and
T.
Yamaguchi
,
Phys. Status Solidi B
259
,
2100414
(
2021
).
17.
R.
Cuscó
,
N.
Domènech-Amador
,
T.
Hatakeyama
,
T.
Yamaguchi
,
T.
Honda
, and
L.
Artús
,
J. Appl. Phys.
117
,
185706
(
2015
).
18.
A.
Segura
,
L.
Artús
,
R.
Cuscó
,
R.
Goldhahn
, and
M.
Feneberg
,
Phys. Rev. Mater.
1
,
024604
(
2017
).
19.
M.
Feneberg
,
J.
Nixdorf
,
M. D.
Neumann
,
N.
Esser
,
L.
Artús
,
R.
Cuscó
,
T.
Yamaguchi
, and
R.
Goldhahn
,
Phys. Rev. Mater.
2
,
044601
(
2018
).
20.
P.
Ning
,
J.
Grümbel
,
J.
Bläsing
,
R.
Goldhahn
,
D.-W.
Jeon
, and
M.
Feneberg
,
Semicond. Sci. Technol.
35
,
095001
(
2020
).
21.
J.
Grümbel
,
R.
Goldhahn
,
D.-W.
Jeon
, and
M.
Feneberg
,
Appl. Phys. Lett.
120
,
022104
(
2022
).
22.
M.
Feneberg
,
J.
Bläsing
,
T.
Sekiyama
,
K.
Ota
,
K.
Akaiwa
,
K.
Ichino
, and
R.
Goldhahn
,
Appl. Phys. Lett.
114
,
142102
(
2019
).
23.
B.
García-Domene
,
J. A.
Sans
,
F. J.
Manjón
,
S. V.
Ovsyannikov
,
L. S.
Dubrovinsky
,
D.
Martinez-Garcia
,
O.
Gomis
,
D.
Errandonea
,
H.
Moutaabbid
,
Y. L.
Godec
 et al.,
J. Phys. Chem. C
119
,
29076
(
2015
).
24.
J. A.
Sans
,
R.
Vilaplana
,
D.
Errandonea
,
V. P.
Cuenca-Gotor
,
B.
García-Domene
,
C.
Popescu
,
F. J.
Manjón
,
A.
Singhal
,
S. N.
Achary
,
D.
Martinez-Garcia
 et al.,
Nanotechnology
28
,
205701
(
2017
).
25.
F.
Fuchs
and
F.
Bechstedt
,
Phys. Rev. B
77
,
155107
(
2008
).
26.
P. D. C.
King
,
T. D.
Veal
,
F.
Fuchs
,
C. Y.
Wang
,
D. J.
Payne
,
A.
Bourlange
,
H.
Zhang
,
G. R.
Bell
,
V.
Cimalla
,
O.
Ambacher
 et al.,
Phys. Rev. B
79
,
205211
(
2009
).
27.
C. Y.
Wang
,
V.
Cimalla
,
H.
Romanus
,
T.
Kups
,
G.
Ecke
,
T.
Stauden
,
M.
Ali
,
V.
Lebedev
,
J.
Pezoldt
, and
O.
Ambacher
,
Appl. Phys. Lett.
89
,
011904
(
2006
).
28.
See http://www.abinit.org for ABINIT is a common project of the Université Catholique de Louvain, Corning Incorporated, and other contributors;
X.
Gonze
,
J.-M.
Beuken
,
R.
Caracas
,
F.
Detraux
,
M.
Fuchs
,
G.-M.
Rignanese
,
L.
Sindic
,
M.
Verstraete
,
G.
Zerah
,
F.
Jollet
 et al.,
Comput. Mater. Sci.
25
,
478
(
2002
).
29.
A.
Cantarero
,
J. Nanophotonics
7
,
071598
(
2013
).
30.
Light Scattering in Solids II
, Topics in Applied Physics Vol.
50
, edited by
M.
Cardona
and
G.
Güntherodt
(
Springer
,
Berlin, Heidelberg
,
1982
).
31.
Z.
Tian
,
J.
Garg
,
K.
Esfarjani
,
T.
Shiga
,
J.
Shiomi
, and
G.
Chen
,
Phys. Rev. B
85
,
184303
(
2012
).
32.
D. A. G.
Bruggeman
,
Ann. Phys.
421
,
160
(
1937
).
33.
S.
Schöche
,
T.
Hofmann
,
V.
Darakchieva
,
N.
Ben Sedrine
,
X.
Wang
,
A.
Yoshikawa
, and
M.
Schubert
,
J. Appl. Phys.
113
,
013502
(
2013
).
34.
P.
Schley
,
R.
Goldhahn
,
C.
Napierala
,
G.
Gobsch
,
J.
Schörmann
,
D. J.
As
,
K.
Lischka
,
M.
Feneberg
, and
K.
Thonke
,
Semicond. Sci. Technol.
23
,
055001
(
2008
).
36.
M.
Kracht
,
A.
Karg
,
M.
Feneberg
,
J.
Bläsing
,
J.
Schörmann
,
R.
Goldhahn
, and
M.
Eickhoff
,
Phys. Rev. Appl.
10
,
024047
(
2018
).
37.
I. H.
Malitson
,
J. Opt. Soc. Am.
52
,
1377
(
1962
).
38.
R. J.
Elliott
,
Phys. Rev.
108
,
1384
(
1957
).
39.
S.
Shokhovets
,
L.
Spieß
, and
G.
Gobsch
,
J. Appl. Phys.
107
,
023509
(
2010
).
40.
R.
Schmidt
,
B.
Rheinländer
,
M.
Schubert
,
D.
Spemann
,
T.
Butz
,
J.
Lenzner
,
E. M.
Kaidashev
,
M.
Lorenz
,
A.
Rahm
,
H. C.
Semmelhack
, and
M.
Grundmann
,
Appl. Phys. Lett.
82
,
2260
(
2003
).
41.
G.
Rossbach
,
M.
Feneberg
,
M.
Röppischer
,
C.
Werner
,
N.
Esser
,
C.
Cobet
,
T.
Meisch
,
K.
Thonke
,
A.
Dadgar
,
J.
Bläsing
,
A.
Krost
, and
R.
Goldhahn
,
Phys. Rev. B
83
,
195202
(
2011
).
42.
H.
Jiang
,
L.
Zhao
,
L.
Gai
,
L.
Ma
,
Y.
Ma
, and
M.
Li
,
CrystEngComm
15
,
7003
(
2013
).
43.
H.
Nishinaka
and
M.
Yoshimoto
,
Cryst. Growth Des.
18
,
4022
(
2018
).
44.
M.
Shubert
,
T. E.
Tiwald
, and
C. M.
Herzinger
,
Phys. Rev. B
61
,
8187
(
2000
).

Supplementary Material