A hydrogen atom, characterized by one unpaired electron and the smallest atomic radius, underlies the operations of various solid-state devices such as transistors, capacitors, solar cells, etc. Given its specific character as donor impurity in oxides, hydrogen may also facilitate efficient electron doping in a wide range of oxide devices. Here, we demonstrate room-temperature electrochemical hydrogenation of an archetypical oxide semiconductor (TiO2) thin film to achieve a 3D-compatible electron doping technique. The hydrogenated region can be precisely defined by photolithography without the influence of polycrystalline grain boundaries. Besides, secondary ion mass spectroscopy with deuterium isotope reveals considerable amount of hydrogen condenses around the TiO2 bottom interface indicating the critical influence of the interface on hydrogen stability. This hydrogen shows excellent stability in contrast to its high diffusivity in bulk TiO2, enabling robust electron doping for oxide thin film devices as well as suggesting stable interface hydrogen reservoir for electrochemical phenomena.

Hydrogen underlies the operations of various solid-state devices, passivating defects in silicon transistors and solar cells,1 causing leakage current in ceramic capacitors and piezoelectric devices,2 mediating charge transfer in fuel cells and atomic switches,3,4etc. Hydrogen can also behave as donor impurity in various oxide semiconductors generating one electron carrier and increasing electrical conductivity.5,6 When considered as donor impurity, hydrogen has a unique character of extremely high permeability in solids,7 inducing various intriguing phenomena at room temperature.8–11 By taking advantage of this characteristics, hydrogen may facilitate efficient electron doping for a range of oxide devices, which is applicable at room temperature in a 3D-compatible way.

Generally speaking, the technique of electron doping is essential for the study of device functionalities and physical phenomena in oxides.12,13 So far, electron doping to oxides has often used in-situ doping during the growth or thermal treatment in the reactive environment.14,15 However, these techniques lack the microscopic controllability of dopant distribution in the device structure, in contrast to the well established ion implantation for conventional semiconductors.16 For this issue, an electrochemical hydrogenation combined with photolithography could provide an alternative pathway, which could control the dopant distribution with high spacial resolution. The exploitation of electrochemistry is reasonable considering a rich variety of electrochemical phenomena in oxides such as anodic oxidization,17,18 electrochromism,19,20 resistance switching,21,22 and so on. While the electrochemical hydrogenation has been used for nanotubes and nanoparticles to optimize the optical, catalytic, and photovoltaic properties,19,23–25 it has rarely been applied to oxide semiconductor thin films in combination with photolithography. In this study, the electrochemical hydrogenation is successfully demonstrated at room temperature in the archetypical wide-gap oxide semiconductor TiO2 thin film with photo-resist mask on the surface [Fig. 1(a)]. The unmasked region of the oxide thin film works as the cathode, inducing the cell reaction with the anode, namely, the external or surface metal electrode. This technique is inherently 3D-compatible with a reasonably high spatial resolution of several tens nanometer, showing the potential to be used for various oxide devices.

FIG. 1.

(a) A schematic illustration of the electrochemical hydrogenation of the TiO2 thin film with a photo-resist mask on the surface. When immersed in the electrolyte, the cell reaction occurs between the unmasked region of the TiO2 thin film (cathode) and the deposited metal electrode (anode). (b) The surface SEM image and (c) the cross section TEM image of the TiO2 thin film on the SiO2/p:Si substrate.

FIG. 1.

(a) A schematic illustration of the electrochemical hydrogenation of the TiO2 thin film with a photo-resist mask on the surface. When immersed in the electrolyte, the cell reaction occurs between the unmasked region of the TiO2 thin film (cathode) and the deposited metal electrode (anode). (b) The surface SEM image and (c) the cross section TEM image of the TiO2 thin film on the SiO2/p:Si substrate.

Close modal

