Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT 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 1,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Electrochemical CO2 Reduction over Metal-/Nitrogen-Doped Graphene Single-Atom Catalysts Modeled Using the Grand-Canonical Density Functional Theory

Renewably driven, electrochemical conversion of carbon dioxide into value-added products is expected to be a critical tool in global decarbonization. However, theoretical studies based on the computational hydrogen electrode largely ignore the nonlinear effects of the applied potential on the calculated results, leading to inaccurate predictions of catalytic behavior or mechanistic pathways. Here, we use grand canonical density functional theory (GC-DFT) to model electrochemical CO2 reduction (CO2R) over metal- and nitrogen-doped graphene catalysts (MNCs) and explicitly include the effects of the applied potential. We used GC-DFT to compute the CO2 to CO reaction intermediate energies at -0.3, -0.7, and -1.2 VSHE catalyzed by MNCs each doped with 1 of the 10 3d block metals coordinated by four pyridinic nitrogen atoms. Our results predict that Sc-, Ti-, Co-, Cu-, and Zn-N4Cs effectively catalyze CO2R at moderate to large reducing potentials (-0.7 to -1.2 VSHE). ZnN4C is a particularly promising electrocatalyst for CO2R to CO both at low and moderate applied potentials based on our thermodynamic analysis. Our findings also explain the observed pH independence of CO production over FeN4C and predict that the rate-determining step of CO2R over FeN4C is not *CO2- formation but rather *CO desorption. Additionally, the GC-DFT-computed density of states analysis illustrates how the electronic states of MNCs and adsorbates change non-uniformly with applied potential, resulting in a significantly increased *CO2- stability relative to other intermediates and demonstrating that the formation of the adsorbed *CO2- anion is critical to CO2R activation. This work demonstrates how GC-DFT paves the way for physically realistic and accurate theoretical simulations of reacting electrochemical systems.

CO2 reduction↗
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 OSTI2020

σ-Aromaticity-Induced Stabilization of Heterometallic Supertetrahedral Clusters [Zn 6 Ge 16 ] 4– and [Cd 6 Ge 16 ] 4–

Here, the largest heterometallic supertetrahedral clusters, [Zn 6 Ge 16 ] 4– and [Cd 6 Ge 16 ] 4– , were directly self-assembled through highly-charged [Ge 4 ] 4– units and transition metal cations, in which 3-center–2-electron σ bonding in Ge 2 Zn or Ge 2 Cd triangles plays a vital role in the stabilization of the whole structure. The cluster structures have an open framework with a large central cavity of diameter 4.6 Å for Zn and 5.0 Å for Cd, respectively. Time-dependent HRESI-MS spectra show that the larger clusters grow from smaller components with a single [Ge 4 ] 4– and ZnMes 2 units. Calculations performed at the DFT level indicate a very large HOMO–LUMO energy gap in [M 6 Ge 16 ] 4– (2.22 eV), suggesting high kinetic stability that may offer opportunities in materials science. These observations offer a new strategy for the assembly of heterometallic clusters with high symmetry.

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

Cylindrical C 96 Fullertubes: A Highly Active Metal‐Free O 2 ‐Reduction Electrocatalyst

Abstract A new isolation protocol was recently reported for highly purified metallic Fullertubes D 5h ‐C 90 , D 3d ‐C 96 , and D 5d ‐C 100, which exhibit unique electronic features. Here, we report the oxygen reduction electrocatalytic behavior of C 60 , C 70 (spheroidal fullerenes), and C 90 , C 96 , and C 100 (tubular fullerenes) using a combination of experimental and theoretical approaches. C 96 (a metal‐free catalyst) displayed remarkable oxygen reduction reaction (ORR) activity, with an onset potential of 0.85 V and a halfway potential of 0.75 V, which are close to the state‐of‐the‐art Pt/C benchmark catalyst values. We achieved an excellent power density of 0.75 W cm −2 using C 96 as a modified cathode in a proton‐exchange membrane fuel cell, comparable to other recently reported efficient metal‐free catalysts. Combined band structure (experimentally calculated) and free‐energy (DFT) investigations show that both favorable energy‐level alignment active catalytic sites on the carbon cage are responsible for the superior activity of C 96 .

