Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusivity”

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 1,747 records · Page 97

Observation of two different cool flame regimes of diethyl ether in a counterflow burner

This short communication reports, for the first time, the existence of two different self-sustaining cool flame regimes of diethyl ether (DEE) in a diffusion counterflow burner: a weaker autoignition-assisted cool flame near the fuel burner and a normal diffusion cool flame near the stagnation plane. Here, the results show that the normal diffusion cool flame extinction limit increases monotonically with the fuel mole fraction, while the autoignition-assisted cool flame approaches a plateau and can exist at a fuel mole fraction below the normal diffusion cool flame. It is shown that both flame regimes are governed by the same low-temperature chain-branching reaction pathway of DEE. By using in situ laser diagnostics, entrainment of unburned fuel stream to oxygen stream at the outer edge of the fuel burner is identified as the governing physical mechanism causing a partially premixed self-sustained hollow cool flame structure. The results reveal that when a fuel with high low-temperature reactivity, two different cool flame regimes can be observed in a counterflow flame experiment. Future studies with high-reactivity fuels in a counterflow burner must ensure to distinguish between the two self-sustaining cool flame regimes. Moreover, the existence of these different cool flame regimes needs to be examined so that they would not trigger an uncontrolled combustion phasing in advanced engines fed with high low-temperature reactivity fuels.

33 ADVANCED PROPULSION SYSTEMS↗

Effect of CO 2 -brine-rock reactions on pore architecture and permeability in dolostone: Implications for CO 2 storage and EOR

Geologic carbon sequestration (GCS) is considered a feasible technology for storing substantive volumes of greenhouse gases in subsurface geological formations. In the reservoir, far from carbon dioxide (CO 2 ) injection wells or in post-injection scenarios, diffusion dominates over advection. This condition conjoins with spatially distributed geochemical reactions to induce heterogeneous changes in pore architecture, i.e. pore body and throat sizes or surface roughness. These changes can affect CO 2 transport properties and storage capacity. In this work, we investigated mineral dissolution and precipitation in dolomite samples saturated with a CO 2 -saturated brine at 93 °C and 34.5 MPa, aged without flow. Two rock types samples, i.e. intergranular- and vuggy-dominant, were selected to investigate changes in pore size, porosity and permeability under reactive conditions. Mineral dissolution and precipitation were characterized using scanning electron microscopy. Changes in pore size were quantified via time-domain nuclear magnetic resonance (TD-NMR) transverse relaxation time (T 2 ) and diffusion coefficient (D) distributions. We show that mineral dissolution likely occurs in highly permeable pathways. These observations are confirmed through analysis of (T 2 ) and diffusion coefficient (D) distributions. In contrast to results during CO 2 -enriched brine continuous injection, mineral precipitation was observed in micropores. The leftward shift of the T 2 peaks, corresponding to micropores, also evidenced mineral precipitation in lowpermeability zones. However, microscale alterations resulted only in a subtle increase in porosity and permeability. Results in this study shed light on effects of geochemical reactions on alteration of rock properties in diffusion-dominated regions during CO 2 storage.

58 GEOSCIENCES↗

Fabrication and thermophysical properties of UO 2 -UB 2 and UO 2 -UB 4 composites sintered via spark plasma sintering

Uranium dioxide (UO 2 ) composites with uranium diboride (UB 2 ) and uranium tetraboride (UB 4 ) have been proposed as advanced fuel candidates due to their high thermal conductivity, high melting point, high fissile density and their ability to incorporate a built-in burnable poison by tailoring the targeted 10 B/ 11 B ratio. As such, it is important to assess the fabrication, and thermal and micromechanical properties of such composites. In this work, UO 2 -UB 2 and UO 2 -UB 4 samples with boride phase fractions of 5, 15 and 30 wt% were fabricated to high densities (above 95 % theoretical density) via spark plasma sintering (SPS). This enabled sintering at relatively low temperatures and short timescales. SPS also aided in maintaining the target phase fractions of the samples as reactions between the constituent phases were suppressed due to the short timescales and reducing environment during sintering. Here, thermal diffusivity measurements from 299 to 1273 K were conducted through laser flash analysis (LFA). The diffusivity increased as a function of boride weight fraction, and UB 2 additions increased the thermal diffusivity of the composites more than UB 4 additions. Assessment of the LFA results indicated in-situ reactions between the UO 2 and boride phases that suppress the thermal diffusivity occur above 800 K for all samples. Oxidation of the boride phase was proposed as the underlying reaction. This was supported by thermodynamic assessments from the literature, as well as microstructural, crystallographic, and nanoindentation characterization performed on these samples.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Atomistic simulations of He bubbles in Beryllium

