Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional calculations”

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

Influence of Hydrogen Isotopes on Vacancy Formation and Antisite Defect Diffusion in Palladium and Vanadium Metals

Density functional theory calculations are performed to study the interaction between hydrogen isotopes, vacancy, and antisite defects in Pd and V. Various defect configurations and defect-defect distances are systematically explored. Binding energies and migration energy barriers are calculated and tabulated. The study provides atomistic data for subsequent mesoscale simulations of tritium, vacancy, and antisite diffusion. In Pd, a strong binding is found between a vacancy and tritium (0.16 eV). This tritium binding increases the thermal concentration of vacancies by a factor of ~10 at 500 ?C. The increase in vacancy concentration enhances V antisite diffusion in Pd by a factor of ~6 at 500 ?C. The influence of tritium is even stronger in V, with a tritium and vacancy binding energy of 0.38 eV. Such a strong binding increases the thermal concentration of vacancies by a factor of ~300 at 500 ?C. The increase in vacancy concentration enhances Pd antisite diffusion in V by a factor of ~640 at 500 ?C. Vanadium and Pd exhibit a strong driving force to intermix with a formation energy of -1.57 eV for V antisite in Pd and -1.05 eV for Pd antisite in V. The results suggest vanadium diffusion into Pd is energetically stronger than the reverse. Zero-point-energy corrections are taken into account and calculations for hydrogen and tritium are presented.

Setyawan, Wahyu↗
DOE OSTI2023

Exploring structural transitions at grain boundaries in Nb using a generalized embedded atom interatomic potential

The advancement in experimental techniques, like the atom probe tomography and high resolution electron microscopy, is fueling interest in studying structural transformations of grain boundaries in metal and alloys to uncover correlations between mechanical properties and solute or impurity segregation to grain boundaries. Atomistic modeling is an important tool that can pinpoint the intricate dynamics of grain boundary phase transitions, but the lack of accurate interatomic potentials needed to simulate the complex dynamics of grain boundary structural transitions and identify different metastable phases has been the bottleneck. To this end, we use niobium as a model body centered cubic (BCC) metal and develop an interatomic potential to study grain boundary phase transitions. The potential for Nb is based on a generalization of the embedded atomic method potential and has sufficient flexibility to learn complex energy landscapes using a small set of training structures. We systematically test and validate the using data from ab initio density functional theory calculations and experiments. Using this potential, we calculate energies of multiple symmetric-tilt grain boundaries spanning a wide range of misorientation angles. Additionally, we explore different metastable structures of the Σ 27(552) [$1\overline{1}0$] grain boundary and use molecular dynamic simulations to study the coexistence of metastable phases and grain boundary transitions at finite temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Machine learning in materials research: Developments over the last decade and challenges for the future

The number of studies that apply machine learning (ML) to materials science has been growing at a rate of approximately 1.67 times per year over the past decade. In this review, I examine this growth in various contexts. First, I present an analysis of the most commonly used tools (software, databases, materials science methods, and ML methods) used within papers that apply ML to materials science. The analysis demonstrates that despite the growth of deep learning techniques, the use of classical machine learning is still dominant as a whole. It also demonstrates how new research can effectively build upon past research, particular in the domain of ML models trained on density functional theory calculation data. Next, I present the progression of best scores as a function of time on the matbench materials science benchmark for formation enthalpy prediction. In particular, a dramatic improvement of 7 times reduction in error is obtained when progressing from feature-based methods that use conventional ML (random forest, support vector regression, etc.) to the use of graph neural network techniques. Finally, I provide views on future challenges and opportunities, focusing on data size and complexity, extrapolation, interpretation, access, and relevance.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The combined effects of Mg 2+ and Sr 2+ incorporation during CaCO 3 precipitation and crystal growth

