Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atom migration”

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

Superatoms as Superior Catalysts: ZrO versus Pd

Abstract Single‐atom catalysts are the focus of studies for over a decade due to their enhanced reactivity at smaller sizes. However, they have limitations as they offer only one active site, which may not be sufficient for reactions requiring the co‐adsorption of multiple reactants. Additionally, atoms can migrate on a substrate and coalesce, resulting in decreased reactivity. Here, an alternate path, a single‐superatom catalyst is provided. Superatoms are clusters of atoms that mimic the chemistry of atoms even if they do not contain a single atom whose chemistry they mimic. Motivated by an experimental paper on the photoelectron‐spectroscopy of negatively charged ions where ZrO is found to mimic properties of a Pd atom, first the reaction of Pd and ZrO with small molecules in the gas‐phase is studied and found that ZrO not only mimics the chemistry of Pd, but is able to activate these molecules more strongly than Pd. A detailed first‐principles study of CO 2 reduction (CO 2 ‐RR) and hydrogen evolution reactions (HER) on Pd and ZrO supported on graphene, Au(111), and Cu(111) surfaces shows that superatoms are indeed superior catalysts. The ability to design numerous superatoms by varying size and composition offers a promising new paradigm for catalyst design and synthesis.

Chemistry↗
DOE OSTI2023

Atomic Dynamics of Multi‐Interfacial Migration and Transformations

Redox-induced interconversions of metal oxidation states typically result in multiple phase boundaries that separate chemically and structurally distinct oxides and suboxides. Directly probing such multi-interfacial reactions is challenging because of the difficulty in simultaneously resolving the multiple reaction fronts at the atomic scale. Using the example of CuO reduction in H 2 gas, a reaction pathway of CuO → monoclinic m-Cu 4 O 3 → Cu 2 O is demonstrated and identifies interfacial reaction fronts at the atomic scale, where the Cu 2 O/m-Cu 4 O 3 interface shows a diffuse-type interfacial transformation; while the lateral flow of interfacial ledges appears to control the m-Cu 4 O 3 /CuO transformation. Together with atomistic modeling, it is shown that such a multi-interface transformation results from the surface-reaction-induced formation of oxygen vacancies that diffuse into deeper atomic layers, thereby resulting in the formation of the lower oxides of Cu 2 O and m-Cu 4 O 3 , and activate the interfacial transformations. In conclusion, these results demonstrate the lively dynamics at the reaction fronts of the multiple interfaces and have substantial implications for controlling the microstructure and interphase boundaries by coupling the interplay between the surface reaction dynamics and the resulting mass transport and phase evolution in the subsurface and bulk.

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

Zirconium Carbide Mediates Coke‐Resistant Methane Dry Reforming on Nickel‐Zirconium Catalysts

Abstract Graphitic deposits anti‐segregate into Ni 0 nanoparticles to provide restored CH 4 adsorption sites and near‐surface/dissolved C atoms, which migrate to the Ni 0 /ZrO 2 interface and induce local Zr x C y formation. The resulting oxygen‐deficient carbidic phase boundary sites assist in the kinetically enhanced CO 2 activation toward CO(g). This interface carbide mechanism allows for enhanced spillover of carbon to the ZrO 2 support, and represents an alternative catalyst regeneration pathway with respect to the reverse oxygen spillover on Ni‐CeZr x O y catalysts. It is therefore rather likely on supports with limited oxygen storage/exchange kinetics but significant carbothermal reducibility.

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

Phosphorus diffusion and deactivation during SiGe oxidation

