Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Cluster diffusion”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Interpretation of Ion Irradiation and Neutron Irradiation Damage in Additively Manufactured 316 Stainless Steel using Multiscale Modeling

The accelerated adoption of nuclear energy necessitates advanced manufacturing technologies, such as additive manufacturing, to meet heightened supply chain requirements and support innovative reactor technologies. Due to the unique microstructural characteristics of additively manufactured materials under distinct solidification conditions, comprehensive evaluation of their performance in reactor environments is essential. The Advanced Materials and Manufacturing Technologies program under the Department of Energy's Office of Nuclear Energy focuses on understanding the irradiation performance and damage evolution of laser powder bed fusion 316 stainless steel, with an emphasis on integrating ion and neutron irradiation data to accelerate the development and qualification of materials for advanced nuclear reactor applications. While ion irradiation is a cost- and time-effective method, modeling and simulation are required to interpret the data for the broader range of irradiation conditions encountered in advanced reactors. In fiscal year 2025, integrated multiscale modeling and simulations were conducted to assess irradiation damage in additively manufactured 316 stainless steel. Key outcomes include predictions of chromium enrichment at grain boundaries, nickel enrichment at dislocation cell walls and void surfaces, and heterogeneous void evolution under ion and neutron irradiation conditions. Cluster dynamics simulations revealed the coarsening of voids at high irradiation temperatures and the suppression of void growth by high network dislocation density, while also demonstrating significant growth and coarsening of voids and self-interstitial atom loops at low dose rates. Machine learning-accelerated atomistic simulations highlighted the impact of the local environment and chromium concentration on vacancy diffusivity, providing key insights on the influence of composition on void swelling and radiation-induced segregation. Additionally, molecular dynamics simulations demonstrated the presence of defect production bias and a significant effect of carbon content on defect cluster behavior. These combined efforts aim to predict the performance of additively manufactured materials under various reactor conditions, supporting their qualification for nuclear reactor applications by interpreting ion irradiation data. This report underscores the potential of integrated multiscale modeling to analyze ion irradiation data in the effort to accelerate the qualification of additively manufactured materials for nuclear reactor components.

316 stainless steel↗
DOE OSTI2022

Diffusion Monte Carlo approaches for studying nuclear quantum effects in fluxional molecules

Abstract Diffusion quantum Monte Carlo (DMC) provides a powerful approach for obtaining the ground state energy and wave function of molecules, ions, and molecular clusters. The approach is uniquely well suited for studies of fluxional molecules, which undergo large amplitude vibrational motions even in their ground state. In contrast to the electronic structure problem, where the wave function must be antisymmetric with respect to exchange of any pair of electrons, the wave function for the ground vibrational state is nodeless. This greatly simplifies the application of DMC for vibrational problems. Because there is not a single potential function that can be used to describe the intramolecular and intermolecular interactions in all molecular systems, most methods that are used to describe nuclear quantum effects rely on a carefully chosen zero‐order description of the molecular vibrations. In contrast, DMC calculations can be performed in Cartesian coordinates, making the DMC algorithm easily transferable between different chemical systems. In this contribution, the theory that underlies DMC will be discussed along with important considerations for performing DMC calculations. Extensions for evaluating vibrationally excited states and molecular properties are also discussed. Insights that can be obtained from DMC calculations are illustrated in the context of the protonated water clusters. This article is categorized under: Molecular and Statistical Mechanics > Molecular Dynamics and Monte‐Carlo Methods Theoretical and Physical Chemistry > Spectroscopy

Chemistry↗
DOE OSTI2022

Development of a creep model informed by lower-length scale simulations to simulate creep in doped UO 2

