Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion process”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Diffusion of Oxygen Isotopes in Thermally Evolving Planetesimals and Size Ranges of Presolar Silicate Grains

Presolar grains are small particles found in meteorites through their isotopic compositions, which are considerably different from those of materials in the solar system. If some isotopes in presolar grains diffused out beyond their grain sizes when they were embedded in parent bodies of meteorites, their isotopic compositions could be washed out, and hence the grains could no longer be identified as presolar grains. We explore this possibility for the first time by self-consistently simulating the thermal evolution of planetesimals and the diffusion length of {sup 18}O in presolar silicate grains. Our results show that presolar silicate grains smaller than ∼0.03 μ m cannot keep their original isotopic compositions even if the host planetesimals experienced a maximum temperature as low as 600 °C. Since this temperature corresponds to that experienced by petrologic type 3 chondrites, isotopic diffusion can constrain the size of presolar silicate grains discovered in such chondrites to be larger than ∼0.03 μ m. We also find that the diffusion length of {sup 18}O reaches ∼0.3–2 μ m in planetesimals that were heated up to 700–800°C. This indicates that, if the original size of presolar grains spans a range from ∼0.001 μ m to ∼0.3 μ m like that in the interstellar medium, then the isotopic records of the presolar grains may be almost completely lost in such highly thermalized parent bodies. We propose that isotopic diffusion could be a key process to control the size distribution and abundance of presolar grains in some types of chondrites.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Modeling of deposit formation in mesoporous substrates via atomic layer deposition: Insights from pore‐scale simulation

Atomic layer deposition (ALD) has been a promising technique in fabricating membranes and tuning their properties with a precision at the atomic level. Fabrication of zeolitic imidazolate framework (ZIF) membranes using the ligand-induced permselectivation (LIPS) method starts with the formation of an oxide in a mesoporous substrate by ALD and is followed by the transformation of this oxide to ZIF using imidazolate vapor treatment. The objective of the ALD step is to block the mesopores with a thin deposit, that is, one with small penetration depth and small thickness on the top surface of the substrate. Unlike typical ALD on nonporous substrates, where all available sites react per ALD cycle, thin deposit formation in a mesoporous substrate requires that only a small fraction of the available deposition sites (i.e., close to the substrate surface) is subjected to ALD. Consequently, reactant dosing and duration of pulses are important process variables which, together with diffusion and reaction kinetics determine the deposit structure. Quantitative understanding of the interplay of these variables and phenomena can enable the rational design of ALD within mesoporous substrates. Here, we extend our earlier modeling effort considering the coexistence of ALD both inside the pores and on the external surface of the substrate. Finite-volume based models were developed and validated to simulate the two distinct modes of deposition cycle by cycle. The total mass uptake of the substrate with ALD cycles can be predicted using the combined surface deposition and pore reaction–diffusion models as affirmed by in situ quartz crystal microbalance experimental data. The ALD reactor model combined with the deposition model can accurately capture the number of ALD cycles needed to block the pores of the substrate. Based on the model, we designed a modified ALD process and examined the performance of the corresponding LIPS membranes. Furthermore, the present modeling work provides a new understanding of the deposit formation via ALD within mesoporous substrates for a variety of membrane applications.

atomic layer deposition↗
DOE OSTI2016

Bond Strength of Plasma Sprayed Zr Diffusion Barrier on U-Mo Fuels

The United States government is continually committed to strengthening nuclear security and nonproliferation for reducing the threat of radioactive materials acquisition by terrorists. To meet this important mission, the US Department of Energy and National Nuclear Security Agency's (NNSA) Office of Material Management and Minimization (M{sup 3}) Reactor Conversion Program (CONVERT) seeks to minimize the use of highly enriched uranium (HEU) in civilian applications through work with governments and facilities around the world to convert research reactors to the use of non-weapons-usable low enriched uranium (LEU) fuels. Each reactor converted not only eliminates the need for HEU at civilian sites, but also reduces the amount of HEU being manufactured, stored, and transported, where the nuclear material is most vulnerable. In instances where suitable LEU fuels do not exist for the conversion of particular reactors, CONVERT contributes to development of new LEU fuels. Currently, the CONVERT Program is committed to development of a high-density U-Mo fuel that will be able to convert high performance research reactors, while establishing an efficient and economic fabrication capability to manufacture the fuel. To that end, Los Alamos National Laboratory (LANL) has been tasked with development of a method to affix zirconium (Zr) metal fuel diffusion barrier via plasma spraying, which will prevent a reaction between the U-Mo fuel and the 6061 aluminum cladding. A successful diffusion barrier requires adequate Zr layer thickness, uniform coverage, sufficient bonding at the Zr/fuel and Zr/aluminum cladding interface, and maintenance of the original nuclear fuel properties i.e. physical dimensions, alloy composition and dispersion, crystallographic phase concentration, neutronics, etc. This work outlines development of the Zr plasma spray processing techniques, parameters affecting the diffusion barrier/nuclear fuel bond strength, and characterization of the Zr diffusion barrier. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Effects of polycrystallinity in nano patterning by ion-beam sputtering