Sanad, Mohamed Fathi↗
DOE OSTI2020

σ‐Aromaticity‐Induced Stabilization of Heterometallic Supertetrahedral Clusters [Zn 6 Ge 16 ] 4− and [Cd 6 Ge 16 ] 4−

Abstract In this work, the largest heterometallic supertetrahedral clusters, [Zn 6 Ge 16 ] 4− and [Cd 6 Ge 16 ] 4− , were directly self‐assembled through highly‐charged [Ge 4 ] 4− units and transition metal cations, in which 3‐center–2‐electron σ bonding in Ge 2 Zn or Ge 2 Cd triangles plays a vital role in the stabilization of the whole structure. The cluster structures have an open framework with a large central cavity of diameter 4.6 Å for Zn and 5.0 Å for Cd, respectively. Time‐dependent HRESI‐MS spectra show that the larger clusters grow from smaller components with a single [Ge 4 ] 4− and ZnMes 2 units. Calculations performed at the DFT level indicate a very large HOMO–LUMO energy gap in [M 6 Ge 16 ] 4− (2.22 eV), suggesting high kinetic stability that may offer opportunities in materials science. These observations offer a new strategy for the assembly of heterometallic clusters with high symmetry.

Xu, Hong‐Lei↗
DOE OSTI2022

Unusual Changes of C−H Bond Lengths in Chiral Zinc Complexes Induced by Noncovalent Interactions

In order to better understand the effect of non-covalent weak interactions on molecules, we have explored a variety of weak interactions, such as improper H-bonding (HB), tetrel bonds (TBs) and halogen bonds, in fluorinated chiral zinc complexes. High resolution neutron diffraction studies revealed a methylene carbon-hydrogen bond elongation and shortening due to TB and improper HB interactions, respectively. To show the accumulative effects of multiple weak interactions on the C-H bond, three types of tetrel bonds have been carefully examined. We have also shown how C-H bond elongation can be easily offset by forming an improper HB with the H atom from this C-H bond. Non-covalent interaction and electrostatic potential analysis investigations have been used to affirm the nature of the interactions based on density functional theory (DFT) and other related calculations.

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

Strong, Yet Split Hydrogen Bonding with Ice Rules in Delafossite (H/D)RhO 2

Despite remaining enigmatic, strong hydrogen bonding provides an advanced design handle for tailoring the properties of functional materials. Here, in this study, 3 R –(H/D)RhO 2 delafossites (prepared by ion exchange of Na + from NaRhO 2 ) contain H/D in linear coordination with O, linking Rh III O 2 layers. Bragg and real-space X-ray and neutron scattering analysis, vibrational and solid-state NMR spectroscopy, and density functional theory (DFT)–based electronic structure calculations have been employed to understand the nature of the hydrogen bonding. Despite short distances between H/D and the two O to which they are bonded, a clear double-minimum corresponding to a shorter and longer (H/D)–O distance is established. The triangular lattices formed by H/D appear to display ice-like disorder, corroborated by low-temperature heat capacity measurements.

NMR spectroscopy↗
DOE OSTI2022

Atomistic and cluster dynamics modeling of fission gas (Xe) diffusivity in TRISO fuel kernels

TRISO fuel particles are candidates for use in next generation reactors including gas reactors, fluoride salt-cooled high temperature reactors, and micro-reactors. The UCO fuel kernel consists of a uranium dioxide (UO) and uranium carbide mixture. The addition of UC helps suppress the formation of carbon monoxide gas, which led to failures during initial TRISO development. The addition of uranium carbide alters the chemistry of the UO kernel, which is known to influence performance parameters such as fission gas diffusivity, although the impact has not been quantified and no models exist that take the change in chemistry into account. Therefore, better understanding and more accurate models of the impact of chemistry on fuel performance are of high priority. In this paper, a first-principles density functional theory (DFT) and empirical potential based multi-scale study has been carried out to model the diffusivity of fission gas xenon (Xe) in UCO TRISO fuel kernels. The focus is on the UO component in the UCO fuel kernels, as that represents the largest volume fraction of the fuel kernels. The study relies on DFT and empirical potential calculations to determine Xe and point defect properties, which are then used in thermodynamic and kinetic models to predict diffusion for intrinsic conditions. In addition, the information is utilized in cluster dynamics simulations using the Centipede code to estimate the impact of irradiation on defect transport. Additionally, the presence of UC or UC in the UCO fuel kernels is shown to have a substantial impact on the UO non-stoichiometry by inducing oxygen vacancies and driving UO sub-stoichiometric, which causes much slower Xe diffusion in UCO compared to light water reactor UO fuel. The application of this model in fuel performance simulations using the Bison code is also demonstrated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Structural, magnetic, optical, dielectric and electronic properties of R 2 NiIrO 6 (R = Pr and Nd): A comprehensive experimental and theoretical investigation