Using molecular dynamics, we predict information at the atomistic scale used to develop a mechanistic UO 2 creep model for use in higher length-scale fuel performance codes. The ultimate objective of the model is to better describe the grain size dependence and therefore impact of doping on creep rates in UO 2 . In a previous NEAMS milestone, we found that Nabarro-Herring (bulk diffusional) creep was too low to capture the experimentally observed creep rates in standard UO 2 . Moreover, in that milestone, other mechanisms were explored, such as Coble (grain boundary) creep and dislocation climb, with each mechanism exhibiting different grain size dependencies. Again, these were orders of magnitude too low to describe the experimental creep rates. In this work, we address the previous assumptions made for the Coble creep mechanism by investigating the diffusivity of various defects at grain boundaries in UO 2 and, critically, to determine if enhanced grain boundary diffusivity allows the model to better reproduce experimental results. The diffusivity as a function of temperature for different concentrations of uranium vacancies and interstitials for bulk UO 2 have also been examined using cluster dynamics. Furthermore, using a concentration dependent segregation model, the concentration of defects at the grain boundary were predicted. This atomistic data was then input into the various creep mechanisms and the creep rates compared to the empirical MATPRO correlation (used in BISON) and experiment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Unexpected Hydroxide Ion Structure and Properties at Low Hydration

Understanding the behavior of hydroxide ions in aqueous and non-aqueous media is fundamental to many chemical, biological, and electrochemical processes. Research has primarily focused on a single fully solvated hydroxide ion, either as an isolated cluster or in bulk. This work presents the first computational study to consider hydroxide under low hydration levels in detail, where the anion may not be fully solvated. Under such conditions, we find that the anions are predominantly present as unique water-bridged hydroxide pair complexes, distinct from previously reported structures under fully hydrated conditions. Although similar hydroxide pair structures were previously reported, we analyze these structures for the first time in the disordered liquid state where they are found to be unusually stable in the presence of bulky quaternary ammonium cations. Our findings help explain the unusual diffusion behavior as well as the higher reactivity of hydroxide anions observed under low hydration conditions.

computational studies↗
DOE OSTI2021

The effect of hydrogen-bonding on flavin’s infrared absorption spectrum

Cluster and continuum solvation computational models are employed to model the effect of hydrogen bonding interactions on the vibrational modes of lumiflavin. Calculated spectra were compared to experimental Fourier-transform infrared (FTIR) spectra in the diagnostic 1450–1800 cm –1 range, where intense V C=C , V C=N , V C 2 =O , and V C 4 =O stretching modes of flavin’s isoalloxazine ring are found. Local mode analysis is used to describe the strength of hydrogen-bonding in cluster models. The computations indicate that V C=C and V C=N mode frequencies are relatively insensitive to intermolecular interactions while the V C 2 =O and V C 4 =O modes are sensitive to direct (and also indirect for V C 2 =O ) hydrogen-bonding interactions. Although flavin is neutral, basis sets without the diffuse functions provide incorrect relative frequencies and intensities. The 6-31+G* basis set is found to be adequate for this system, and there is limited benefit to considering larger basis sets. Calculated vibrational mode frequencies agree with experimentally determined frequencies in solution when cluster models with multiple water molecules are used. Accurate simulation of relative FTIR band intensities, on the other hand, requires a continuum (or possibly quantum mechanical/molecular mechanical) model that accounts for long-range electrostatic effects. Lastly, an experimental peak at ca. 1624 cm –1 that is typically assigned to the V C 2 =O vibrational stretching mode has a complicated shape that suggests multiple underlying contributions. Our calculations show that this band has contributions from both the C 6 –C 7 and C 2 = O stretching vibrations.

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

Benchmarking Correlation-Consistent Basis Sets for Frequency-Dependent Polarizabilities with Multiresolution Analysis

This paper presents the first converged frequency-dependent HF polarizability results for general molecules, on a set of 89 closed-shell atoms and molecules. The solver employs multiresolution analysis (MRA) in a multiwavelet basis to compute both ground and response states to a guaranteed precision, which are validated against independent numerical grid calculations on atoms and linear molecules. The MRA ground-state energies and response properties are used to evaluate results in correlation-consistent basis sets up to 5Z augmented with either single or double diffuse functions and core-polarization functions. Systematic trends are revealed through consideration of chemical composition as well as the use of machine learning to cluster convergence trends, the latter suggesting the possibility of learning and correcting basis-set error.

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