Formation of He bubbles can have a significant effect on the microstructural evolution and properties of irradiated materials. In this report, we use atomistic simulations based on machine learning potentials to investigate the fundamental behavior of He bubbles in Be, with a specific focus on the shape, stability, and diffusivity of bubbles. Stability of He bubbles is quantified in terms of formation energies, which are determined as a function of the ratio of He/V. We find that He bubbles become unstable with respect to plastic deformation through punch-out dislocations around the bubble when the He/V ratio is larger than ~1.25, and the punch-out process induces the change of the regular bubble shape. In general, the bubble shape of He in Be is found to be ellipsoid-like. It is also found that for a fixed He/V ratio, the bubble attracts vacancies to become larger in size. If the bubble size is constant, the bubble attracts additional He atoms until the punch-out reaction occurs. The dominant diffusion mechanism of He bubbles changes from surface diffusion to volume diffusion as the temperature is increased, with a crossover occurring at about 900 K.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Mobility assessment of the BCC and carbide phases in the C-Nb, C-U and Nb-U systems

Uranium carbides with refractory metal additions are considered for Gen IV nuclear reactors and nuclear thermal propulsion as fuels for their high-temperature and corrosion resistant properties. Understanding kinetic effects that dictate microstructural evolution during fabrication and operating conditions is essential to advance technological development of these fuels. This work presents the development of an atomic mobility database for C-Nb-U systems based off available experimental data supported with ab-initio methods. The mobility assessments and uncertainty quantification (using Markov chain Monte Carlo) were conducted in the Kawin software. Carbon diffusion is considered dominant, as metal diffusion is much slower, with niobium diffusion being even slower and rate limiting than uranium metal. We provide a comprehensive and self-consistent thermo-kinetic database that is validated by diffusion couple simulations through Kawin. In conclusion, this enables prediction of microstructural and phase evolution critical for the development and lifetime assessment of next generation nuclear fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Coating developments towards enabling aluminum as a bipolar plate material for PEM fuel cells

Bipolar plates represent a significant portion of the cost and weight of a proton exchange membrane fuel cell stack. As a result, there has been significant interest in using low weight and cost materials, such as aluminum. Aluminum has good mechanical and physical properties; however, it performs poorly in the corrosive environment within the fuel cell stack. To overcome the corrosion performance issue, there have been efforts towards developing coating processes to enable aluminum to attain corrosion targets. In an effort to enable aluminum as a viable bipolar plate material, a solid phase processing (diffusion bonding) approach for bonding a titanium foil to aluminum is investigated in this study. The microstructure, corrosion performance, mechanical properties and electrical resistance are investigated. To avoid the known reduction in contact resistance performance of titanium after exposure to bipolar plate media, the effect of addition of highly conductive gold particles after diffusion bonding, on contact resistance and corrosion behavior, is also studied and compared with diffusion bonded Ti to Al. In conclusion, the results indicate that diffusion bonded Ti to Al provides a viable alternative material combination for bipolar plate that avoids the micro pores and crevasses that are associated with vapor deposition or electroplating.

25 ENERGY STORAGE↗

On the thermal coarsening and transformation of nanoscale oxide inclusions in 316L stainless steel manufactured by laser powder bed fusion and its influence on impact toughness

The thermal evolution of nanoscale oxide inclusions in 316L stainless steel (SS) manufactured by laser powder bed fusion additive manufacturing (AM) was explored. The size, chemical composition, morphology, and distribution of the oxides were characterized as the function of heat treatment conditions. The study revealed the mechanistic driving force of the rapid oxide coarsening during recrystallization. Ostwald ripening governs oxide coarsening. The active grain boundary-oxide interaction at the early stage of recrystallization accelerated oxide coarsening via enhanced solute diffusion along grain boundaries. Pipe diffusion along dislocation cellular boundaries has a negligible contribution to oxide coarsening. At high temperatures (T> 1065 °C), although lattice diffusion primarily controlled the oxide growth, the contribution from the grain-boundary diffusion was necessary. The transformation from MnSiO 3 to CrMn 2 O 4 took place in the un-recrystallized grains but was not observed when recrystallization started. The interaction of grain boundary and oxides during recrystallization resulted in a high fraction of oxides accumulated at grain boundaries. While oxide coarsening does not significantly alter the toughness value, grain-boundary oxides promote microvoid formation and intergranular fracture under Charpy impact in the recrystallized AM 316L SS.

