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,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Development of molecular cluster models to probe pyrite surface reactivity

Abstract The recent discovery that anaerobic methanogens can reductively dissolve pyrite and utilize dissolution products as a source of iron and sulfur to meet their biosynthetic demands for these elements prompted the development of atomic‐scale nanoparticle models, as maquettes of reactive surface sites, for describing the fundamental redox steps that take place at the mineral surface during reduction. The given report describes our computational approach for modeling n (FeS 2 ) nanoparticles originated from mineral bulk structure. These maquettes contain a comprehensive set of coordinatively unsaturated Fe (II) sites that are connected via a range of persulfide (S 2 2− ) ligation. In addition to the specific maquettes with n = 8, 18, and 32 FeS 2 units, we established guidelines for obtaining low‐energy structures by considering the pattern of ionic, covalent, and magnetic interactions among the metal and ligand sites. The developed models serve as computational nano‐reactors that can be used to describe the reductive dissolution mechanism of pyrite to better understand the reactive sites on the mineral, where microbial extracellular electron‐transfer reactions can occur.

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

Interstitial Atoms and the Frustrated and Allowed Structural Transitions Principle: Tunability in the Electronic Structure of AuCu 3 ‐type Frameworks in Dy 4 T 1− x Ga 12 (T = Ag, Ir)

In this Article, we explore how the chemical pressure (CP) features of an intermetallic phase may provide opportunities to couple perturbations in electron count with the stabilization of the underlying geometrical structure. AuCu 3 ‐type LnGa 3 (Ln = lanthanide or group 3 metal) phases contain octahedral cavities of negative CP held open by overly compressed Ln–Ga interactions, leading to a series of transition metal‐stuffed derivatives. We present new additions to this family with the synthesis and crystal structures of Dy 4 T 1−x Ga 12 with (T, x) = (Ag, 0.29) and (Ir, 0.15), adopting Y 4 PdGa 12 ‐type superstructures of the AuCu 3 ‐type. density functional Ttheory (DFT)‐CP calculations, when adjusted to avoid dipolar CP features, affirm that T atom incorporation provides a mechanism for the relief of packing tensions, while electronic density of states distributions illustrate that the T atoms serve largely as electron or hole donors to the band structure, as needed for them to attain d 10 configurations. The maximum obtainable value for x may be limited by a mismatch between the Fermi energy and pseudogap, in line with the balance of factors envisioned by the frustrated and allowed structural transitions principle. Furthermore, trends in resistivity measurements on T = Ir, Pd, and Ag compounds are interpretable in terms of the varying degrees of disorder arising from x < 1.0.

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

Single-phase Ru 1- x-y Mn x Co y O 2 nanoparticles as highly effective oxygen reduction electrocatalysts in alkaline media with enhanced stability and fuel-tolerance

Here, a series of single-phase binary and ternary ruthenium and transition metal oxide nanoparticles supported on Vulcan - Ru 1-x MnxO2/C, Ru1-xCoxO 2 /C, Ru 1-x-y MnxCo y O 2 /C, Ru 1-x Fe x O 2 /C, Ru 1-x Ni x O 2 /C and Ru 1-x V x O 2 /C, as well as RuO 2 /C, MnO 2 /C, Co 3 O 4 /C and Mn 3-x Co x O 4 /C - were synthesized via a facile method and studied for their oxygen reduction reaction (ORR) activity in alkaline media. Single-phase Ru 1-x Mn x O 2 /C, Ru 1-x-y Mn x Co y O 2 /C and Ru 1-x Co x O 2 /C catalysts significantly boosted ORR kinetics in alkaline media. Moreover, they were H 2 and methanol tolerant, and exhibited long term stability due to their single-phase structure. Ru 0.85 Mn 0.15 O 2 /C stood out as the most active catalyst among them. A volcano relationship of the ORR activity of these binary and ternary metal oxide catalysts vs. the O adsorption energy was found, with the maximum activity observed for the Ru 1-x Mn x O 2 /C, Ru 1-x-y Mn x Co y O 2 /C and Ru 1-x Co x O 2 /C catalysts due to their optimal O binding energy. These materials have potential applications as highly active ORR catalysts in alkaline fuel cells and metal-air batteries.

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

