Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion coefficients”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

On the Importance of Using Event-Specific Wave Diffusion Rates in Modeling Diffuse Electron Precipitation

A few to tens of keV electron precipitation that carries substantial energy source down to the upper atmosphere to create aurora is manifested as an important magnetosphere-ionosphere coupling process. The precipitation is usually caused by scattering processes associated with plasma waves in the magnetosphere. The scattering process is often quantified by wave diffusion rates that indicate how fast an electron is scattered. Global models commonly use diffusion coefficients that are derived from statistical wave models. However, due to the statistical nature, many localized, transient features could be smeared out. In this study, we investigate electron precipitation using event-specific diffusion coefficients that are obtained based on simultaneous in-situ measured/inferred, rather than statistical, chorus wave dynamics. We find that the application of the event-specific diffusion coefficients associated with a more dynamic and intense chorus wave model leads more electrons, particularly at several to tens of keV in the dawn-to-noon sector at L > 3, to precipitate than using statistical coefficients. Here, the new simulation roughly captures both the intensity and variability of the precipitating flux as detected by the NOAA/POES satellites. Ionospheric electron density in the lower E region (100–120 km) observed by the mid-latitude Millstone Hill radar is also much better reproduced, while the case using statistical diffusion coefficients underestimates the ionization rate. This study implies the importance of using event-specific diffusion rates in simulating the diffuse electron precipitation and understanding the magnetosphere-ionosphere coupling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Accelerated kinetic Monte Carlo: A case study; vacancy and dumbbell interstitial diffusion traps in concentrated solid solution alloys

Vacancy and self-interstitial atomic diffusion coefficients in concentrated solid solution alloys can have a non-monotonic concentration dependence. Here, the kinetics of monovacancies and $\langle$100$\rangle$ dumbbell interstitials in Ni–Fe alloys are assessed using lattice kinetic Monte Carlo (kMC). The non-monotonicity is associated with superbasins, which impels using accelerated kMC methods. Detailed implementation prescriptions for first passage time analysis kMC (FPTA-kMC), mean rate method kMC (MRM-kMC), and accelerated superbasin kMC (AS-kMC) are given. The accelerated methods are benchmarked in the context of diffusion coefficient calculations. The benchmarks indicate that MRM-kMC underestimates diffusion coefficients, while AS-kMC overestimates them. In this application, MRM-kMC and AS-kMC are computationally more efficient than the more accurate FPTA-kMC. Our calculations indicate that composition dependence of migration energies is at the origin of the vacancy’s non-monotonic behavior. In contrast, the difference between formation energies of Ni–Ni, Ni–Fe, and Fe–Fe dumbbell interstitials is at the origin of their non-monotonic diffusion behavior. Additionally, the migration barrier crossover composition—based on the situation where Ni or Fe atom jumps have lower energy barrier than the other one—is introduced. KMC simulations indicate that the interplay between composition dependent crossover of migration energy and geometrical site percolation explains the non-monotonic concentration-dependence of atomic diffusion coefficients.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Diffusion and phase formation in the γ-uranium-technetium system

Phase formation in the U-Tc binary system at 800 °C was investigated using a diffusion couple experiment. Scanning electron microscopy (SEM) and energy dispersive X-ray spectroscopy (EDS) identified four novel potential intermetallic phases - U 7 Tc 3 , U 13 Tc 12 , U 3 Tc 5 , and UTc 4 . Diffusion coefficients were calculated for the intermetallic phases using the Boltzmann-Matano method and were respectively found to be – 120, 38.2, 15.6, and 1.51 × 10 −18 m 2 /s. Tc also exhibits a solid-solution phase with high penetration into the U with a diffusion coefficient of ∼ 10 −14 m 2 /s. Furthermore, these findings expand the number of known U-Tc phases and provide the first diffusion coefficients for the U-Tc system, and contribute valuable data to the broader field of actinide metallurgy.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Analysis of Silver Release from Furnace-Tested TRISO Particles