36 MATERIALS SCIENCE↗

First in operando transient grating spectroscopy measurements on tungsten during high flux plasma operation in PISCES-RF

This work details the construction and first in operando transient grating spectroscopy measurements conducted on PISCES-RF during plasma operation. A preliminary study on a tungsten sample with varied plasma species (Ar and D 2 ), ion flux (4 × 10 21 to 6 × 10 22 ion/m 2 /s), ion energy (5 to 80 eV), and surface temperature (20 to 200 °C) is presented. Prior to plasma exposure, the thermal diffusivity and surface acoustic wave (SAW) speed, (6.9 ± 0.2) × 10 −5 m 2 /s and (2.66 ± 0.02) km/s, were measured in situ with a thermal grating period of 12.5 μm. During Ar plasma exposure, these two quantities vary according to the sample surface temperature. At low flux and varied ion energy, D 2 plasma exposure also results in thermal diffusivity and SAW speed varying with surface temperature alone. At high flux, D 2 plasma exposure results in no convincing change to the thermal diffusivity, but the SAW speed is reduced. After returning to 20 °C, post plasma exposures, the thermal diffusivity and SAW speed almost recover to the initial value prior to exposure.

Surface acoustic wave SAW↗

Investigation of desorption of hydrogen gas from polymer matrix using thermal desorption analysis and finite element modeling

Polymer materials are used widely as sealing materials and liners in hydrogen storage and transportation infrastructure. During operating conditions, polymers are exposed to high pressure hydrogen gas periodically. When in contact, the hydrogen gas diffuses through the polymer material with ease due to its smaller sized molecules compared to the polymer molecules. In the event of rapid decompression of outside hydrogen gas, the hydrogen molecules within the polymer diffuse out. This process is known as the desorption of polymer and can cause blisters or cracking causing permanent damage to the polymer. Here, we used a combined experimental and modeling approach to investigate the desorption of the hydrogen gas in the polymer. Thermal desorption analysis (TDA) experiments were performed on the polymer sample exposed to the high-pressure hydrogen gas to measure the diffusion coefficient and equilibrium hydrogen concentration. Experiments also provided the quantity of hydrogen within the polymer during the desorption process which followed Fick’s law. However, for higher pressurization values, some anomaly was observed in the quantity of hydrogen remaining in the polymer sample. A continuum mechanics based coupled diffusion–deformation-damage model was deployed to simulate the gas desorption process in the material using the finite element method to understand its effect on the material. It was found that the anomaly in TDA results were due to the blister formation on the polymer surface which resulted in longer hydrogen retention. Simulations also showed that blister formation can be affected by multiple parameters such as the addition of filler particles, and pre-existing cavity size.

36 MATERIALS SCIENCE↗

Ebullition dominates methane fluxes from the water surface across different ecohydrological patches in a temperate freshwater marsh at the end of the growing season

Measurements of the spatial heterogeneity of methane fluxes in wetlands are critical to better understand and predict methane emissions at the ecosystem scale. However, the within-wetland spatial heterogeneity of fluxes is rarely assessed. Here, we use a spatially balanced rapid chamber-based survey of methane at different ecohydrological patches within a temperate freshwater marsh. We measured fluxes exclusively from the water surface without including vegetation. We further used the data from chamber measurements to partition diffusive and ebullitive fluxes. Three ecohydrological patches were distinguishable in the marsh, defined by the type and presence/absence of vegetation. Here, these patches were emergent vegetation, floating-leaved, and open water. Net methane fluxes from the water surface (diffusion plus ebullition) in emergent vegetation patches were larger than in the floating-leaved vegetation and open water patches (p < 0.05). Diffusive fluxes, representing a sizable smaller fraction of net fluxes, were larger in vegetated than in unvegetated patches (p < 0.05), while ebullitive fluxes mirrored the magnitude and differences observed in the net fluxes.Moreover, pooled net and ebullitive fluxes across patches (but not diffusive fluxes) were negatively correlated with water levels, the primary variable affecting patch distribution. Altogether, our results indicate that the differences among ecohydrological patches are driven by ebullition, ultimately highlighting challenges faced by scientists and practitioners in the field and modelers seeking to improve the predictability and resolution of wetland biogeochemical models.