Double perovskites are highly promising materials capable of exhibiting a wide variety of phenomena. In this work, we perform a comprehensive experimental and theoretical study of polycrystalline R 2 NiIrO 6 (R = Pr and Nd) compounds. Both compounds were synthesised using the solid-state reaction method. Rietveld refinement confirmed a monoclinic structure with the P2 1 /n space group for both compounds. The scanning electron images showed the average grain sizes of 0.55 μm for R = Pr and 0.46 μm for R = Nd. Fourier transform infrared ra- diation spectra of the two compounds presented two intense bands at 470 cm -1 and 540 cm -1 . The optical measurements revealed that the band gaps of the compounds were in the visible absorption range. The field- cooled magnetisation - field hysteresis measurements indicated exchange bias properties in the synthesised compounds at low temperatures. Both temperature and frequency variation of dielectric constant and loss tangent measurements were conducted. The frequency-dependent ac conductivity measurements indicated that the conductivity increases with the increase of frequency as well as temperature. The Nd 2 NiIrO 6 compound showed lower ac conductivities compared to its isostructural Pr 2 NiIrO 6 compound. The atomic and electronic structures of Nd 2 NiIrO 6 and Pr 2 NiIrO 6 were explored using the spin-polarised calculations performed within the DFT+U method. Our results suggested that the inclusion of on-site correlations and repulsions for the d-states of atoms was necessary in order to obtain finite band gaps of Nd 2 NiIrO 6 and Pr 2 NiIrO 6 systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Wet to Dry Controls Lanthanide Scandate Synthesis

The choice of temperature and gas conditions used in a water pressure-controlled reactor is guided by density functional theory (DFT) to synthesize nearly phase pure lanthanide scandate nanoparticles (LnScO 3 , Ln = La, Nd, Sm, Gd). In this synthetic method, low water vapor partial pressures, well below water’s gas liquidus, inhibit particle growth, while an excess of water vapor results in undesired rare earth hydroxide and oxyhydroxide secondary phases. The optimal humidity for high-purity LnScO 3 particle synthesis is shown to vary with the lanthanide; DFT is used to calculate the thermodynamics of secondary phase formation for each lanthanide tested, such that the role of water vapor may be quantified and used to maintain phase purity (greater than 96 mol%) across the series. Furthermore, the combination of thermodynamic calculation and experimental confirmation with this pressure-controlled reactor provides an opportunity to explore analogous syntheses of other inorganic perovskite nanoparticles.

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

Passivating Nucleobases Bring Charge Transfer Character to Optically Active Transitions in Small Silver Nanoclusters

DNA-wrapped silver nanoclusters (DNA–AgNCs) are known for their efficient luminescence. However, their emission is highly sensitive to the DNA sequence, the cluster size, and its charge state. To get better insights into photophysics of these hybrid systems, simulations based on density functional theory (DFT) are performed. Our calculations elucidate the effect of the structural conformations, charges, solvent polarity, and passivating bases on optical spectra of DNA–AgNCs containing five and six Ag atoms. It is found that inclusion of water in calculations as a polar solvent media results in stabilization of nonplanar conformations of base-passivated clusters, while their planar conformations are more stable in vacuum, similar to the bare Ag5 and Ag6 clusters. Cytosines and guanines interact with the cluster twice stronger than thymines, due to their larger dipole moments. In addition to the base–cluster interactions, hydrogen bonds between bases notably contribute to the structure stabilization. While the relative intensity, line width, and the energy of absorption peaks are slightly changing depending on the cluster charge, conformations, and base types, the overall spectral shape with five well-resolved bands at 2.5–5.5 eV is consistent for all structures. Independent of the passivating bases and the cluster size and charge, the low energy optical transitions at 2.5–3.5 eV exhibit a metal to ligand charge transfer (MLCT) character with the main contribution emerging from Ag-core to the bases. Cytosines facilitate the MLCT character to a larger degree comparing to the other bases. However, the doublet transitions in clusters with the open shell electronic structure (Ag5 and Ag6 + ) result in appearance of additional red-shifted (<2.5 eV) and optically weak band with negligible MLCT character. The passivated clusters with the closed shell electronic structure (Ag5 + and Ag6) exhibit higher optical intensity of their lowest transitions with much higher MLCT contribution, thus having better potential for emission, than their open shell counterparts.

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

