skip to main content


Title: First-principles description of oxygen self-diffusion in rutile TiO 2 : assessment of uncertainties due to enthalpy and entropy contributions
Properties related to transport such as self-diffusion coefficients are relevant to fuel cells, electrolysis cells, and chemical/gas sensors. Prediction of self-diffusion coefficients from first-principles involves precise determination of both enthalpy and entropy contributions for point defect formation and migration. We use first-principles density functional theory to estimate the self-diffusion coefficient for neutral O0i and doubly ionized O i 2− interstitial oxygen in rutile TiO 2 and compare the results to prior isotope diffusion experiments. In addition to formation and migration energy, detailed estimates of formation and migration entropy incorporating both vibrational and ionization components are included. Distinct migration pathways, both based on an interstitialcy mechanism, are identified for O0i and O i 2− . These result in self-diffusion coefficients that differ by several orders of magnitude, sufficient to resolve the charge state of the diffusing species to be O i 2− in experiment. The main sources of error when comparing computed parameters to those obtained from experiment are considered, demonstrating that uncertainties due to computed defect formation and migration entropies are comparable in magnitude to those due to computed defect formation and migration energies. Even so, the composite uncertainty seems to limit the accuracy of first-principles calculations to within a factor of ±10 3 , demonstrating that direct connections between computation and experiment are now increasingly possible.  more » « less
Award ID(s):
1709327
NSF-PAR ID:
10090772
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Physical Chemistry Chemical Physics
Volume:
20
Issue:
25
ISSN:
1463-9076
Page Range / eLocation ID:
17448 to 17457
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Hematite (α-Fe 2 O 3 ) is a promising transition metal oxide for various energy conversion and storage applications due to its advantages of low cost, high abundance, and good chemical stability. However, its low carrier mobility and electrical conductivity have hindered the wide application of hematite-based devices. Fundamentally, this is mainly caused by the formation of small polarons, which show conduction through thermally activated hopping. Atomic doping is one of the most promising approaches for improving the electrical conductivity in hematite. However, its impact on the carrier mobility and electrical conductivity of hematite at the atomic level remains to be illusive. In this work, through a kinetic Monte-Carlo sampling approach for diffusion coefficients combined with carrier concentrations computed under charge neutrality conditions, we obtained the electrical conductivity of the doped hematite. We considered the contributions from individual Fe–O layers, given that the in-plane carrier transport dominates. We then studied how different dopants impact the carrier mobility in hematite using Sn, Ti, and Nb as prototypical examples. We found that the carrier mobility change is closely correlated with the local distortion of Fe–Fe pairs, i.e. the more stretched the Fe–Fe pairs are compared to the pristine systems, the lower the carrier mobility will be. Therefore, elements which limit the distortion of Fe–Fe pair distances from pristine are more desired for higher carrier mobility in hematite. The calculated local structure and pair distribution functions of the doped systems have remarkable agreement with the experimental EXAFS measurements on hematite nanowires, which further validates our first-principles predictions. Our work revealed how dopants impact the carrier mobility and electrical conductivity of hematite and provided practical guidelines to experimentalists on the choice of dopants for the optimal electrical conductivity of hematite and the performance of hematite-based devices. 
    more » « less
  2. Resonant tunneling diodes (RTDs) have come full-circle in the past 10 years after their demonstration in the early 1990s as the fastest room-temperature semiconductor oscillator, displaying experimental results up to 712 GHz and fmax values exceeding 1.0 THz [1]. Now the RTD is once again the preeminent electronic oscillator above 1.0 THz and is being implemented as a coherent source [2] and a self-oscillating mixer [3], amongst other applications. This paper concerns RTD electroluminescence – an effect that has been studied very little in the past 30+ years of RTD development, and not at room temperature. We present experiments and modeling of an n-type In0.53Ga0.47As/AlAs double-barrier RTD operating as a cross-gap light emitter at ~300K. The MBE-growth stack is shown in Fig. 1(a). A 15-μm-diam-mesa device was defined by standard planar processing including a top annular ohmic contact with a 5-μm-diam pinhole in the center to couple out enough of the internal emission for accurate free-space power measurements [4]. The emission spectra have the behavior displayed in Fig. 1(b), parameterized by bias voltage (VB). The long wavelength emission edge is at  = 1684 nm - close to the In0.53Ga0.47As bandgap energy of Ug ≈ 0.75 eV at 300 K. The spectral peaks for VB = 2.8 and 3.0 V both occur around  = 1550 nm (h = 0.75 eV), so blue-shifted relative to the peak of the “ideal”, bulk InGaAs emission spectrum shown in Fig. 1(b) [5]. These results are consistent with the model displayed in Fig. 1(c), whereby the broad emission peak is attributed to the radiative recombination between electrons accumulated on the emitter side, and holes generated on the emitter side by interband tunneling with current density Jinter. The blue-shifted main peak is attributed to the quantum-size effect on the emitter side, which creates a radiative recombination rate RN,2 comparable to the band-edge cross-gap rate RN,1. Further support for this model is provided by the shorter wavelength and weaker emission peak shown in Fig. 1(b) around = 1148 nm. Our quantum mechanical calculations attribute this to radiative recombination RR,3 in the RTD quantum well between the electron ground-state level E1,e, and the hole level E1,h. To further test the model and estimate quantum efficiencies, we conducted optical power measurements using a large-area Ge photodiode located ≈3 mm away from the RTD pinhole, and having spectral response between 800 and 1800 nm with a peak responsivity of ≈0.85 A/W at  =1550 nm. Simultaneous I-V and L-V plots were obtained and are plotted in Fig. 2(a) with positive bias on the top contact (emitter on the bottom). The I-V curve displays a pronounced NDR region having a current peak-to-valley current ratio of 10.7 (typical for In0.53Ga0.47As RTDs). The external quantum efficiency (EQE) was calculated from EQE = e∙IP/(∙IE∙h) where IP is the photodiode dc current and IE the RTD current. The plot of EQE is shown in Fig. 2(b) where we see a very rapid rise with VB, but a maximum value (at VB= 3.0 V) of only ≈2×10-5. To extract the internal quantum efficiency (IQE), we use the expression EQE= c ∙i ∙r ≡ c∙IQE where ci, and r are the optical-coupling, electrical-injection, and radiative recombination efficiencies, respectively [6]. Our separate optical calculations yield c≈3.4×10-4 (limited primarily by the small pinhole) from which we obtain the curve of IQE plotted in Fig. 2(b) (right-hand scale). The maximum value of IQE (again at VB = 3.0 V) is 6.0%. From the implicit definition of IQE in terms of i and r given above, and the fact that the recombination efficiency in In0.53Ga0.47As is likely limited by Auger scattering, this result for IQE suggests that i might be significantly high. To estimate i, we have used the experimental total current of Fig. 2(a), the Kane two-band model of interband tunneling [7] computed in conjunction with a solution to Poisson’s equation across the entire structure, and a rate-equation model of Auger recombination on the emitter side [6] assuming a free-electron density of 2×1018 cm3. We focus on the high-bias regime above VB = 2.5 V of Fig. 2(a) where most of the interband tunneling should occur in the depletion region on the collector side [Jinter,2 in Fig. 1(c)]. And because of the high-quality of the InGaAs/AlAs heterostructure (very few traps or deep levels), most of the holes should reach the emitter side by some combination of drift, diffusion, and tunneling through the valence-band double barriers (Type-I offset) between InGaAs and AlAs. The computed interband current density Jinter is shown in Fig. 3(a) along with the total current density Jtot. At the maximum Jinter (at VB=3.0 V) of 7.4×102 A/cm2, we get i = Jinter/Jtot = 0.18, which is surprisingly high considering there is no p-type doping in the device. When combined with the Auger-limited r of 0.41 and c ≈ 3.4×10-4, we find a model value of IQE = 7.4% in good agreement with experiment. This leads to the model values for EQE plotted in Fig. 2(b) - also in good agreement with experiment. Finally, we address the high Jinter and consider a possible universal nature of the light-emission mechanism. Fig. 3(b) shows the tunneling probability T according to the Kane two-band model in the three materials, In0.53Ga0.47As, GaAs, and GaN, following our observation of a similar electroluminescence mechanism in GaN/AlN RTDs (due to strong polarization field of wurtzite structures) [8]. The expression is Tinter = (2/9)∙exp[(-2 ∙Ug 2 ∙me)/(2h∙P∙E)], where Ug is the bandgap energy, P is the valence-to-conduction-band momentum matrix element, and E is the electric field. Values for the highest calculated internal E fields for the InGaAs and GaN are also shown, indicating that Tinter in those structures approaches values of ~10-5. As shown, a GaAs RTD would require an internal field of ~6×105 V/cm, which is rarely realized in standard GaAs RTDs, perhaps explaining why there have been few if any reports of room-temperature electroluminescence in the GaAs devices. [1] E.R. Brown,et al., Appl. Phys. Lett., vol. 58, 2291, 1991. [5] S. Sze, Physics of Semiconductor Devices, 2nd Ed. 12.2.1 (Wiley, 1981). [2] M. Feiginov et al., Appl. Phys. Lett., 99, 233506, 2011. [6] L. Coldren, Diode Lasers and Photonic Integrated Circuits, (Wiley, 1995). [3] Y. Nishida et al., Nature Sci. Reports, 9, 18125, 2019. [7] E.O. Kane, J. of Appl. Phy 32, 83 (1961). [4] P. Fakhimi, et al., 2019 DRC Conference Digest. [8] T. Growden, et al., Nature Light: Science & Applications 7, 17150 (2018). [5] S. Sze, Physics of Semiconductor Devices, 2nd Ed. 12.2.1 (Wiley, 1981). [6] L. Coldren, Diode Lasers and Photonic Integrated Circuits, (Wiley, 1995). [7] E.O. Kane, J. of Appl. Phy 32, 83 (1961). [8] T. Growden, et al., Nature Light: Science & Applications 7, 17150 (2018). 
    more » « less
  3. Abstract

    Many of the studies on the entropy‐stabilized oxide (Mg0.2Co0.2Ni0.2Cu0.2Zn0.2)O have been heavily application‐based. Previous works have studied effects of cation stoichiometry on the entropy‐driven reaction to form a single phase, but a fundamental exploration of the effects of anion stoichiometry and/or redox chemistry on electrical properties is lacking. Using near‐edge X‐ray absorption fine structure (NEXAFS) and electrical measurements, we show that oxidizing thin film samples of (Mg0.2Co0.2Ni0.2Cu0.2Zn0.2)O affects primarily the valence of Co, leaving the other cations in this high‐entropy system unchanged. This oxidation increases electrical conduction in these thin films, which occurs via small polaron hopping mediated by the Co valence shift from 2+ to a mixed 2+/3+ state. In parallel, we show that bulk samples sintered in an oxygen‐rich atmosphere have a lower activation energy for electrical conduction than those equilibrated in a nitrogen (reducing) atmosphere. Combining feasible defect compensation scenarios with electrical impedance measurements and NEXAFS data, we propose a self‐consistent interpretation of Co redox‐mediated small polaron conduction as the dominant method of charge transfer in this system.

     
    more » « less
  4. Abstract

    Entropy‐stabilized oxide (ESO) research has primarily focused on discovering unprecedented structures, chemistries, and properties in the single‐phase state. However, few studies discuss the impacts of entropy stabilization and secondary phases on functionality and in particular, electrical conductivity. To address this gap, electrical transport mechanisms in the canonical ESO rocksalt (Co,Cu,Mg,Ni,Zn)O are assessed as a function of secondary phase content. When single‐phase, the oxide conducts electrons via Cu+/Cu2+small polarons. After 2 h of heat treatment, Cu‐rich tenorite secondary phases form at some grain boundaries (GBs), enhancing grain interior electronic conductivity by tuning defect chemistry toward higher Cu+carrier concentrations. 24 h of heat treatment yields Cu‐rich tenorite at all GBs, followed by the formation of anisotropic Cu‐rich tenorite and equiaxed Co‐rich spinel secondary phases in grains, further enhancing grain interior electronic conductivity but slowing electronic transport across the tenorite‐rich GBs. Across all samples, the total electrical conductivity increases (and decreases reversibly) by four orders of magnitude with heat‐treatment‐induced phase transformation by tuning the grains’ defect chemistry toward higher carrier concentration and lower migration activation energy. This work demonstrates the potential to selectively grow secondary phases in ESO grains and at GBs, thereby tuning the electrical properties using microstructure design, nanoscale engineering, and heat treatment, paving the way to develop many novel materials.

     
    more » « less
  5. Abstract Knowledge of oxygen diffusion in garnet is crucial for a correct interpretation of oxygen isotope signatures in natural samples. A series of experiments was undertaken to determine the diffusivity of oxygen in garnet, which remains poorly constrained. The first suite included high-pressure (HP), nominally dry experiments performed in piston-cylinder apparatus at: (1) T = 1050–1600 °C and P = 1.5 GPa and (2) T = 1500 °C and P = 2.5 GPa using yttrium aluminum garnet (YAG; Y3Al5O12) cubes. Second, HP H2O-saturated experiments were conducted at T = 900 °C and P = 1.0–1.5 GPa, wherein YAG crystals were packed into a YAG + Corundum powder, along with 18O-enriched H2O. Third, 1 atm experiments with YAG cubes were performed in a gas-mixing furnace at T = 1500–1600 °C under Ar flux. Finally, an experiment at T = 900 °C and P = 1.0 GPa was done using a pyrope cube embedded into pyrope powder and 18O-enriched H2O. Experiments using grossular were not successful. Profiles of 18O/(18O+16O) in the experimental charges were analyzed with three different secondary ion mass spectrometers (SIMS): sensitive high-resolution ion microprobe (SHRIMP II and SI), CAMECA IMS-1280, and NanoSIMS. Considering only the measured length of 18O diffusion profiles, similar results were obtained for YAG and pyrope annealed at 900 °C, suggesting limited effects of chemical composition on oxygen diffusivity. However, in both garnet types, several profiles deviate from the error function geometry, suggesting that the behavior of O in garnet cannot be fully described as simple concentration-independent diffusion, certainly in YAG and likely in natural pyrope as well. The experimental results are better described by invoking O diffusion via two distinct pathways with an inter-site reaction allowing O to move between these pathways. Modeling this process yields two diffusion coefficients (D values) for O, one of which is approximately two orders of magnitude higher than the other. Taken together, Arrhenius relationships are:log⁡Dm2s-1=-7.2(±1.3)+(-321(±32)kJmol-12.303RT) for the slow pathway, andlog⁡Dm2s-1=-5.4(±0.7)+(-321(±20)kJmol-12.303RT) for the fast pathway. We interpret the two pathways as representing diffusion following vacancy and inter-stitial mechanisms, respectively. Regardless, our new data suggest that the slow mechanism is prevalent in garnet with natural compositions, and thus is likely to control the retentivity of oxygen isotopic signatures in natural samples. The diffusivity of oxygen is similar to Fe-Mn diffusivity in garnet at 1000–1100 °C and Ca diffusivity at 850 °C. However, the activation energy for O diffusion is larger, leading to lower diffusivities at P-T conditions characterizing crustal metamorphism. Therefore, original O isotopic signatures can be retained in garnets showing major element zoning partially re-equilibrated by diffusion, with the uncertainty caveat of extrapolating the experimental data to lower temperature conditions. 
    more » « less