Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy density functional”

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

Hourglass-type bulk Ni 3 d band and Ce 4 f Kondo resonance states in the potential topological Kondo semimetal CeNiSn via angle-resolved photoemission spectroscopy

The electronic structure of CeNiSn, which is a potential topological Kondo insulator and a Dirac nodal-loop semimetal, has been investigated by employing temperature (T) dependent angle-resolved photoemission spectroscopy (ARPES). The Fermi surfaces (FSs) and the band structures of CeNiSn for three orthogonal crystallographic planes are measured, in which both the very dispersive bands and the flat bands are observed, having mainly the Ni 3d character and the Ce 4f character, respectively. The measured FSs and ARPES bands agree reasonably well with the density functional theory (DFT) calculations. The Fermi-edge (E F ) photon energy (hν) map along k b (=k (010) ) shows that the metallic E F -crossing states on the (010) surface have the three-dimensional character, suggesting that the observed E F -crossing metallic states do not correspond to the topological surface states of the two-dimensional character. On the other hand, albeit weak, the features of the hourglass-type bulk band crossings are observed along SXS, with the energies and the slopes being similar to those predicted by the DFT calculations, supporting the Dirac semimetallic nature of CeNiSn. In T-dependent ARPES, the Ce 4f Kondo resonance states are clearly revealed at low T, which become much suppressed above ~80 K. This feature is consistent with the Kondo temperature of CeNiSn, estimated from its ρ(T) data. This work demonstrates the importance of the coherent Kondo states in determining the topological properties of CeNiSn.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Accelerating Multiscale Materials Modeling with Machine Learning

The focus of this project is to accelerate and transform the workflow of multiscale materials modeling by developing an integrated toolchain seamlessly combining DFT, SNAP, LAMMPS, (shown in Figure 1-1) and a machine-learning (ML) model that will more efficiently extract information from a smaller set of first-principles calculations. Our ML model enables us to accelerate first-principles data generation by interpolating existing high fidelity data, and extend the simulation scale by extrapolating high fidelity data (10 2 atoms) to the mesoscale (10 4 atoms). It encodes the underlying physics of atomic interactions on the microscopic scale by adapting a variety of ML techniques such as deep neural networks (DNNs), and graph neural networks (GNNs). We developed a new surrogate model for density functional theory using deep neural networks. The developed ML surrogate is demonstrated in a workflow to generate accurate band energies, total energies, and density of the 298K and 933K Aluminum systems. Furthermore, the models can be used to predict the quantities of interest for systems with more number of atoms than the training data set. We have demonstrated that the ML model can be used to compute the quantities of interest for systems with 100,000 Al atoms. When compared with 2000 Al system the new surrogate model is as accurate as DFT, but three orders of magnitude faster. We also explored optimal experimental design techniques to choose the training data and novel Graph Neural Networks to train on smaller data sets. These are promising methods that need to be explored in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Solution–Doped Donor–Acceptor Copolymers Based on Diketopyrrolopyrrole and 3, 3'–Bis (2–(2–(2–Methoxyethoxy) Ethoxy) ethoxy)–2, 2'–Bithiophene Exhibiting Outstanding Thermoelectric Power Factors with p –Dopants

The design of polymeric semiconductors exhibiting high electrical conductivity (σ) and thermoelectric power factor (PF) will be vital for flexible large-area electronics. In this work, four polymers based on diketopyrrolopyrrole (DPP), 2,3-dihydrothieno[3,4-b][1,4]dioxine (EDOT), thieno[3,2-b]thiophene (TT), and 3, 3'-bis (2-(2-(2-methoxyethoxy) ethoxy) ethoxy)-2, 2'-bithiophene (MEET) are investigated as side-chains, with the MEET polymers newly synthesized for this study. These polymers are systematically doped with tetrafluorotetracyanoquinodimethane ( F 4 TCNQ), CF3SO3H, and the synthesized dopant Cp(CN) 3 -(COOMe) 3 , differing in geometry and electron affinity. The DPP-EDOT-based polymer containing MEET as side-chains exhibits the highest conductivity (σ) ≈700 S cm–1 in this series with the acidic dopant (CF 3 SO 3 H). This polymer also shows the lowest oxidation potential by cyclic voltammetry (CV), the strongest intermolecular interactions evidenced by differential scanning calorimetry (DSC), and has the most oxygen-based functionality for possible hydrogen bonding and ionic screening. Other polymers exhibit high σ ≈300–500 S cm–1 and power factor up to 300 µW m –1 K –2 . The mechanism of conductivity is predominantly electronic, as validated by time-dependent conductance studies and transient thermo voltage monitoring over time, including for those doped with the acid. Furthermore, these materials maintain significant thermal stability and air stability over ≈6 weeks. Density functional theory calculations reveal molecular geometries and inform about frontier energy levels. Raman spectroscopy, in conjunction with scanning electron microscopy (SEM-EDS) and x-ray diffraction, provides insight into the solid-state microstructure and degree of phase separation of the doped polymer films. Infrared spectroscopy enables this study to further quantify the degree of charge transfer from polymer to dopant.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Radionuclide incorporation in negative thermal expansion α-Zr(WO 4 ) 2 : A density functional theory study