High temperature testing of intact TRISO particles previously irradiated in the AGR-5/6/7 experiment was performed in the Furnace for Irradiated TRISO Testing (FITT) to directly confirm diffusive release of silver and europium from intact TRISO particles. Testing was conducted from 1,100–1,600°C for exposure times up to 100 h to directly confirm silver through layer release below safety testing temperatures. The FITT analysis showed highest levels of silver release at 1,300–1,400°C which confirms athermal release behaviors previous observed in step-wise and transient safety tests. Additionally, release was non-uniform with some particles releasing a majority of their inventory while others did not appear to release silver under identical testing conditions, which was consistent with historic observations. An assessment of the effective silver diffusion coefficient in the SiC layer was conducted and indicated maximum values in the 1,300–1,400°C range. The magnitude of the calculated diffusion coefficients also exceeded currently accepted diffusion coefficients. The release behavior of europium from intact particles was also analyzed in FITT for at 1,450–1,550°C for 500 h to 750 h. Direct confirmation of europium release below safety testing temperatures was confirmed absent contributions from matrix release. The analysis indicated europium release follows a general Arrhenius behavior and suggests general uniform release and indicates irradiation conditions influence observed release response. Calculated diffusion kinetics agree well with historic experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Reduced diffusion and enhanced retention of multiple radionuclides from pore structure characterization of barrier materials for enhanced repository performance

Fluid flow and chemical transport in porous media are the macroscopic consequences of pore structure, which integrates geometry (e.g., pore size and surface area, pore-size distribution) and topology (e.g., pore connectivity). Low-permeability geological media whose pores are poorly interconnected will exhibit the characteristics of anomalous diffusion and sample size-dependent effective porosity, which will strongly impact long-term net diffusion and retention of radionuclides in geological repository settings involving different host rocks and barrier materials. A suite of innovative and complementary experimental approaches is utilized to study the microscopic pore structure and macroscopic fluid flow & chemical transport for a range of host rocks and barrier materials, in addition to standard clay minerals and reference rocks. With a particular focus on quantifying the presence and magnitude of “isolated” pores for a reduced effective porosity in low-permeability geomedia, the integrated methodologies for basic properties and pore structure characterization of these geomedia include X-ray diffraction, thin section petrography, grain size distribution, water immersion porosimetry after vacuum-pulling for full saturation, mercury intrusion porosimetry, nitrogen physisorption, scanning electron microscopy, X-ray computed tomography, and (ultra-)small angle neutron (X-ray) scattering. In addition, custom-designed gas diffusion, tracer recipe involving a range of anionic and cationic chemicals with subsequent analyses by laser ablation and inductively coupled plasma-mass spectrometry, along with batch sorption, column transport, and imbibition tests were conducted for coupled effects of pore structure and chemical retention/transport. From the perspectives of pore structure in conjunction with multiple and complementary approaches to examining a range of sample sizes under different observational scales, we find that the poor pore connectivity is prevalent in low-permeability media (mudstone and crystalline rock) that is related to geological processes (e.g., compaction, diagenesis and thermal maturation). For example, the deep and organic matter-rich mudstones have a much smaller effective porosity than the total porosity (as a result of poor pore connectivity) and associated diffusion coefficient, and the effective porosity & diffusion coefficients are also dependent upon the sample sizes used in the measurement. Similarly, most of the pore space in the shallow mudstone is also controlled by pore-throat diameters in the 5-50 nm range of intergranular pore types from its fine-grained nature, but with an overall good pore connectivity. However, the nm-sized pore space (physically pore-network architecture) and strong sorption capacities (chemical retention from clay minerals) of both shallow and deep mudstones lead to the synergistic retention of cationic radionuclides and their utilities as effective host rocks and barrier materials. Our unique approaches of studying how the micro-scale pore structure affect macro-scale fluid flow, diffusion & retention, and chemical transport produce improved mechanistic understanding, and realistic quantification, of diffusion and retention of typical radionuclides in a range of generic host rocks and barrier materials (clay/shale, salt, crystalline rock, and tuff), with the overall results leading to scientifically-based understanding of enhanced isolation (from both diffusion and retention) of radionuclides and improved confidence on the long-term performance of geological repository to store high-level radioactive wastes. In addition to the training of 25 undergraduates, graduates, and postdocs of UTA, the scientists (organizations) involved in performing this work (e.g., discussion, sample sharing, and operation of SANS and SAXS instruments) include Ed Matteo, Yifeng Wang, and Kristopher Kuhlman (Sandia National Laboratories), Jens Birkholzer, Liange Zheng, Tim Kneafsey, and Sharon Borglin (Lawrence Berkeley National Laboratory), Mavrik Zavarin (Lawrence Livermore National Laboratory), Yukio Tachi and Yuta Fukatsu (Japan Atomic Energy Agency), Mieke de Craen (Euridice, Belgium), Markus Bleuel (NIST), Wei-Ren Chen, Gergely Nagy, Changwoo Do, William Heller, Larry Anovitz, and Kenneth Littrell (ORNL), as well as Jan Illvsky, Ivan Kuzmenko, Ju-Sang Park and Jon Almers (ANL). Key deliverables include a total of 13 peer-reviewed journal articles (nine published and three under review), 23 presentations at scientific conferences (AAPG, AAPG Southwest Section, AGU, Asian Clay Conference, GSA, GSA South-Central Section, IHLRWM, InterPore, International Conference on Chemistry and Migration Behavior of Actinides and Fission Products in the Geosphere, International Conference on Coupled Processes in Fractured Geological Media: Observation, Modeling and Application), and academic institutions (UTA, New Mexico State University; University of Poitiers, France; University of Helsinki, Finland; Uppsala University, Sweden; Istanbul Technical University, Turkey) and other organizations (Andra, France; Posiva Oy, Finland).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Unravelling chloride transport/microstructure relationships for blended-cement pastes with the mini-migration method