Dopant profiles near the semiconductor–oxide interface are critical for microelectronic device performance. As the incorporation of Si 1−x Ge x into transistors continues to increase, it is necessary to understand the behavior of dopants in Si 1−x Ge x . In this paper, the diffusion and electrical activation of phosphorus within a strained, single-crystal Si 0.7 Ge 0.3 layer on Si during oxidation are reported. Both layers were uniformly doped, in situ, with an average phosphorus concentration of 4 × 10 19 atoms/cm 3 . After high-temperature oxidation, secondary ion mass spectrometry measurements revealed that the bulk of the phosphorus diffuses out of only the SiGe layer and segregates at the oxidizing SiGe–SiO 2 interface. Hall effect measurements corroborate the observed phosphorus loss and show that the phosphorus diffusing to the oxidizing interface is electrically inactive. Through density functional theory (DFT) calculations, it is shown that phosphorus interstitials prefer sites near the SiGe–SiO 2 interface. Finally, based on a combination of experimental data and DFT calculations, we propose that the phosphorus atoms are displaced from their lattice sites by Ge interstitials that are generated during SiGe oxidation. The phosphorus atoms then migrate toward the SiGe–SiO 2 interface through a novel mechanism of hopping between Ge sites as P–Ge split interstitials. Once they reach the interface, they are electrically inactive, potentially in the form of interstitial clusters or as part of the reconstructed interface or oxide.

Thornton, Chappel S.↗
DOE OSTI2024

Initial Development of Fusion Magnet Simulation Capabilities for Performance and Safety Evaluation Using the MOOSE Framework

Fusion energy holds the promise of being a transformative technology as a carbon-neutral, sustainable source of energy. Whole device modeling and the development of fusion digital twins will be increasingly important for emerging fusion device concepts at both national laboratories and within the commercial fusion industry. However, meeting the challenge of whole device modeling of fusion energy devices requires robust, multiphysics, multiscale modeling and simulation technologies capable of running on large-scale supercomputers. Detailed analysis of individual systems at-scale is also required to ensure safe and efficient operation as well as provide the safety basis for future device designs and licensing activities. In a tokamak, toroidal and poloidal magnets confine and shape the fusion plasma to promote the fusion reaction. High plasma temperatures and high magnetic field requirements in modern design concepts (leading to high amounts of energy stored within each magnet) impose electrical, thermal, and mechanical loads on the magnet components, which in turn impacts the safety considerations of the magnet and their supporting systems. Idaho National Laboratory (INL) has a history of working in this space, including development and benchmarking of the Magnetic System Circuitry Analysis Program (MSCAP) and Magnet Arcing (MAGARC) codes to study magnet quench events; notably, MAGARC was used to study quenching during the ITER Engineering Design Activity. However, these legacy codes and capabilities are not parallel and scalable, and new tools are required for future advances in this area, which leads to the INL-developed Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Developed originally for fission reactor systems under United States Department of Energy, Office of Nuclear Energy modeling and simulation programs, the MOOSE framework has been well-suited to multiscale, multiphysics modeling and simulation needs for nuclear systems. The framework is open-source, well-tested, under continuous development and deployment, and developed to a Nuclear Quality Assurance, Level 1 software quality standard. MOOSE has also been used in the fusion space previously in several projects: INL’s Tritium Migration Analysis Program, Version 8 (TMAP8) for tritium migration, UK Atomic Energy Authority’s A Unified Resource for OpenMC (fusion) Reactor Applications (AURORA) code for fusion thermo-mechanical and neutronics analysis, and Argonne National Laboratory’s Cardinal for high-fidelity computational fluid dynamics and neutronics. However, to model superconducting magnets, several MOOSE enhancements are required: additions to the current MOOSE electromagnetic capabilities, new material libraries for superconductors of interest (such as YBCO), as well as fusion-specific models for thermo-mechanics. This talk will discuss initial development activities to build these capabilities in MOOSE, focusing on initial validation and benchmarking activities. Proposed coupling workflows and future work to support the simulation of fusion magnets and magnet structural assemblies for performance and safety evaluation in MOOSE will also be discussed.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Shear strain alters the structure and migration mechanism of self-interstitial atoms in copper

Here, we use atomistic modeling to show that externally applied shear strain causes the lowest energy self-interstitial atom (SIA) structure in copper (Cu) to change from a <100>-type dumbbell to a <110>-type dumbbell. Concurrently, SIA migration switches from the 3-D random walk characteristic of <100>-type dumbbells to a 1-D mechanism analogous to that of crowdion SIAs. Furthermore, the relative energies of these two dumbbell structures as a function of strain are well predicted using elastic dipole tensors computed at zero strain, indicating that examination of these tensors may be used to assess the likelihood of strain-induced SIA structure transitions in other materials. Changes in lowest energy SIA structures and associated migration mechanisms stand to impact predictions of SIA behavior in irradiated solids.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Thermal gradient effect on helium and self-interstitial transport in tungsten