The incorporation of uranium, plutonium and technetium in the negative thermal expansion (NTE) α-Zr (WO 4 ) 2 has been investigated within the framework of density functional theory (DFT). It is found that the vacancy formation energies of the charged vacancies are overall larger than that of its counterpart neutral Frenkel defects and Schottky defects. DFT calculations suggest that U and Pu substitutions for the Zr site are preferred in α-Zr (WO 4 ) 2 . In case of Tc substitution, both Tc(IV) for the Zr site and Tc(VII) for the W site are considered under oxygen-poor and oxygen-rich conditions, while Tc(VII) substitution can be improved significantly by including Y 2 O 3 (charge compensation).

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

Pseudoatom molecular dynamics plasma microfields

Spectral lines are powerful diagnostic tools for both laboratory and astrophysical plasmas, as their shape is sensitive to the plasma environment. The low-frequency component of the electric microfield is an important input for semi-analytic line broadening codes. Here, in this paper, we detail a new method of calculating plasma microfields using configuration-resolved pseudoatom molecular dynamics. This approach accounts for both quantum atomic structure and N-body effects, similar to density functional theory molecular dynamics, but with less computational cost. We present pseudoatom microfields at plasma conditions relevant for recent high energy density laboratory astrophysics experiments conducted at the Sandia Z-Machine, National Ignition Facility, and Linac Coherent Light Source. Compared to established microfield codes we find moderate deviations at solid density conditions and strong agreement at lower plasma densities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2024

Surface Studies of β-1,3,5,7-Tetranitro-1,3,5,7-Tetrazoctane and Pentaerythritol Tetranitrate from Density Functional Tight-Binding Calculations and Implications on Crystal Shape

Here, we use density functional tight-binding (DFTB) theory to calculate the surface energies of two energetic crystals: monoclinic β-1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (β-HMX) and tetragonal pentaerythritol tetranitrate (PETN). The results are then employed to determine crystal shapes using the Bravais–Friedel–Donnay–Harker, attachment energy, and surface energy models. We find that energy-based models yield predictions in good agreement with experimental observations. Additionally, we propose a simple model that reframes surface energy as a measure of the lost intermolecular interactions during the formation of a surface from the bulk. The model accurately captures the results from the DFTB calculations and enables us to explain and predict surface energies as a function of the local molecular environment.

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

Single-Molecule Kinetics of Styrene Hydrogenation on Silica-Supported Vanadium: The Role of Disorder for Single-Atom Catalysts

We report a theoretical approach for the study of supported atom catalysis is developed based on recent advances in the study of single-molecule kinetics. This view is particularly useful in exhibiting the role of disorder in single-atom and single-site catalysts on amorphous supports. The distribution of passage times (or waiting times) through a complex catalytic network originating from a set of coupled active sites is described by a probability distribution function (PDF), f(t), that reflects the local environment of the reaction center. An efficient algorithm is developed based on the linear algebra of the Markov transition matrix that produces f(t) or its moments. The kinetics of the hydrogenation reaction of styrene on an organovanadium(III) catalyst supported on amorphous silica is studied. A kinetic model consisting of three intertwined catalytic cycles emanating from three chemically distinct active sites is proposed to describe the chemistry. Density functional theory (DFT) calculations are employed to determine the free energy barriers of the reactions, which are used to construct the rate coefficient matrix. The disorder induced by the amorphous support material is divided into a low-dimensional short-range component reflecting the covalent structures near the reaction center and a weaker long-range component modeling the bulk randomness. The results are computed and analyzed for a wide range of concentration values and disorder scenarios. The unusual structure in the f(t) PDF is found to occur for certain cases that reveal the contribution of multiple catalytic pathways acting in concert.

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

Density Functional Theory Study of Iron–Oxygen Divacancies in Magnetite (Fe 3 O 4 ) and Hematite (Fe 2 O 3 )

