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

PFOA catalytic reaction mechanisms on zerovalent iron (Fe 0 ): A DFT investigation

In this study, density functional theory (DFT) calculations with dispersion corrections are employed to investigate catalytic reactions of perfluorooctanoic acid (PFOA) onto zerovalent iron (Fe 0 ). The main goal of this investigation is to explain the reaction mechanisms, including bond dissociation energies, activation barriers, and rate-determining steps (RDSs), of PFOA degradation. Additionally, we aim to reveal the catalytic effects of Fe 0 by comparing the reaction energy profile for the degradation of PFOA both in solution (i.e., isolated) and after adsorption on Fe 0 . Along with investigating previously reported pathways, this study proposes a new PFOA degradation pathway on Fe0. Overall, the results revealed that Fe 0 plays a key role in the catalytic degradation of PFOA, significantly affecting the RDSs and increasing defluorination rates. For instance, the activation barrier for the decarboxylation step significantly decreased from 3.08 eV to 1.11 eV over the Fe0 surface, while the activation barrier for F dissociation decreased from ∼5 eV to 0.54 eV. Additionally, the proposed pathway in this study (i.e., successive defluorination) was found to be the most favorable pathway for PFOA degradation on Fe 0 . Taken together, the DFT results are consistent with reported experimental findings and provide critical insights that will help advance the design of catalysts for PFAS degradation.

Catalysis↗
DOE OSTI2023

Electronic Structure and Mechanical Properties of Solvated Montmorillonite Clay Using Large-Scale DFT Method

Montmorillonite clay (MMT) has been widely used in engineering and environmental applications as a landfill barrier and toxic waste repository due to its unique property as an expandable clay mineral that can absorb water easily. This absorption process rendered MMT to be highly exothermic due to electrostatic interactions among molecules and hydrogen bonds between surface atoms. A detailed study of a large supercell model of structural clay enables us to predict long-term nuclear waste storage. Herein, a large solvent MMT model with 4071 atoms is studied using ab initio density functional theory. The DFT calculation and analysis clarify the important issues, such as bond strength, solvation effect, elasticity, and seismic wave velocities. These results are compared to our previous study on crystalline MMT (dry). The solvated MMT has reduced shear modulus (G), bulk modulus (K), and Young’s modulus (E). We observe that the conduction band (CB) in the density of states (DOS) of solvated MMT model has a single, conspicuous peak at –8.5 eV. Moreover, the atom-resolved partial density of states (PDOS) summarizes the roles played by each atom in the DOS. These findings illuminate numerous potential sophisticated applications of MMT clay.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Thermal Stability of Schottky Contacts and Rearrangement of Defects in β ‐Ga 2 O 3 Crystals

Abstract The thermal stability of different Schottky contacts (Au, Pt, and Ni) on (100) β ‐Ga 2 O 3 single crystals grown by the Czochralski method is investigated. Besides the examination of the Schottky barrier parameters, contact‐dependent defect levels are investigated by deep‐level transient spectroscopy (DLTS) in a 100–650 K (ramp‐up) and 650–100 K (ramp‐down) temperature cycle. Several defect levels are detected below the conduction band minimum at 0.41, 0.60, 0.77, 0.96, and 1.17 eV. In the temperature ramp‐down DLTS, the 1.17 eV level disappears, and the 0.60 eV level appears for all Schottky contacts. DFT calculations suggest that rearrangement and dissociation of a single hydrogen from a doubly‐hydrogenated Ga─O divacancy complex occurs during the temperature sweep under bias. The trap level at 0.96 eV only appears after the thermal load for the Ni contact, in contrast to Au and Pt, where it is present without a thermal budget. Temperature‐dependent leakage current (at −4 V) measurements indicate oxidation of Ni, and further thermodynamic analysis suggests alloying of Au‐Ga atoms at the Au/ β ‐Ga 2 O 3 interface. These studies provide insight into the behavior induced by these common Schottky contacts and the alteration associated with temperature cycling.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Thermoelectric Behavior of Structurally Disordered Lanthanide and Actinide Silicides

