Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “charge clusters”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

A Simple Mechanism for Fractionating Oxygen Isotopes in the Solar Nebula

Lightning in the Solar Nebula is caused by the tribo-electric charging of dust grains carried by massive turbulent flows and driven by the accretion energy in the disk: it has long been one agent assumed responsible for the formation of chondrules. The degree to which charge separation can occur is dependent upon a number of factors, including the concentration of radioactive sources and the total level of ionization in the nebula, and these factors determine the maximum energy likely to be released by a single bolt. While chondrule formation requires a massive discharge, even a small lightning bolt can vaporize grains in the ionized discharge channel. Experimental studies have shown that silica, iron silicate and iron oxide grains formed from a high voltage discharge in hydrogen rich gas containing some oxygen produces solids that are enriched in O-17 and O-18 relative to the composition of the starting gas. Vaporization of silicates produces SiO, metal and free oxygen atoms in each discharge and these species will immediately begin to recondense from the hot plasma. Freshly condensed grains are incrementally enriched in heavy oxygen while the gas is enriched in O-16. Repeated evaporation and condensation of silicates in continuously occurring lightning discharges will monotonically increase the fractionation of oxygen isotopes between the O-17 and O-18 rich dust and the O-16 rich gas. The first mass independently fractionated refractory oxide particles were produced in the lab following the condensation of a flowing gas mixture containing variable amounts of hydrogen, silane, pentacarbonyl iron and oxygen that passed through a high voltage discharge powered by a Tesla coil. While the exact chemical pathway is still uncertain, the most probable reaction mechanisms involve oxidation of the growing refractory clusters by O3, OH or O atoms. This model has some interesting consequences for chemical processes in the early solar nebula. Chemical fractionation of recondensed dust evaporated via lightning discharges should be strongly time dependent. At earlier times, the accretion rate is maximal, thus driving strong turbulence, energetic grain-grain collisions, tribo-electric charging and charge separation, leading to frequent, powerful lightning discharges. As the accretion rate diminishes, turbulence decreases and lightning discharges will become both less powerful and less frequent, thus decreasing the rate of dust-gas fractionation. The most rapid increase in the formation of O-16 poor dust will occur early in nebular history. Generation of fractionated dust should be distributed throughout the inner disk. Once condensed, grain dispersal would average out any significant isotopic anomalies within the inner disk.

Nuth, Joseph A., III↗
DOE OSTI2024

Importance of Electron Correlation on the Geometry and Electronic Structure of [2Fe–2S] Systems: A Benchmark Study of the [Fe 2 S 2 (SCH 3 ) 4 ] 2–,3–,4– , [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S- p -tol) 4 ] 2– , and [Fe 2 S 2 (S- o -xyl) 4 ] 2– Complexes

Iron–sulfur clusters are crucial for biological electron transport and catalysis. Obtaining accurate geometries, energetics, manifolds of their excited electronic states, and reduction energies is important to understand their role in these processes. Here, using a [2Fe–2S] model complex with Fe II and Fe III oxidation states, which leads to different charges, i.e., [Fe 2 S 2 (SMe) 4 ] 2–,3–,4– , we benchmarked a variety of computational methodologies ranging from density functional theory (DFT) to post-Hartree–Fock methods, including complete active space self-consistent field (CASSCF), multireference configuration interaction, the second-order N-electron valence state perturbation theory (NEVPT2), and the linearized integrand approximation of adiabatic connection (AC0) approaches. Additionally, we studied three experimentally well-characterized complexes, [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S-o-tol) 4 ] 2– , and [Fe 2 S 2 (S-o-xyl) 4 ] 2– , via DFT methods. We conclude that the dynamic electron correlation is important for accurately predicting the geometry of these complexes. Broken symmetry (BS) DFT correctly predicts experimental geometries of low-spin multiplicity, while CASSCF does not. However, BS-DFT significantly overestimates the difference between the low- and high-spin electronic states for a given oxidation state. At the same time, CASSCF underestimates it but provides relative energies closer to the reference NEVPT2 results. Finally, AC0 provides energetics of NEVPT2 quality with the additional advantage of being able to use large CASSCF sizes. NEVPT2 gives the best estimates of the Fe III /Fe III → Fe II /Fe III (4.27 eV) and Fe II /F III → Fe II /F II (7.72 eV) reduction energies. The results provide insight into the electronic structure of these complexes and assist in the understanding of their physical properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Microclusters of Kinked Silicon Nanowires Synthesized by a Recyclable Iodide Process for High-Performance Lithium-Ion Battery Anodes