Employing graphites with distinctly different mean grain sizes, we study the effects of polycrystallinity on the pattern formation by ion-beam sputtering. The grains influence the growth of the ripples in a highly anisotropic fashion; both the mean uninterrupted ripple length along the ridges and the surface width depend on the mean size of the grains, which is attributed to the large sputter yield at the grain boundary compared with that on the terrace. In contrast, the ripple wavelength does not depend on the mean size of the grains, indicating that the mass transport across the grain boundaries should efficiently proceed by both thermal diffusion and ion-induced processes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Spatially Mapping the CO 2 Alkaline Sorbent Diffuse Microenvironment Using Operando Raman Spectroscopy

When designing a chemical process, the local balance of transport and kinetics, collectively referred to as the diffuse microenvironment, plays a critical role in performance but is difficult to directly observe. This work demonstrates a method of two-dimensional spatial chemical mapping of the diffuse microenvironment in the context of alkali metal hydroxide direct air capture of carbon dioxide using a custom operando gas-absorption flow cell along with confocal Raman spectroscopy. Notably, we observe the concentration boundary layer near the gas–liquid interface and elucidate the interplay of carbonate and bicarbonate ions within it while inferring local hydroxide depletion through continuum modeling. These first of their kind observations provide a technique to compare the performance of direct air capture solvents based on diffuse microenvironment dynamics while also providing metrics important for air contactor design such as boundary layer thickness. Overall, this work showcases a new experimental platform to study interfacial diffuse microenvironments in and outside of the field of direct air capture of carbon dioxide.

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

Review: Induced Seismicity During Geoenergy Development—A Hydromechanical Perspective

The basic triggering mechanism underlying induced seismicity traces back to the mid-1960s that relied on the process of pore-fluid pressure diffusion. The last decade has experienced a renaissance of induced seismicity research and data proliferation. An unprecedent opportunity is presented to us to synthesize the robust growth in knowledge. The objective of this article is to provide a concise review of the triggering mechanisms of induced earthquakes with a focus on hydro-mechanical processes. Four mechanisms are reviewed: pore-fluid pressure diffusion, poroelastic stress, Coulomb static stress transfer, and aseismic slip. For each, an introduction of the concept is presented, followed by case studies. Diving into these mechanisms sheds light on several outstanding questions. For example, why did some earthquakes occur far from fluid injection or after injection stopped? Our review converges on the following conclusions: (a) Pore-fluid pressure diffusion remains a basic mechanism for initiating inducing seismicity in the near-field. (b) Poroelastic stresses and aseismic slip play an important role in inducing seismicity in regions beyond the influence of pore-fluid pressure diffusion. (c) Coulomb static stress transfer from earlier seismicity is shown to be a viable mechanism for increasing stresses on mainshock faults. (d) Multiple mechanisms have operated concurrently or consecutively at most induced seismicity sites. (e) Carbon dioxide injection is succeeding without inducing earthquakes and much can be learned from its success. Future research opportunities exist in deepening the understanding of physical and chemical processes in the nexus of geoenergy development and fluid motion in the Earth’s crust.

58 GEOSCIENCES↗
DOE OSTI2022

Advection versus diffusion in Richtmyer-Meshkov mixing