Calcite (CaCO 3 ) composition and properties are defined by the chemical environment in which CaCO 3 forms. However, a complete understanding of the relationship between aqueous chemistry during calcite precipitation and resulting chemical and physical CaCO 3 properties remains elusive; therefore, we present an investigation into the coupled effects of divalent cations Sr 2+ and Mg 2+ on CaCO 3 precipitation and subsequent crystal growth. Through chemical analysis of the aqueous phases and microscopy of the resulting calcite phases in compliment with density functional theory calculations, we elucidate the relationship between crystal growth and the resulting composition (elemental and isotopic) of calcite. The results of this experimental and modeling work suggest that Mg 2+ and Sr 2+ have cation-specific impacts that inhibit calcite crystal growth, including: (1) Sr 2+ incorporates more readily into calcite than Mg 2+ (D Sr > D Mg ), and increasing [Sr 2+ ] t or [Mg 2+ ] t increases D Sr ; (2) the inclusion of Mg 2+ into structure leads to a reduction in the calcite unit cell volume, whereas Sr 2+ leads to an expansion; (3) the inclusion of both Mg 2+ and Sr 2+ results in a distribution of unit cell impacts based on the relative positions of the Sr 2+ and Mg 2+ in the lattice. These experiments were conducted at saturation indices of CaCO 3 of ~4.1, favoring rapid precipitation. This rapid precipitation resulted in observed Sr isotope fractionation confirming Sr isotopic fractionation is dependent upon the precipitation rate. We further note that the precipitation and growth of calcite favors the incorporation of the lighter 86 Sr isotope over the heavier 87 Sr isotope, regardless of the initial solution conditions, and the degree of fractionation increases with D Sr . In sum, these results demonstrate the impact of solution environment to influence the incorporation behavior and crystal growth behavior of calcite. These factors are important to understand in order to effectively use geochemical signatures resulting from calcite precipitation or dissolution to gain specific information.

58 GEOSCIENCES↗
DOE OSTI2022

Pressure dependent magnetic properties on bulk CrBr 3 single crystals

The van der Waals class of materials offer an approach to two-dimensional magnetism as their spin fluctuations can be tuned upon exfoliation of layers. Moreover, it has recently been shown that spin-lattice coupling and long-range magnetic ordering can be modified with pressure in van der Waals materials. In this work, the magnetic properties of quasi two-dimensional CrBr 3 are reported applying hydrostatic pressure. The application of pressure up to 0.844 GPa shows a 1.77% decrease in saturation magnetization with a decrease in the Curie temperature from 33.05 to 30.41 K. Density functional theory calculations with pressure up to 1 GPa show a reduction in volume and interplanar distance as pressure increases. To further understand the magnetic properties with applied pressure, the magnetocrystalline anisotropy energy (MAE) and exchange coupling parameter (J) are calculated. There is small decrease in MAE and the first nearest neighbor interaction (J 1 ) (U = 2.7 eV and J = 0.7 eV) is increasing with respect to increasing pressure. Overall, CrBr 3 displays ferromagnetic interlayer coupling and the calculated exchange coupling and MAE parameters match well with the observations from the experimental work.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Magnetic ordering in Eu 2 In and Eu 2 Sn

Eu 2 In and Eu 2 Sn crystallize in the orthorhombic Co 2 Si-type structure (oP12, Pnma, No. 62) with In and Sn atoms occupying one 4c site and the Eu atoms filling two other 4c sites. Eu 2 In has a nearly ideal first-order magnetostructural transition (FOMT) at 55 K with a hysteresis of less than 0.1 K, a large entropy change and an adiabatic temperature change of 5.0 K in a field of 2 T. The anhysteretic nature of the FOMT is likely due to there being no change in cell symmetry and relatively small changes in the lattice parameters. There is no magnetostructural transition in Eu2Sn. In this work we present the results of powder neutron diffraction, magnetization, and Eu Mössbauer spectroscopy aimed to investigate the nature of magnetic order for both Eu 2 In and Eu 2 Sn. The Eu Mössbauer spectrum of Eu 2 In at 5 K shows two equal area components, consistent with Eu occupying two equal multiplicity crystallographic sites. However, the different hyperfine fields (B) of 27 T and 17 T suggest that the magnetic environments of the Eu moments on the two 4c sites are different. Neutron diffraction data at 2.5 K show that in Eu 2 In the order is ferromagnetic, with Eu moments on both Eu sites oriented parallel to the a-axis; moment values of 6.8 μB and 6.5 μB were found. For Eu 2 Sn measurements find two antiferromagnetic transitions, which are corroborated by neutron diffraction. Analysis of density-functional theory calculations shows negligible energy difference between differing magnetic configurations, indirectly supporting stability of multiple magnetic structures observed experimentally. While the transition at T N1 = 30 K corresponds to the formation of a simple k 1 = 0 antiferromagnetic structure with Eu-moments pointing along the b-axis, at T N2 = 13 K a coexisting second magnetic order with k 2 = [0, ½, ½] appears.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Metal-organic framework supported single-site nickel catalysts for butene dimerization