First-wall materials in a fusion reactor are expected to withstand harsh conditions, with high heat and particle fluxes that modify the materials microstructure. These fluxes will create strong gradients of temperature and concentration of diverse species. Besides the He ash and the hydrogenic species, neutron particles generated in the fusion reaction will collide with the material creating intrinsic defects, such as vacancies, self-interstitials atoms (SIAs), and clusters of such point defects. These defects and the He atoms will then migrate in the presence of the aforementioned gradients. In this study, we use nonequilibrium molecular dynamics to analyze the transport of He and SIAs in the presence of a thermal gradient in tungsten. We observe that, in all cases, the defects and impurity atoms tend to migrate toward the hot regions of the tungsten sample. The resulting species concentration profiles are exponential distributions, rising toward the hot regions of the sample, in agreement with irreversible thermodynamics analysis. For both He atoms and SIAs, we find that the resulting species flux is directed opposite to the heat flux, indicating that species transport is governed by a Soret effect (thermal-gradient-driven diffusion) characterized by a negative heat of transport that drives species diffusion uphill (from the cooler to the hot regions of the sample). Here, we demonstrate that the steady-state species profiles obtained accounting for the Soret effect vary significantly from those where temperature-gradient-driven transport is not considered and discuss the implications of such a Soret effect on the response to plasma exposure of plasma-facing tungsten.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The effect of composition on helium bubble nucleation in complex concentrated alloys

Molecular dynamics simulations of helium implantation have been performed in the MoNbTaTi complex concentrated alloy (CCA) to determine the impact of composition on helium clustering and bubble nucleation. It was found that compositions with lower interstitial atomic volume resulted in smaller cluster sizes, higher helium-to-vacancy ratios, and larger migration barriers. As the atomic volume decreases, it is harder to accommodate additional helium in the lattice and for helium to diffuse between interstitial sites, resulting in smaller cluster sizes. In particular, increasing molybdenum directly correlated with a decreased interstitial atomic volume and increased migration barrier. Niobium exhibited the opposite trend, where more niobium in the material resulted in larger interstitial atomic volumes and cluster sizes but lower migration barriers. The atomic volume was determined to be an indicator of the susceptibility of the material to helium bubble growth and could potentially be used as a metric to screen compositions for resistance to helium damage. In this way, properties such as available interstitial volumes could be used to screen materials more rapidly given the vast compositional space of CCAs.