The Richtmyer-Meshkov (RM) instability is one of the most severe degradation mechanisms for inertial confinement fusion (ICF), and mitigating it has been a priority for the global ICF effort. In this Letter, the instability's ability to atomically mix is linked to its background decay of residual turbulent energy. Here, we show how recently derived inequalities from the mathematical theory of PDEs constrain the evolution. A model RM process at leading order may diffusively mix or retain imprints of its initial structures indefinitely, depending on initial conditions, and there exists a theoretical range of zero-mixing for certain values of parameters. The results may apply to other systems resembling scalar transport in decaying turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Do TFSA anions slither? Pressure exposes the role of TFSA conformational exchange in self-diffusion

Multi-nuclear ( 1 H, 2 H, and 19 F) magnetic resonance spectroscopy techniques as functions of temperature and pressure were applied to the study of selectively deuterated 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)amide (EMIM TFSA) ionic liquid isotopologues and related ionic liquids. For EMIM TFSA, temperature-dependent 2 H T 1 data indicate stronger electric field gradients in the alkyl chain region compared to the imidazolium ring. Most significantly, the pressure dependences of the EMIM and TFSA self-diffusion coefficients revealed that the displacements of the cations and anions are independent, with diffusion of the TFSA anions being slowed much more by increasing pressure than for the EMIM cations, as shown by their respective activation volumes (28.8 ± 2.5 cm³/mol for TFSA vs. 14.6 ± 1.3 cm³/mol for EMIM). Increasing pressure may lower the mobility of the TFSA anion by hindering its interconversion between trans and cis conformers, a process that is coupled to diffusion according to published molecular dynamics simulations. Measured activation volumes (ΔV ‡ ) for ion self-diffusion in EMIM bis(fluoromethylsulfonyl)amide and EMIM tetrafluoroborate support this hypothesis.

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

Mass Transport Limitations and Kinetic Consequences of Corn Stover Deacetylation

Alkaline pretreatment of herbaceous feedstocks such as corn stover prior to mechanical refining and enzymatic saccharification improves downstream sugar yields by removing acetyl moieties from hemicellulose. However, the relationship between transport phenomena and deacetylation kinetics is virtually unknown for such feedstocks and this pretreatment process. Here, we report the development of an experimentally validated reaction–diffusion model for the deacetylation of corn stover. A tissue-specific transport model is used to estimate transport-independent kinetic rate constants for the reactive extraction of acetate, hemicellulose and lignin from corn stover under representative alkaline conditions (5–7 g L -1 NaOH, 10 wt% solids loadings) and at low to mild temperatures (4–70°C) selected to attenuate individual component extraction rates under differential kinetic regimes. The underlying transport model is based on microstructural characteristics of corn stover derived from statistically meaningful geometric particle and pore measurements. These physical descriptors are incorporated into distinct particle models of the three major anatomical fractions (cobs, husks and stalks) alongside an unsorted, aggregate corn stover particle, capturing average Feret lengths of 917–1239 μm and length-to-width aspect ratios of 1.8–2.9 for this highly heterogeneous feedstock. Individual reaction–diffusion models and their resulting particle model ensembles are used to validate and predict anatomically-specific and bulk feedstock performance under kinetic-controlled vs. diffusion-controlled regimes. In general, deacetylation kinetics and mass transfer processes are predicted to compete on similar time and length scales, emphasizing the significance of intraparticle transport phenomena. Critically, we predict that typical corn stover particles as small as ~2.3 mm in length are entirely diffusion-limited for acetate extraction, with experimental effectiveness factors calculated to be 0.50 for such processes. Debilitatingly low effectiveness factors of 0.021–0.054 are uncovered for cobs—implying that intraparticle mass transfer resistances may impair observable kinetic measurements of this anatomical fraction by up to 98%. These first-reported quantitative maps of reaction vs. diffusion control link fundamental insights into corn stover anatomy, biopolymer composition, practical size reduction thresholds and their kinetic consequences. These results offer a guidepost for industrial deacetylation reactor design, scale-up and feedstock selection, further establishing deacetylation as a viable biorefinery pretreatment for the conversion of lignocellulosics into value-added fuels and chemicals.

09 BIOMASS FUELS↗
DOE OSTI2019

Beyond the effectiveness factor: Multi-step reactions with intraparticle diffusion limitations

