We present the growth and characterization of layered heterostructures comprised of LaTiO3 and SrTiO3 epitaxially grown on Si (001). Magnetotransport measurements show that the sheet carrier densities of the heterostructures scale with the number of LaTiO3/SrTiO3 interfaces, consistent with the presence of an interfacial 2-dimensional electron gas (2DEG) at each interface. Sheet carrier densities of 8.9 × 1014 cm−2 per interface are observed. Integration of such high density oxide 2DEGs on silicon provides a bridge between the exceptional properties and functionalities of oxide 2DEGs and microelectronic technologies.

Thin film crystalline oxides present a rich and diverse field of study, stimulating research in both fundamental condensed matter physics and the development of applications based on their material properties. Advancements in thin film growth and characterization have enabled artificial oxide heterostructures to be engineered with atomic-scale precision. These advancements have led to the discovery of interesting physical phenomena that arise at oxide interfaces, including switchable magnetism,1 superconductivity,2 and metallic conduction.3 Much of this rich behavior has been observed in 2-dimensional electron gases (2DEGs) that form at oxide interfaces in which a discontinuity in polarization exists. The RTiO3/SrTiO3 (RTO/STO) system is a quintessential example of such a polar interface, where R denotes a rare earth element such as La,4 Gd,5 Sm,6 or Nd.7 Along the {100} directions, RTiO3 has alternating planes of RO and TiO2 that are nominally +1 and −1 in ionic charge, respectively. In contrast, the alternating planes of SrO and TiO2 along the {100} directions of SrTiO3 are charge neutral. The formation of a 2DEG compensates or screens the discontinuity in polarization at the interface. In theory, 2DEGs that have a charge density of ∼3.3 × 1014 charges/cm2 (0.5 electrons per unit cell (u.c.)) fully compensate or screen the discontinuity in polarization. The densities of such 2DEGs are an order of magnitude higher than those of 2DEGs found in state-of-the-art AlGaN/GaN heterostructures8 and have been experimentally confirmed at interfaces between GdTiO3/STO.9 The high carrier densities offered by oxide 2DEGs could lead to the development of technologies for infrared and optical plasmonics,10,11 “charge gain” devices,12,13 and highly scaled devices in general.14–16 

In order to explore the potential application of oxide 2DEGs in optical and microelectronic devices, the integration of oxide 2DEGs on Si is required. Integration of crystalline oxides on a Si platform circumvents the problem of growing larger variations of ceramic substrates typically used to grow oxide heterostructures. Furthermore, the integration of oxide 2DEGs on silicon offers significant advantages for heat dissipation for potential power applications, since the thermal conductivity of silicon is an order of magnitude higher than that of typical oxide substrates. Here, we utilize oxide molecular beam epitaxy (MBE) to grow heterostructures comprised of LaTiO3 (LTO) and STO on Si (001). Structural characterization indicates the LTO/STO heterostructures are single crystalline and epitaxial to the Si substrate. The sheet resistances and carrier densities of the heterostructures are strongly dependent on the number of LTO/STO interfaces, consistent with the presence of high density 2DEGs at the interfaces between LTO and STO.

LTO/STO heterostructures are grown using reactive oxide MBE in a custom chamber operating at a base pressure of 1 × 10−10 Torr. The structures are epitaxially grown on undoped (>3000 Ω cm), (001) oriented silicon wafers that are grown using the float-zone technique (Virginia Semiconductor, Inc.). Prior to film deposition, the wafers are cleaned in piranha solution (3:1 H2SO4:H2O2, 130 °C), etched in dilute hydrofluoric acid (3:70 with water), and subsequently oxidized for 30 s under an ultraviolet lamp to form a layer of SiOx. The resulting protective layer of SiOx is desorbed from the surface by heating the substrate to 900 °C (thermocouple temperature) in ultra-high vacuum (UHV) until a 2 × 1 reconstruction is observed in the reflection high energy electron diffraction (RHEED) pattern, indicating the presence of a clean dimerized surface,17 as shown in Fig. 1(a). The substrate is then cooled to 670 °C, after which a 1/2 monolayer (ML) of Sr (99.99%, Sigma-Aldrich) is deposited to form a template for the subsequent deposition and crystallization of STO.18 The substrate is then cooled to room temperature, after which an additional 1/8 ML Sr is deposited and oxidized in a background pressure of 5 × 10−7 Torr of O2. Finally, deposition of 7/8 ML of SrO followed by codeposition of 3/2 ML of SrO and 2 ML of TiO2 (99.99% Ti, Alfa Aesar) is performed in a background pressure of 5 × 10−7 Torr of O2. The amorphous SrO and TiO2 layers are then annealed in UHV at 720 °C to form 2.5 u.c. of crystalline STO. Subsequent layers of STO and LTO (99.9% La, ESPI Metals) of various thicknesses and layering sequences, discussed below, are deposited at 600 °C. The STO layers are grown in an oxygen background pressure of 5 × 10−7 Torr, whereas LTO is grown at a reduced pressure of 1 × 10−7 Torr to avoid the formation of a pyrochlore phase (La2Ti2O7).19,20