Facile electrocatalytic proton reduction by a [Fe–Fe]-hydrogenase bio-inspired synthetic model bearing a terminal CN — ligand

An azadithiolate bridged CN — bound pentacarbonyl bis-iron complex, mimicking the active site of [Fe–Fe] H 2 ase is synthesized. The geometric and electronic structure of this complex is elucidated using a combination of EXAFS analysis, infrared and Mössbauer spectroscopy and DFT calculations. The electrochemical investigations show that complex 1 effectively reduces H + to H 2 between pH 0–3 at diffusion-controlled rates (10 11 M —1 s —1 ) i.e. 10 8 s —1 at pH 3 with an overpotential of 140 mV. Electrochemical analysis and DFT calculations suggests that a CN— ligand increases the p K a of the cluster enabling hydrogen production from its Fe(I)–Fe(0) state at pHs much higher and overpotential much lower than its precursor bis-iron hexacarbonyl model which is active in its Fe(0)–Fe(0) state. The formation of a terminal Fe–H species, evidenced by spectroelectrochemistry in organic solvent, via a rate determining proton coupled electron transfer step and protonation of the adjacent azadithiolate, lowers the kinetic barrier leading to diffusion controlled rates of H 2 evolution. The stereo-electronic factors enhance its catalytic rate by 3 order of magnitude relative to a bis-iron hexacarbonyl precursor at the same pH and potential.

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

Argyrodite-type advanced lithium conductors and transport mechanisms beyond paddle-wheel effect

Abstract Development of next-generation solid-state Li-ion batteries requires not only electrolytes with high room-temperature (RT) ionic conductivities but also a fundamental understanding of the ionic transport in solids. In spite of considerable work, only a few lithium conductors are known with the highest RT ionic conductivities ~ 0.01 S/cm and the lowest activation energies ~0.2 eV. New design strategy and novel ionic conduction mechanism are needed to expand the pool of high-performance lithium conductors as well as achieve even higher RT ionic conductivities. Here, we theoretically show that lithium conductors with RT ionic conductivity over 0.1 S/cm and low activation energies ~ 0.1 eV can be achieved by incorporating cluster-dynamics into an argyrodite structure. The extraordinary superionic metrics are supported by conduction mechanism characterized as a relay between local and long-range ionic diffusions, as well as correlational dynamics beyond the paddle-wheel effect.

25 ENERGY STORAGE↗
DOE OSTI2020

Grain detachment and transport clogging during mineral dissolution in carbonate rocks with permeable grain boundaries

In this paper, we explore the effect of permeable grain boundaries on carbonate rock dissolution in order to improve our understanding of grain detachment and migration during reactive flow. To do so we investigated the effects of grain size and enhanced reactivity along grain boundaries on global dissolution kinetics. Variations of permeability and porosity were calculated, and exponential relationships were observed in fractured rocks. Our model employed a reactive transport framework based on the Darcy-Brinkman equation to simulate calcite dissolution in carbonate rocks composed of microporous grains and, included fluid transport along grain boundaries. The model includes fluid flow, solute transport by advection–diffusion, heterogeneous reaction between fluid and minerals and grain detachment with subsequent grain transport in macropores. The migration of solid particles due to dissolution was based on cluster analysis and local movement, and the results show that grain detachment can lead to significant decreases in permeability due to the clogging of transport pathways. Microporous media with smaller grain sizes (fine grains in this study) showed a higher average reaction rate than those with coarser grains. In addition, as the overall rates of geochemical processes are commonly affected by the presence of texture heterogeneities such as fractures, a single fracture (macropore) introduced into a microporous matrix composed of permeable grains and grain boundaries was modelled as a large channel connecting the inlet and outlet of the simulation domain. It was found that macropore clogging by grain detachment temporary decreased permeability, lowering the long-term global reaction rates. Lastly, we observed that increases in flow rate can reduce detachment and fracture clogging by reducing local dissolution along grain boundaries.

58 GEOSCIENCES↗
DOE OSTI2020