The electrochemical hydrogenation is demonstrated in polycrystalline anatase TiO2 thin films on SiO2/p:Si substrates. Polycrystalline anatase TiO2 thin films with the thickness of 30-40 nm are fabricated by amorphous deposition at room temperature and crystallization during the thermal treatment.26 Amorphous deposition is performed on the SiO2(100-500 nm)/Si substrate by pulsed laser deposition (KrF) at the 1 Pa oxygen pressure. Crystallization is induced by thermal treatment in O2 at 500-600 oC for 30-60 min. Scanning electron microscopy (SEM) in Fig. 1(b) shows anisotropic polycrystalline grain with a micrometer length scale in the in-plane direction. Transmission electron microscopy (TEM) in Fig. 1(c) shows a well-defined plate-like anatase grain, which helps us elucidate the role of grain boundary in the hydrogenation process as shown later.

Electrochemical hydrogenation is performed with the three-terminal potentiostat. The working electrode is the TiO2 thin film with a small gold pad on the surface for the electrical contact. The counter electrode is platinum and the reference electrode is Ag|AgCl in the saturated KCl solution. The working electrode and the counter electrode are immersed in the potassium phosphate buffer solution with pH = 6.3, which is bridged to the reference electrode via the KCl salt bridge. After keeping the target potential for 1 min, the working electrode is rinsed by water, and the electrical measurement is performed. Here, the electrical potential is defined by the external voltage applied on the working electrode with respect to the reference electrode, from which the electrode potential of the reference electrode (0.199 V vs. SHE) is also subtracted.

As shown in Fig. 2(a), the sheet resistance of the dried TiO2 thin film decreases from ∼109 Ω (insulating) at -0.6 V vs. standard hydrogen electrode (SHE), and saturates at ∼104 Ω (conductive) below -0.8 V vs. SHE (filled diamonds). This decrease in the sheet resistance indicates the electrochemical reduction of TiO2, in other words, the electron doping. The saturation of the sheet resistance implies the reduction reaction of TiO2 is impeded by the side reaction such as hydrogen generation. This reduced TiO2 thin film is further exposed to various electrical potentials subsequently, and the sheet resistance is found to be modulated within one order reversibly as shown by open diamonds and square in Fig. 2(a). However, the sheet resistance never recovers to the initial value ∼109 Ω, possibly because the positive electrical potential on TiO2 forms a blocking contact (depletion region) between TiO2 and electrolyte.

FIG. 2.

(a) Sheet resistance of the TiO2 thin film as a function of its electrical potential in the electrolyte. The sheet resistance is measured in the air after rinsing and drying of the sample. The empty diamonds are the data after the complete electrochemical doping at -0.9V vs. SHE for 1 min and the subsequent exposure to each electrical potential for 1 min. The empty square is the data after the complete electrochemical doping at -0.9V vs. SHE, the subsequent oxidization at 1 V for 1 min, and the subsequent exposure to the target potential for 1 min. (b,c) In-plane X-ray diffraction of the TiO2 thin film before doping (“Initial”) and after doping (“Doped”) (b) in a wide range and (c) in a narrow range around the anatase TiO2 (101) peak. The indices in (b) are for anatase TiO2.

FIG. 2.

(a) Sheet resistance of the TiO2 thin film as a function of its electrical potential in the electrolyte. The sheet resistance is measured in the air after rinsing and drying of the sample. The empty diamonds are the data after the complete electrochemical doping at -0.9V vs. SHE for 1 min and the subsequent exposure to each electrical potential for 1 min. The empty square is the data after the complete electrochemical doping at -0.9V vs. SHE, the subsequent oxidization at 1 V for 1 min, and the subsequent exposure to the target potential for 1 min. (b,c) In-plane X-ray diffraction of the TiO2 thin film before doping (“Initial”) and after doping (“Doped”) (b) in a wide range and (c) in a narrow range around the anatase TiO2 (101) peak. The indices in (b) are for anatase TiO2.

Close modal