54 ENVIRONMENTAL SCIENCES↗

Ultrafast Oxygen Conduction in Sillén Oxychlorides

Oxygen ion conductors are crucial for enhancing the efficiency of various clean energy technologies, including fuel cells, solid oxide air batteries, electrolyzers, membranes, sensors, and more. In this study, a structure-similarity analysis of ≈62k oxygen-containing compounds identified the MBi 2 O 4 X (M = rare-earth element, X = halogen element) family as promising candidates for fast oxygen transport. Among these, LaBi 2 O 4 Cl is found as an ultrafast oxygen conductor with an ultralow migration barrier of 0.1 eV based on ab initio studies. Its 2D layered structure, featuring a “triple fluorite” layer, supports diffusion of both oxygen vacancies and interstitials. In addition to vacancy diffusion with a 0.1 eV barrier, ab initio studies show interstitial diffusion exhibits a modest barrier of 0.6–0.8 eV. Frenkel pairs are found to be the dominant defects in intrinsic LaBi 2 O 4 Cl, facilitating significant vacancy-mediated oxygen diffusion at elevated temperatures. With 2.8% oxygen vacancies, LaBi 2 O 4 Cl is predicted to achieve a conductivity of 0.3 S/cm at 25 °C in a single crystal. Experimental synthesis and characterization of polycrystalline LaBi 2 O 4 Cl and Sr-doped LaBi 2 O 4 Cl revealed conductivity exceeding that of YSZ and LSGM below 400 °C, with lower activation energies, achieving a total conductivity of 0.1−0.2 mS/cm at 300 °C. Here, while these results confirm its potential of fast oxygen transport, we suggest further experimental optimization of LaBi 2 O 4 Cl, including aliovalent doping and microstructure refinement, could significantly enhance its performance, facilitating fast oxygen conduction approaching room temperature.

Defects↗

Solid 3D Li–S Battery Design via Stacking 2D Conductive Microporous Coordination Polymers and Amorphous Li–S Layers

To make a lithium-sulfur (Li-S) battery practical, not only high gravimetric energy capacity is important, but also high volumetric energy capacity will be required. The currently explored Li-S cathode designs often deploy systems with liquid electrolyte infiltration, hence with relatively low volumetric capacity. In the current study, we theoretically test a compact solid three-dimensional (3D) design (more like a Li-ion battery cathode than a conventional Li-S cathode) consisted of a sandwich structure alternating between the two-dimensional (2D) Mn-hexaaminobenzene-based coordination polymer (2D Mn-HAB-CP) layer and the amorphous Li-S layer. We study the theoretical limits for both its gravimetric and volumetric energy capacity, as well as its structural stability and Li diffusion within the cathode system. To study the Li diffusion within an amorphous system, we also develop a pull-atom molecular dynamics (PA-MD) to calculate the barrier heights of such disordered systems. We reveal the mechanism that determines the Li diffusion in the amorphous layer of the system. Overall, we find such a 3D solid Li-S cathode can be practical, with sufficient large gravimetric and volumetric energy capacity, as well as the Li diffusion constant. It also solves many other common Li-S cathode problems, from Li polysulfide dissolution to electrical insulating, and structure instabilities.

25 ENERGY STORAGE↗

Molecular Dynamics Simulation Study of the Protonation State Dependence of Glutamic Acid Transport through a Cyclic Peptide Nanotube