Lanthanide and actinide intermetallic materials exhibit intriguing electronic, structural, and magnetic properties arising from their partially filled f orbitals. Strong f-electron correlations can significantly modify the electronic density of states, carrier scattering, and transport behavior, making these systems promising candidates for unconventional thermoelectric materials. In previous studies, single crystals of A1.33T4Al8Si2 (A = Ce, Th, U, Np; T = Co, Ni) were synthesized. This structure consists of ordered TAl2 corrugated double layers separated by disordered An–Si monolayers. The inherent disorder within the An–Si layers has been linked to unusual magnetic behavior, such as spin-glass phenomena, and is expected to strongly influence phonon scattering and charge-carrier transport relevant to thermoelectric performance. However, the thermoelectric properties of this material family have not previously been explored. In this work, DFT calculations were performed to evaluate the electronic density of states, band structures, and thermoelectric properties, including the effects of doping. To validate these predictions, Ce1.33Co4Al8Si2 single crystals were grown and experimentally characterized. Additionally, a new Yb1.33Co4Al8Si2 phase was synthesized and characterized using single-crystal XRD, SEM-EDS, magnetization measurements, and thermoelectric measurements. These results provide insight into structure–property relationships governing thermoelectric behavior in f-electron intermetallics and offer guidance for future thermoelectric materials.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2021

Mechanistic Investigation of Dimethylmercury Formation Mediated by a Sulfide Mineral Surface

Mercury (Hg) pollution is a global environmental problem, but many key transformations among 29 various Hg species are not well understood. Although anaerobic microorganisms are the primary 30 producers of monomethylmercury (MMeHg), the abiotic formation of dimethylmercury (DMeHg) 31 from MMeHg have been suggested to account for a large portion of DMeHg in the ocean. It has 32 been shown that abiotic formation of DMeHg from MMeHg can be facilitated by reduced sulfur 33 groups on mineral surfaces. However, a detailed mechanistic understanding of this reaction is 34 lacking. Here we perform density functional theory (DFT) calculations to explore mechanisms of 35 DMeHg formation on the 110 surface of a CdS(s) (hawleyite) nanoparticle. We show that 36 coordination of MMeHg substituents to adjacent reduced sulfur groups protruding from the surface 37 facilitate DMeHg formation and that the reaction proceeds through direct transmethylation from 38 one MMeHg substituent to another. Solvation effects were found to play an important role in the 39 surface reconstruction of the nanoparticle and in decreasing the energetic barrier for DMeHg 40 formation. We also investigate DMeHg formation from MMeHg mediated by organic dithiols and 41 compare the results to those mediated by the mineral surface. In contrast to the surface reaction, 42 the calculations suggest that dithiol-mediated formation of DMeHg instead proceeds through a 43 mechanism in which Hg–Hg interactions activate a methyl group from one MMeHg substituent 44 for intramolecular transfer to the adjacent MMeHg substituent. These findings refine previously 45 proposed mechanisms of abiotic DMeHg formation and fill gaps in our understanding of global 46 Hg cycling in the environment.

Lian, Peng↗
DOE OSTI2020

Selective acetylene hydrogenation over single metal atoms supported on Fe 3 O 4 (001): A first-principle study

Supported single-atom catalysts (SACs) have gained increasing attention for improved catalytic activity and selectivity for industrially-relevant reactions. In this study, we explore the hydrogenation of acetylene over single Pt, Ru, Rh, Pd, and Ir atoms supported on the Fe 3 O 4 (001) surface using density functional theory (DFT) calculations. The thermodynamic profile of H diffusion is significantly changed with the type of single metal atoms used and suggests that H spillover from the single atom dopant to the Fe 3 O 4 (001) surface is favored and will likely lead to high H coverages in the functioning catalyst. Correspondingly, as the surface H coverage increases, the important desorption step of ethylene becomes energetically competitive against the detrimental hydrogenation steps of ethylene to ethane. Overall, we show herein that the selectivity of SACs towards ethylene production can be tuned by considering the proper type of metal and controlling the redox state of the support.

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