Point and extended defects in heteroepitaxial 𝛽−G⁡a 2 ⁢O 3 films

G⁡a 2 ⁢O 3 is emerging as an excellent potential semiconductor for high power and optoelectronic devices. However, the successful development of G⁡a 2⁢ O 3 in a wide range of applications requires a full understanding of the role and nature of its point and extended defects. Here, in this work, high quality epitaxial G⁡a 2 ⁢O 3 films were grown on sapphire substrates by metal-organic chemical vapor deposition and fully characterized in terms of structural, optical, and electrical properties. Then defects in the films were investigated by a combination of depth-resolved Doppler broadening and lifetime of positron annihilation spectroscopies and thermally stimulated emission (TSE). Positron annihilation techniques can provide information about the nature and concentration of defects in the films, while TSE reveals the energy level of defects in the bandgap. Despite very good structural properties, the films exhibit short positron diffusion length, which is an indication of high defect density and long positron lifetime, a sign for the formation of Ga vacancy related defects and large vacancy clusters. These defects act as deep and shallow traps for charge carriers as revealed from TSE, which explains the reason behind the difficulty of developing conductive G⁡a 2 ⁢O 3 films on non-native substrates. Positron lifetime measurements also show nonuniform distribution of vacancy clusters throughout the film depth. Further, the work investigates the modification of defect nature and properties through thermal treatment in various environments. It demonstrates the sensitivity of G⁡a 2 ⁢O 3 microstructures to the growth and thermal treatment environments and the significant effect of modifying defect structure on the bandgap and optical and electrical properties of G⁡a 2 ⁢O 3 .

defects↗
DOE OSTI2023

Spatially Distributed Ramp Reversal Memory in VO 2

Ramp‐reversal memory has recently been discovered in several insulator‐to‐metal transition materials where a non‐volatile resistance change can be set by repeatedly driving the material partway through the transition. This study uses optical microscopy to track the location and internal structure of accumulated memory as a thin film of VO 2 is temperature cycled through multiple training subloops. These measurements reveal that the gain of insulator phase fraction between consecutive subloops occurs primarily through front propagation at the insulator‐metal boundaries. By analyzing transition temperature maps, it is found, surprisingly, that the memory is also stored deep inside both insulating and metallic clusters throughout the entire sample, making the metal‐insulator coexistence landscape more rugged. This non‐volatile memory is reset after heating the sample to higher temperatures, as expected. Diffusion of point defects is proposed to account for the observed memory writing and subsequent erasing over the entire sample surface. By spatially mapping the location and character of non‐volatile memory encoding in VO 2 , this study results enable the targeting of specific local regions in the film where the full insulator‐to‐metal resistivity change can be harnessed in order to maximize the working range of memory elements for conventional and neuromorphic computing applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Atom Efficiency of Pd Sites for Methane Combustion: Single Atom Catalysts Versus Nanocatalysts

Methane combustion is an important reaction for energy production and methane removal from the atmosphere. This reaction highly relies on the use of noble metal Pd-based catalysts, which therefore drives the pursuit of catalysts with high atomic dispersion and activity. In this work, Pd/ceria catalysts dominated with Pd single atoms or nanosized Pd clusters (∼1 nm) are prepared and characterized by combining high-resolution high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM), in situ diffuse reflectance infrared Fourier-transform spectroscopy (DRIFTS), and Raman and X-ray absorption spectroscopy (XAS) techniques. By comparing the turnover frequencies (TOF; per every Pd atom) of Pd/ceria single atom catalysts and nanocatalysts, it is found that the atom efficiency of Pd is increased by 10 ∼30 times from single atom catalysts to nanocatalysts. For Pd single atom catalysts, although their activity can be tuned by changing the local structures, the intrinsic activity and number of active sites need to be further improved by engineering the surfaces of supports. For nanosized Pd species, despite the high TOF, the Pd atoms in the bulk structure are not directly participating in the catalytic reaction. Furthermore, this work highlights the importance of increasing the intrinsic activity of individual noble atoms, as well as the homogeneity of their local structures. For Pd/ceria systems reported in this work, our results indicate that from the application point of view, at the current stage, it is not practical to replace Pd nanocatalysts with single atom catalysts for methane combustion.

