Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ATOMIC MODELS”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Irradiation-enhanced diffusion and diffusion-limited creep in U 3 Si 2

U 3 Si 2 is an advanced fuel candidate due to its relatively high fissile density and attractive thermal properties. Compared to standard UO 2 fuel, there are significant data gaps for the thermophysical and thermomechanical properties of U 3 Si 2 . Point defect concentrations and mobilities under irradiation govern a number of important fuel performance properties, such as creep and fission gas release. In this work, we utilized density functional theory (DFT) data to inform a cluster dynamics framework to predict point defect concentrations in U 3 Si 2 under irradiation. Molecular dynamics (MD) simulations were used to examine the contribution of atomic mixing during ballistic cascades to diffusion, as well as the diffusivity of U and Si at grain boundaries. These atomic scale models for diffusivity were then used to inform a creep model based on bulk (Nabarro-Herring) and grain boundary (Coble) diffusional creep, and climb-limited dislocation creep. The model compares well against available experimental data and has been implemented in the BISON fuel performance code. Finally, a demonstration case using simple power profiles has been carried out, showing that negligible creep occurs due to the low temperatures experienced by U 3 Si 2 in-reactor, a consequence of its high thermal conductivity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Development and application of a uranium mononitride (UN) potential: Thermomechanical properties and Xe diffusion

We report that atomic-scale modeling of thermophysical and defect properties of uranium mononitride (UN) plays an important role in establishing a better understanding and improved models of UN fuel performance. Having an accurate interatomic potential is crucial for generating reliable data at finite temperatures using molecular dynamic simulations. We report a new interatomic potential for UN, based on a combination of many-body and pairwise interactions, a simple form that we later show could be easily adapted to include Xe-U and Xe-N interactions, i.e., generating a U-N-Xe interatomic potential. The potential was fitted to experimental thermal expansion and single crystal elastic constants, as well as Frenkel, Schottky, anti-Schottky, and antisite pair reaction energies from density functional theory (DFT) calculations. Using the potential, we successfully reproduced experimental lattice parameters, thermal expansion, single crystal elastic constants, and temperature dependent heat capacity. The potential also performs reasonably well in reproducing the energy of the aforementioned stoichiometric defect reactions and defect migration barriers calculated using DFT. However, the potential underestimates the energy difference between the tetrahedral and dumbbell uranium interstitials, and a more complex potential form might be needed to overcome this issue. The potential was also used to predict UN single crystal elastic constants and elastic properties at different temperatures, showing that UN becomes softer and more compressible with increasing temperature. We also compare our potential against literature data from previous empirical potentials, demonstrating similar or better behavior depending on the property of interest. To enable the simulation of Xe in UN a Buckingham potential has been fitted to DFT-derived Xe incorporation energies. The potential was then used to determine the activation energy for Xe diffusion due to various Xe-containing defects, with {XeU : VU} exhibiting the lowest activation energy.

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

Local atomic structure and lattice defect analysis in heavily Co-doped ZnS thin films using X-ray absorption fine structure spectroscopy

Iin this work, X-ray absorption fine structure (XAFS) spectroscopy was used to address both the local atomic structure and the local electronic structure around Zn, Co, and S in Co-doped ZnS thin films. X-ray absorption near-edge spectroscopy (XANES) revealed changes in the strong pre-edge feature of the S K-edge spectra that depend on the Co concentration and are related to atomic bound-state transitions through hybridization between $3d$ transition metal and $sp$ host semiconductor electronic states. These changes reveal intrinsic effects, such as S vacancies and dopant effects. From the Co K-edge XANES spectra, substitution of Co 2+ into ZnS is observed. The local atomic structure around the dopant, obtained by XAFS spectroscopy, indicates a reduction of the Co–S bond length. Optical analysis shows deep absorption centers or intermediate states in the bandgap, ascribed to $d$-$d$ transitions in the Co atoms that modify the electronic band structure and cause a decrease in the bandgap as a function of Co concentration. This work provides direct evidence of dopant effects in heavily Co-doped ZnS, including bond distances and Debye-Waller factors, providing key elements to generate a realistic atomic structural model, which is crucial to understand the observed electronic properties of this material.

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