Density functional theory (DFT) calculations are employed to investigate the formation energies, charge redistribution, and binding energies of iron–oxygen divacancies in magnetite (Fe 3 O 4 ) and hematite (Fe 2 O 3 ). For magnetite, we focus on the low-temperature phase to explore variations with local environments. Building on previous DFT calculations of the variations in formation energies for oxygen vacancies with local charge and spin order in magnetite, we extend this analysis to include octahedral iron vacancies before analyzing the iron–oxygen divacancies. We also assessed the relative stability of iron–oxygen divacancies by comparing their formation energies with those of individual vacancies. Our findings reveal a significant energetic driving force for the formation of divacancy clusters, particularly in magnetite, where divacancies in the +1 charge state exhibit formation energies comparable to those of neutral iron vacancies under oxidizing conditions. In hematite, the results indicate a strong tendency for oxygen vacancies to bind to iron vacancies. These results highlight the significance of iron–oxygen vacancy complexes in the transport properties of iron oxides, with particular relevance to diffusion mechanisms under irradiation conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High-Accuracy Semiempirical Quantum Models Based on a Minimal Training Set

A great need exists for computationally efficient quantum simulation approaches that can achieve an accuracy similar to high-level theories at a fraction of the computational cost. In this regard, we have leveraged a machine-learned interaction potential based on Chebyshev polynomials to improve density functional tight binding (DFTB) models for organic materials. The benefit of our approach is two-fold: (1) many-body interactions can be corrected for in a systematic and rapidly tunable process, and (2) high-level quantum accuracy for a broad range of compounds can be achieved with ~0.3% of data required for one advanced deep learning potential. Our model exhibits both transferability and extensibility through comparison to quantum chemical results for organic clusters, solid carbon phases, and molecular crystal phase stability rankings. Overall, our efforts thus allow for high-throughput physical and chemical predictions with up to coupled-cluster accuracy for systems that are computationally intractable with standard approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Acetic Acid Adsorption and Reactions on Ni(110)

Acetic acid adsorption and reactions at multiple surface coverage values on Ni(110) were studied with temperature-programmed desorption (TPD) and infrared reflection absorption spectroscopy (IRAS) at 90–500 K. The experimental measurements were interpreted with density functional theory (DFT) calculations that provided information on adsorbate geometries, energies, and vibrational modes. Below the monolayer saturation coverage of 0.36 ML at 90 K, acetic acid adsorbs mostly molecularly. Above this coverage, a physisorbed layer is formed with dimers and catemers, without detectable monomers. Dimers and catemers desorb as molecular acetic acid at 157 and 172 K, respectively. Between 90 and 200 K, the O–H bond in acetic acid breaks to form bridge-bonded bidentate acetate that becomes the dominant surface species. Desorption-limited hydrogen evolution is observed at 265 K. However, even after the acetate formation, acetic acid desorbs molecularly at 200–300 K due to recombination. Minor surface species observed at 200 K, acetyls or acetates with a carbonyl group, decompose below 350 K and generate adsorbed carbon monoxide. At 350 K, the surface likely undergoes restructuring, the extent of which increases with acetic acid coverage. The initial dominant bridge-bonded bidentate acetate species formed below 200 K remain on the surface, but they now mostly adsorb on the restructured sites. Here, the acetates and all other remaining hydrocarbon species decompose simultaneously at 425 K in a narrow temperature range with concurrent evolution of hydrogen, carbon monoxide, and carbon dioxide. Above 425 K, only carbon remains on the surface.

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

Active Learning Guided Computational Discovery of Plant-Based Redoxmers for Organic Nonaqueous Redox Flow Batteries

Organic nonaqueous redox flow batteries (ONRFBs) are promising energy storage devices due to their scalability and reliance on sourceable materials. However, finding suitable redox-active organic molecules (redoxmers) for these batteries remains a challenge. Using plant-based compounds as precursors for these redoxmers can decrease their costs and environmental toxicity. In this computational study, flavonoid molecules have been examined as potential redoxmers for ONRFBs. Flavone and isoflavone derivatives were selected as catholyte (positive charge carrier) and anolyte (negative charge carrier) molecules, respectively. To drive their redox potentials to the opposite extremes, in silico derivatization was performed using a novel algorithm to generate a library of > 40000 candidate molecules that penalizes overly complex structures. A multiobjective Bayesian optimization based active learning algorithm was then used to identify best redoxmer candidates in these search spaces. Furthermore, our study provides methodologies for molecular design and optimization of natural scaffolds and highlights the need of incorporating expert chemistry awareness of the natural products and the basic rules of synthetic chemistry in machine learning.

25 ENERGY STORAGE↗
DOE OSTI2021

Nonintuitive Surface Self-Assembly of Functionalized Molecules on Ag(111)