Coverage-dependent structures and thermodynamic stability of intercalated Gd layers beneath buffer-layer graphene on SiC(0001)

Electronic properties of two-dimensional (2D) materials are strongly influenced by their atomic arrangements, making the theoretically-aided characterization of experimentally-synthesized 2D structures crucial. Using first-principles density functional theory, we analyze nearly 200 configurations of intercalated Gd layers beneath buffer-layer graphene on SiC(0001) over a Gd coverage range of 0.01 < θ < 1.2. By fully relaxing selectively-constructed configurations at each coverage within a large, low-strain supercell, we determine the coverage dependence of the chemical potential for intercalated Gd structures. Thermodynamically-preferred configurations below θ ≈ 0.8 form single-atom-thick monolayers, while 3D-like or multilayer structures emerge beyond θ ≈ 0.9. Most structures are amorphous-like, including the configuration at the chemical potential minimum around θ ≈ 0.4. In contrast, a strongly stretched Gd(0001)-like monolayer at θ = 1/3 and a nearly perfect Gd(0001) monolayer at θ = 1 are significantly less favorable with 0.16 eV and 0.82 eV higher chemical potentials above the minimum, respectively. Furthermore, the graphene layer decoupled by intercalated Gd near the chemical potential minimum is significantly flatter compared to its morphology above intercalated 3D structures at higher coverages and nearly isolated Gd atoms in the lowest coverage region. In conclusion, these findings align with our experimental results and underscore the need for further research on this unique intercalated system, which holds significant potential for diverse applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mechanistic understanding of support effect on the activity and selectivity of indium oxide catalysts for CO 2 hydrogenation

Herein we present a mechanistic study on the support effect (ZrO 2 and CeO 2 ) of In 2 O 3 catalysts in CO 2 hydrogenation by a combined experimental and computational approach. Kinetic experiments and surface characterization suggested that the activity of In 2 O 3 catalysts cannot be simply correlated with the abundance of surface oxygen vacancies (O v ) formed by either H 2 -reduction or thermal treatment, which has been frequently invoked in previous studies. The support effect should originate from the electronic interactions between In 2 O 3 and the support oxide, rather than geometric factors or the difference in the particle size of In 2 O 3 . Theoretical modelling revealed that surface O v facilitate the formation and stabilization of the formate (HCOO*) intermediate. While a carbonate-like structure is favored for CO 2 adsorption on CeO 2 -supported or unsupported In 2 O 3 catalysts, CO 2 tends to bind strongly in a bent configuration on the O v site at the In 2 O 3 -ZrO 2 interface. The distinct CO 2 adsorption structures on different supported In 2 O 3 catalysts may account for the different reaction energy profiles in the subsequent hydrogenation reactions, especially the rate-limiting step, i.e., hydrogenation of HCOO* to CH 2 O* and methoxy (CH 3 O*). The relatively higher methanol selectivity of In 2 O 3 catalyst supported on ZrO 2 with respect to that on CeO 2 are suggested to stem from the greater energy difference (Δ$E_a$) between the parallel hydrogenation and C-O bond cleavage of HCOO*, which leads to the formation of methanol and CO, respectively. We report this study underlines the important role of metal-oxide-interface in determining the catalytic behavior of oxide-supported In 2 O 3 catalysts in CO 2 conversion.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

An automated cluster surface scanning method for exploring reaction paths on metal-cluster surfaces