Though common in many industrial chemical processes, multi-step reactions occurring within porous catalyst particles provide significant challenges to previously established modeling frameworks. In this work, a system of reaction-diffusion equations is solved analytically via eigenanalysis. The intraparticle concentration profiles are integrated over the particle volume to obtain effectiveness factors for each reacting species, referred to herein as a multi-step effectiveness vector (MEV). The MEV solution is tested for the case of uniform flow past a stationary catalyst particle with a state-of-the-art, multi-step reaction scheme for fluid catalytic cracking (FCC). Outputs from particle-resolved direct numerical simulation (DNS) are compared to outputs from ordinary differentials with effective reaction closures coming from the MEV solution, the classic single-step effectiveness factor (SEF) (R. Aris, 'The mathematical theory of diffusion and reaction in permeable catalysts: The theory of the steady state,' Oxford University Press, 1975), and utilizing the intrinsic rates as-is. When compared to DNS, the MEV solution shows excellent agreement for the formation of all species and quantitatively captures the intraparticle concentration profiles. By contrast, the intrinsic treatment shows an over-prediction for product formation and the SEF treatment shows mixed levels of accuracy that depend upon the reaction coupling. Specifically, the SEF treatment accurately characterizes product formation if the kinetic scheme is weakly coupled but shows significant departure if stronger coupling is present. The computational cost associated with the MEV solution is associated with eigendecomposition of the Thiele matrix, an N x N matrix where N is the number of reacting components. Overall, the method is shown to be relatively easy to implement and may significantly enhance the fidelity of chemical reactor simulations that employ multi-step reaction schemes.

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

Scaling Limits and Generic Bounds for Exploration Processes

We consider exploration algorithms of the random sequential adsorption type both for homogeneous random graphs and random geometric graphs based on spatial Poisson processes. At each step, a vertex of the graph becomes active and its neighboring nodes become blocked. Given an initial number of vertices N growing to infinity, we study statistical properties of the proportion of explored (active or blocked) nodes in time using scaling limits. We obtain exact limits for homogeneous graphs and prove an explicit central limit theorem for the final proportion of active nodes, known as the jamming constant, through a diffusion approximation for the exploration process which can be described as a unidimensional process. We then focus on bounding the trajectories of such exploration processes on random geometric graphs, i.e., random sequential adsorption. As opposed to exploration processes on homogeneous random graphs, these do not allow for such a dimensional reduction. Instead we derive a fundamental relationship between the number of explored nodes and the discovered volume in the spatial process, and we obtain generic bounds for the fluid limit and jamming constant: bounds that are independent of the dimension of space and the detailed shape of the volume associated to the discovered node. Lastly, using coupling techinques, we give trajectorial interpretations of the generic bounds.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Valuation of exotic options in the framework of Levy processes

In this paper we explore a straightforward procedure to price derivatives by using the Monte Carlo approach when the underlying process is a jump-diffusion. We have compared the Black-Scholes model with one of its extensions that is the Merton model. The latter model is better in capturing the market’s phenomena and is comparative to stochastic volatility models in terms of pricing accuracy. We have presented simulations of asset paths and pricing of barrier options for both Geometric Brownian motion and exponential Levy processes as it is the concrete case of the Merton model. A desired level of accuracy is obtained with simple computer operations in MATLAB for efficient computational time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Hybrid finite element and Brownian dynamics method for charged particles

Diffusion is often the rate-determining step in many biological processes. Currently, the two main computational methods for studying diffusion are stochastic methods, such as Brownian dynamics, and continuum methods, such as the finite element method. A previous study introduced a new hybrid diffusion method that couples the strengths of each of these two methods, but was limited by the lack of interactions among the particles; the force on each particle had to be from an external field. This study further develops the method to allow charged particles. The method is derived for a general multidimensional system and is presented using a basic test case for a one-dimensional linear system with one charged species and a radially symmetric system with three charged species.

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

Polymer-Stabilized Liquid Metal Nanoparticles as a Scalable Current Collector Engineering Approach Enabling Lithium Metal Anodes