McCarthy, Megan Jeanne [Sandia National Laboratori↗
DOE OSTI2023

Mechanism of Nitrogen-Carbon Bond Formation from Iron(IV) Disilylhydrazido Intermediates during N 2 Reduction

We recently reported a reaction sequence that activates C–H bonds in simple arenes as well as the N–N triple bond in N 2 , delivering the aryl group to N 2 to form a new N–C bond ( Nature 2020 , 584, 221). This enables the transformation of abundant feedstocks (arenes and N 2 ) into N-containing organic compounds. The key N–C bond forming step occurs upon partial silylation of N 2 . However, the pathway through which reduction, silylation, and migration occurred was unknown. Here, in this work, we describe synthetic, structural, magnetic, spectroscopic, kinetic, and computational studies that elucidate the steps of this transformation. N 2 must be silylated twice at the distal N atom before aryl migration can occur, and sequential silyl radical and silyl cation addition is a kinetically competent pathway to a formally iron(IV)–NN(SiMe 3 ) 2 intermediate that can be isolated at low temperature. Kinetic studies show its first-order conversion to the migrated product, and DFT calculations indicate a concerted transition state for migration. The electronic structure of the formally iron(IV) intermediate is examined using DFT and CASSCF calculations, which reveal contributions from iron(II) and iron(III) resonance forms with oxidized NNSi 2 ligands. The depletion of electron density from the Fe-coordinated N atom makes it electrophilic enough to accept the incoming aryl group. This new pathway for the N–C bond formation offers a method for functionalizing N 2 using organometallic chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1981

Cluster adsorption on amorphous and crystalline surfaces - A molecular dynamics study of model Pt on Cu and model Pd on Pt

Molecular dynamics was used to study the structure, dispersion and short-time behavior of ten-atom clusters adsorbed onto amorphous and crystalline substrates, in which the cluster atoms differed from the substrate atoms. Two adatom-substrate model systems were chosen; one, in which the interaction energy between adatom pairs was greater than that between substrate pairs, and the other, in which the reverse was true. At relatively low temperature ranges, increased dispersion of cluster atoms occurred: (a) on the amorphous substrate as compared to the FCC(100) surface, (b) with increasing reduced temperature, and (c) with adatom-substrate interaction energy stronger than adatom-adatom interaction. Two-dimensional clusters (rafts) on the FCC(100) surface displayed migration of edge atoms only, indicating a mechanism for the cluster rotation and shape changes found in experimental studies.

Garofalini, S. H.↗
DOE OSTI2025

Mo Atom Rearrangement Drives Layer-Dependent Reactivity in Two-Dimensional MoS 2

Two-dimensional (2D) materials offer a valuable platform for manipulating and studying chemical reactions at the atomic level, owing to the ease of controlling their microscopic structure at the nanometer scale. While extensive research has been conducted on the structure-dependent chemical activity of 2D materials, the influence of structural transformation during the reaction has remained largely unexplored. In this work, we report the layer-dependent chemical reactivity of MoS 2 during a nitridation atomic substitution reaction and attribute it to the rearrangement of Mo atoms. Our results show that the chemical reactivity of MoS 2 decreases as the number of layers is reduced in the few-layer regime. In particular, monolayer MoS 2 exhibits significantly lower reactivity compared with its few-layer and multilayer counterparts. Atomic-resolution transmission electron microscopy (TEM) reveals that MoN nanonetworks form as reaction products from monolayer and bilayer MoS 2 , with the continuity of the MoN crystals increasing with layer number, consistent with the local conductivity mapping data. The layer-dependent reactivity is attributed to the relative stability of the hypothetically formed MoN phase, which retains the number of Mo atomic layers present in the precursor. Specifically, the low chemical reactivity of monolayer MoS 2 is attributed to the high energy cost associated with Mo atom diffusion and migration necessary to form multilayer Mo lattices in the thermodynamically stable MoN phase. In conclusion, this study underscores the critical role of lattice rearrangement in governing chemical reactivity and highlights the potential of 2D materials as versatile platforms for advancing the understanding of materials chemistry at the atomic scale.

Chemical reactivity↗
NASA NTRS1998

Computational Studies of the Interaction of H/H2 with Diamond and Silicon Surfaces

The interaction of hydrogen atoms and molecules with diamond and silicon surfaces is important in several important applications. Two areas that we are interested are: 1) tribology (molecular level friction) and 2) the role of H atoms in silicon chemical vapor deposition (CVD). In the tribology area, H atoms can be used to tie off dangling bonds, which otherwise form bonds between adjacent surfaces, and lead to resistance to sliding the surfaces by each other. Processes which are important in understanding molecular level friction include barriers to addition of H/H2 to the surface and barriers to migration of H atoms on the surface. In the silicon CVD area, we have studied the process of H2 elimination from the 100 surface of silicon. Cluster models for the dime surfaces of diamond are presented. The unrelaxed 100 surface has carbene like surface carbon atoms; however, for the relaxed surface these dimerize to give rows of surface dimers and there is a significant amount of p bonding between the radical orbitals of the dimer. The 110 surface has zig-zag rows of carbon atoms with a dangling bond on each carbon atom. These dangling bonds are hybridized away from each other and thus interact less strongly than for the 100 surface. Finally, the 111 surface has surface C atoms arranged in a triangular pattern and the surface dangling bonds are well separated from each other (second nearest neighbor distance) leading to almost no interaction between adjacent dangling bonds. These qualitative features may be quantified by computing the overlap of adjacent dangling bonds in a GVB(pp) calculation. The overlaps are 0.462, 0.292, and 0.016 for the diamond 100, 110, and 111 surfaces, respectively.