Enhanced Oxide Reduction by Hydrogen at Cuprous Oxide–Copper Interfaces near Ascending Step Edges

We report that understanding the dynamic processes involved in the interaction of hydrogen with oxides is of fundamental importance in catalysis. This paper probes the reduction of Cu 2 O-‘29’ surfaces by hydrogen at room temperature combining in situ ambient pressure scanning tunneling microscopy (STM), X-ray photoelectron spectroscopy (XPS), and density functional theory (DFT) calculations. Reduction of the atomic layer thin Cu 2 O film is observed to be preferentially initiated2 at step edges and terrace defects, where perfect Cu 2 O(111) terraces are found to be stable towards hydrogenation under the same conditions. After a long induction period, regions of partially reduced Cu 2 O-‘29’ and metallic Cu co-exist before the surface is fully reduced to Cu(111). The reduction rate strongly depends on the nature of nearby Cu step edges. We propose a mechanism for the reduction of Cu 2 O-‘29’ by hydrogen where free copper atoms from ascending metallic step edges facilitate the formation of active ensembles for H 2 dissociation and transfer H to the edges of Cu2O regions

08 HYDROGEN↗
DOE OSTI2024

Rationalizing Acidic Oxygen Evolution Reaction over IrO 2 : Essential Role of Hydronium Cation

The development of active, stable, and more affordable electrocatalysts for acidic oxygen evolution reaction (OER) is of great importance for the practical application of electrolyzers and the advancement of renewable energy conversion technologies. Currently, IrO 2 is the only catalyst with high stability and activity, but a high cost. Further optimization of the catalyst is limited by the lack of understanding of catalytic behaviors at the acid-IrO 2 interface. Here, in strong interaction with the experiment, we develop an explicit model based on grand-canonical density function theory (GC-DFT) calculations to describe acidic OER over IrO 2 . Compared to the explicit models reported previously, hydronium cations (H 3 O + ) are introduced at the electrochemical interface in the current model. As a result, a variation in stable IrO 2 surface configuration under the OER operating condition from previously proposed complete *O-coverage to a mixture coverage of *OH and *O is revealed, which is well supported by in situ Raman measurements. In addition, the accuracy of predicted overpotential is increased in comparison with the experimentally measured. More importantly, in this study, an alteration of the potential limiting step from previously identified *O→*OOH to *OH→*O is observed, which opens new opportunities to advance the IrO 2 -based catalysts for acidic OER.

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

Nonprecious Single Atom Catalyst for Methane Pyrolysis

The development of a suitable catalytic system for methane pyrolysis reactions requires a detailed investigation of the activation energy of C-H bonds on catalysts, as well as their stability against sintering and coke formation. In this work, both single-metal Ni atoms and small clusters of Ni atoms deposited on titanium nitride (TiN) plasmonic nanoparticles were characterized for the C-H bond activation of a methane pyrolysis reaction using ab initio spin-polarized density functional theory (DFT) calculations. The present work shows the complete reaction pathway, including energy barriers for C-H bond activation and dehydrogenated fragments, during the methane pyrolysis reaction on catalytic systems. Interestingly, the C-H bond activation barriers were low for both Ni single-atom and Ni-clusters, showing the energy barriers of ~1.10 eV and ~0.88 eV, respectively. Additionally, single-atom Ni-TiN showed weaker binding to adsorbates, and a net endothermic reaction pathway indicated that the single-atom Ni-TiN was expected to resist coke formation on its surface. However, these Ni single-atom catalysts can sinter, aggregate into a small cluster, and form a coke layer from the highly exothermic reaction pathway that the cluster takes despite the facile reaction pathway.