Dendrites and dead lithium formation over prolonged cycling have long been challenges that hinder the safe implementation of metallic Li anodes. In this study, we employ polymer-stabilized liquid metal nanoparticles (LM-P NPs) of eutectic gallium indium (EGaIn) to create uniform Li nucleation sites enabling homogeneous lithium electrodeposition. Block copolymers of poly(ethylene oxide) and poly(acrylic acid) (PEO-b-PAA) were grafted onto the EGaIn surface, forming stabilized, well-dispersed NPs. Using a scalable spray coating approach, LM-P NPs were fabricated on copper current collectors, providing lithiophilic PEO sites and interactive carboxyl groups to guide Li deposition. The Li-EGaIn alloying process greatly reduced the Li + diffusion barrier, enabling fast Li transport through the coating layer, resulting in decreased nucleation overpotential. Therefore, about five times lower Li nucleation overpotential was obtained on the LM-P modified Cu with an optimal composition of the polymers than the bare Cu substrates. DFT computations was used to reveal the binding properties between the LM-P layer and Li. Due to the regulated Li plating/stripping process, as-obtained 30 μm Li anodes paired with LiNi 0.8 Co 0.1 Mn 0.1 O 2 (NCM8 11 ) with a negative/positive electrode capacity (N/P) ratio ~ 10 exhibited stable cycling performance at 0.5C for over 250 cycles, with an average Coulombic efficiency of 99.55%. Ultrathin Li (1 μm) anodes with an N/P ratio ~ 0.6 were also demonstrated in Li|LiFePO 4 cells, which examined the stabilization of Li by LM-P NPs and monitored practical loadings of Li anodes that are close to anode-free systems.

25 ENERGY STORAGE↗
DOE OSTI2018

Sinter-Bonding of AISI 316L and 17-4 PH Stainless Steels

The sinter-bonding method was used to produce AISI 316L/17-4 PH components by conventional powder metallurgy technology. The production process consisted of following steps: the powder layering, co-pressing and co-sintering processes. The co-sintering process was carried out at different sintering temperatures (1200, 1250 and 1300 °C), and times (60, 90 and 120 min) in hydrogen atmosphere. The effect of co-pressing and co-sintering conditions on the ability to create a permanent connection between AISI 316L and 17-4 PH stainless steels was investigated. The densification, dimensional change and microstructural evolution of sintered AISI 316L/17-4 PH components were studied. The increase in co-pressing pressure as well as using of the longer time and the higher temperature of co-sintering process contributed to increase in AISI 316L/17-4 PH component density. The dilatometric study as well as dimensional changes measurements revealed an appropriate sintering compatibility for AISI 316L and 17-4PH steels when the sintering process was carried out at the temperature of 1200 °C for time ranging from 60 to 120 min and at 1250 °C for shorter sintering times. The sintering at temperature of 1300 °C (regardless of other manufacturing conditions) led to a significant deviation from the cylindrical shape of samples, as a result of difference of shrinkage of AISI 316L and 17-4 PH steels. Microstructural study indicated that AISI 316L and 17-4 PH steels creates a good diffusion joint (without any cracks and warpage) during the high-temperature sintering process. Nickel diffusion takes place from the AISI 316L to 17-4 PH stainless steel, and copper and molybdenum diffuse in opposite direction of nickel diffusion. It is possible to achieve a solid connection between two materials exhibiting difference in chemical composition by using sinter-bonding and conventional powder metallurgy technology.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Challenges and Recent Developments of Neutron Non-Destructive Assay with UO{sub 2}F{sub 2}-Bearing Materials - 17336