The effect of the protonation state of glutamic acid on its translocation through cyclic peptide nanotubes (CPNs) was assessed by using molecular dynamics (MD) simulations. Anionic (GLU–), neutral zwitterionic (GLU0), and cationic (GLU+) forms of glutamic acid were selected as three different protonation states for an analysis of energetics and diffusivity for acid transport across a cyclic decapeptide nanotube. Based on the solubility-diffusion model, permeability coefficients for the three protonation states of the acid were calculated and compared with experimental results for CPN-mediated glutamate transport through CPNs. Potential of mean force (PMF) calculations reveal that, due to the cation-selective nature of the lumen of CPNs, GLU–, so-called glutamate, shows significantly high free energy barriers, while GLU+ displays deep energy wells and GLU0 has mild free energy barriers and wells inside the CPN. The considerable energy barriers for GLU– inside CPNs are mainly attributed to unfavorable interactions with DMPC bilayers and CPNs and are reduced by favorable interactions with channel water molecules through attractive electrostatic interactions and hydrogen bonding. Unlike the distinct PMF curves, position-dependent diffusion coefficient profiles exhibit comparable frictional behaviors regardless of the charge status of three protonation states due to similar confined environments imposed by the lumen of the CPN. The calculated permeability coefficients for the three protonation states clearly demonstrate that glutamic acid has a strong protonation state dependence for its transport through CPNs, as determined by the energetics rather than the diffusivity of the protonation state. In addition, the permeability coefficients also imply that GLU– is unlikely to pass through a CPN due to the high energy barriers inside the CPN, which is in disagreement with experimental measurements, where a considerable amount of glutamate permeating through the CPN was detected. To resolve the discrepancy between this work and the experimental observations, several possibilities are proposed, including a large concentration gradient of glutamate between the inside and outside of lipid vesicles and bilayers in the experiments, the glutamate activity difference between our MD simulations and experiments, an overestimation of energy barriers due to the artifacts imposed in MD simulations, and/or finally a transformation of the protonation state from GLU– to GLU0 to reduce the energy barriers. Altogether, our study demonstrates that the protonation state of glutamic acid has a strong effect on the transport of the acid and suggests a possible protonation state change for glutamate permeating through CPNs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Insights into Rotational and Translational Dynamics in Mixtures of Ethylene Glycol and Choline Chloride Using Nuclear Magnetic Resonance Techniques

This work examines molecular dynamics and interactions in ethylene glycol–choline chloride (EG–ChCl) mixtures across 0–33 mol % ChCl, spanning the true eutectic region near 17–20 mol % and the commonly used 1:2 formulation. We combine pulsed-field-gradient (PFG) diffusion, fast-field-cycling (FFC) relaxometry, temperature-dependent 13 C T 1 , and nuclear Overhauser effect spectroscopy (NOESY) to disentangle local from macroscopic dynamics. PFG and FFC show that both translational and average rotational motions largely track the strong increase in viscosity with ChCl content, with ethylene glycol consistently diffusing faster than the choline cation and no global dynamical anomaly at the eutectic composition. More subtle, site-specific composition effects nevertheless emerge. The ratio of the diffusion coefficient of the hydroxyl group of choline to the diffusion coefficient of the methyl group of choline displays a shallow minimum in the 17–25 mol % region, indicating a modest change in how the hydroxyl-bearing end of choline samples the underlying translational motion relative to the methyl groups. 13 C T 1 analysis shows that rotational correlation times at 25 °C generally increase with ChCl, reflecting viscosity-coupled slowing, while the CH 2 –N α site exhibits a small but reproducible deviation from this monotonic trend near the eutectic. NOESY spectra at similar compositions reveal enhanced cross-relaxation between EG and choline protons, consistent with increased headgroup–solvent contact density rather than a wholesale structural rearrangement. Overall, our multitechnique study demonstrates that EG–ChCl dynamics are predominantly viscosity-dominated, with the eutectic region acting as a subtle dynamical crossover where specific choline segments become maximally coupled to the hydrogen-bond network. These insights refine the structure–dynamics picture of choline-chloride DESs and provide practical guidance for tuning composition in electrochemical, separation, and catalytic applications.

diffusion↗

UiO-66 Metal–Organic Framework as an Anode for a Potassium-Ion Battery: Quantum Mechanical Analysis

The natural abundance of potassium in the earth’s crust is 1000 times higher than that of lithium, so energy technologies built on potassium are more sustainable. Potassium-ion batteries have attracted considerable attention because of their relatively low cost and high operating potential, but questions remain about the best anode material for such batteries. Here, in this paper, we report first-principles computations based on density functional theory to investigate the performance of the UiO-66 metal–organic framework as an anode material for potassium-ion batteries; the goal is to provide a fundamental understanding of metal–organic framework (MOF)-based electrodes to guide the design and development of high-performance potassium-ion batteries. Our study includes the stability and electronic properties of potassiated structures and the mechanisms of potassium intercalation and diffusion in the framework. The results indicate that UiO-66 has a maximum specific capacity of 644 mAh/g as the anode of a potassium-ion battery. During potassiation, we observe charge transfer from potassium to carbon or oxygen of UiO-66 near the intercalated K. During K diffusion, the K migrates along the UiO-66 framework with a maximal migration energy barrier of 0.377 eV in the optimal pathway, which is much lower than the barriers for Li and Na diffusion in UiO-66. The diffusion coefficient of K in the anode is several orders of magnitude larger than those of Li and Na. This favors potassium ions over lithium ions or sodium ions when UiO-66 is the anode.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomic-Scale Insights into Carbon Dissolution in α-, γ-, and θ-Al 2 O 3 : Phase-dependent Transport Dynamics from First-Principles Calculations