08 HYDROGEN↗
DOE OSTI2021

Canted antiferromagnetic order and spin dynamics in the honeycomb-lattice compound Tb 2 Ir 3 Ga 9

Here, single-crystal neutron diffraction, inelastic neutron scattering, bulk magnetization measurements, and first-principles calculations are used to investigate the magnetic properties of the honeycomb lattice Tb 2 Ir 3 Ga 9 . While the R ln 2 magnetic contribution to the low-temperature entropy indicates a J eff = 1 / 2 moment for the lowest-energy crystal-field doublet, the Tb 3 + ions form a canted antiferromagnetic structure below 12.5 K. Due to the Dzyaloshinskii-Moriya interactions, the Tb moments in the a b plane are slightly canted towards b with a canted moment of 1.22 μ B per formula unit. A minimal x x z spin Hamiltonian is used to simultaneously fit the spin-wave frequencies along the high-symmetry directions and the field dependence of the magnetization along the three crystallographic axes. Long-range magnetic interactions for both in-plane and out-of-plane couplings up to the second nearest neighbors are needed to account for the observed static and dynamic properties. The z component of the exchange interactions between Tb moments is larger than the x and y components. This compound also exhibits bond-dependent exchange with negligible nearest-neighbor exchange coupling between moments parallel and perpendicular to the 4 f orbitals. Despite the J eff = 1 / 2 moments, the spin Hamiltonian is denominated by a large in-plane anisotropy K z ~ – 1 meV . DFT calculations confirm the antiferromagnetic ground state and the substantial interplane coupling at larger Tb-Tb distances.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Elucidating the Chemical Pathways Responsible for the Sooting Tendency of 1 and 2-phenylethanol

Yield Sooting Index (YSI) measurements have shown that oxygenated aromatic compounds (OACs) tend to have lower YSI than aromatic hydrocarbon (AHC) compounds. For example, typical AHCs such as toluene and ethyl benzene have a YSI of 170 and 216, respectively, in contrast, OACs such as phenol and anisole have a YSI of 81 and 111, respectively. However, this trend is not always true as was observed for the structural isomers 1-phenylethanol (1PE, YSI=142) and 2-phenylethanol (2PE, YSI=207), where 2PE contains a YSI more representative of AHCs than OACs. We applied flow reactor experiments and density functional theory (DFT) calculations to examine how oxygen functionality present in 1PE and 2PE alters the reaction pathways leading to the observed difference in soot formation. It was determined that the proximity of the oxygen functional group to the aromatic ring determines whether the oxygen remains attached to the primary reacting species (for 1PE) or is eliminated early in the combustion sequence (for 2PE). For these alcohols, preservation of the oxygen in the molecule leads to further OACs, while loss of the oxygen leads to AHCs and benzyl radical. The direct pathways to AHCs and benzyl radical result in the higher YSI observed for 2PE.

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

Understanding of Ag Nanocatalysts for Electrocatalytic CO2 Conversion: Effects of Particle Size and Carbon Support

In this talk, we combined ultrahigh vacuum (UHV) surface science techniques, electrochemical measurements, and computational modeling to investigate Ag based electrocatalysts for CO2 reduction reaction (CO2RR). Our goal is to understand the critical characteristics governing the activity and selectivity of Ag electrocatalysts. Ag electrocatalysts were grown on highly oriented pyrolytic graphite (HOPG) in the UHV chamber, characterized with X-ray photoelectron spectroscopy (XPS) and scanning tunneling microscopy (STM), and then tested in a custom-built gastight H-cell. Supported by computational modeling based on density functional theory (DFT) calculations and microkinetic modeling (MKM), our studies revealed a strong size-dependent electrocatalytic CO2-to-CO conversion of the Ag nanoparticle electrocatalysts with average particle diameter between 2 to 6 nm. Smaller diameter (< 3 nm) particles favored H2 evolution reaction (HER) due to a high population of Ag edge sites, whereas larger diameter particles favored CO2RR as the population of Ag(100) surface sites grew. We further discovered that electronic interactions between small diameter Ag particles and highly defective carbon supports could break the size-dependent CO2RR reactivity, resulting in highly selective (CO Faradaic Efficiency > 90%) and active Ag nanoparticle electrocatalysts with sizes < 2 nm diameter. This knowledge is key to understand electrocatalysts performance and to ultimately guide electrocatalyst design