Constitutive Model of Radiation Aging Effects in Filled Silicone Elastomers under Strain

Filled silicone elastomers, an essential component in many technological applications, are often subjected to controlled or unintended radiation for a variety of reasons. Radiation exposure can lead to permanent mechanical and structural changes in the material, which is manifested as altered mechanical response, and in some cases, a permanent set. For unfilled elastomers, network theories developed and refined over decades can explain these effects in terms of chain-scission and cross-link formation and a hypothesis involving independent networks formed at different strain levels of the material. In this work, we expose a filled silicone rubber to gamma radiation while being under finite elongational strain and show that the observed mechanical and structural changes can be quantitatively modeled within the same theoretical framework developed for unfilled elastomers as long as nuances associated with the Mullins effect are accounted for in a consistent manner. In this work, we employ Ogden’s incompressible hyperelastic model within the framework of Tobolsky’s two-network scheme to describe the observed permanent set and mechanical modulus changes as a function of radiation dosage. In the process, we conclude that gamma radiation induces both direct cross-linking at chain crossings (H-links) and main-chain-scission followed by cross-linking (Y-links). We provide an estimate of the ratio of chain-scission to cross-linking rates, which is in reasonable agreement with previous experimental estimate from Charlesby–Pinner analysis. We use density functional theory (DFT)-based quantum mechanical calculations to explore the stability of –Si and –SiO radicals that form upon a radiation-induced chain-scission event, which sheds light on the relative rates of Y-linking and H-linking processes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Database of Nonaqueous Proton-Conducting Materials

This work presents the assembly of 48 papers, representing 74 different compounds and blends, into a machine-readable database of nonaqueous proton-conducting materials. SMILES was used to encode the chemical structures of the molecules, and we tabulated the reported proton conductivity, proton diffusion coefficient, and material composition for a total of 3152 data points. The data spans a broad range of temperatures ranging from -70 to 260 °C. To explore this landscape of nonaqueous proton conductors, DFT was used to calculate the proton affinity of 18 unique proton carriers. The results were then compared to the activation energy derived from fitting experimental data to the Arrhenius equation. It was found that while the widely recognized positive correlation between the activation energy and proton affinity may hold among closely related molecules, this correlation does not necessarily apply across a broader range of molecules. This work serves as an example of the potential analyses that can be conducted using literature data combined with emerging research tools in computation and data science to address specific materials design problems.

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

Ultrafast conformational dynamics of Rydberg-excited N -methyl piperidine

We have observed the ultrafast conformational dynamics of electronically excited N-methyl piperidine (NMP) using time-resolved Rydberg fingerprint spectroscopy. Optical excitation at various wavelengths ranging from 212 nm to 229 nm leads to the 3s or 3p Rydberg states and induces coherent oscillatory motions with periods of about 700 fs. Furthermore, these coherent motions survive the internal conversion from 3p to 3s but then dephase on a time scale of a few oscillations. Intramolecular vibrational energy redistribution on a picosecond time scale leads to an equilibrium between two conformeric structures that are separated in binding energy by 0.09 eV. Model calculations using the DFT-SIC method are in excellent agreement with the experiments and identify the conformers as the chair and twist structures of NMP. The analysis of the equilibrium parameters at long time delays as a function of excitation wavelength allows for the extraction of thermodynamic parameters for the conformeric transformation. We derive an enthalpy of the chair to twist reaction in the 3s excited state of 62 meV with an entropy of 19.70 J mol –1 K –1 . An activation energy of 276 meV is also obtained with a kinetic model.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Atomistic mechanisms for catalytic transformations of NO to NH3, N2O, and N2 by Pd