Metal-cluster surfaces present a wide variety of unique coordination environments. This complexity makes it difficult to manually probe the surface reactivity of such clusters. Here, we present a simple and automated method to systematically discover reaction pathways on cluster surfaces, based on the automated cluster surface scanning (ACSS) technique for mapping out potential energy surfaces. We showcase our method on 55-atom icosahedral Cu and Ag clusters, where we determine the activation energies of four elementary steps common in heterogeneous catalysis – hydrogen recombination (H* + H* → H 2 * + *), oxygen recombination (O* + O* → O 2 * + *), water formation (OH* + H* → H 2 O* + *), and CO oxidation (CO* + O* → CO 2 * + *) – with density functional theory calculations (DFT-PBE + D3). We show that the ACSS method requires significantly less human effort than the established manually performed climbing-image nudged elastic band (MP + CI-NEB) technique and locates transition states with comparable accuracy (root-mean-squared error of 0.10 eV) and similar computational cost. Rigorous sampling of the potential energy surface with the ACSS method allows one to locate all lowest-energy reaction pathways obtained via the MP+CI-NEB approach, as well as alternative pathways that one may have missed with the MP+CI-NEB approach due to the many possible pathways available on these clusters. The accuracy and efficiency afforded by the ACSS method could enable high-throughput exploration of the diverse reactivity of metal clusters.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electron-phonon coupling from GW perturbation theory: Practical workflow combining BerkeleyGW, ABINIT, and EPW

Here we present a workflow of practical calculations of electron-phonon (e-ph) coupling with many-electron correlation effects included using the GW perturbation theory (GWPT). This workflow combines BerkeleyGW, ABINIT, and EPW software packages to enable accurate e-ph calculations at the GW self-energy level, going beyond standard calculations based on density functional theory (DFT) and density-functional perturbation theory (DFPT). This workflow begins with DFT and DFPT calculations (ABINIT) as starting point, followed by GW and GWPT calculations (BerkeleyGW) for the quasiparticle band structures and e-ph matrix elements on coarse electron k- and phonon q-grids, which are then interpolated to finer grids through Wannier interpolation (EPW) for computations of various e-ph coupling determined physical quantities such as the electron self-energies or solutions of anisotropic Eliashberg equations, among others. A gauge-recovering symmetry unfolding technique is developed to reduce the computational cost of GWPT (as well as DFPT) while fulfilling the gauge consistency requirement for Wannier interpolation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Low rotational barriers for the most dynamically active methyl groups in the proposed antiviral drugs for treatment of SARS-CoV-2, apilimod and tetrandrine

A recent screening study highlighted a molecular compound, apilimod, for its efficacy against the SARS-CoV-2 virus, while another compound, tetrandrine, demonstrated a remarkable synergy with the benchmark antiviral drug, remdesivir. In this work, we find that because of significantly reduced potential energy barriers, which also give rise to pronounced quantum effects, the rotational dynamics of the most dynamically active methyl groups in apilimod and tetrandrine are much faster than those in remdesivir. Because dynamics of methyl groups are essential for biochemical activity, screening studies based on the computed potential energy profiles may help identify promising candidates within a given class of drugs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Crystal growth, physical and optical properties of TlSr 2 Cl 5 and TlSr 2 Br 5

Small diameter (Ø 16 mm) TlSr 2 Cl 5 and TlSr 2 Br 5 crystals were grown by the Vertical Bridgman method. X-ray diffraction measurements show that both have the monoclinic crystal structure with space group P21/c. TlSr 2 Cl 5 and TlSr 2 Br 5 have a density of 4.14 g/cm 3 and 5.03 g/cm 3 , respectively. The effective Z of TlSr 2 Cl 5 and TlSr 2 Br 5 is 63.7 and 58.6, respectively. Radioluminescence spectra of TlSr 2 Cl 5 and TlSr 2 Br 5 feature a broad emission band peaking at 440 and 445 nm, respectively. As a result the light yield of TlSr 2 Cl 5 and TlSr 2 Br 5 was estimated to be 17,000 and 45,000 ph/MeV, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2020

TlSr 2 I 5 :Eu 2+ - A new high density scintillator for gamma-ray detection