In spite of the large decrease in sheet resistance, the TiO2 thin film shows negligible change in its crystalline structure by electrochemical hydrogenation. X-ray diffraction before (“Initial”) and after (“Doped”) hydrogenation shows no difference neither in the peak position, the peak intensity, nor the peak width [Figs. 2(b) and 2(c)], indicating negligible change in the lattice constant and crystallinity. Here, X-ray diffraction is performed in the plane by using grazing incident X-ray in order to maximize the diffraction intensity. The angle between the X-ray source and the sample surface is 0.3° (below the total reflection angle), and the detector is parallel to the sample surface (2θ = ω = 0.3°). In the plane, the angle between the X-ray source and the detector (2θχ) is kept two times larger than the angle of in-plane sample rotation (Φ). All the diffraction peaks in Fig. 2(b) correspond to the anatase phase of TiO2. The absence of structural change in TiO2 is also confirmed by Raman spectroscopy (not shown), indicating a great advantage of this technique for damage-less electron doping.

In order to elucidate the mechanism of hydrogenation and also the ability of restricting the hydrogenated region, the electrochemical process is performed through the line-and-space photo resist mask on the TiO2 surface. In this case, because the sheet resistance of the TiO2 thin film saturates below -0.8 V vs. SHE as shown in Fig. 2(a), the same electrochemical reduction occurs simply by putting aluminum anode (-1.66 V vs. SHE) on the TiO2 surface as shown in Fig. 1(a), and immersing the sample in diluted hydrochloric acid or diluted tetramethylammonium hydroxide solution for 1 min without external voltage. Then, the photo resist mask is removed by acetone, and the conductive area on the TiO2 thin film is measured by conductive atomic force microscopy (cAFM). Here, the AFM and cAFM is performed by using Cypher S of Oxford Instruments. The iridium coated silicon tip is grounded, and the hydrogenated TiO2 thin film is electrically biased by 3 V during the measurement. The tip curvature radius of 28 nm enables the spatial resolution of a few tens nm.

The cAFM image clearly shows the line-and-space pattern corresponding to the area which is uncovered by the photo-resist mask [Fig. 3(a)]. This line-and-space conductivity pattern is in contrast to the simultaneously measured TiO2 surface morphology [Fig. 3(b)], which only shows a flat surface with a root mean square of ∼0.3 nm. In this electrochemical process, there are three possible mechanisms of TiO2 reduction: the reduction simply through the top surface, the reduction mediated by the grain boundary, and the reduction via the underlying SiO2 layer.27,28 However, the last possibility can be denied by the experimentally observed conductivity pattern because the reduction from the underneath would not be affected by the surface mask and lead to a homogeneous conductivity pattern. In order to further clarify the reduction mechanism, the magnified scan is performed for the square region in Figs. 3(a) and 3(b). The surface morphology in this region shows several grain boundaries as shown by arrows in Figure 3(d). On the other hand, the corresponding cAFM image in Figure 3(c) shows nontrivial conductivity even in the area far from the grain boundaries, indicating the reduction occurs relatively homogeneously inside the TiO2 grain. Besides, the edge line of the doped region is straight and almost unaffected by the position of each grain boundary. This negligible influence of the grain boundary indicates the reduction reaction is not mediated by the grain boundary but mainly occurs through the grain surface. This result is also consistent with the ability of electrochemical reduction even in the anatase single crystal, which has no grain boundary.10 It should be noted the apparent high conductivity along the grain boundary in Fig. 3(c) does not necessarily mean the higher conductive at the boundary, but includes crosstalks from the surface morphology signal.

FIG. 3.

(a,c) The local electrical conductivity and (b,d) the surface morphology of the TiO2 thin film which is electrochemically hydrogenated through the line-and-space resist mask. (a,b) are the wide scans and (c,d) are the magnified scans inside the squares in (a,b). The white arrows in (d) indicate the grain boundaries in the TiO2 thin film.

FIG. 3.

(a,c) The local electrical conductivity and (b,d) the surface morphology of the TiO2 thin film which is electrochemically hydrogenated through the line-and-space resist mask. (a,b) are the wide scans and (c,d) are the magnified scans inside the squares in (a,b). The white arrows in (d) indicate the grain boundaries in the TiO2 thin film.