Theoretical study of CH 3 Cl-N 2 line shapes in the ν 1 band. Line mixing effects in QR doublets and QQ k sub-branches

Line coupling and line mixing effects of CH 3 Cl lines in the ν 1 band perturbed by N 2 have been investigated. We have taken into account the non-diagonality of the exp(–S 2 ) operator within specified line spaces as well as the k-degeneracy of the transitions (due to the double degeneracy of the j,k levels with k ≠ 0). These transitions should be considered as doublets coupled by the line mixing process. A new problem appears in the calculation when the atom-atom potential model is introduced. In order to overcome this difficulty, a pragmatic approach is proposed. Comparisons of theoretically calculated matrix elements of W with measurements of QR(j,k) doublets as well as some QQ k sub-branches, which are strongly affected by line mixing, have been made. Finally, the results show that the formalism improved in this way leads to rather accurate predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Energy and mechanical properties predictions in Fe-Ni binary system by ab initio calculations

Here, the fcc random solid solutions in Fe-Ni binary system were investigated using ab initio calculations. Density functional theory (DFT) was applied along with supercells generated by the Special Quasi-Random Structure (SQS) method to predict the enthalpy of mixing and magnetic and elastic properties. The results were compared with relevant experimental and computational data produced by other groups. The effect of SQS supercell size and relaxation procedure on the predicted properties of the alloys were also studied. Our calculations show that 16-atoms SQS model provides adequate approach to the mixing energy of the alloy. However, calculation of mechanical properties requires larger supercell. Excellent agreement between experimental and calculated elastic properties was found in fcc region of Fe-Ni system.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Investigating radioactive negative ion production via double electron capture

The relative cross sections for radioactive negative ion production via double electron capture have been measured for collisions between a 40 keV projectile beam of uranium-238 and potassium vapor. This was performed at the collinear resonance ionization spectroscopy (CRIS) experiment at CERN-ISOLDE and is a step towards measuring the electron affinities (EAs) of elements that cannot be efficiently produced in negative ion sources at radioactive ion beam (RIB) facilities. This includes short-lived radioactive isotopes that have low production quantities and heavy and superheavy elements that systematically have smaller EAs than work functions of available ion source materials. Negative ions are particularly sensitive to electron–electron correlation effects, which make such studies ideal for benchmarking atomic structure models that go beyond the independent particle model. While the EAs of most light elements have been measured, experimental investigations on heavier elements, namely the actinides, remain scarce due to their radioactive nature and production difficulty. By developing negative ion production by charge exchange, we aim to make these studies feasible at RIB facilities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Validation of the single-event method for low-energy electron transport via stopping power calculations with $\mathrm{MCNP}$

Monte Carlo simulations of low-energy ( <50 keV) electron transport in matter are essential for a broad range of application fields. Several Monte Carlo codes have developed specialized treatments for this case, but a comprehensive validation of low-energy electron transport for general-purpose simulations remains lacking in the literature. One approach to accomplish this validation is calculation of stopping power using low-energy electron transport physics, as stopping power is a fundamental radiation transport quantity which must be simulated accurately for nearly any application. Here in this work, we use the Monte Carlo N-Particle (MCNP) radiation transport code with the single-event method for electron transport to calculate stopping powers of low-energy electrons (50 eV to 30 keV) in 41 elemental solids, 14 compound solids, and five rare gas solids, comparing simulation results to published semi-empirical stopping power calculations from optical measurements. In general, the simulations give good agreement (typically within ±10%) with semi-empirical stopping power values at higher energies: 300 eV and above for most elemental solids, 1 keV and above for compound solids, and 400 eV and above for rare gas solids. Agreement between MCNP and semi-empirical values is worse below these energies. The most significant source of error is the EPRDATA14 cross section data, which does not account for changes in electronic structure due to solid-state bonding, particularly in compound materials. The simplistic model of atomic excitation used to generate the EPRDATA14 cross sections is another key source of error. Additionally, the breakdown of the continuous slowing-down approximation introduces significant uncertainty at low energies, although this is a limitation of the calculation method and not of the simulation procedure. Accounting for these and other uncertainty sources, the single event method in MCNP is robust and able to give good accuracy for a variety of low-energy electron transport problems through diverse kinds of materials.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2016