FIG. 1.

Sequential in situ RHEED images of (a) a clean dimerized surface of Si showing a 2 × 1 reconstruction off of the 〈110〉 azimuth, (b) 4.5 u.c. of STO along 〈100〉, (c) 10 u.c. of LTO along 〈100〉, and (d) the final 15 u.c. of STO along 〈100〉. The LTO films have rougher RHEED patterns and show weak streaks in between the principal streaks, possibly indicating distortions of the oxygen octahedra in the perovskite structure.

FIG. 1.

Sequential in situ RHEED images of (a) a clean dimerized surface of Si showing a 2 × 1 reconstruction off of the 〈110〉 azimuth, (b) 4.5 u.c. of STO along 〈100〉, (c) 10 u.c. of LTO along 〈100〉, and (d) the final 15 u.c. of STO along 〈100〉. The LTO films have rougher RHEED patterns and show weak streaks in between the principal streaks, possibly indicating distortions of the oxygen octahedra in the perovskite structure.

Close modal

RHEED and x-ray diffraction (XRD) measurements indicate that the heterostructures are single crystalline and epitaxial to the underlying Si substrate. Sequential RHEED images for a tri-layer heterostructure comprised of 15 u.c. STO/10 u.c. LTO/4.5 u.c. STO/Si are shown in Figs. 1(b)1(d). The streaky RHEED patterns and high contrast of the streaks from the background indicate the growth mode is 2-dimensional and that the epitaxial layers are of excellent crystallinity. XRD 2θω spectra (Rigaku Smartlab diffractometer) are shown in Fig. 2(a). Data from all of the films are fit numerically and out-of-plane lattice constants (aoop) for both the STO and LTO layers are summarized in Table I. In comparison, bulk lattice parameters are 3.905 Å and 3.97 Å for STO and LTO,21 respectively, with the latter being represented by a pseudocubic unit cell. The slightly larger (smaller) measured aoop for LTO (STO) suggests that the LTO and STO are coherently strained to one another with the STO tensiley strained and the LTO compressively strained. After five unit cells, the LTO/STO layers are expected to be relaxed with respect to the silicon.22 

FIG. 2.

Structure of a 15 u.c. STO/10 u.c. LTO/4.5 u.c. STO/Si film. (a) 2θ–ω XRD spectrum around the 002 STO Bragg reflection. The data are shown in blue, with a fit to the nominal thicknesses overlaid as an orange line. The inset shows the surface morphology of a 2 × 2 μm2 region of the film scanned by tapping mode AFM. The film surface is smooth and continuous, with a root mean square roughness of ∼0.3 nm. (b) HAADF STEM image taken along the silicon 〈110〉 direction, showing the spatial separation between the different layers of the film.

FIG. 2.

Structure of a 15 u.c. STO/10 u.c. LTO/4.5 u.c. STO/Si film. (a) 2θ–ω XRD spectrum around the 002 STO Bragg reflection. The data are shown in blue, with a fit to the nominal thicknesses overlaid as an orange line. The inset shows the surface morphology of a 2 × 2 μm2 region of the film scanned by tapping mode AFM. The film surface is smooth and continuous, with a root mean square roughness of ∼0.3 nm. (b) HAADF STEM image taken along the silicon 〈110〉 direction, showing the spatial separation between the different layers of the film.

Close modal
TABLE I.

Number of interfaces and measured carrier concentrations at T = 140 K.