Close modal

The specific reduction reaction in this electrochemical process would be the hydrogenation according to the past studies.9,10,25 The hydrogenation of TiO2 is actually confirmed by secondary ion mass spectroscopy (SIMS). The SIMS is measured by NANO SCIENCE CORPORATION with 2 keV cesium ion beam in ADEPT-1010 (ULVAC-PHI, Inc.). The backside SIMS (shown later) is measured by Toray Research Center, Inc. with 3 keV cesium ion beam in ADEPT-1010 (ULVAC-PHI, Inc.). When SIMS is applied to hydrogen, it should be distinguished from the hydrogen inside the surface adsorbed water. To address this issue, the TiO2 thin film is electrochemically reduced in the electrolyte using D2O, and is subsequently rinsed in H2O in order to exchange the adsorbed water. In the experiment, this “Doped” sample is prepared by immersing the TiO2 thin film with a small aluminum pad in 0.5 mol L-1 deuterium chloride solution for 1 min, and then by rinsing it in H2O for 30 sec. The “Non-doped” sample (for reference) is prepared by immersing the TiO2 thin film without an aluminum pad in 0.5 mol L-1 deuterium chloride solution for 1 min, and then by rinsing it in H2O for 30 sec.

The SIMS spectra in Fig. 4(a) show a large amount of deuterium in this “Doped” sample while showing almost no deuterium in the reference “Non-doped” sample. The deuterium concentration integrated throughout the spectrum (from the SiO2 layer to the TiO2 surface) is around 1 × 1014 cm-2, approximately consistent with the sheet electron density from the Hall measurement (N2D = 6 × 1013 cm-2). It should be noted Al and Cl are also measured (using 2 keV oxygen ion beam only for Al), but the obtained concentration is below ∼3 × 1016 cm-3 for Al, and below ∼3 × 1017 cm-3 for Cl.

FIG. 4.

(a) SIMS spectra of deuterium for the electrochemically hydrogenated sample (“Doped”) and the reference sample (“Non-doped”). The backside SIMS spectrum is also shown for the hydrogenated sample (“Doped, back”). The dashed line indicates the TiO2/SiO2 interface. (b) The SIMS spectra taken from the top side of the TiO2 surface and (c) the back-side SIMS spectra taken from the SiO2 side, for H (open diamond) and D (filled diamond) in the electrochemically hydrogenated sample. Both in (b,c), the large H concentration may originate from the adsorbed water in the grain boundary as well as the preparation process for the SIMS measurement. (d) TiO2 thickness dependence of the sheet electron density (N2D) from the Hall measurement. The dashed line is the least square fitting of the data.

FIG. 4.

(a) SIMS spectra of deuterium for the electrochemically hydrogenated sample (“Doped”) and the reference sample (“Non-doped”). The backside SIMS spectrum is also shown for the hydrogenated sample (“Doped, back”). The dashed line indicates the TiO2/SiO2 interface. (b) The SIMS spectra taken from the top side of the TiO2 surface and (c) the back-side SIMS spectra taken from the SiO2 side, for H (open diamond) and D (filled diamond) in the electrochemically hydrogenated sample. Both in (b,c), the large H concentration may originate from the adsorbed water in the grain boundary as well as the preparation process for the SIMS measurement. (d) TiO2 thickness dependence of the sheet electron density (N2D) from the Hall measurement. The dashed line is the least square fitting of the data.

Close modal