Pd↗
DOE OSTI2022

Reaction–drift–diffusion models from master equations: application to material defects

We present a general method to produce well-conditioned continuum reaction–drift–diffusion equations directly from master equations on a discrete, periodic state space. We assume the underlying data to be kinetic Monte Carlo models (i.e. continuous-time Markov chains) produced from atomic sampling of point defects in locally periodic environments, such as perfect lattices, ordered surface structures or dislocation cores, possibly under the influence of a slowly varying external field. Our approach also applies to any discrete, periodic Markov chain. Here, the analysis identifies a previously omitted non-equilibrium drift term, present even in the absence of external forces, which can compete in magnitude with the reaction rates, thus being essential to correctly capture the kinetics. To remove fast modes which hinder time integration, we use a generalized Bloch relation to efficiently calculate the eigenspectrum of the master equation. A well conditioned continuum equation then emerges by searching for spectral gaps in the long wavelength limit, using an established kinetic clustering algorithm to define a proper reduced, Markovian state space.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The K2 M67 Study: Precise Mass for a Turnoff Star in the Old Open Cluster M67

We present a study of the bright detached eclipsing main-sequence binary WOCS 11028 (Sanders 617) in the open cluster M67. Although the binary has only one eclipse per orbital cycle, we show that the masses of the stars can be derived very precisely thanks to a strong constraint on the orbital inclination: M {sub A} = 1.222 ± 0.006M {sub ⊙} and M {sub B} = 0.909 ± 0.004M {sub ⊙}. We use a spectral energy distribution fitting method to derive characteristics of the component stars in lieu of the precise radii that would normally be derived from a doubly eclipsing binary. The deconvolution of the SEDs reveals that the brighter component of the binary is at the faint turnoff point for the cluster—a distinct evolutionary point that occurs after the convective core has been established and while the star is in the middle of its movement toward lower surface temperatures, before the so-called hook at the end of the main sequence. The measurements are in distinct disagreement with evolution models at solar metallicity: higher metal abundances are needed to reproduce the characteristics of WOCS 11028 A. We discuss the changes to model physics that are likely to be needed to address the discrepancies. The clearest conclusions are that diffusion is probably necessary to reconcile spectroscopic abundances of M67 stars with the need for higher metallicity models and that reduced strength convective overshooting is occurring for stars at the turnoff. At super-solar bulk metallicity, various indicators agree on a cluster age between about 3.5 and 4.0 Gyr.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Cluster Frustration in the Breathing Pyrochlore Magnet LiGaCr 4 S 8

We present a comprehensive neutron scattering study of the breathing pyrochlore magnet LiGaCr 4 S 8 . We observe an unconventional magnetic excitation spectrum with a separation of high and low-energy spin dynamics in the correlated paramagnetic regime above a spin-freezing transition at 12(2) K. By fitting to magnetic diffuse-scattering data, we parameterize the spin Hamiltonian. We find that interactions are ferromagnetic within the large and small tetrahedra of the breathing pyrochlore lattice, but antiferromagnetic further-neighbor interactions are also essential to explain our data, in qualitative agreement with density-functional theory predictions [Ghosh et al., npj Quantum Mater. 4, 63 (2019)]. We explain the origin of geometrical frustration in LiGaCr 4 S 8 in terms of net antiferromagnetic coupling between emergent tetrahedral spin clusters that occupy a face-centered lattice. Furthermore, our results provide insight into the emergence of frustration in the presence of strong further-neighbor couplings, and a blueprint for the determination of magnetic interactions in classical spin liquids.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Gamma radiation-induced defects in KCl, MgCl 2 , and ZnCl 2 salts at room temperature