Walch, Stephen P.↗
DOE OSTI2025

Nanocluster Active Sites Formed on Heterogeneous Thermal Catalysts and Electrocatalysts by Operando Reactive Environments

Here, in this Viewpoint, we summarize our recent studies in this endeavor, ending with a forward-looking perspective into the remaining challenges and open questions for consideration by the catalysis community. Through a computationally efficient framework of systematically investigating metal atom ejection from and migration onto metal surfaces, leading to the formation of metal atom clusters, we show that the surface of a metal catalyst is dynamic and much less rigid than previously thought. We suggest that the operando formed metal atom clusters are novel active sites, which can dominate the activity of steps or kinks (let aside terraces) from where their constituent atoms were ejected. Through more readily responding to reaction conditions by forming catalytically active clusters, solid metal catalysts approach not only homogeneous catalyst motifs, thus establishing a bridge between homogeneous and heterogeneous catalysis, but also liquid metals which have been suggested to be better than traditional heterogeneous catalysts.

adsorption↗
DOE OSTI2025

Selenium Migration and Local Structures in Cu‐Doped CdSeTe Solar Cells after Aging

Selenium grading plays a critical role in state-of-the-art Cadmium Telluride photovoltaic cells by enhancing long-wavelength absorption and extending minority carrier lifetimes —key to enabling the current performance record of 23.08%. However, very little is understood about selenium motion. In this study, a comprehensive, multimodal, and multiscale approach is employed to investigate Se migration and local structural changes in copper (Cu)-doped CdSeTe solar cells subjected to accelerated stress. X-ray fluorescence (XRF) microscopy shows unexpected levels of Se diffusion after 500 h under heat (75°C) and light (0.8 suns, 80 mW/cm 2 ), suggesting the coexistence of fast and slow diffusion channels even at low temperatures, with unexpectedly low activation energies (<0.85 eV). X-ray Absorption Near Edge Structure (XANES) analysis indicates a preferential migration of Se atoms to anionic lattice sites and a reduction in Se-Cl co-passivation at Te-terminated dislocation cores. Furthermore, these findings point to a reconfiguration of Se local environments and highlight the potential role of extended structural defects in enabling Se transport at low temperatures. Additionally, XANES results suggest that the presence of metallic Cu across the absorber layer may contribute to back-contact degradation and reduced hole density in both fresh and aged devices.

14 SOLAR ENERGY↗
DOE OSTI2022

Reactivity of Iridium Complexes of a Triphosphorus-Pincer Ligand Based on a Secondary Phosphine. Catalytic Alkane Dehydrogenation and the Origin of Extremely High Activity

The selective functionalization of alkanes and alkyl groups is a major goal of chemical catalysis. Toward this end, a bulky triphosphine with a central secondary phosphino group, bis(2-di-t-butyl-phosphinophenyl)phosphine ( tBu P H PP), has been synthesized. When complexed to iridium, it adopts a meridional (“pincer”) configuration. The secondary phosphino H atom can undergo migration to iridium to give an anionic phosphido-based–pincer ( tBu PPP) complex. Stoichiometric reactions of the ( tBu PPP)Ir complexes reflect a distribution of steric bulk around the iridium center in which the coordination site trans to the phosphido group is quite crowded; one coordination site cis to the phosphido is even more crowded; and the remaining site is particularly open. The ( tBu PPP)Ir precursors are the most active catalysts reported to date for dehydrogenation of n-alkanes, by about 2 orders of magnitude. The electronic properties of the iridium center are similar to that of well-known analogous ( R PCP)Ir catalysts. Accordingly, DFT calculations predict that ( tBu PPP)Ir and ( tBu PCP)Ir are, intrinsically, comparably active for alkane dehydrogenation. While dehydrogenation by ( R PCP)Ir proceeds through an intermediate trans-(PCP)IrH 2 (alkene), ( tBu PPP)Ir follows a pathway proceeding via cis-(PPP)IrH 2 (alkene), thereby circumventing unfavorable placement of the alkene at the bulky site trans to phosphorus. ( tBu PPP)Ir and ( tBu PCP)Ir, however, have analogous resting states: square planar (pincer)Ir(alkene). Alkene coordination at the crowded trans site is therefore unavoidable in the resting states. Furthermore, the resting state of the ( tBu PPP)Ir catalyst is destabilized by the architecture of the ligand, and this is largely responsible for its unusually high catalytic activity.