The industrial pollutant NO is a potential threat to the environment and to human health. Thus, selective catalytic reduction of NO into harmless N2, NH3, and/or N2O gas is of great interest. Among many catalysts, metal Pd has been demonstrated to be most efficient for selectivity of reducing NO to N2. However, the reduction mechanism of NO on Pd, especially the route of N−N bond formation, remains unclear, impeding the development of new, improved catalysts. We report here the elementary reaction steps in the reaction pathway of reducing NO to NH3, N2O, and N2, based on density functional theory (DFT)-based quantum mechanics calculations. We show that the formation of N2O proceeds through an Eley-Rideal (E−R) reaction pathway that couples one adsorbed NO* with one non−adsorbed NO from the solvent or gas phase. This reaction requires high NO* surface coverage, leading first to the formation of the trans-(NO)2* intermediate with a low N−N coupling barrier (0.58 eV). Notably, trans-(NO)2* will continue to react with NO in the solvent to form N2O, that has not been reported. With the consumption of NO and the formation of N2O* in the solvent, the Langmuir-Hinshelwood (L-H) mechanism will dominate at this time, and N2O* will be reduced by hydrogenation at a low chemical barrier (0.42 eV) to form N2. In contrast, NH3 is completely formed by the L-H reaction, which has a higher chemical barrier (0.87 eV). Our predicted E-R reaction has not previously been reported, but it explains some existing experimental observations. In addition, we examine how catalyst activity might be improved by doping a single metal atom (M) at the NO* adsorption site to form M/Pd and show its influence on the barrier for forming the N−N bond to provide control over the product distribution.

Physics↗
DOE OSTI2024

Segregation of Al and its effect on coercivity in Nd-Fe-B

Doping a small amount of Al can effectively enhance coercivity in Nd-Fe-B magnets. Here we investigated the partitioning behavior of Al and its effect on coercivity in Nd-Fe-B using first principles DFT (density functional theory) calculation and micromagnetic simulation. The calculated substitution energies of Fe by Al are negative at the crystallographic sites of 4c and 8j 2 while they are positive values at the other sites in Nd 2 Fe 14 B (2:14:1), implying a small solubility of Al in 2:14:1. Further, Al prefers to segregate at grain boundary (GB) and stabilize the Nd-Fe-Al phase with a Nd 6 Fe 13 Si-type tetragonal structure (6:13:1). The formation of the antiferromagnetic or weak ferrimagnetic 6:13:1-like phase depletes Fe and reduces the amount of ferromagnetic Nd–Fe type grain boundary phase (GBP), which weakens the inter-grain magnetic interaction. Micromagnetic simulations indicate that the 6:13:1-like GBP increases the pinning field of magnetic domain wall at GB and suppresses the nucleation of reversal magnetic domain on the grain surface of 2:14:1 during demagnetization process. The formation of Al-rich shell on 2:14:1 grain surface can further moderately increase the domain pinning field at GB and the nucleation field of 2:14:1 grain. Developing novel processing method to tailor Al segregation and promote formation of 6:13:1-like phase at GB can be a promising approach to improve coercivity in Nd-Fe-B magnet.

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

Electronic properties of single-layer CoO 2 /Au(111)

Here we report direct measurements via angle-resolved photoemission spectroscopy (ARPES) of the electronic dispersion of single-layer (SL) CoO 2 . The Fermi contour consists of a large hole pocket centered at the $\overline{\Gamma}$ point. To interpret the ARPES results, we use density functional theory (DFT) in combination with the multi-orbital Gutzwiller Approximation (DFT+GA), basing our calculations on crystalline structure parameters derived from x-ray photoelectron diffraction and low-energy electron diffraction. Our calculations are in good agreement with the measured dispersion. We conclude that the material is a moderately correlated metal. We also discuss substrate effects, and the influence of hydroxylation on the CoO 2 SL electronic structure.

36 MATERIALS SCIENCE↗