Silicon shows great promise as a high-capacity anode material for lithium-ion batteries. Nanostructuring silicon minimizes the impact of fracturing during charging and discharging. However, synthesizing nanostructured silicon typically requires complicated procedures with high manufacturing costs. Additionally, these complicated procedures typically have poor secondary particle formation, a requirement to achieve a high tap density. Here, a cost-effective synthesis procedure which generates nearly ideal secondary particle clusters of nanostructured silicon is developed. The cost-effectiveness is a result of the in operando generation of silicon iodide from iodine gas and low-grade silicon microparticle. Decomposition of silicon iodide into crystalline silicon and iodine gas enables recycling of the iodine gas, allowing for near full reutilization of iodine in the following cycles. Additionally, the optimal nanostructures and microstructures of silicon synthesized by the recyclable iodide decomposition reaction enables 83.6% capacity retention over 1000 cycles. The good performance is a result of well-maintained morphology during cycling, enabling reutilization of the solid electrolyte interphase.

25 ENERGY STORAGE↗
DOE OSTI2025

Structure–property relations of binary ferrite melts

Molten ferrite systems are used in the smelting and refining processes in steelmaking, to reduce the loss of metals in slags and to accelerate reaction rates. Here, high-energy x-ray diffraction experiments have been performed on aerodynamically levitated molten spheres of 43BaO–57FeO X and 43SrO–57FeO X at 1873 K using laser beam heating. The composition was varied within the range of x = 1–1.5 by changing the oxygen partial pressure of the levitation gas. The corresponding x-ray pair distribution functions have been interpreted using empirical potential structure refinement (EPSR) modeling. In oxygen-rich melts (x = 1.5), our EPSR models indicate very similar structures for the different alkaline-earth liquids, with both the Ba–O and Sr–O coordination numbers to be ∼8.4 and the total Fe–O coordination numbers ∼5.7. However, our models show that in reducing environments, the Fe 3+ and Fe 2+ ions exhibit very different behaviors in the Ba- and Sr-ferrite liquids. In the Ba-ferrite melt, the Fe 3+ –O coordination number decreases from 5.7 (at x = 1.5) to 5.2 (at x = 1.07), whereas Fe 2+ –O remains constant at ∼5.0 across the same compositional range. In the Sr melts, both the Fe 2+ –O and Fe 3+ –O coordination numbers rise from ∼5.7 (at x = 1.5) to 6.3 (at x = 1.07). All models show the structures to be heterogeneous with intertwined nanometer sized clusters or channels of Ba/Sr–O and Fe–O polyhedra that grow as oxygen content is reduced. Changes in the viscosity and electrical properties are interpreted in terms of the number of bridging and non-bridging oxygens associated with FeO 4 tetrahedra and concentration of charge carriers, respectively.