A chloride mini-migration method is proposed to estimate effective diffusion coefficients at the scale of the cement paste and to investigate mechanisms with complementary microstructure analyses on the same material. Blended-cement pastes with a wide range of properties were investigated in this study: systems at w/b = 0.3–0.5 including Portland cement, white Portland cement, slag-Portland cement, fly ash, glass powder, and/or limestone and calcined clay. Comparisons showed the relative and combined importance of three main parameters on the effective diffusion coefficient: the porosity, the pore connectivity parameter and the conductivity of the pore solution (low values of the latter two are generally key aspects of the high resistance of blended-cement systems against chloride ingress). Notably, a general correlation for all the investigated systems was established between the effective diffusion coefficient and the bulk conductivity (whereas no simple correlation was observed between the effective diffusion coefficient and the formation factor or the pore connectivity parameter).

36 MATERIALS SCIENCE↗
DOE OSTI2025

Cosmic-ray transport in inhomogeneous media

A theory of cosmic-ray transport in multiphase diffusive media is developed, with the specific application to cases in which the cosmic-ray diffusion coefficient has large spatial fluctuations that may be inherently multiscale. We demonstrate that the resulting transport of cosmic rays is diffusive in the long-time limit, with an average diffusion coefficient equal to the harmonic mean of the spatially varying diffusion coefficient. Thus, cosmic-ray transport is dominated by areas of low diffusion even if these areas occupy a relatively small, but not infinitesimal, fraction of the volume. On intermediate time-scales, the cosmic rays experience transient effective subdiffusion, as a result of low-diffusion regions interrupting long flights through high-diffusion regions. In the simplified case of a two-phase medium, we show that the extent and extremity of the subdiffusivity of cosmic-ray transport is controlled by the spectral exponent of the distribution of patch sizes of each of the phases. We finally show that, despite strongly influencing the confinement times, the multiphase medium is only capable of altering the energy dependence of cosmic-ray transport when there is a moderate (but not excessive) level of perpendicular diffusion across magnetic-field lines.

Astronomy and Astrophysics↗
DOE OSTI2024

On-the-fly machine learned force fields for the study of warm dense matter: Application to diffusion and viscosity of CH