In the enrichment process of an operating Gaseous Diffusion Plant (GDP), such as the Portsmouth and Paducah facilities, UF6 gas will immediately hydrolyze forming the chemical compound UO{sub 2}F{sub 2} (Uranyl Fluoride). This mostly happens in Alumina (Al{sub 2}O{sub 3}) traps where the UF6 adsorbs to moistened Alumina pellets forming a UO{sub 2}F{sub 2} layer around the pellet. Similarly, NaF.UF{sub 6} absorption traps produce the complex Na{sub 2}UF{sub 8} which also hydrolyzes when exposed to moist air. UO{sub 2}F{sub 2} can also be formed in many other processes where UF6 gas comes into contact with moisture within structures and objects such as pipes, converters compressors and coolers. When the structures are refurbished, or designated for D and D, UO{sub 2}F{sub 2}-bearing objects are generated both in situ and ex situ requiring fissile (U-235) characterization. In addition, ex situ waste streams are generated when UO{sub 2}F{sub 2}-bearing materials are placed in drummed or boxed waste containers. Gamma spectroscopy Non-Destructive Assay (NDA) is commonly utilized to characterize UO{sub 2}F{sub 2}-bearing waste but generally has high associated Total Measurement Uncertainties (TMU's) and, in some cases, Data Quality Objectives (DQO's) cannot be met. In addition, high UO{sub 2}F{sub 2} mass loadings, such as hold-up materials, require gamma modeling when infinitely thick to the Uranium energy lines. A common neutron NDA technique, or method, is to exploit the neutrons generated from the F-19(a, n)Na-22 reaction which has an appreciable reaction cross-section. The isotopes of Uranium are a-emitters, and, in particular, U-234 has a significantly higher a-decay reaction cross-section leading to a significant relative neutron Specific Activity (nSA). Neutron NDA is appealing since, in some cases, it has potentially lower TMU's than gamma NDA and neutron detector slabs facilitate surrounding large, complex-shaped, objects. Characterizing UO{sub 2}F{sub 2}-bearing waste for fissile content using neutrons generated from U-234, requires the U-235/U-234 ratio as well as the enrichment as well as the corresponding nSA's. It is also noted that U-238 has a considerable, relative, spontaneous fission nSA and this is a higher concern for low-enrichments where the mass of U-238 dominates. Additional neutron interfering elements include cosmic-ray induced spallation within the objects and background from other, in situ, UO{sub 2}F{sub 2}-bearing objects. In this paper we present the challenges and recently developed hardware and algorithmic strategies for quantifying fissile mass both in situ and ex situ using neutron NDA that works in conjunction to gamma spectroscopy NDA, or process knowledge, for objects and drummed/boxed waste containing UO{sub 2}F{sub 2}-bearing materials leading to minimal TMU's and achievable DQO's. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Investigation of Cu(In,Ga)Se{sub 2} polycrystalline growth: Ga diffusion and surface morphology evolution

Graphical abstract: - Highlights: • Ga diffusion in CIGS absorption layer after annealing treatment. • Phenomenon of surface reconstruction after annealing treatment. • Understand selenium effect on CIGS annealing process. • Explain the kinetic of Ga diffusion and MoSe{sub 2} formation. - Abstract: We report a study of selenization and annealing treatment of copper indium gallium selenide (CIGS) film. Morphologies and composition of surface and cross section were observed by scanning electron microscopy (SEM) equipped with Energy Dispersive Spectroscopy (EDS). X-ray diffraction (XRD) and Raman spectra were used to investigate film structure. Depth profiles of element distributions were detected by Auger electron spectroscopy (AES). A double-layer structure was formed in the film by selenizing metallic precursor at 450 °C. Further annealing at 600 °C in pure argon enhanced gallium diffusion from the bottom to the top of the film, while additional selenium in the annealing had a negative effect. A MoSe{sub 2} layer was detected between CIGS and Mo layers with annealing in additional Se. The annealing treatment also significantly modified the film surface morphology. A large amount of triangular and polygon shaped islands were observed by SEM. That might be due to different nucleation kinetics for different crystal facets.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Temperature- and composition-dependent hydrogen diffusivity in palladium from statistically-averaged molecular dynamics

Solid-state hydrogen storage materials undergo complex phase transformations whose kinetics is often limited by hydrogen diffusion. Among metal hydrides, palladium hydride undergoes a diffusional phase transformation upon hydrogen uptake, during which the hydrogen diffusivity varies with hydrogen composition and temperature. Here we perform robust statistically-averaged molecular dynamics simulations to obtain a well-converged analytical expression for hydrogen diffusivity in bulk palladium that is valid throughout all stages of the reaction. Our studies confirm significant dependence of the diffusivity on composition and temperature that elucidate key trends in the available experimental measurements. Whereas at low hydrogen compositions, a single process dominates, at high hydrogen compositions, diffusion is found to exhibit behavior consistent with multiple hopping barriers. Further analysis, supported by nudged elastic band computations, suggests that the multi-barrier diffusion can be interpreted as two distinct mechanisms corresponding to hydrogen-rich and hydrogen-poor local environments.

36 MATERIALS SCIENCE↗