α-Al 2 O 3 exhibits superior carburizing corrosion resistance compared to metastable γ-Al 2 O 3 and θ-Al 2 O 3 phases in high-temperature CO 2 environments, yet its atomic-scale origins remain unclear. Using first-principles density functional theory, we systematically investigate carbon dissolution and diffusion in α-Al 2 O 3 , γ-Al 2 O 3 , and θ-Al 2 O 3 , including the effects of oxygen (O) and aluminum (Al) vacancies. Our results show that α-Al 2 O 3 consistently exhibits higher carbon solution enthalpies than γ-Al 2 O 3 and θ-Al 2 O 3 in both pristine and defective structures, indicating lower intrinsic carbon solubility in α-Al 2 O 3 . Vacancies significantly enhance carbon incorporation: O vacancies reduce solution enthalpy, while Al vacancies further amplify this effect, with a strong preference for carbon at Al vacancy sites. Carbon diffusion barriers are also highest in α-Al 2 O 3 , reflecting slower carbon mobility. Al vacancies increase diffusion barriers across all phases, while O vacancies raise barriers in α- and γ-Al 2 O 3 but slightly lower them in θ-Al 2 O 3 . These results reveal a dual mechanism behind the carburizing resistance of α-Al 2 O 3 : reduced carbon solubility and elevated diffusion barriers. Furthermore, this work provides atomic-scale insights to guide the design of alumina-based materials with improved carburizing resistance through phase selection and defect engineering.

36 MATERIALS SCIENCE↗

Strong Enhancement of Nanoconfined Water Mobility by a Structure Breaking Salt

For the majority of the water present on earth, the two most important factors influencing its behavior are confinement, in either inorganic or organic matrixes, and the presence of solutes. Here, we investigate the effect of confinement in 3 nm pores on water diffusivity in aqueous solutions with archetypical solutes, a structure making (kosmotrope) NaCl and a structure breaking (chaotrope) KCl, up to 1.0 M in concentration. The water diffusivity in bulk aqueous solutions in such a concentration range is known to decrease very slightly in the presence of NaCl and increase very slightly in the presence of KCl. However, here we observe the water diffusivity in confined H 2 O–KCl increases by a factor of 2 compared to the pure water diffusivity in the same confinement. This unusually strong cumulative effect of confinement and a structure breaking additive may have profound implications for the mobility and transport of aqueous species in nature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Adatom-Driven Oxygen Intermixing during the Deposition of Oxide Thin Films by Molecular Beam Epitaxy

Thin film deposition from the vapor phase is a complex process involving adatom adsorption, movement, and incorporation into the growing film. While simulations have been able to capture key aspects of these processes to reproduce film morphology, these models typically ignore any subsurface diffusion that occurs. The results of these diffusion and intermixing processes are often observed experimentally, but cannot be explained by invoking bulk lattice diffusion at the deposition temperature. In this work, we present quantitative experimental data that reveals anion and cation intermixing over long length scales during the deposition of epitaxial Fe 2 O 3 and Cr 2 O 3 films and heterostructures by oxygen-plasma-assisted molecular beam epitaxy. We track this diffusion by incorporating well-defined tracer layers containing 18 O and/or 57 Fe, and measure their redistribution on the nanometer scale with atom probe tomography. Molecular dynamics simulations suggest potential intermixing events, which are then examined via nudged elastic band calculations. These reveal that adatoms on the film surface act to “pull up” subsurface O and Fe. Subsequent ring-like rotation mechanisms involving both adatom and subsurface anions then facilitate their mixing. A simple model incorporating these events qualitatively agrees with the experimental intermixing results. In addition to film deposition, these intermixing mechanisms may to be operant during other surface-mediated processes such as heterogeneous catalysis and corrosion.

77 NANOSCIENCE AND NANOTECHNOLOGY↗