We develop a framework for on-the-fly machine learned force field (MLFF) molecular dynamics (MD) simulations of warm dense matter (WDM). In particular, we employ an MLFF scheme based on the kernel method and Bayesian linear regression, with the training data generated from the Kohn–Sham density functional theory (DFT) using the Gauss spectral quadrature method, within which we calculate energies, atomic forces, and stresses. We verify the accuracy of the formalism by comparing the predicted properties of warm dense carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that ab initio MD simulations of WDM can be accelerated by up to three orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the diffusion coefficients and shear viscosity of CH at a density of 1 g/cm3 and temperatures in the range of 75 000–750 000 K. We find that the self- and inter-diffusion coefficients and the viscosity obey a power law with temperature, and that the diffusion coefficient results suggest a weak coupling between C and H in CH. In addition, we find agreement within standard deviation with previous results for C and CH but disagreement for H, demonstrating the need for ab initio calculations as presented here.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

The nanocaterpillar's random walk: diffusion with ligand–receptor contacts

Particles with ligand–receptor contacts bind and unbind fluctuating “legs” to surfaces, whose fluctuations cause the particle to diffuse. Quantifying the diffusion of such “nanoscale caterpillars” is a challenge, since binding events often occur on very short time and length scales. Here, in this work, we derive an analytical formula, validated by simulations, for the long time translational diffusion coefficient of an overdamped nanocaterpillar, under a range of modeling assumptions. We demonstrate that the effective diffusion coefficient, which depends on the microscopic parameters governing the legs, can be orders of magnitude smaller than the background diffusion coefficient. Furthermore it varies rapidly with temperature, and reproduces the striking variations seen in existing data and our own measurements of the diffusion of DNA-coated colloids. Our model gives insight into the mechanism of motion, and allows us to ask: when does a nanocaterpillar prefer to move by sliding, where one leg is always linked to the surface, and when does it prefer to move by hopping, which requires all legs to unbind simultaneously? We compare a range of systems (viruses, molecular motors, white blood cells, protein cargos in the nuclear pore complex, bacteria such as Escherichia coli, and DNA-coated colloids) and present guidelines to control the mode of motion for materials design.

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

Hydrogen diffusion induced dislocation transformations in a nickel superalloy

The diffusion of hydrogen in metals and alloys induces embrittlement that can adversely affect the structural properties. We examine the adsorption and diffusion of hydrogen in Inconel-718 (IN-718), and scrutinize the ensuing effects on the dislocation behavior in the alloy to elucidate the fundamental mechanisms of hydrogen-microstructure interactions from classical molecular simulations. Hydrogen adsorption increases with time until the surface saturates, while hydrogen diffusion exhibits strong temperature dependence, with diffusion coefficients converging above 1300 K regardless of the initial hydrogen concentration in the alloy. The diffusion in IN-718 is significantly sluggish than in pure Ni, Fe, or Cr, and is strongly impacted by hydrogen concentrations, resulting in an order of magnitude higher diffusion coefficient for hydrogen (10-14 m2/s relative to 10-15 m2/s) at high concentrations, especially below 600 K. Hydrogen diffusion coefficient varies from 10-12 to 10-15 m2/s in IN-718 depending on temperature (500–1400 K). More critically, our results reveal that increasing hydrogen concentration induces microstructural changes in the alloy, transforming perfect dislocations into stair-rods and Shockley partials, with higher temperatures favoring the latter. The results are significant for hydrogen fuel applications to gain insights into the materials chemistry for designing safer and more efficient propulsion systems, particularly in high-performance environments related to controlled hydrogen combustion applications.

Banerjee, Tanumoy↗
DOE OSTI2020

Data on the comprehensive first-principles diffusion study of the aluminum-magnesium system

First-principles calculation of diffusion coefficients between Mg and Al is investigated comprehensively using density functional theory (DFT). The effect of different uncertainty sources arising from first principles calculations has been investigated systematically. These sources include the diffusion model, energetic, entropic and attempt frequency calculations. Variation in self and impurity diffusion coefficients of Mg and Al in stable phases are quantified using different DFT settings and compared with the experiments. Using the optimal DFT settings, diffusion coefficients in metastable phases of Al and Mg are predicted. The dataset refers to “An integrated experimental and computational study of diffusion and atomic mobility of the aluminum-magnesium system” [1].

36 MATERIALS SCIENCE↗
DOE OSTI2020

Generalized hydrodynamic analysis of transport through a finite open nanopore for two-component single-file systems