Homotopic sites in a well-controlled environment are not only ideal systems for mechanistic studies, but also allow optimal control of catalytic transformations. Sites having only a single metal cation and sites consisting of metal oxo complexes with few nickel (Ni) cations supported on the nodes of UiO-66 metal-organic framework (Ni-UiO-66) are studied for 1-butene dimerization. Monomeric Ni sites, which bind to the Zr 6 node via two Zr-OH(µ3) linkages, are active and selective for the dimerization of 1-butene to linear and mono-branched C 8 isomers. Ni oxo complexes with few Ni cations show lower activity and promote the oligomerization of transiently formed C 8 isomers. In conclusion, Kohn-Sham density function theory calculations combined with spectroscopic measurements and kinetic analyses indicate that dimerization follows a Cossee-Arlman reaction mechanism.

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

Hydrogen activation by rhodium under the cover of a copper oxide thin film

Here, the activation of reactants by catalytically active metal sites at metal-oxide interfaces is important for understanding the effect of metal-support interactions on nanoparticle catalysts and for tuning activity and selectivity. Using a combined experimental and theoretical approach, we studied the activation of H 2 and the effect of CO poisoning on isolated Rh atoms completely or partially covered by a copper oxide (Cu 2 O) thin film. Temperature-programmed desorption (TPD) experiments conducted in ultra-high vacuum (UHV) show that neither a partially nor a fully oxidized Cu 2 O layer grown on a Rh/Cu(111) single-atom alloy can activate hydrogen in UHV. However, in situ ambient pressure X-ray photoelectron spectroscopy (AP-XPS) experiments performed at elevated H 2 pressures reveal that Rh significantly accelerates the reduction of these Cu 2 O thin films by hydrogen. Remarkably, the fastest reduction rate is observed for the fully oxidized sample with all Rh sites covered by Cu 2 O. Both TPD and AP-XPS data demonstrate that these covered Rh sites are inaccessible to CO, indicating that Rh under Cu 2 O is active for H 2 dissociation but cannot be poisoned by CO. In contrast, an incomplete oxide film leaves some of the Rh sites exposed and accessible to CO, and hence prone to CO poisoning. Density functional theory calculations demonstrate that unlike many reactions in which hydrogen activation is rate limiting, the rate-determining step in the dissociation of H 2 on thin-film Cu 2 O with Rh underneath is the adsorption of H 2 on the buried Rh site, and once adsorbed, the dissociation of H 2 is barrierless. These calculations also explain why H 2 can only be activated at higher pressures. Together, these results highlight how different the reactivity of atomically dispersed Rh in Cu can be depending on its accessibility through the oxide layer, providing a way to engineer Rh sites that are active for hydrogen activation but resilient to CO poisoning.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Cross-aldol condensation on missing linker sites of metal–organic framework UiO-66

Missing linkers induce coordinatively unsaturated Lewis acid-base pair sites in metal–organic frameworks such as UiO-66 (Zr). The Lewis acid-base pair is induced by replacing the benzene dicarboxylic acid linker with formic acid, which subsequently is thermally decomposed. The concentration of coordinatively unsaturated Zr 4+ cations was titrated with acetone. In parallel to their concentration, the catalytic activity increased while the selectivity for aldol condensation between furfural and acetone remained constant. The missing linkers generate a complex ensemble consisting of an (at least partly) accessible Zr 4+ cation surrounded by oxygen, catalyzing the cross-aldol condensation. As a result, density functional theory calculations, matching the experimental findings provide a more detailed understanding of the elementary steps during the condensation reaction.

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

Thermodynamic evaluation and modeling of the CeF 3 -LaF 3 -LiF-NaF system

The thermodynamic properties of the LiF-NaF-CeF 3 -LaF 3 pseudo-quaternary system have been evaluated, and models for the relevant phases have been revised or newly developed. The solid solubility of NaF in LaF 3 and CeF 3 was described in the model based on reported phase-equilibria observations and results from density functional theory calculations. The Gibbs energy functions for the NaCeF 4 and NaLaF 4 phases have been updated with recently reported thermochemical measurements. The Gibbs energy functions and model relations for the LiF-NaF-LaF 3 and LiF-NaF-CeF 3 pseudo-ternary systems were successfully optimized based on thermochemical properties and phase equilibria observations.

CALPHAD↗
DOE OSTI2023

Synergistic double-atom catalysts of metal-boron anchored on g-C 2 N for electrochemical nitrogen reduction: Mechanistic insight and catalyst screening