Experimental and theoretical studies of laser-induced breakdown spectroscopy emission from iron oxide: Studies of atmospheric effects

In this paper, we report on a comprehensive study of the emission spectra from laser-induced breakdown spectroscopy (LIBS) measurements on iron oxide. Measurements have been made of the emission from Fe 2 O 3 under atmospheres of air, He, and Ar, and at different atmospheric pressures. The effect of varying the time delay of the measurement is also explored. Theoretical calculations were performed to analyze the plasma conditions and find that a reasonably consistent picture of the change in plasma temperature and density for different atmospheric conditions can be reached. We also investigate the sensitivity of the OI 777 nm emission lines to the plasma conditions, something that has not been explored in detail in the previous work. In conclusion, we also show that LIBS can be used to differentiate between FeO and Fe 2 O 3 by examining the ratio of the intensities of selected Fe emission to O emission lines.

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

Role of Cr-rich carbide precipitates in the intergranular oxidation of Ni-Cr alloys

The influence of grain boundary Cr carbides on the intergranular oxidation behavior of a Ni-16Cr alloy is considered using atomistic-to-mesoscale model and atom probe tomography. Cr carbide strongly perturbs the collective reactive dynamics of oxidizing species and alloy elements in the intergranular region. Strong attractive interactions between oxygen and carbide create a driving force for Cr and Ni accumulation in the grain boundary adjacent to the carbide, and for the depletion of Cr and Ni ahead of the oxidation front beyond the carbide. Furthermore, the results shed light on the mechanism of carbide-assisted protection of the alloy against grain boundary corrosion/oxidation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Advancing the Performance of Lithium-Rich Oxides in Cooperation with Inherent Complexities: Engineering Domain Structure

The nanodomain structure of lithium- and manganese-rich, composite cathode materials has been systematically altered through synthesis conditions while keeping larger-scale morphological differences to a minimum. Clear changes in electrochemical performance across the samples studied are observed, especially with respect to the anomalous impedance at low states-of-charge. Atomic-scale modeling, coupled to electrochemical measurements and physical characterization, reveals the local consequences of the high-voltage activation charge process as a function of specific domain structures. Furthermore, the results explain the influence of different postactivation structures on the insertion of Li-ions and the associated impedance during discharge. This work adds to our series of studies on lithium- and manganese-rich oxides, demonstrating that the inherent performance of LMRs can be greatly enhanced through rational, highly controlled synthetic strategies based on an understanding of synthesis–structure–property relationships across length scales.

25 ENERGY STORAGE↗
DOE OSTI2017

Stabilization of Lithium Transition Metal Silicates in the Olivine Structure

While olivine LiFePO 4 shows amongst the best electrochemical properties of Li-ion positive electrodes with respect to rate behavior owing to facile Li + migration pathways in the framework, replacing the [PO 4 ] 3- polyanion with a silicate [SO 4 ] 4- moitie in olivine is desirable. This would allow additional balancing alkali content and hence electron transfer, and increase the capacity. We demonstrate the first stabilization of a lithium transition-metal silicate (as a pure silicate) in the olivine structure type. Using LiInSiO 4 and LiScSiO 4 as the parent materials, transition metal (Mn, Fe, Co) substitutions on the In/Sc site were investigated by computational modelling via atomic scale simulation. Transition metal substitution was found to be only favourable for Co, a finding confirmed by the successful solid state synthesis of olivine LixInyCo 2-x-y SiO 4 . Finally, the stabilization of the structure was achieved by entropy provided by cation disorder.

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

A Multiple Filter Based Neural Network Approach to the Extrapolation of Adsorption Energies on Metal Surfaces for Catalysis Applications