In this paper we report on the scintillation properties of TlSr 2 I 5 doped with Eu 2+ , a novel thallium-containing high-resolution scintillator for gamma-ray spectroscopy. Here, small diameter, good quality crystals of TlSr 2 I 5 :Eu (TSI) with different Eu 2+ concentrations were grown by the vertical Bridgman method. X-ray diffraction measurements show that single crystals of TSI belong to the monoclinic system with space group P2 1 /c. TlSr 2 I 5 has a density of 5.32 g/cm 3 and effective Z of 60. The X-ray excited emission of Eu 2+ doped TlSr 2 I 5 features a broad emission band peaking between 460–470 nm. The light yield of TlSr 2 I 5 crystal doped with 1% Eu 2+ is measured to be ~72,000 ph/MeV with an energy resolution of 2.8% at 662 keV. The scintillation decay time which is characteristic of Eu 2+ shows two components, with 90% of the light in a ~500 ns component and the rest in a longer component of ~3 μs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Effects of co-adsorbed water on different bond cleavages involved in acetic acid decomposition on Pt (111)

Acetic acid decomposition on Pt (111) in the presence of co-adsorbed water is a good model system for oxygenate decomposition on Pt (111) in aqueous phase, with application in hydrogen production and biomass conversion. Here, in this study, we present a density functional theory (DFT) theory calculation of how co-adsorbed water affects different bond cleavages of acetic acid decomposition on Pt (111). The presence of co-adsorbed water generally enhances O$-$H bond cleavage while inhibiting OC$-$O and OC$-$OH bond cleavage. The influence of co-adsorbed water on C$-$H bond cleavage varies the most and depends on the nature of the transition state and how co-adsorbed water stabilizes the initial and final state. Although these trends are useful as general guidance, they are not sufficient to predict the effect on a complex reaction network such as acetic acid decomposition on Pt (111) which has several parallel reaction paths with similar energies. In the absence of co-adsorbed water, the two lowest energy pathways are decarboxylation (DCX) and decarbonylation (DCN) pathways through a common CH 2 COO intermediate, in which the DCX pathway (C$-$H bond cleavage of CH 2 COO) is more favorable than DCN pathway (OC$-$O bond cleavage of CH 2 COO). In the presence of co-adsorbed water, the energy difference between Csingle bondH bond cleavage and OC$-$O bond cleavage of CH 2 COO increases, suggesting an increase in the favorability of acetic acid decarboxylation (formation of carbon dioxide) over acetic acid decarbonylation (formation of carbon monoxide) on Pt (111).

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

Exploring Ca–Ce–M–O (M = 3d Transition Metal) Oxide Perovskites for Solar Thermochemical Applications

Solar thermochemical (STC) processes hold promise as efficient ways to generate renewable fuels, fuel precursors, or chemical feedstocks using concentrated sunlight. Specifically, one actively researched approach is the two-step STC cycle, which uses a redox-active, off-stoichiometric, transition-metal oxide material to split water and/or CO 2 , generating H 2 and/or CO, respectively, or syngas (a combination of H 2 and CO). Identifying novel metal oxides that yield larger reduction extents (practically achievable off-stoichiometries) than the state-of-the-art CeO 2 is critical. Here, we explore the chemical space of Ca–Ce–M–O (M = 3d transition metal, except Cu and Zn) metal oxide perovskites, with Ca and/or Ce occupying the A site and M occupying the B site within an ABO 3 framework, as potential STC candidates. We use density functional theory (DFT)-based calculations and systematically evaluate the oxygen vacancy (VaO) formation energy (≈ enthalpy of reduction in an STC cycle), electronic properties, thermodynamic stability of CaMO 3 , CeMO 3 , and Ca 0.5 Ce 0.5 MO 3 perovskites, and the VaO formation energy within Ca 0.5 Ce 0.5 Ti 0.5 Mg 0.5 O 3 perovskite. We consider only Ca and/or Ce on the A site because of their similar size and the potential redox activity of Ce 4+ . If both Ce and M exhibit simultaneous reduction with Va O formation, the resulting perovskite could exhibit a larger entropy of reduction than a single cation reduction. The increased entropy produces increased reduction for fixed temperature, partial pressure of oxygen, and reduction enthalpy, and therefore increased STC efficiency. Importantly, we identify Ca 0.5 Ce 0.5 MnO 3 , Ca 0.5 Ce 0.5 FeO 3 , and Ca 0.5 Ce 0.5 VO 3 to be promising candidates based on their Va O formation energy and thermodynamic (meta)stability. Moreover, based on our calculated on-site magnetic moments, electron density of states, and electron density differences between pristine and defective structures, we find Ca 0.5 Ce 0.5 MnO 3 to exhibit simultaneous reduction of both Ce 4+ (A-site) and Mn 3+ (B-site), highlighting a particularly promising candidate for STC applications with a predicted higher entropy of reduction than CeO 2 . Lastly, we extract metrics that govern the trends in Va O formation energies, such as standard reduction potentials, and provide pointers for further experimental and theoretical studies, which will enable the design of improved materials for the STC cycle.