ia Structureb ns (cm−2)c ns/i (cm−2) n3D (cm−3)d n3D,max (cm−3)e rint (nm)f aoop,STO (nm)g aoop,LTO (nm)h
(15,10)1  1.7 × 1015  8.5 × 1014  1.5 × 1021  6.7 × 1021  0.4  0.3904  0.3995 
(7.5,5)2  3.9 × 1015  9.8 × 1014  3.4 × 1021  6.7 × 1021  0.3  0.3883  0.3996 
(4,4)3i  5.4 × 1015  9.0 × 1014  4.7 × 1021  8.4 × 1021  0.4  0.3864  0.3985 
(3,3)4i  6.7 × 1015  8.4 × 1014  5.8 × 1021  8.4 × 1021  0.6  0.3861  0.3995 
ia Structureb ns (cm−2)c ns/i (cm−2) n3D (cm−3)d n3D,max (cm−3)e rint (nm)f aoop,STO (nm)g aoop,LTO (nm)h
(15,10)1  1.7 × 1015  8.5 × 1014  1.5 × 1021  6.7 × 1021  0.4  0.3904  0.3995 
(7.5,5)2  3.9 × 1015  9.8 × 1014  3.4 × 1021  6.7 × 1021  0.3  0.3883  0.3996 
(4,4)3i  5.4 × 1015  9.0 × 1014  4.7 × 1021  8.4 × 1021  0.4  0.3864  0.3985 
(3,3)4i  6.7 × 1015  8.4 × 1014  5.8 × 1021  8.4 × 1021  0.6  0.3861  0.3995 
a

i is the number of LTO/STO interfaces in the film.

b

The film structure is expressed in the format (m,n)x, where x is the number of repeats of a bilayer stack containing m u.c. of STO on n u.c. of LTO. All films are grown on a 4.5 u.c. STO buffer layer on silicon.

c

ns is measured sheet carrier density of the entire film.

d

n3D is the 3-dimensional carrier density of the film under the assumption that conduction occurs throughout the entire film.

e

n3D,max is the theoretical maximum 3-dimensional carrier density given complete intermixing of the LTO and STO in the film, i.e., assuming the entire film is LaxSr1−xTiO3, x = 0.5 for i = 2, 4 and x = 0.4 for i = 6, 8.

f

rint is the average interface roughness (±0.1 nm) extracted from fits from x-ray reflectivity data.

g

aoop,STO is the average out-of-plane lattice parameter (±0.0005 nm) extracted for the STO layers in the film.

h

aoop,LTO is the average out-of-plane lattice parameter (±0.0005 nm) extracted for the LTO layers in the film.

i

The 6- and 8-interface samples have one additional unit cell of STO at the top in order to keep the total film constant at 29.5 u.c.

One result of relaxation during film growth is increased interface roughness. Fits to x-ray reflectivity (XRR) measurements indicate that the average interface roughness is on the order of 1 u.c. (see Table I), with the roughness lowest in the 4.5 u.c. STO buffer layer on Si and highest in the topmost layer. The film roughness persists to the surface of the heterostructure stack, as indicated by atomic force microscopy (AFM, Digital Instruments Nanoscope III) shown in the right inset of Fig. 2. A high-angle annular dark-field (HAADF) image shown in Fig. 2(b) obtained using a scanning transmission electron microscope operating at 200 kV (STEM, Tecnai Osiris, FEI) indicates that the interfaces between the STO/Si and LTO/STO are abrupt at atomic length scales, although the interface between STO/LTO appears less sharp, which is consistent with the extracted roughness values from the XRR data.

Electrical transport measurements indicate that the conductivity within our LTO/STO/Si heterostructures arises from the formation of a 2DEG at the interface between LTO and STO. Four-point electrical transport measurements (Quantum Design, Physical Properties Measurement System) are performed as a function of temperature in the Van der Pauw geometry. Electrical contact to the LTO/STO/Si heterostructures is established by wire bonding Al wires to Au contacts that are deposited on the corners of 5 × 5 mm2 samples. Figure 3 shows sheet resistance measurements of tri-layer heterostructures of the format 15 u.c. STO/x u.c. LTO/4.5 u.c. STO/Si as a function of the number of LTO unit cells, x. For x = 0, i.e., no LTO layer, the STO film exhibits semiconducting behavior, which we attribute to the presence of residual oxygen vacancies that act as n-type dopants. The insertion of just 2 u.c. (x = 2) of LTO into the STO film induces a significant reduction in the sheet resistance, as shown in Fig. 3. However, this decrease in sheet resistance is independent of the thickness of the LTO layer; namely, the sheet resistance does not scale with LTO thickness as shown in Fig. 3 for x = 10. From this, we can rule out conduction through the bulk of LTO, which has been shown when the films are grown oxygen-rich20 or when LTO is compressively strained, where films compressively strained above 1.6% yield conductive behavior.21 Measurements of the lattice parameter for the LTO/STO heterostructures discussed above indicate a strain less than that required for conduction in the LTO. In addition, we calculate an upper bound for compressive strain in the LTO due to the STO in the heterostructure to be 0.7%, significantly smaller than the amount observed to yield a metal-insulator transition. From both the measured strain and the LTO thickness dependence, we conclude that conduction within our LTO/STO heterostructures occurs at the LTO/STO interface.