Ag nanoparticles↗
DOE OSTI2025

Screening of transition metals for PFAS adsorption: A comparative DFT investigation

Transition metals are promising candidates for catalytic destruction of per- and polyfluoroalkyl substances (PFAS). This study employed density functional theory (DFT) calculations to systematically investigate the adsorption of PFAS, which is a critical step in the catalytic process, onto a group of transition metals, namely Cu, Pd, Pt, and Rh. The investigation considered the influence of different PFAS characteristics, including chain length, functional group (i.e., sulfonic vs. carboxylic), and protonation state (i.e., protonated, anionic, and deprotonated radical) on the adsorption mechanism. Overall, the findings indicated that the adsorption of PFAS on transition metals exhibited thermodynamically favorable energetics, with Rh showing the most favorable adsorption for PFAS. However, the adsorption strength was significantly dependent on the protonation state of PFAS. Cu was the highest electron donor, while Pt was the highest electron acceptor. Results were compared to experimental studies in the literature for qualitative validation.

Adsorption↗
DOE OSTI2022

Structural and Bonding Analysis in Monomeric Actinide(IV) Oxalate from Th(IV) to Pu(IV): Comparison with the An(IV) Nitrate Series

In this work, single-crystal X-ray diffraction (SC-XRD) structures and Raman spectra of a series of new isomorphous molecular An(IV)-oxalate compounds (Th, U, Np, and Pu) are reported. These complexes are crystallized with cobalt(III) hexamine ([Co(NH 3 ) 6 ] 3+ ) as the counter cations, [Co(NH 3 ) 6 ] 2 [An(C 2 O 4 ) 5 ]·4H 2 O, revealing five bidentate nonbridging oxalate ligands in the first coordination sphere (CN = 10). The nonbridging oxalate is rather uncommon for An(IV)-oxalate systems, which are widely characterized as polymeric compounds. Density functional theory (DFT) calculations were performed to examine the bonding between An(IV) cations and oxalate ligands. For comparison, we also report results obtained for the An(IV)-hexanitrate series, [(C 2 H 5 ) 4 N] 2 [An(NO 3 ) 6 ] (with An = Th, U, Np, Pu, and Ce), which consists of O-donor ligands as well but with a larger coordination number (CN = 12). The bonding analysis confirms that the actinide–oxygen bond is predominantly ionic with a minor increase in covalency from Th to U and slight variations from U to Pu. Further comparison showed that the charge transfer increases slightly when increasing the number of anions in the coordination sphere (C 2 O$_4^{2–}$: CN = 10; NO$_3^–$: CN = 12), but covalent effects as indicated by the amount of internuclear electron density accumulation are small and similar for oxalate and nitrate.

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

Effect of impurities on hydrogen defect stability and migration barrier in yttrium dihydride crystal