Computational catalyst discovery involves the development of microkinetic reactor models based on estimated parameters determined from density functional theory (DFT). For complex surface chemistries, the number of reaction intermediates can be very large, and the cost of calculating the adsorption energies by DFT for all surface intermediates even for one active site model can become prohibitive. In this paper, we have identified appropriate descriptors and machine learning models that can be used to predict a significant part of these adsorption energies given data on the rest of them. Furthermore, our investigations also included the case when the species data used to train the predictive model are of different size relative to the species the model tries to predict—this is an extrapolation in the data space which is typically difficult with regular machine learning models. Due to the relative size of the available data sets, we have attempted to extrapolate from the larger species to the smaller ones in the current work. Here, we have developed a neural network based predictive model that combines an established additive atomic contribution based model with the concepts of a convolutional neural network that, when extrapolating, achieves a statistically significant improvement over the previous models.

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

Bond Dissociation Energies of the Actinide Halides AnX, An = Ac–Lr and X = F–I, Utilizing Relativistic Composite Coupled Cluster Approaches

Bond dissociation energies (BDEs) have been calculated for the set of actinide halides AnX with An=Ac, Pa, and Np-Lr and X=F-I. Two composite thermochemistry methods based on the Feller-Peterson-Dixon (FPD) approach have been utilized, one involving spinor-based relativistic CCSD(T) calculations where spin-orbit (SO) was included at the orbital level and another using scalar relativistic CCSD(T) with a posteriori SO contributions based on 2-component multireference configuration interaction calculations. The method that was chosen for a given actinide halide was based on which representation yielded the best single determinant reference determinant for the coupled cluster calculation. The spinor-based method was chosen for all cases except for AmX, CmX, and BkX. Both composite approaches included contributions accounting for basis set truncation, outer-core correlation, the Gaunt interaction, and QED. The scalar FPD results, as well as the spinor-based calculations for AcF, also included higher order electron correlation up through CCSDT(Q). In addition to BDEs, CCSD(T) equilibrium bond lengths, harmonic frequencies, and vibrational anharmonicity constants are reported for all species. Last, the FPD BDEs for the fluorides were used to confirm the trend across the actinide series previously predicted by Gibson using bonding models based atomic promotion energies that provide a single 6d electron for bonding. In particular the local minimum in the BDEs at AmF is confirmed in the present calculations. Furthermore, the BDEs for LrX are predicted to be slightly larger than those of AcX, making them the largest in the actinide halide series.

Actinides↗
DOE OSTI2025

First-Principles Insights into the Thermocatalytic Cracking of Ammonia-Hydrogen Blends on Fe(110). 2. Kinetics

Ammonia (NH 3 ) is an energy-rich molecule that is routinely synthesized from nitrogen (N 2 ) and hydrogen (H 2 ). NH 3 ’s more favorable physical properties compared to H 2 suggests it may offer a way to more conveniently store, transport, and, when needed, extract H 2 via thermal decomposition. However, the high kinetic barrier and endoergicity to decompose to H 2 and N 2 require high temperatures. The standard reaction free energy indicates nearly 100% thermodynamic conversion to the diatomic molecules only at ~673 K and higher. However, even at these temperatures, a catalyst, e.g., iron (Fe), is needed for favorable kinetic conversion. Here, in this study, we explore via density functional theory the kinetics of NH 3 decomposition on the most stable facet of body-centered cubic Fe, namely, (110), under typical high-temperature and finite-pressure operando conditions. We predict coverage-dependent energetics of elementary surface reactions, often neglected in atomic-scale modeling. From these models, we find the recombinative desorption of adsorbed N as N 2 is rate-determining at 573.15–773.15 K and even at an extreme case of 1173.15 K. From microkinetic modeling, we find that the steady-state turnover frequencies (TOFs) for N 2 and H 2 generation rates (r$_{H_2}$) depend exponentially on temperature. The catalyst achieves a steady-state TOF of 36.4 s –1 and an r$_{H_2}$ of 0.107 μmol cm –2 s –1 for a feed of 1.8 bar NH 3 with 0.2 bar H 2 at 1173.15 K. However, at 773.15 K, with the same feed composition and velocity, the steady-state TOF and r$_{H_2}$ decrease to 0.14 s –1 and 4.10 × 10 –4 μmol cm –2 s –1 , respectively, as the process is significantly hindered by slow N 2 desorption. Although at first glance counterintuitive, our simulations suggest that surface modifications that reduce Fe’s reactivity toward NH x species should enhance its overall NH 3 decomposition activity.

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