FIG. 3.

Sheet resistance plotted against temperature for films of 15 u.c. STO/x u.c. LTO/4.5 u.c. STO grown on undoped float zone silicon, whose sheet resistance is also shown by itself for comparison (black squares). The 19.5 u.c. STO/Si sample (dark cyan circles) has semiconducting behavior, but once LTO is inserted into the structure, the sheet resistance drops and is unaffected by the thickness of LTO (blue diamonds and red triangles).

FIG. 3.

Sheet resistance plotted against temperature for films of 15 u.c. STO/x u.c. LTO/4.5 u.c. STO grown on undoped float zone silicon, whose sheet resistance is also shown by itself for comparison (black squares). The 19.5 u.c. STO/Si sample (dark cyan circles) has semiconducting behavior, but once LTO is inserted into the structure, the sheet resistance drops and is unaffected by the thickness of LTO (blue diamonds and red triangles).

Close modal

To further establish the presence of interfacial conductivity, the number of interfaces within our heterostructures is varied while keeping the overall thickness of the heterostructures constant. Heterostructures with 2, 4, 6, and 8 LTO/STO interfaces are grown and electrically characterized, as summarized in Table I. All samples have an identical total film thickness of 29.5 u.c., consisting of 25 u.c. of heterostructures on a 4.5 u.c. STO buffer deposited on silicon. Figure 4(a) shows that the sheet resistance is strongly dependent on the number of LTO/STO interfaces, decreasing monotonically with increasing number of interfaces.

FIG. 4.

Effect of the number of interfaces on conductivity. All films shown have the same total thickness of 29.5 u.c., but with a differing number of LTO/STO interfaces. (a) shows the sheet resistance monotonically drops as the number of interfaces is increased from 2 to 8. (b) shows the scaling of total sheet carrier density with the number of interfaces, as extracted from Hall measurements at 140 K. For conduction to be an interface effect, one expects a linear increase of the total sheet carrier density as the number of interfaces is increased. (c) shows the sheet carrier density per interface as a function of the number of interfaces. One should again expect the per-interface density to remain constant as the total number of interfaces increases. The red dashed lines for (b) and (c) are guides to the eye and show a linear fit to the data (black circles).

FIG. 4.

Effect of the number of interfaces on conductivity. All films shown have the same total thickness of 29.5 u.c., but with a differing number of LTO/STO interfaces. (a) shows the sheet resistance monotonically drops as the number of interfaces is increased from 2 to 8. (b) shows the scaling of total sheet carrier density with the number of interfaces, as extracted from Hall measurements at 140 K. For conduction to be an interface effect, one expects a linear increase of the total sheet carrier density as the number of interfaces is increased. (c) shows the sheet carrier density per interface as a function of the number of interfaces. One should again expect the per-interface density to remain constant as the total number of interfaces increases. The red dashed lines for (b) and (c) are guides to the eye and show a linear fit to the data (black circles).

Close modal