Room temperature post-irradiation measurements of diffuse reflectance and electron paramagnetic resonance spectroscopies were made to characterize the long-lived radiation-induced species formed from the gamma irradiation of solid KCl, MgCl 2 , and ZnCl 2 salts up to 100 kGy. The method used showed results consistent with those reported for electron and gamma irradiation of KCl in single crystals. Thermal bleaching of irradiated KCl demonstrated accelerated disaggregation of defect clusters above 400 K, due to decomposition of Cl 3 - . The defects formed in irradiated MgCl 2 comprised a mixture of Cl 3 - , F-centers, and Mg + associated as M-centers. Further, Mg metal cluster formation was also observed at 100 kGy, in addition to accelerated destruction of F-centers above 20 kGy. Irradiated ZnCl 2 afforded the formation of Cl 2 - due to its high ionization potential and crystalline structure, which decreases recombination. The presence of aggregates in all cases indicates the high diffusion of radicals and the predominance of secondary processes at 295 K. Additionally, thermal bleaching studies showed that chloride aggregates’ stability increases with the ionization potential of the cation present. The characterization of long-lived radiolytic transients of pure salts provides important information for the understanding of complex salt mixtures under the action of gamma radiation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Insights into palladium poisoning of Cu/SSZ-13 selective catalytic reduction catalysts

Here, the impacts of Pd poisoning to the activity, selectivity, and hydrothermal stability of Cu/SSZ-13 selective catalytic reduction (SCR) catalysts are reported. Pd lowers DeNOx efficiency of Cu/SSZ-13 via two mechanisms: (1) displacing SCR active sites in the form of isolated Pd-ions, and (2) catalyzing non-selective NH 3 oxidation in the form of PdO. The first mechanism works by the diffusion of isolated Pd-ions into chabazite cages to displace ZCu II OH SCR active species, and it occurs on Cu/SSZ-13 catalysts with and without external surface CuO clusters in similar fashion. In contrast, the second mechanism works differently with and without external surface CuO clusters. For Cu/SSZ-13 catalyst without external CuO clusters (i.e., primarily isolated Cu-ions), PdO leads to significant decrease of DeNOx efficiency at reaction temperatures above ~400°C due to its non-selective NH 3 oxidation activity at high temperatures. However, this poisoning effect becomes much less impactful on Cu/SSZ-13 catalyst containing CuO clusters. This is due to the formation of CuO-PdO solid solution via interactions between CuO and PdO, which reduces non-selective NH 3 oxidation potential. Furthermore, this solid solution formation even mitigates adverse effects caused by hydrothermal aging. Hence, the poisoning effects of Pd are closely related to Cu speciation and spatial distribution of a Cu/SSZ-13 catalyst. Finally, this study suggests a strategy in eliminating Pd poisoning, that is, the introduction of an oxide phase that effectively traps PdO but does not adversely influence SCR.

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

Magneto-structural studies of an unusual [Mn III Mn II Gd III (OR) 4 ] 4– partial cubane from 2,2'-bis- p- t Bu-calix[4]arene

Reaction of 2,2'-bis- p - t Bu-calix[4]arene (H 8 L) with MnCl 2 ·4H 2 O, GdCl 3 ·6H 2 O and 2,6-pyridinedimethanol (H 2 pdm) affords [Mn III Mn II Gd III (H 3 L)(pdmH)(pdm)(MeOH) 2 (dmf)]·3MeCN·dmf ( 3 ·3MeCN·dmf) upon vapour diffusion of MeCN into the basic dmf/MeOH mother liquor. 3 crystallises in the tetragonal space group P 4 1 2 1 2 with the asymmetric unit comprising the entire cluster. The highly unusual core contains a triangular arrangement of Mn III Mn II Gd III ions housed within a [Mn III Mn II Gd III (OR) 4 ] 4– partial cubane. Magnetic susceptibility and magnetisation data reveal best fit parameters J Mn(II)–Mn(III) = +0.415 cm –1 , J Mn(III)–Gd(III) = +0.221 cm –1 , J Mn(II)–Gd(III) = –0.258 cm –1 and D Mn(III) = –4.139 cm –1 . Theoretically derived magnetic exchange interactions, anisotropy parameters, and magneto-structural correlations for 3 are in excellent agreement with the experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