Importantly, the deuterium profile in Fig. 4(a) is strongly inhomogeneous in the plane-normal direction, and shows a larger concentration on the side of the TiO2/SiO2 bottom interface. This inhomogeneous deuterium profile is not attributed to the deuterium diffusion during the measurement because the back-side SIMS, which is measured in the opposite direction from the SiO2 side, also shows the similar inhomogeneous profile (“Doped, back” in Fig. 4(a)). In order to further investigate the influence of the interface, the N2D in the hydrogenated TiO2 thin films is plotted as a function of the TiO2 film thicknesses in Fig. 4(d). A nontrivial amount of N2D remains at the zero thickness limit (∼4 × 1013 cm-2), indicating a considerable amount of donor hydrogen is located at the interface rather than inside the TiO2 thin film. The existence of the interfacial hydrogen is also consistent with the large interfacial deuterium concentration for the “Doped” sample in Fig. 4(a). This interfacial deuterium concentration is unaccountable only by the matrix effect when it is compared with the SIMS spectra of hydrogen (H) as shown in Figs. 4(b) and 4(c). In both cases, the interfacial concentration of D is much larger than that of H although the matrix effects for H and D have approximately similar trend to each other.29 It should be noted the interfacial deuterium concentration becomes smaller in the back-side SIMS because the thermal process during the preparation process for this measurement (60 °C for 6 hours) may partially release deuterium.

One of the concerns about the electrochemical hydrogenation is the low stability of doped hydrogen, which is previously reported.9 Indeed, most electrochemically doped hydrogen in bulk TiO2 does not form -OH bonds, and hence, diffuses out only in a few hours. The formation of rigid -OH bonds requires the thermal treatment in the hydrogen atmosphere. Surprisingly, the electrochemically doped hydrogen in the TiO2 thin film at room temperature is stable as indicated by the time dependence of N2D [Fig. 5(a)]. Although N2D decreases by 20 % in the first day, the decay becomes slower, and N2D almost saturates after several tens days. Fig. 5(b) shows the log-log plot of N2D as a function of time, where the simple extrapolation of the experimental data implies N2D can be maintained above 1 × 1013 cm-2 for as long as 106 days (∼2,700 years). Generally, the hydrogen can be stabilized by dangling bonds and interstitial sites, for example, the dangling bonds in amorphous silicon,30 and the interstitial sites in LnNi5-type (Ln: lanthanoid) hydrogen storage alloy.31 Considering a lot of dangling bonds and interstitial sites exist at the interface, they may account for the experimentally observed high stability of hydrogen. It should be noted that the strong interaction of the interface with hydrogen has also been recognized in other materials systems such as the SiO2/Si interface in silicon transistors,32 the Pd/oxide interface in hydrogen gas sensors,33 and the metal/ceramics interface in the mechanical component.34 

FIG. 5.

(a) The sheet resistance (Rsh), the sheet electron density (N2D), and the Hall electron mobility (μHall) of the hydrogenated TiO2 thin film as functions of the exposure time to the air. (b) The log-log plot of N2Dvs. time. The dashed line is the linear extrapolation of the data after two days.

FIG. 5.

(a) The sheet resistance (Rsh), the sheet electron density (N2D), and the Hall electron mobility (μHall) of the hydrogenated TiO2 thin film as functions of the exposure time to the air. (b) The log-log plot of N2Dvs. time. The dashed line is the linear extrapolation of the data after two days.

Close modal

One of the unique characteristics of this electrochemical hydrogenation with respect to other doping technique is the inhomogeneous dopant distribution perpendicular to the interface; the hydrogen concentration is largest at the interface and becomes smaller far from the interface [Fig. 4(a)]. Obviously, this distribution has no problem when we just impart electrical conductivity to the TiO2 thin films or nanoparticles. In such cases, the inhomogeneous distribution would be rather advantageous due to the high hydrogen stability as shown in Fig. 5. On the other hand, the inhomogeneous distribution may be problematic for fundamental research, which often needs to control the absolute dopant concentration with perfect spatial homogeneity. In such cases, the thermal treatment in hydrogen atmosphere is more preferable because it may achieve homogeneous hydrogen distribution by stabilizing hydrogen in the bulk part of TiO2.9 Thus, the inhomogeneity of electrochemical hydrogenation makes this technique more suitable for simply imparting electrical conductivity to the TiO2 thin films and nanoparticles than for precisely controlling hydrogen dopant concentration.