Single-file diffusion (SFD) in finite open nanopores is characterized by nonzero spatially varying tracer diffusion coefficients within a generalized hydrodynamic description. This contrasts with infinite SFD systems where tracer diffusivity vanishes. In standard tracer counterpermeation (TCP) analysis, two reservoirs, each containing a different species, are connected to opposite ends of a finite pore. Here, we implement an extended TCP analysis to allow the two reservoirs to contain slightly different mixtures of the two species. Then, determination of diffusion fluxes through the pore allows extraction of diffusion coefficients for near-constant partial concentrations of the two species. This analysis is applied for a lattice-gas model describing two-component SFD through a finite linear pore represented by a one-dimensional array of cells. Two types of particles, A and B, can hop only to adjacent empty cells with generally different rates, $h_A$ and $h_B$. Particles are noninteracting other than exclusion of multiple cell occupancy. Results reveal generalized hydrodynamic tracer diffusion coefficients which adopt small values inversely proportional to pore length in the pore center, but which are strongly enhanced near pore openings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Investigating Morphology and Diffusion in Simulations of Precise Anion-Conducting Polymers

Using atomistic molecular dynamics simulations, we investigate the morphology and transport properties of a new class of polymers which are functionalized with quaternary ammonium groups for use as anion exchange membranes. The polymers are precision polyolefins with either a trimethylammonium (p5CNMe3) or a dimethyl-hexyl ammonium (p5CNMe2Hx) pendant group at every fifth carbon along a polyethylene backbone. Simulations are performed at hydration levels of 5, 10, 15, and 20 water molecules per ammonium group. The hydrated polymers form nanoscale, percolated hydrophilic domains (water channels) in the hydrophobic polymer matrix that become wider with increasing water content. Water and hydroxide anion diffusion coefficients also increase with increasing water content. The morphology of the water domains is similar in both polymers, while the diffusion coefficients are somewhat lower in p5CNMe2Hx at fixed water content. Furthermore, the diffusion coefficients in both polymers fall on the same curve as a function of the fractal dimension of the percolated water channels, which appears to be a useful scalar measure of the effects of the nanoscale morphology on water and hydroxide anion transport.

Anion exchange membrane↗
DOE OSTI2025

Vacancy-Dependent Diffusion Mechanism in Oxygen-Defective SrFeO 3 Perovskite Materials: First-Principles Density Functional Theory and Experimental Approach

Understanding oxygen diffusion at the atomic scale in SrFeO 3−δ perovskites is crucial for developing oxygen storage materials with optimal performance. Such materials are required to have high stability, corrosion resistance, and acceptable oxygen storage capacity at moderate operating temperatures and pressures. Here, in this study, we used first-principles density functional theory and thermogravimetric analysis to study the vacancy-dependent oxygen diffusion in oxygen-deficient SrFeO 3−δ (δ = 0, 0.065, 0.125, 0.25, 0.5) perovskites. The electronic structures, including the partial- and spin-resolved density of states, for different SrFeO 3−δ phases were calculated and compared with available experimental and theoretical results. By mapping the migration pathways, we investigated diffusion mechanisms and calculated the energy barriers for oxygen diffusion in cubic, orthorhombic, and brownmillerite phases of SrFeO 3−δ perovskites. Using the calculated energy barriers, we deduced the diffusion time scales and diffusion coefficients within SrFeO 3−δ . A diffusion coefficient on the order of 10 –8 m 2 /s was obtained for SrFeO 2.875 . We experimentally investigated the roles of temperature and oxygen partial pressures on the redox kinetics and deduced the kinetics rate and diffusion density, which agreed well with the calculated values for the density of diffusing oxygen vacancy in the lattice. Our results showed that the energy barrier tends to reduce at higher oxygen concentrations. Our results serve as an important guideline for designing oxygen storage materials with optimal redox kinetics.

chemical looping with oxygen uncoupling (CLOU)↗
DOE OSTI2026

Nanodomain Formation and Temperature-Dependent Diffusion in Deep Eutectic Solvents Revealed by Single-Molecule Tracking