Benmore, Chris J. [Argonne National Laboratory (AN↗
DOE OSTI2023

Nonlocal features of the spin-orbit exciton in Kitaev materials

A comparative resonant inelastic x-ray scattering (RIXS) study of three well-known Kitaev materials is presented: 𝛼−Li 2 ⁢IrO 3 , Na 2 ⁢IrO 3 , and 𝛼−RuCl 3 . Despite similar low-energy physics, these materials show distinct electronic properties, such as the large difference in the size of the charge gap. The RIXS spectra of the spin-orbit exciton for these materials show remarkably similar three-peak features, including sharp low energy peak (peak A) as well as transitions between 𝑗 eff = 1/2 and 𝑗 eff = 3/2 states. Comparison of experimental spectra with cluster calculations reveals that the observed three-peak structure reflects the significant role that nonlocal physics plays in the electronic structure of these materials. Here, in particular, the low-energy peak A arises from a holon-doublon pair rather than a conventional particle-hole exciton as proposed earlier. Our study suggests that while spin-orbit assisted Mott insulator is still the best description for these materials, electron itinerancy cannot be ignored when formulating low-energy Hamiltonian of these materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1984

The nature of the metal-CO interaction and bonding

The adsorption of CO on metal surfaces is represented by molecular orbital cluster models of CO at an on top site and adsorbed normal to the surface carbon end down. Ab initio SCF and MCSCF calculations are performed for several clusters. The new constrained space orbital variation CSOV approach is used to analyze the bonding and to compare CO adsorption on Al, representative of sp metals, with that on Cu, representative of transition metals. There is a large repulsion between the superposed free CO and metal charge distributions which is considerably smaller for Cu than for Al because there are fewer valence sigma electrons for Cu than for Al. The CSOV analysis shows that the metal to CO pi donation is much more important than the CO to metal sigma donation. It is also shown that for Cu, the d pi contribution to the metal pi donation is larger than the valence 4p pi contribution. The d pi donation is compared between Fe, Ni, and Cu and this donation and the metal-CO interaction are found to be different in the order Fe greater than Ni greater than Cu.

Bagus, P. S.↗
NASA NTRS2002

VIZARD: analysis of Affymetrix Arabidopsis GeneChip data

SUMMARY: The Affymetrix GeneChip Arabidopsis genome array has proved to be a very powerful tool for the analysis of gene expression in Arabidopsis thaliana, the most commonly studied plant model organism. VIZARD is a Java program created at the University of California, Berkeley, to facilitate analysis of Arabidopsis GeneChip data. It includes several integrated tools for filtering, sorting, clustering and visualization of gene expression data as well as tools for the discovery of regulatory motifs in upstream sequences. VIZARD also includes annotation and upstream sequence databases for the majority of genes represented on the Affymetrix Arabidopsis GeneChip array. AVAILABILITY: VIZARD is available free of charge for educational, research, and not-for-profit purposes, and can be downloaded at http://www.anm.f2s.com/research/vizard/ CONTACT: moseyko@uclink4.berkeley.edu.

Non-NASA Center↗
DOE OSTI2022

Domain architecture and catalysis of the Staphylococcus aureus fatty acid kinase

Fatty acid kinase (Fak) is a two-component enzyme that generates acyl-phosphate for phospholipid synthesis. Fak consists of a kinase domain protein (FakA) that phosphorylates a fatty acid enveloped by a fatty acid binding protein (FakB). The structural basis for FakB function has been established, but little is known about FakA. Here, we used limited proteolysis to define three separate FakA domains: the amino terminal FakA_N, the central FakA_L, and the carboxy terminal FakA_C. The isolated domains lack kinase activity, but activity is restored when FakA_N and FakA_L are present individually or connected as FakA_NL. The X-ray structure of the monomeric FakA_N captures the product complex with ADP and two Mg 2+ ions bound at the nucleotide site. The FakA_L domain encodes the dimerization interface along with conserved catalytic residues Cys240, His282, and His284. AlphaFold analysis of FakA_L predicts the catalytic residues are spatially clustered and pointing away from the dimerization surface. Furthermore, the X-ray structure of FakA_C shows that it consists of two subdomains that are structurally related to FakB. Analytical ultracentrifugation demonstrates that FakA_C binds FakB, and site-directed mutagenesis confirms that a positively charged wedge on FakB meshes with a negatively charged groove on FakA_C. Finally, small angle X-ray scattering analysis is consistent with freely rotating FakA_N and FakA_C domains tethered by flexible linkers to FakA_L. These data reveal specific roles for the three independently folded FakA protein domains in substrate binding and catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Searching for metastable particles using graph computing

The reconstruction of charged particle trajectories at the Large Hadron Collider and future colliders relies on energy depositions in sensors placed at distances ranging from a centimeter to a meter from the colliding beams. We propose a method of detecting charged particles that decay invisibly after traversing a short distance of about 25 cm inside the experimental apparatus. One of the decay products may constitute the dark matter known to be 84% of all matter at galactic and cosmological distance scales. Our method uses graph computing to cluster spacepoints recorded by two-dimensional silicon pixel sensors into mathematically-defined patterns. The algorithm may be implemented on silicon-based integrated circuits using field-programmable gate array technology to augment or replace traditional computing platforms.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Fast Proton Transfer and Hydrogen Evolution Reactivity Mediated by [Co 13 C 2 (CO) 24 ] 4–

A common approach to speeding up proton transfer (PT) by molecular catalysts is manipulation of the secondary coordina-tion sphere with proton relays and these enhance overall reaction rates by orders of magnitude. In contrast, heterogenous electrocatalysts have band structures that promote facile PT concerted with electron transfer (ET), known as the Volmer mechanism. Here in this paper, we show that [Co 13 C 2 (CO) 24 ] 4– , containing multiple Co-Co bonds to statistically enhance observed rates of PT, promotes PT on the order of 2.3 × 10 9 M -1 s -1 which suggests a diffusion-limited rate. The fast ET and PT chemistry is at-tributed to the delocalized electronic structure of [Co 13 C 2 (CO) 24 ] 4– . Electrochemical characterization of [Co 13 C 2 (CO) 24 ] 4– in the presence and absence of protons reveals ET kinetics and diffusion behavior similar to other small clusters such as nano-materials and fullerenes.

08 HYDROGEN↗
DOE OSTI2020

Genesis of the periodic lattice distortions in the charge density wave state of 2 H- TaSe 2

Using high-energy x-ray diffraction and large-scale three-dimensonal structure modeling, we reveal the evolution of local lattice distortions in the archetypal charge density wave (CDW) system 2 H- TaSe 2 as it is cooled from its normal state above room temperature down to temperature where Ta atoms form a long-range ordered superstructure. In particular, we find that the structure is formed via a gradual in-plane clustering of Ta atoms exhibiting unusually short bonding distances already near room temperature, and not via a cooperative distortion of Ta atomic planes taking place at the critical ordering temperature. Our findings clarify the debated local symmetry and magnitude of the periodic lattice distortions in the CDW state of 2 H- TaSe 2 , and emphasize the role of locally correlated lattice distortions as a precursor of CDWs in low-dimensional solids. Further, we demonstrate an efficient experimental approach to study them.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Measurement of $\nu_\mu$ CC Interactions With Two-Proton Final State in MINERvA

This dissertation presents a measurement of charged–current (CC) muon–neutrino interactions with exactly two protons and no pions in the final state (CC~$2p\,0\pi$), using data collected by the MINERvA detector in the NuMI medium–energy beam at Fermilab. Such two–proton topologies are a sensitive probe of nuclear dynamics in the few–GeV regime, including multi–nucleon correlations (npnh, notably $2p2h$) and intranuclear final–state interactions (FSI) such as pion absorption and nucleon rescattering. A precise experimental characterization of these processes is essential both for neutrino–interaction theory and for reducing systematic uncertainties in oscillation experiments that rely on accurate modeling of neutrino–nucleus interactions. Events are selected by requiring a $\nu_\mu$ CC interaction with a reconstructed $\mu^-$ and two proton tracks originating from a common vertex in MINERvA’s finely segmented scintillator tracker, with no reconstructed mesons. Muon charge and momentum are constrained by matching to the MINOS Near Detector, while proton identification exploits energy–loss profiles and stopping–proton features. Backgrounds from pion–producing channels that enter the signal region through FSI or reconstruction effects are constrained with data–driven sidebands (Michel–electron and isolated–cluster “blob” samples) and tuned via a simultaneous fit across signal and sideband regions. To correct detector resolution and acceptance effects, the analysis employs iterative Bayesian unfolding with extensive validation: statistical pseudo–experiments, and robustness checks against generator systematic “universes” and additional strong shape warps. Single–differential cross sections are reported for three observables tailored to the two–proton final state: the opening–angle cosine $\cos\!\left(\theta_{pp}\right)$, the leading–proton momentum, and the subleading–proton momentum. Systematic uncertainties include contributions from neutrino flux, interaction modeling (e.g., npnh and resonance parameters, pion FSI), and detector response (calibration, reconstruction efficiencies). The resulting distributions provide targeted constraints on the interplay of multi–nucleon dynamics and FSI that shape CC~$2p\,0\pi$ final states on hydrocarbon. Comparisons to modern GENIE–based simulations highlight kinematic regions where model components require refinement. These measurements thus inform generator tuning and improve the reliability of neutrino–energy reconstruction strategies for current and future long–baseline oscillation programs.

Syrotenko, Vladyslav S. [Tufts U.]↗
DOE OSTI2021

Role of Fe Doping on Local Structure and Electrical and Magnetic Properties of PbTiO 3

In this work, the local structure and multiferroic properties of iron-doped lead titanate (PbTi 1-x Fe x O 3-δ ) samples was investigated over the entire composition range (x = 0-1). Inherent polarization in PbTiO 3 decreases due to Fe 3+ incorporation up to the solubility limit (x similar to 0.3), although homogeneous doping persists only up to x = 0.1. Ti prefers highly distorted oxygen octahedra for any x value, while Fe prefers more symmetric O-deficient polyhedra (Fe-O n ). The charge compensating oxygen vacancies induce local tilting of the Fe-O n polyhedra beyond a critical x value (x ≥ 0.2), promoting magnetic interaction between two adjacent Fe atoms. The strain induced by local heterogeneity could act as a coupling force between magnetic and ferroelectric properties. Fe-rich clusters evolve into ferromagnetic PbFe 12 O 19 with increased Fe doping. PbTi 1-x Fe x O 3-δ (x ≥ 0.3) samples therefore have separate origins for the ferroelectric (PbTi 1-x Fe x O 3-δ ) and magnetic (PbFe 12 O 19 ) phases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Determination of Intermolecular Distances in Dilute Solutions of Macroions and Their Direct Correlations with Self-Assembly and Macrophase Transitions

Here, this work demonstrates new attempts to determine the intermolecular distances between charged macroionic solutes in their dilute solutions, which regulate the solute’s microphase (self-assembly) and macrophase transitions. Small-angle X-ray scattering (SAXS) and analytical ultracentrifuge (AUC) techniques were applied to determine the intermolecular distances for the 2.42 nm-sized, spherical uranyl peroxide molecular cluster {U 60 } in their self-assembled states in dilute aqueous solution and concentrated phases, respectively. The counterion-mediated attraction among {U 60 } leads to characteristic, inter-{U 60 } distances in solutions, which increase gradually with decreasing the strength of introduced counterions (e.g., lower valency). The gradual increment of inter-{U 60 } distance demonstrates a nice correlation with the macroscopic phase transitions of {U 60 }, from single crystals to concentrated fluids containing rigid 2-D sheets, then to dilute solutions containing a small amount of standalone, floating 2-D sheets of {U 60 }, and finally to dilute solutions containing a limited amount of self-assembled single-layered, spherical blackberry structures of different sizes from the bending of more flexible 2-D sheets.

36 MATERIALS SCIENCE↗
NASA NTRS2006

NASA Radiation Track Image GUI for Assessing Space Radiation Biological Effects

The high-charge high-energy (HZE) ion components of the galactic cosmic rays when compared to terrestrial forms of radiations present unique challenges to biological systems. In this paper we present a deoxyribonucleic acid (DNA) breakage model to visualize and analyze the impact of chromatin domains and DNA loops on clustering of DNA damage from X rays, protons, and HZE ions. Our model of DNA breakage is based on a stochastic process of DNA double-strand break (DSB) formulation that includes the amorphous model of the radiation track and a polymer model of DNA packed in the cell nucleus. Our model is a Monte-Carlo simulation based on a randomly located DSB cluster formulation that accomodates both high- and low-linear energy transfer radiations. We demonstrate that HZE ions have a strong impact on DSB clustering, both along the chromosome length and in the nucleus volume. The effects of chromosomal domains and DNA loops on the DSB fragment-size distribution and the spatial distribution of DSB in the nucleus were studied. We compare our model predictions with the spatial distribution of DSB obtained from experiments. The implications of our model predictions for radiation protection are discussed.

Ponomarev, Artem L.↗
DOE OSTI2025

Impact of solvation on the electronic resonances in uracil

Interactions of low-energy electrons with the DNA and RNA nucleobases are known to form metastable states, known as electronic resonances. In this work, we study electron attachment to solvated uracil, an RNA nucleobase, using the orbital stabilization method at the Equation of Motion-Coupled Cluster for Electron Affinities with Singles and Doubles (EOM-EA-CCSD) level of theory with the Effective Fragment Potential (EFP) solvation method. We benchmarked the approach using multireference methods, as well as by comparing EFP and full quantum calculations. The impact of solvation on the first one particle (1p) shape resonance, formed by electron attachment to the π* LUMO orbital, as well as the first two particle one hole (2p1h) resonance, formed by electron attachment to neutral uracil's π–π* excited state, was investigated. We used molecular dynamics simulations for solvent configurations and applied charge stabilization technique-based biased sampling to procure configurations adequate to cover the entire range of the electron attachment energy distribution. The electron attachment energy in solution is found to be distributed over a wide range of energies, between 4.6 eV to 6.8 eV for the 2p1h resonance, and between −0.1 eV to 2 eV for the 1p resonance. The solvent effects were similar for the two resonances, indicating that the exact electron density of the state is not as important as the solvent configurations. Multireference calculations extended the findings showing that solvation effects are similar for the lowest four resonances, further indicating that the specific solute electron density is not as important, but rather the water configurations play the most important role in solvation effects. Lastly, by comparing bulk solvation to clusters of uracil with a few water molecules around it, we find that the impact of microsolvation is very different from that of bulk solvation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Exploring the Role of Inorganic and Organic Interfaces on CO 2 and CH 4 Partitioning: Case Study of Silica, Illite, Calcite, and Kerogen Nanopores on Gas Adsorption and Nanoscale Transport Behaviors

The adsorption, partitioning, and diffusion of CO 2 and CH 4 in organic–inorganic pores composed of kerogen with silica, illite, or calcite are studied using density functional theory and classical molecular dynamics simulations. The adsorption and partitioning behavior of CO 2 and CH 4 molecules is found to be a function of chemistry of the solid interface, pore size, and surface area. CO 2 molecules are preferentially adsorbed compared to CH 4 molecules on silica, illite, calcite, and kerogen surfaces due to significant contributions of electrostatic interactions with the atoms of the associated surfaces. CO 2 and CH 4 molecules display higher adsorption energy on calcite and silica compared to those on illite and kerogen because of enhanced interactions with the positively charged calcium ions on the calcite surface and the hydroxyl group (-OH) on the silica surface. Enhanced internal nanopores and the availability of adsorption active sites in illite and kerogen matrices aid in higher partitioning of the gas molecules into these pores. CO 2 has a significant influence on the clustering and swelling of kerogen fragments compared to CH 4 . Higher CO 2 adsorption on inorganic and kerogen surfaces results in lower overall self-diffusivities compared to CH 4 . This study illustrates surface and pore size controls on the adsorption, partitioning, and self-diffusivities of CO 2 and CH 4 at conditions relevant for subsurface energy applications. The chemistry of the surfaces along with the availability of internal pores and active adsorption sites influences the adsorption behavior of gases. Our results suggest that CO 2 stored in depleted gas reservoirs may preferentially inhabit nanoconfined pores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

The role of water vapor during the synthesis of hydrogen doped In 2 O 3

Hydrogen-doped In 2 O 3 synthesized using magnetron sputtering or atomic layer deposition in the presence of water vapor has high transparency in the near infra-red region and mobility values that are almost three to four times of those in commercially used In 2 O 3 :Sn at a much lower carrier concentration. However, simple questions like what happens when water molecules enter In 2 O 3 or why the dominant charge carrier changes from an oxygen vacancy in In 2 O 3 :Sn to H + in In 2 O 3 are not clear. Using hybrid functional based density functional theory calculations, in this study we show that water molecules spontaneously split into H + and OH – . The H + interstitial bonds with an oxygen in In 2 O 3 , and the OH – interstitial occupies a vacant lattice site or is trapped by an oxygen vacancy. Thus, in agreement with experimental findings, our results suggest that the concentration of oxygen vacancies decreases as the partial pressure of water increases. The charge transition level, ϵ ( – / + ) , of a hydrogen interstitial and the Fermi level at which the formation energies of OH – and H + are equal lie above the conduction band minimum, suggesting that the H + interstitial has a lower formation energy. In addition, water molecules do not exhibit any propensity to form clusters, meaning that H + ions can diffuse and remain spatially separated from each other.

36 MATERIALS SCIENCE↗