In conclusion, we demonstrate the electrochemical hydrogenation of the polycrystalline anatase TiO2 thin film at room temperature, transforming the insulating TiO2 to the highly conductive n-type semiconductor. The hydrogenation occurs almost homogeneously in the plane and little influenced by the grain boundary. As a result, the hydrogenated area can be precisely controlled by photo-resist mask on the surface, providing a powerful electron doping technique for oxide devices. In the plane-normal direction, it is indicated a considerable amount of hydrogen exists around the bottom interface in the hydrogenated TiO2 thin film. This interfacial hydrogen reservoir will affect electrochemical phenomena in TiO2, as well as facilitate novel devices such as ion-electron transducers and neuromorphic devices.35,36

This work was supported by MEXT/JSPS KAKENHI Grant Number JP17K18869.

1.
J. I.
Pankove
and
N. M.
Johnson
,
Semiconductors and Semimetals
34
,
1
(
1991
).
2.
W. P.
Chen
,
Y.
Wang
,
J. Y.
Dai
,
S. G.
Lu
,
X. X.
Wang
,
P. F.
Lee
,
H. L. W.
Chan
, and
C. L.
Choy
,
Appl. Phys. Lett.
84
,
103
(
2004
).
3.
G. N. Q.
Minh
,
J. Am. Ceram. Soc.
76
,
563
(
1993
).
4.
T.
Tsuruoka
,
I.
Valov
,
C.
Mannequin
,
T.
Hasegawa
,
R.
Waser
, and
M.
Ando
,
Jpn. J. Appl. Phys.
55
,
06GJ09
(
2016
).
5.
C. G.
Van de Walle
,
Phys. Rev. Lett.
85
,
1012
(
2000
).
6.
C.
Kilic
and
A.
Zunger
,
Appl. Phys. Lett.
81
,
73
(
2002
).
7.
H.
Iwahara
,
Solid State Ionics
52
,
99
(
1992
).
8.
P. F.
Chester
and
D. H.
Bradhurst
,
Nature
199
,
1056
(
1963
).
9.
L. A.
Harris
and
R.
Schumacher
,
J. Electrochem. Soc.
127
,
1186
(
1980
).
10.
H.
Pelouchova
,
P.
Janda
,
J.
Weber
, and
I.
Kavan
,
J. Electroanal. Chem
566
,
73
(
2004
).
11.
Z.
Zhang
,
D.
Schwanz
,
B.
Narayanan
,
M.
Kotiuga
,
J. A.
Dura
,
M.
Cherukara
,
H.
Zhou
,
J. W.
Freeland
,
J.
Li
,
R.
Sutarto
,
F.
He
,
C.
Wu
,
J.
Zhu
,
Y.
Sun
,
K.
Ramadoss
,
S. S.
Nonnenmann
,
N.
Yu
,
R.
Comin
,
K. M.
Rabe
,
S. K. R. S.
Sankaranarayanan
, and
S.
Ramanathan
,
Nature
553
,
68
(
2017
).
12.
T.
Minami
,
Semicond. Sci. Technol
20
,
S35
(
2005
).
13.
M.
Imada
,
A.
Fujimori
, and
Y.
Tokura
,
Rev. Mod. Phys.
70
,
1039
(
1998
).
14.
Y.
Cong
,
J.
Zhang
,
F.
Chen
, and
M.
Anpo
,
J. Phys. Chem. C
111
,
6976
(
2007
).
15.
X.
Chen
,
L.
Liu
,
P. Y.
Yu
, and
S. S.
Mao
,
Science
331
,
746
(
2011
).
16.
S. M.
Sze
and
M. K.
Lee
,
Semiconductor Devices: Physics and Technology
, 3rd edn (
John Wiley & Sons
,
Hoboken
,
2012
).
17.
F.
Keller
,
M. S.
Hunter
, and
D. L.
Robinson
,
J. Electrochem. Soc.
100
,
411
(
1953
).
18.
G. K.
Mor
,
O. K.
Varghese
,
M.
Paulose
,
K.
Shankar
, and
C. A.
Grimes
,
Solar Energy Mater. Solar Cells
90
,
2011
(
2006
).
19.
S.
Berger
,
A.
Ghicov
,
Y.-C.
Nah
, and
P.
Schmuki
,
Langmuir
25
,
4841
(
2009
).
20.
C. G.
Granqvist
,
Solar Energy Mater. Solar Cells
60
,
201
(
2000
).
21.
R.
Waser
,
R.
Dittmann
,
G.
Staikov
, and
K.
Szot
,
Adv. Mater.
21
,
2632
(
2009
).
22.
D.-H.
Kwon
,
K. M.
Kim
,
J. H.
Jang
,
J. M.
Jeon
,
M. H.
Lee
,
G. H.
Kim
,
X.-S.
Li
,
G.-S.
Park
,
B.
Lee
,
S.
Han
,
M.
Kim
, and
C. S.
Hwang
,
Nat. Nano.
5
,
148
(
2009
).
23.
Y.
Alivov
,
V.
Singh
,
Y.
Ding
,
L. J.
Cerkovkik
, and
P.
Nagpal
,
Nanoscale
6
,
10839
(
2014
).
24.
T.
Close
,
G.
Tulsyan
,
C. A.
Diaz
,
S. J.
Weinstein
, and
C.
Richter
,
Nature Nano
10
,
418
(
2015
).
25.
J. M.
Jiménez
,
G. R.
Bourret
,
T.
Berger
, and
K. P.
McKenna
,
J. Am. Chem. Soc.
138
,
15956
(
2016
).
26.
T.
Yajima
,
G.
Oike
,
T.
Nishimura
, and
A.
Toriumi
,
Phys. Stat. Sol. A
213
,
2196
(
2016
).
27.
S. R.
Hofstein
,
IEEE Trans. Elec. Dev.
ED-14
,
749
(
1967
).
28.
J.-Y.
Lee
,
J. G.
Lee
,
S.-H.
Lee
,
M.
Seo
,
L.
Piao
,
J. H.
Bae
,
S. Y.
Lim
,
Y. J.
Park
, and
T. D.
Chung
,
Nature Commun.
4
,
2766
(
2013
).
29.
E. H.
Hauri
,
A. M.
Shaw
,
J.
Wang
,
J. E.
Dixon
,
P. L.
King
, and
C.
Mandeville
,
Chem. Geol.
235
,
352
(
2006
).
30.
T.
Tiedje
,
C. R.
Wronski
,
B.
Abeles
, and
J. M.
Cebulka
,
Solar Cells
2
,
301
(
1980
).
31.
H.
Senoh
,
N.
Takeichi
,
T.
Kiyobayashi
,
H.
Tanaka
,
H. T.
Takeshita
,
T.
Oishi
, and
N.
Kuriyama
,
J. Alloys Comp.
404-406
,
47
(
2005
).
32.
E.
Cartier
,
J. H.
Stathis
, and
D. A.
Buchanan
,
Appl. Phys. Lett.
63
,
1510
(
1993
).
33.
L.-G.
Ekedahl
,
M.
Eriksson
, and
I.
Lundström
,
Acc. Chem. Res.
31
,
249
(
1998
).
34.
L. M.
Liu
,
S. Q.
Wang
, and
H. Q.
Ye
,
J. Phys.: Condens. Matter
17
,
5335
(
2005
).
35.
D. J.
Carrad
,
A. B.
Mostert
,
A. R.
Ullah
,
A. M.
Burke
,
H. J.
Joyce
,
H. H.
Tan
,
C.
Jagadish
,
P.
Krogstrup
,
J.
Nygard
,
P.
Meredith
, and
A. P.
Micolich
,
Nano Lett
17
,
827
(
2017
).
36.
C. J.
Wan
,
L. Q.
Zhu
,
Y. H.
Liu
,
P.
Feng
,
Z. P.
Liu
,
H. L.
Cao
,
P.
Xiao
,
Y.
Shi
, and
Q.
Wan
,
Adv. Mater.
28
,
3557
(
2016
).