02 PETROLEUM↗
DOE OSTI2022

Cost-effective carbon fiber precursor selections of polyacrylonitrile-derived blend polymers: carbonization chemistry and structural characterizations

Blending polyacrylonitrile (PAN) with plastic wastes and bio-based polymers provides a convenient and inexpensive method to realize cost-effective carbon fiber (CF) precursors. In this work, PAN-based blend precursors are investigated using ReaxFF reactive molecular dynamics simulations with respect to the formation of all-carbon rings, the evolutions of oxygen-containing and nitrogen-containing species, and the migration of carbon atoms to form turbostratic graphene. From these simulations, we identify that PAN/cellulose (CL) blend manifests the highest carbon yield and the most substantial all-carbon ring formation. This ReaxFF-based finding is confirmed by Raman and TEM experiments indicating high crystallinity for PAN/CL-derived blend CFs. We trace the pathway of gasification and carbonization of PAN/CL to elaborate the mechanism of the formation of all-carbon ring networks. We discover that the acetals of CL can catalyze the cyclization of the blend precursor, allowing for the search for CL derivatives or the other kinds of bio-based polymers with similar functionalities as alternative blends. In addition, we examine the structural characteristics using the carbon–carbon (C–C) radial distribution functions, C–C bond length distributions, and sp2 C atom ratios for the four representative precursors, i.e., PAN, oxidized PAN, PAN/nylon 6,6, and PAN/CL. Our simulation results show the most extensive all-carbon ring cluster and graphitic structure growths for PAN/CL. Here, we propose PAN/CL as a cost-effective alternative CF precursor, since (a) CL is naturally abundant and eco-friendly for production, (b) the blend precursor PAN/CL does not require oxidation treatment, (c) PAN/CL has a high carbon yield with substantial all-carbon ring formation, and (d) PAN/CL based CFs potentially provide a mechanical property enhancement.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Carboxylic acid induced restructuring of the Fe 3 O 4 (001) surface

The redox properties of Fe 3 O 4 surfaces are central to many catalytic processes and enable dynamic, reduction-induced morphological restructuring during reactions. Here, we investigate the structural changes of the Fe 3 O 4 (001) surface during the decomposition of formic and acetic acids using scanning tunneling microscopy, X-ray photoelectron spectroscopy, and density functional theory (DFT) calculations. Both acids readily deprotonate, forming ordered carboxylate overlayers on the surface. Product formation pathways involve the removal of lattice oxygen, resulting in extensive surface restructuring. For formic acid, only a modest level of surface oxygen removal (∼3%) is observed, resulting in elongated pits along the octahedral Fe rows and exhibiting an aspect ratio of ∼3. In contrast, acetic acid induces more extensive reduction, with the removal of ∼20% of surface oxygen, yielding significantly larger pits while maintaining a similar aspect ratio. Repeated exposure to acetic acid further enlarges the pits, indicating preferential etching at step edges. DFT calculations reveal a mechanistic sequence in which lattice oxygen removal destabilizes adjacent Fe atoms, promoting their migration into the bulk and subsequent pit propagation and step edge formation. Together, these findings provide atomistic insights into the coupling between carboxylic acid conversion and oxide surface restructuring, underscoring the strong interplay between redox chemistry and morphological changes on catalytically active Fe 3 O 4 surfaces.

Ortiz-Garcia, José J. [Pacific Northwest National ↗