Hall measurements show that the sheet carrier density scales with the number of interfaces within the heterostructures, indicating that the conductivity occurs predominantly at LTO/STO interfaces. We present Hall data measured at an intermediate temperature of ∼140 K, where the Hall resistance is linear over the measured field, indicating that the effects of thermally generated carriers in the Si are minimized. Scaling of the sheet carrier density with the number of interfaces within the LTO/STO/Si heterostructures is observed, as shown and summarized in Fig. 4(b) and Table I, respectively. The average sheet carrier density per LTO/STO interface, shown in Fig. 4(c), is ∼8.9 × 1014 cm−2, which is nearly a factor of 3 higher than the theoretical limit of 3.3 × 1014 cm−2 needed to compensate for the polar discontinuity at the interface. While oxygen vacancies in the STO generated during growth of LTO under low O2 pressures have been cited as a potential source of additional carriers,23,24 we believe the higher measured densities may arise from intermixing at the LTO/STO interfaces. Because the increase of carrier density does not scale with the total number of STO layers, the additional carriers cannot be attributed to oxygen vacancies in the bulk. Fig. 4(b) shows an extrapolated zero-interface intercept reasonably close to the origin, suggesting that other non-interfacial means of conduction are not significant. However, we cannot rule out oxygen vacancies that may be localized to the LTO/STO interfaces. Alternatively, intermixing of La and Sr at each LTO/STO interface could generate additional carriers. Table I shows the theoretical maximum to the 3D carrier concentration, assuming the LTO/STO heterostructures are completely intermixed to form solid solution LaxSr1−xTiO3. The measured 3D carrier density of each LTO/STO heterostructure is lower than the corresponding solid-solution equivalent, indicating the conductivity is predominantly an interfacial effect.25 These conclusions are also consistent with the measured interface roughness for the films. Therefore, La/Sr intermixing at interfaces may explain the higher than expected measured carrier concentration.

In summary, we have demonstrated epitaxial growth of the LTO/STO oxide system onto silicon by molecular beam epitaxy. We show that the heterostructures are single crystalline and epitaxial to the Si substrate. The LTO/STO interface generates a 2DEG, which we measure to have a carrier concentration of 8.9 × 1014 cm−2. We attribute the higher than theoretical density to additional carriers originating from intermixing at the LTO/STO interfaces, but also note that other mechanisms such as interface-bound oxygen vacancies may play a role as well. The integration of a high density oxide 2DEG on Si provides a pathway to explore the exceptional properties of this electronic system in microelectronic and photonic devices.

This work was funded by the Office of Naval Research Multidisciplinary University Research Initiative to support the EXtreme Electron DEvices (EXEDE) program, along with partial support for by the National Science Foundation through NSF DMR-1309868 and NSF MRSEC DMR-1119826. Marvin Wint, Carol Jenkins, and Alexis Ernst are acknowledged for valuable technical assistance and TEM sample preparation.