The impurity or alloying atoms in YH 2 can alter the local electronic structure and so the hydrogen defect stability, as well as the H migration barrier energy. Thus, DFT calculations were employed to determine the effect of foreign elements from alkali and alkaline earth metals to transition metals and one critical impurity element, O, on H vacancy stability and retention characteristics in YH2. Results revealed that alloying elements act as hydrogen vacancy sinks by reducing the vacancy formation energy at neighboring sites. The implantation of non-magnetic foreign elements (s1, s2, and d10 valence electrons) in hydrogen energy landscape was calculated to be minor; while the hydrogen vacancy formation energy was reduced from 1.37 eV to 1.00 eV, the migration energy barrier of hydrogen was increased from 0.87 eV to 1.15 eV for non-magnetic foreign elements. The migration energy barrier monotonically decreased with increasing d-shell occupancy, reaching as low as 0.4 eV for Cr, Mo(d4), and Fe (d4). Alloying with late transition metals (d8 and d9) moderately impacted the hydrogen vacancy formation. Finally, it was found to be O addition into the YH 2- lattice did not alter the energy landscape of hydrogen vacancies. Since alloyed YH 2 has not been studied extensively, this study provides an atomistic understanding how alloying elements and impurities trap vacancies and affects hydrogen mobility YH 2 . Meanwhile, the main findings of this study may serve as guidelines for introducing alloying elements in ZrH 2 as well.

08 HYDROGEN↗
DOE OSTI2025

Mechanistic Insights into Acetate Selectivity on Intermetallic CuPd(110) in CO Reduction

Experimental studies demonstrate that CuPd(110) uniquely favors acetate formation during CO reduction (CORR), contrasting with the preference for ethylene on Cu surfaces. To elucidate this selectivity, we employed explicit solvation density functional theory (DFT) calculations to investigate the reaction mechanism from both thermodynamic and kinetic angles. Here, our findings reveal that on CuPd(110), the acetate-pathway intermediate H 2 CCO is thermodynamically favored at experimental conditions, while CHCHO─a precursor to ethylene─is preferred on Cu(111). Beyond thermodynamics, we find that H 2 CCO is kinetically accessible under the experimental conditions on CuPd(110), aiding acetate formation. Electron density difference analyses further corroborate distinct protonation preferences supporting this mechanism. We propose a thermodynamic screening parameter based on the Gibbs free energy, G H 2 CCO < G CHCHO , as a guide for designing Cu-based catalysts with enhanced acetate selectivity. These results offer critical mechanistic insights into the CORR product distribution and a rational framework for future catalyst design.

Acetate↗
DOE OSTI2023

Single and dual metal atom catalysts for enhanced singlet oxygen generation and oxygen reduction reaction

We demonstrate rational design of graphene-supported single and dual metal atom catalysts (SACs and DACs) for photocatalytic applications, such as singlet oxygen ( 1 O 2 ) sensitization and H 2 O 2 production. Here we combine density functional theory (DFT) and time-dependent DFT (TD-DFT) calculations with experimental verifications. We found a synergistic effect between triplet sensitization and triplet–triplet (Dexter) energy transfer; both play a role in the photocatalytic activity through the volcano plot of 3d transition metal SACs. More specifically, FeN 4 -SAC exhibits a low ISC energy gap (ΔE ISC ) of 0.039 eV, compared with 0.108 eV for FeNiN8-DACs, both possessing a high Bader charge transfer of 0.366 e- and 0.405 e-, respectively. Guided by these computational results, we synthesized a series of SACs and a DAC and confirmed their structures with scanning transmission electron microscopy (STEM) along with the X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS). We then confirm their band structures with low-energy inverse photoemission spectroscopy (LEIPS) and UV-vis-NIR. Subsequently, we synthesized the catalysts for the photooxygenation of anthracene and two-electron oxygen reduction reaction (ORR) to measure their photocatalytic activity. We found that H 2 O 2 production through the two-electron ORR competes with the 1 O 2 generation through Dexter energy transfer. FeN 4 -SAC demonstrates a high photooxygenation conversion of 86% and a high 1 O 2 quantum yield of 1.04, obtained from electron spin resonance (ESR) spectroscopy, with low H 2 O 2 production. In contrast, NiN 4 -SAC exhibits a low 1 O 2 generation and a high H 2 O 2 production mainly because of the high Gibbs free energy of the OOH* intermediate. In conclusion, this work proposes an effective DFT-guided strategy for designing SACs and DACs for various photocatalytic applications.

36 MATERIALS SCIENCE↗