Segregation and Migration of the Oxygen Vacancies in the 3 (111) Tilt Grain Boundaries of Ceria

In nanocrystalline materials, defect-grain boundary (GB) interaction plays a key role in determining the structure stability, as well as size-dependent ionic, electronic, magnetic and chemical properties. In this study, we systematically investigated using density functional theory segregation and migration of oxygen vacancies at the Σ3 [110] / (111) grain boundary of ceria. Three oxygen layers near the GB are predicted to be segregation sites for oxygen vacancies. Moreover, the presence of oxygen vacancies stabilizes this tilt GB at a low Fermi level and/or oxygen poor conditions. An atomic strain model was proposed to rationalize layer dependency of the relaxation energy for +2 charged oxygen vacancy. The structural origin of large relaxation energies at layers 1 and 2 was determined to be free-volume space that induces ion relaxation towards the GB. Our results not only pave the way for improving the oxygen transport near GBs of ceria, but also provide important insights into engineering the GB structure for better ionic, magnetic and chemical properties of nanocrystalline ceria.

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

Atomic-Scale Surface Structure and CO Tolerance of NiCu Single-Atom Alloys

Ni is one of the most extensively utilized metals in industrial catalysis. For example, Ni is the catalyst of choice for the steam reforming of hydrocarbons. However, pure Ni also detrimentally catalyzes the formation of graphitic carbon, which in turn leads to coking and deactivation of the catalyst. It has been shown that alloying small amounts of a less reactive metal like Au into Ni can alleviate this issue by breaking up the larger Ni ensembles that promote coke formation. We are taking the opposite of this approach by alloying very small amounts of Ni into Cu, another catalytically less active host metal to create single Ni atom sites. In this way our single-atom alloy approach has the potential to greatly enhance catalytic selectivity and reduce poisoning, analogous to other single-atom alloys such as PtCu and PdCu. Herein we report the atomic-scale surface structure and local geometry of low coverages of Ni deposited on a Cu(111) single crystal with scanning tunneling microscopy. At 433 K, low concentrations of Ni alloy in the Cu host as a single-atom alloy in Ni-rich brims along ascending step edges. To support our STM assignments of the single atom dispersion of Ni, reflection absorption infrared spectroscopy of CO on NiCu was performed. To access the binding strength of CO to isolated Ni sites, we use temperature-programmed desorption studies which reveal that CO binds more weakly to single Ni atoms in Cu compared to Ni(111), indicating that NiCu single-atom alloys are promising for catalytic applications in which CO poisoning is an issue. Altogether, these results provide a guide for the preparation of NiCu single-atom alloy model catalysts that are predicted by theory to be promising for a number of reactions.

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

Insight into induced charges at metal surfaces and biointerfaces using a polarizable Lennard–Jones potential

Metallic nanostructures have become popular for applications in therapeutics, catalysts, imaging, and gene delivery. Molecular dynamics simulations are gaining influence to predict nanostructure assembly and performance; however, instantaneous polarization effects due to induced charges in the free electron gas are not routinely included. Here we present a simple, compatible, and accurate polarizable potential for gold that consists of a Lennard–Jones potential and a harmonically coupled core-shell charge pair for every metal atom. The model reproduces the classical image potential of adsorbed ions as well as surface, bulk, and aqueous interfacial properties in excellent agreement with experiment. Induced charges affect the adsorption of ions onto gold surfaces in the gas phase at a strength similar to chemical bonds while ions and charged peptides in solution are influenced at a strength similar to intermolecular bonds. The proposed model can be applied to complex gold interfaces, electrode processes, and extended to other metals.

77 NANOSCIENCE AND NANOTECHNOLOGY↗