1.
Y.-H.
Chu
,
L. W.
Martin
,
M. B.
Holcomb
,
M.
Gajek
,
S.-J.
Han
,
Q.
He
,
N.
Balke
,
C.-H.
Yang
,
D.
Lee
,
W.
Hu
,
Q.
Zhan
,
P.-L.
Yang
,
A.
Fraile-Rodriguez
,
A.
Scholl
,
S. X.
Wang
, and
R.
Ramesh
,
Nat. Mater.
7
,
478
(
2008
).
2.
C. H.
Ahn
,
S.
Gariglio
,
P.
Paruch
,
T.
Tybell
,
L.
Antognazza
, and
J.-M.
Triscone
,
Science
284
,
1152
(
1999
).
3.
A.
Ohtomo
and
H. Y.
Hwang
,
Nature
427
,
423
(
2004
).
4.
K.
Shibuya
,
T.
Ohnishi
,
M.
Kawasaki
,
H.
Koinuma
, and
M.
Lippmaa
,
Jpn. J. Appl. Phys., Part 2
43
,
L1178
(
2004
).
5.
P.
Moetakef
,
J. Y.
Zhang
,
A.
Kozhanov
,
B.
Jalan
,
R.
Seshadri
,
S. J.
Allen
, and
S.
Stemmer
,
Appl. Phys. Lett.
98
,
112110
(
2011
).
6.
J. Y.
Zhang
,
C. A.
Jackson
,
R.
Chen
,
S.
Raghavan
,
P.
Moetakef
,
L.
Balents
, and
S.
Stemmer
,
Phys. Rev. B: Condens. Matter Mater. Phys.
89
,
075140
(
2014
).
7.
P.
Xu
,
D.
Phelan
,
J. S.
Jeong
,
K. A.
Mkhoyan
, and
B.
Jalan
,
Appl. Phys. Lett.
104
,
082109
(
2014
).
8.
O.
Ambacher
,
J.
Smart
,
J. R.
Shealy
,
N. G.
Weimann
,
K.
Chu
,
M.
Murphy
,
W. J.
Schaff
,
L. F.
Eastman
,
R.
Dimitrov
,
L.
Wittmer
,
M.
Stutzmann
,
W.
Rieger
, and
J.
Hilsenbeck
,
J. Appl. Phys.
85
,
3222
(
1999
).
9.
P.
Moetakef
,
T. A.
Cain
,
D. G.
Ouellette
,
J. Y.
Zhang
,
D. O.
Klenov
,
A.
Janotti
,
C. G.
Van de Walle
,
S.
Rajan
,
S. J.
Allen
, and
S.
Stemmer
,
Appl. Phys. Lett.
99
,
232116
(
2011
).
10.
J. A.
Schuller
,
E. S.
Barnard
,
W. S.
Cai
,
Y. C.
Jun
,
J. S.
White
, and
M. L.
Brongersma
,
Nat. Mater.
9
,
193
(
2010
).
11.
M. L.
Brongersma
and
V.
Shalaev
,
Science
328
,
440
(
2010
).
12.
J.
Son
,
S.
Rajan
,
S.
Stemmer
, and
S. J.
Allen
,
J. Appl. Phys.
110
,
084503
(
2011
).
13.
M.
Boucherit
,
O.
Shoron
,
C. A.
Jackson
,
T. A.
Cain
,
M. L. C.
Buffon
,
C.
Polchinski
,
S.
Stemmer
, and
S.
Rajan
,
Appl. Phys. Lett.
104
,
182904
(
2014
).
14.
S.
Rajan
,
A.
Chini
,
M. H.
Wong
,
J. S.
Speck
, and
U. K.
Mishra
,
J. Appl. Phys.
102
,
044501
(
2007
).
15.
P. S.
Park
and
S.
Rajan
,
IEEE Trans. Electron Devices
58
,
704
(
2011
).
16.
D. N.
Nath
,
S.
Keller
,
E.
Hsieh
,
S. P.
DenBaars
,
U. K.
Mishra
, and
S.
Rajan
,
Appl. Phys. Lett.
97
,
162106
(
2010
).
17.
T.
Sakamoto
,
T.
Kawamura
, and
G.
Hashiguchi
,
Appl. Phys. Lett.
48
,
1612
(
1986
).
18.
R. A.
McKee
,
F. J.
Walker
, and
M. F.
Chisholm
,
Phys. Rev. Lett.
81
,
3014
(
1998
).
19.
A.
Ohtomo
,
D. A.
Muller
,
J. L.
Grazul
, and
H. Y.
Hwang
,
Appl. Phys. Lett.
80
,
3922
(
2002
).
20.
F.
Lichtenberg
,
A.
Herrnberger
,
K.
Wiedenmann
, and
J.
Mannhart
,
Prog. Solid State Chem.
29
,
1
(
2001
).
21.
F. J.
Wong
,
S.-H.
Baek
,
R. V.
Chopdekar
,
V. V.
Mehta
,
H.-W.
Jang
,
C.-B.
Eom
, and
Y.
Suzuki
,
Phys. Rev. B: Condens. Matter Mater. Phys.
81
,
161101
(R) (
2010
).
22.
D. P.
Kumah
,
J. W.
Reiner
,
Y.
Segal
,
A. M.
Kolpak
,
Z.
Zhang
,
D.
Su
,
Y.
Zhu
,
M. S.
Sawicki
,
C. C.
Broadbridge
,
C. H.
Ahn
, and
F. J.
Walker
,
Appl. Phys. Lett.
97
,
251902
(
2010
).
23.
W.
Siemons
,
G.
Koster
,
H.
Yamamoto
,
W. A.
Harrison
,
G.
Lucovsky
,
T. H.
Geballe
,
D. H. A.
Blank
, and
M. R.
Beasley
,
Phys. Rev. Lett.
98
,
196802
(
2007
).
24.
A.
Kalabukhov
,
R.
Gunnarsson
,
J.
Borjesson
,
E.
Olsson
,
T.
Claeson
, and
D.
Winkler
,
Phys. Rev. B: Condens. Matter Mater. Phys.
75
,
121404
(R) (
2007
).
25.
Y.
Aiura
,
H.
Bando
,
I.
Hase
,
Y.
Nishihara
,
Y.
Haruyama
, and
H.
Suzuki
,
Superlattices Microstruct.
21
,
321
(
1997
).