14 SOLAR ENERGY↗
DOE OSTI2025

Structural Distortions and Uniaxial Negative Thermal Expansion in the Polar Dion–Jacobson Oxide RbNdTa 2 O 7

We provide deeper insight into the crystal structures, sequential structural phase transitions (I2cm → Cmce → I4/mcm → P4/mmm), thermal expansion, and electronic properties of the n = 2 Dion–Jacobson polar oxide RbNdTa 2 O 7 , through X-ray powder diffraction, neutron powder diffraction, Raman studies, and density functional theory calculations. We observed a uniaxial negative thermal expansion (NTE) across the first-order transition, I2cm → Cmce, where the unit cell contracts along the c-axis, which is driven by a contraction of the NdTa 2 O 6 layer. Here, this NTE occurs within the temperature range of the first-order phase transition and contrasts with the corkscrew mechanism typically observed in Ruddlesden–Popper phases. In RbNdTa 2 O 7 , the I2cm (hybrid improper ferroelectric) → Cmce (antipolar) transition involves crucial changes in the bond lengths of Nd and Ta polyhedra, coupled with polar to antipolar displacement of the Nd ions, leading to a net contraction in the NdTa 2 O 6 layer along the c-axis, while preserving the overall octahedral tilting magnitude. This transition highlights the intricate interplay between the Nd and Ta coordination and the associated TaO 6 distortions. Temperature-dependent Raman spectra analysis further confirms the first-order structural transition and associated NTE, providing evidence for increased bond stiffness across this transition. Additionally, using neutron powder diffraction, we have determined that the transition I4/mcm → P4/mmm occurs at approximately 1150 K. Finally, we have calculated from DFT + U, the partial density of states, the energy bandgaps, and effective masses of the charge carriers of the polar ground structure.

Chemical structure↗
DOE OSTI2025

Toward High-Voltage Cathodes for Zinc-Ion Batteries: Discovery Pipeline and Material Design Rules

Efficient energy storage systems are crucial to address the intermittency of renewable energy sources. As multivalent batteries, Zn-ion batteries (ZIBs), while inherently low voltage, offer a promising low-cost alternative to Li-ion batteries due to the viable use of zinc as the anode. However, to maximize the potential impact of ZIBs, rechargeable cathodes with improved Zn diffusion are needed. To better understand the chemical and structural factors influencing Zn-ion mobility within battery electrode materials, we employ a high-throughput computational screening approach to systematically evaluate candidate intercalation hosts for ZIB cathodes, expanding the chemical search space on empty intercalation hosts that do not contain Zn. We leverage a high-throughput screening funnel to identify promising cathodes in ZIBs, integrating screening criteria with density functional theory (DFT)-based calculations of Zn2+ intercalation and diffusion inside the host materials. Using these data, we identify the design principles that favor Zn-ion mobility in candidate cathode materials. Building on previous work on divalent-ion cathodes, this study broadens the chemical space for next-generation multivalent energy storage systems.

electrodes↗