The rational design of a novel catalytic center with a sound basis remains both challenging and rewarding for the electrochemical reduction of N 2 (eNRR), which has provided a feasible route for achieving clean and sustainable NH 3 production under ambient conditions. Herein, using density functional theory calculations, we demonstrate that hybrid metal (M)-boron (B) double-atom catalysts (DACs) embedded in g-C 2 N substrate (M-B@C 2 N, M = 3d, 4d and 5d transition metals) can achieve both high catalytic activity and high selectivity in eNRR. The proposed M-B@C 2 N DACs have exhibited impressive feasibility and stability thanks to the resilient and robust C 2 N substrate with abundant pyridinic N atoms distributed among right-sized pore structures. Our results reveal that like the metal center, the embedded B atom can actively involve in Ntriple bondN bond activation via π*-backdonation mechanism concomitant with the substantial charge transfer to adsorbed *N 2 , leading to sizable Nsingle bondN bond elongation. Accordingly, both adsorption energy and Nsingle bondN bond length of *N 2 can be employed as catalytic descriptors for predicting eNRR activity in terms of the limiting potentials (UL). Using high-throughput screening method, we found that six M-B@C 2 N candidates have stood out as the outstanding electrocatalysts for driving eNRR, namely, M = Ti (U L = 0 V), Mo (U L = 0 V), Nb (U L = -0.04 V), W (U L = -0.23 V), Zr (U L = -0.26 V), V (U L = -0.28 V). The underlying origin is attributed to the balanced and constrained N-affinity of M-B dual site working in synergy, which can thus be used as one important guide of catalyst design.

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

In situ preparation of visible-light-driven carbon quantum dots/NaBiO 3 hybrid materials for the photoreduction of Cr(VI)

Here, different carbon quantum dots (CQDs)/NaBiO 3 hybrid materials were synthesized as photocatalysts to effectively utilize visible light for the photocatalytic degradation of contaminants effectively. These hybrid materials exhibit an enhanced photocatalytic reduction of hexavalent chromium (Cr(VI)) in the aqueous medium. Zero-dimensional nanoparticles of CQDs were embedded within the two-dimensional NaBiO 3 nanosheets by the hydrothermal process. Compared with that of the pure NaBiO 3 nanosheets, the photocatalytic performance of the hybrid catalysts was significantly high and 6 wt.% CQDs/NaBiO 3 catalyst exhibited better photocatalytic performance. We performed the first-principles density functional theory calculations to study the interfacial properties of pure NaBiO 3 nanosheets and hybrid photocatalysts, and confirmed the CQDs played an important role in the CQDs/NaBiO 3 composites. The experimental results indicated that the enhanced reduction of Cr(VI) was probably due to the high loading of CQDs (electron acceptor) on NaBiO 3 , which made NaBiO 3 nanomaterials to respond in visible light and significantly improved their electron-hole separation efficiency.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The correlation between N deficiency and the mechanical properties of the Ti 2 AlNy MAX phase

The role of X deficiency on the mechanical properties of MAX phases was studied by synthesizing Ti 2 AlN through powder metallurgy in stoichiometric and sub/extra-stoichiometric nitrogen compositions. XRD analyses and ab initio calculations indicate that nitrogen vacancies result in a lattice contraction predominantly along the c-axis. The elastic moduli and intrinsic hardness of substoichiometric Ti 2 AlN 0.9 measured from nanoindentation tests are shown to be slightly smaller than that of Ti 2 AlN. The key mechanical indexes, bulk (B), shear (G) and Young’s (E) moduli as well as the hardness variation are calculated in density functional theory, and show different responses depending on the concentration of N vacancies. Overall, this joint experimental and theoretical study provides a full understanding of the energetics, chemical bonding, electronic structure, and mechanics of the N deficient MAX phases which would increase the application of nitride ceramics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Impact of fission product inclusion on phase development in U 3 Si 2 fuel