Deep eutectic solvents (DESs) are typically regarded as homogeneous liquids; however, recent work shows that many exhibit nanoscale structural heterogeneity. Most studies attribute these nanoscale features to short-range chemical interactions. It is still unclear whether a long-range physical mechanism also plays a role. Here, in this study, we examined the nanoscale structure in two hydrophobic DESs, 1:3 tetrabutylammonium bromide: l-menthol (DES-butyl) and 1:3 tetraoctylammonium bromide: l-menthol (DES-octyl). The notation 1:3 represents the molar ratio of the hydrogen bond acceptors to hydrogen bond donors used in the synthesis of the DESs. Single-molecule tracking (SMT) coupled with maximum entropy method (MEM) analysis was used to measure the number of diffusion populations of a dilute concentration of an added fluorescent probe. The presence of more than one population of diffusion coefficients indicates the existence of multiple local environments for the fluorescent probe (i.e., nanoscale structures in the DES). DES-butyl showed a relatively narrow diffusion coefficient distribution centered at 0.55 μm 2 /s, whereas DES-octyl displayed two distinct diffusing populations at 20 °C, with diffusion coefficients of 0.12 μm 2 /s and 0.53 μm 2 /s for the slow and fast populations, respectively. As DES-octyl was heated, the slow-diffusing population steadily diminished and disappeared above ∼30 °C, indicating that the nanodomains present at lower temperatures collapse as the liquid becomes more thermodynamically mixed. This temperature-dependent homogenization is consistent with a physical mechanism of nanostructure formation, for example, liquid–liquid phase separation (LLPS), wherein the structure is not driven solely by specific chemical interactions. The SMT-MEM results suggest that a long-range physical mechanism is the most plausible origin of the measured nanoscale structure in DES-octyl.

Opare-Addo, Jemima [Ames Laboratory (AMES), Ames, ↗
DOE OSTI2020

Stratification of polymer mixtures in drying droplets: Hydrodynamics and diffusion

In this work, we study the evaporation-induced stratification of a mixture of short and long polymer chains in a drying droplet using molecular simulations. We systematically investigate the effects of hydrodynamic interactions (HI) on this process by comparing hybrid simulations accounting for HI between polymers through the multiparticle collision dynamics technique with free-draining Langevin dynamics simulations neglecting the same. We find that the dried supraparticle morphologies are homogeneous when HI are included but are stratified in core–shell structures (with the short polymers forming the shell) when HI are neglected. The simulation methodology unambiguously attributes this difference to the treatment of the solvent in the two models. We rationalize the presence (or absence) of stratification by measuring phenomenological multicomponent diffusion coefficients for the polymer mixtures. The diffusion coefficients show the importance of not only solvent backflow but also HI between polymers in controlling the dried supraparticle morphology.

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

Theoretical Analysis of Experimental Data of Sodium Diffusion in Oxidized Molybdenum Thin Films

In this work, the diffusion process of sodium (Na) in molybdenum (Mo) thin films while it was deposited on soda lime glass (SLG) was studied. A small amount of oxygen was present in the chamber while the direct-current (DC) magnetron sputtering was used for the deposition. The substrate temperatures were varied to observe its effect. Such molybdenum films, with or without oxidations, are often used in thin film solar cells, either as back contact or as hole transport layers. Secondary ion mass spectrometry (SIMS) was used to quantify the concentration of the species. A grain diffusion mechanistic model incorporating the effect of grain and grain boundary geometrical shape and size was developed. The model was used to provide an in-depth theoretical analysis of the sodium diffusion in molybdenum thin films that lead to the measured SIMS data. It was observed that not only diffusion coefficients should be considered when analyzing diffusion processes in thin films but also the ratio of grain boundary size to grain size. Both depend on substrate temperature and directly affect the amount of diffused species in the film. The data were analyzed under the light of the film growth speed versus diffusion front speed, the effect of oxygen content, and the effect of substrate temperature on the overall diffusion process. The temperature inversely affects the ratio of grain boundary size and grain size and directly affects the diffusion coefficient, which leads to a preferable temperature at which the highest amount of alkali can be found in the film.

36 MATERIALS SCIENCE↗