The fabrication of nanomaterials involves selfordering processes of functional molecules on inorganic surfaces. To obtain specific molecular arrangements, a common strategy is to equip molecules with functional groups. However, focusing on the functional groups alone does not provide a comprehensive picture. Especially at interfaces, processes that govern self-ordering are complex and involve various physical and chemical effects, often leading to unexpected structures, as we showcase here on the example of a homologous series of quinones on Ag(111). Naively, one could expect that such quinones, which all bear the same functionalization, form similar motifs. In salient contrast, our joint theoretical and experimental study shows that profoundly different structures are formed. Using a machine-learning-based structure search algorithm, we find that this is due to a shift of the balance of three antagonizing driving forces: adsorbate-substrate interactions governing adsorption sites, adsorbate-adsorbate interactions favoring close packing, and steric hindrance inhibiting certain otherwise energetically beneficial molecular arrangements. The theoretical structures show excellent agreement with our experimental characterizations of the organic/ inorganic interfaces, both for the unit cell sizes and the orientations of the molecules within. The nonintuitive interplay of similarly important interaction mechanisms will continue to be a challenging aspect for the design of functional interfaces. With a detailed examination of all driving forces, we are, however, still able to devise a design principle for self-assembly of functionalized molecules.

structure prediction↗
DOE OSTI2023

Proton Relay for the Rate Enhancement of Electrochemical Hydrogen Reactions at Heterogeneous Interfaces

Proton transfer is critically important to many electrocatalytic reactions, and directed proton delivery could open new avenues for the design of electrocatalysts. However, although this approach has been successful in molecular electrocatalysis, proton transfer has not received the same attention in heterogeneous electrocatalyst design. Here, in this work, we report that a metal oxide proton relay can be built within heterogeneous electrocatalyst architectures and improves the kinetics of electrochemical hydrogen evolution and oxidation reactions. The volcano-type relationship between activity enhancement and pK a of amine additives confirms this improvement; we observe maximum rate enhancement when the pK a of a proton relay matches the pH of the electrolyte solution. Density-functional-theory-based reactivity studies reveal a decreased proton transfer energy barrier with a metal oxide proton relay. These findings demonstrate the possibility of controlling the proton delivery and enhancing the reaction kinetics by tuning the chemical properties and structures at heterogeneous interfaces.

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

Electronically-coupled redox centers in trimetallic cobalt complexes

Synthesis and isolation of molecular building blocks of metal–organic frameworks (MOFs) can provide unique opportunities for characterization that would otherwise be inaccessible due to the heterogeneous nature of MOFs. Herein, we report a series of trinuclear cobalt complexes incorporating dithiolene ligands, triphenylene-2,3,6,7,10,11-hexathiolate (THT) (1 3+ ), and benzene hexathiolate (BHT) (2 3+ ), with 1,1,1,-tris(diphenylphosphinomethyl)ethane (triphos) employed as the capping ligand. Single crystal X-ray analyses of 1 3+ and 2 3+ display three five-coordinate cobalt centers bound to the triphos and dithiolene ligands in a distorted square pyramidal geometry. Cyclic voltammetry studies of 1 3+ and 2 3+ reveal three redox features associated with the formation of mixed valence states due to the sequential reduction of the redox-active metal centers (Co III/II ). Using this electrochemical data, the comproportionality values were determined for 1 and 2 (log K c = 1.4 and 1.5 for 1, and 4.7 and 5.8 for 2), suggesting strong resonance-stabilized coupling of the metal centers, with stronger electronic coupling observed for complex 2 compared to that for complex 1. Cyclic voltammetry studies were also performed in solvents of varying polarity, whereupon the difference in the standard potentials (ΔE 1/2 ) for 1 and 2 was found to shift as a function of the polarity of the solvent, indicating a negative correlation between the dielectric constant of the electrochemical medium and the stability of the mixed valence species. Spectroelectrochemical studies of in situ generated multi-valent (MV) states of complexes 1 and 2 display characteristic NIR intervalence charge transfer (IVCT) bands, and analysis of the IVCT transitions for complex 2 suggests a weakly coupled class II multi-valent species and relatively large electronic coupling factors (1700 cm –1 for the first multi-valent state of 2 2+ , and 1400 and 4000 cm –1 for the second multi-valent state of 2 + ). Here, density functional theory (DFT) calculations indicate a significant deviation in relative energies of the frontier orbitals of complexes 1 3+ , 2 3+ , and 3 + that contrasts those calculated for the analogous trinuclear cobalt dithiolene complexes employing pentamethylcyclopentadienyl (Cp*) as the capping ligand (Co 3 Cp* 3 THT and Co 3 Cp* 3 BHT, respectively), and may be a result of the cationic nature of complexes 1 3+ , 2 3+ , and 3 + .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