Due to its high thermal conductivity and uranium density, U 3 S i2 has been considered as a candidate for use as an accident tolerant fuel (ATF). In order to fully assess its suitability and performance as a fuel, the impact of fission products (FPs) on the stability and performance of U 3 Si 2 must be investigated. The interactions of FPs and U 3 Si 2 have had relatively little study until now and require experimental and computational examination. U 3 Si 2 was doped with individual FPs to explore U-Si-FP interactions and phase equilibria that may impact the performance of the ATF during irradiation. Elemental Ce, Mo, Y, or Zr were used to individually dope U 3 Si 2 at a concentration of 5 wt% FP. A diffusion couple of a 1:1 Mo:Zr alloy and U 3 Si 2 was heated to 1200 °C in order to consider the impacts of multiple FPs on the stability and structure of the fuel. Samples were characterized for FP solubility and secondary phase formation using electron microscopy, energy dispersive spectroscopy, and x-ray diffraction. First principles density functional theory calculations complemented the experimental effort to understand FP behavior. Experimental and computational findings were used in the development of a thermodynamic database containing 8 major FPs and their associated silicide phases. Finally, fuel compositions generated from depletion calculations were used to thermodynamically model the equilibrium phases of the fuel undergoing burnup.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

First-principles study of the surface properties of uranium carbides

Uranium carbides have attracted renewed interest as advanced nuclear fuels for Generation IV reactors. As an important property required for gas bubble modeling in nuclear fuels, the surface energy of uranium carbides is scarce in literature. In this work, we study the surface properties of uranium carbides by first-principles density functional theory calculations. Surface orientations with maximum Miller index up to 3, 2 and 2 are investigated for UC, U 2 C 3 and α-UC 2 , respectively. By studying the effects of surface termination and chemical potential on surface energy, we identify the factors that determines the surface stability. From the calculated surface energies, the surface properties of uranium carbide single crystals are obtained from Wulff construction, including equilibrium morphology, dominant surface orientation and area weighted surface energy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Insights into the thermal decomposition of plutonium(IV) oxalate – a DFT study of the intermediate structures

The thermal decomposition of plutonium oxalate to oxide is one of the most studied reactions in actinide chemistry but the intermediates have been the subject of debate for decades. Recent experimental data suggest that the decomposition of Pu(IV) oxalate in air undergoes dehydration first, then reduction to Pu(III) oxalate. The precise structural modifications that take place are unknown as experiments have not been able to fully characterize the intermediates at the microscopic level. To rectify this, we employed solid state density functional theory calculations at the PBE-D3 level with a Hubbard U correction to model the structures and energetics of potential dehydrated Pu(IV) and Pu(III) oxalate intermediate compounds. Based on the theoretical study presented here, the anhydrous analogues of the known hydrated Pu(IV) and Pu(III) oxalates are the preferred crystal structures formed through an overall exothermic reaction process. However, decomposition could proceed through the formation of a higher energy, more complicated 3D lattice structure with frustrated oxalate binding. It is expected that the intermediates presented here could be identified using spectroscopic techniques to enable further insight into the reaction mechanism.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Helium interaction with solutes and impurities in neutron-irradiated nanostructured ferritic alloys: A first principles study

Density functional theory calculations are performed to explore the binding between He and alloying solutes, impurities, and transmutation products expected in neutron irradiated nanostructured ferritic alloys (NFAs), here 14YWT is taken as an example. Elements that exhibit significant binding (attraction) with an interstitial He are Y (binding energy = 0.46 eV), Mg (0.32), O (0.33), Ti (0.16), and C (0.15). Those that provide significant binding to a substitutional He are O (1.44), Y (1.24), N (0.73), H (0.56), Mg (0.52), Ti (0.34), Si (0.34), C (0.33), Al (0.32), Ni (0.26), Ta (0.23), and Mn (0.16). The presence of these elements in Fe matrix could reduce the transport of He towards oxide particles, dislocations, and internal boundaries, and could promote He bubble nucleation in the matrix. For convenience, we compile existing binding energy data of He with He n and He n V (He-vacancy) clusters. Dissociation pathway analysis reveals that, in general, the most likely dissociation of a He n V cluster is by a sequential emission of individual He atoms. Furthermore, larger bubbles are more prone to dissociation than smaller ones. In addition, we estimate the binding energy (segregation energy) of He in bulk Y 2 Ti 2 O 7 (YTO) single crystal, YTO/Fe interface, and YTO particle embedded in Fe, with respect to interstitial He in Fe, from existing formation energies of He in these structures. We also compile available data of He binding with Fe self-interstitial atom (SIA), SIA clusters, and edge and screw dislocations. Note that given the absence of DFT data, the binding with SIA clusters and dislocations are gathered from simulations with empirical potentials. Finally, the data presented in this paper is important to inform multiscale simulations of He bubble accumulation.

36 MATERIALS SCIENCE↗