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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Epitaxial growth, magnetoresistance, and electronic band structure of GdSb magnetic semimetal films

Motivated by observations of extreme magnetoresistance (XMR) in bulk crystals of rare-earth monopnictide (RE-V) compounds and emerging applications in novel spintronic and plasmonic devices based on thin-film semimetals, we have investigated the electronic band structure and transport behavior of epitaxial GdSb thin films grown on III-V semiconductor surfaces. The Gd 3 + ion in GdSb has a high spin S=7/2 and no orbital angular momentum, serving as a model system for studying the effects of antiferromagnetic order and strong exchange coupling on the resulting Fermi surface and magnetotransport properties of RE-Vs. Here, we present a surface and structural characterization study mapping the optimal synthesis window of thin epitaxial GdSb films grown on III-V lattice-matched buffer layers via molecular-beam epitaxy. To determine the factors limiting XMR in RE-V thin films and provide a benchmark for band-structure predictions of topological phases of RE-Vs, the electronic band structure of GdSb thin films is studied, comparing carrier densities extracted from magnetotransport, angle-resolved photoemission spectroscopy (ARPES), and density-functional theory (DFT) calculations. ARPES shows a hole-carrier rich, topologically trivial, semimetallic band structure close to complete electron-hole compensation, with quantum confinement effects in the thin films observed through the presence of quantum-well states. DFT-predicted Fermi wave vectors are in excellent agreement with values obtained from quantum oscillations observed in magnetic field-dependent resistivity measurements. An electron-rich Hall coefficient is measured despite the higher hole-carrier density, attributed to the higher electron Hall mobility. The carrier mobilities are limited by surface and interface scattering, resulting in lower magnetoresistance than that measured for bulk crystals.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Correlating Experimentally Determined Hydrogen Binding Energy with Hydrogen Evolution Activity over Metal Monolayers on Molybdenum Nitride

It is well established that hydrogen binding energy (HBE) is a key descriptor for hydrogen evolution reaction (HER) activity, and such a relationship is a useful tool for searching efficient and cost-effective HER catalysts. However, in almost all cases, the HBE values are obtained from density functional theory (DFT) calculations. Here in this study, temperature programmed desorption (TPD) was used to experimentally determine the HBE values of metal monolayers supported on molybdenum nitride (Mo 2 N), and electrochemical measurements were performed on the same surfaces. Combined DFT and kinetic Monte Carlo (kMC) simulations were used to validate the trend observed with TPD and the electrochemical HER activity. Depositing one monolayer of Pt on Mo 2 N led to similar HBE values seen for bulk Pt, and electrochemical measurements showed that monolayer Pt on Mo 2 N had similar HER activity to that of bulk Pt. Similar studies were also performed for monolayer Pd on Mo 2 N to confirm the correlation. This work also demonstrates that Mo 2 N is a promising support to reduce precious metal loading in HER catalysts.

08 HYDROGEN↗
DOE OSTI2025

Revisiting point defect thermodynamics in group IVB and VB transition metal carbides

We present a comprehensive re-examination of point defect thermodynamics in group IVB and VB transition metal carbides (TMCs) with the rocksalt structure using a combination of density functional theory (DFT) calculations and a statistical mechanical Wagner-Schottky model within the canonical ensemble. The most stable configurations of point defects were discovered using basin-hopping global optimization, driven by either a machine learning interatomic potential (MLIP) or DFT. A key finding is the identification of previously unreported dicarbon antisites—a C–C dimer occupying a metal site—as the structural (constitutional) defects on the carbon-rich side of stoichiometry in all group IVB and VB TMCs except TaC. Furthermore, dicarbon antisite-containing thermal defect complexes, such as quadruple and interbranch defects, can dominate in TMCs under specific stoichiometric and temperature conditions. In conclusion, by incorporating dicarbon antisites into the defect landscape, this work provides a revised understanding of the thermodynamics of point defects in TMCs.

Carbides↗
DOE OSTI2021

Kohn-Sham Solver (KSSOLV) v2.0

KSSOLV is a MATLAB toolbox for solving Kohn-Sham density functional theory based electronic structure eigenvalue problems. It uses an object oriented features of MATLAB to represent atom, molecules, wavefunctions and Hamiltonians and their operations. It is designed to make it easier for users to prototype and test new algorithms for solving the Kohn-Sham problem. KSSOLV2.0 contains significant improvement over the original KSSOLV described in a paper published in ACM Transaction on Mathematical Software (attached). In addition to performing ground state calculation for small molecules, it can also perform geometry optimization for both molecules and solids. It uses standard pseudopotentials and implements local density approximation, generalized gradient approximation and hybrid functionals. Future releases will also include time-dependent DFT and post DFT calculations such as the GW quasi-particle energy calculation and Bethe-Salpeter equation solver for optical absorption.

Yang, Chao↗
DOE OSTI2025

Machine Learning-Driven Solvent Screening for Biobased 2,3-Butanediol Extraction

Biobased 2,3-butanediol (2,3-BDO) is a valuable biomass-derived chemical due to its versatility in being transformed into a wide variety of products. However, the separation and purification of 2,3-BDO from fermentation broth remain a significant challenge owing to its high boiling point and hydrophilic nature. Herein, we developed a machine learning (ML)-based screening workflow that uses molecular calculations as training data and requires only a small number of experimental measurements for validation to identify alternative solvent candidates for the liquid–liquid extraction (LLE) of 2,3-BDO from aqueous solution. In particular, 130 density functional theory (DFT) calculations with the implicit solvation method not only built a correlation between the computational partition coefficient and the experimental distribution coefficient of 2,3-BDO but also parameterized an Extra-Trees ML model to screen the distribution coefficient for a wider range of 6717 organic solvents. The experimental measurements of only 24 solvents were needed to validate the computational results. A list of 50 prioritized solvents was proposed for 2,3-BDO LLE, and seven additional experimental measurements were conducted to further verify our selected solvents. The impact of the extraction temperature and solvent-to-feed ratio was also investigated for selected solvents in experiments. Furthermore, this work suggested alternative solvents for 2,3-BDO LLE and proposed a versatile workflow that requires fewer experiments and can be applied to a broader range of LLE studies.

Extraction↗
DOE OSTI2020

Local Strain and Polarization Mapping in Ferrielectric Materials

CuInP 2 S 6 (CIPS) is a van der Waals material that has attracted attention because of its unusual properties. Recently, a combination of density functional theory (DFT) calculations and piezoresponse force microscopy (PFM) showed that CIPS is a uniaxial quadruple-well ferrielectric featuring two polar phases and a total of four polarization states that can be controlled by external strain. In this study, we combine DFT and PFM to investigate the stress-dependent piezoelectric properties of CIPS, which have so far remained unexplored. The two different polarization phases are predicted to differ in their mechanical properties and the stress sensitivity of their piezoelectric constants. This knowledge is applied to the interpretation of ferroelectric domain images, which enables investigation of local strain and stress distributions. The interplay of theory and experiment produces polarization maps and layer spacings which we compare to macroscopic X-ray measurements. We found that the sample contains only the low-polarization phase and that domains of one polarization orientation are strained, whereas domains of the opposite polarization direction are fully relaxed. The described nanoscale imaging methodology is applicable to any material for which the relationship between electromechanical and mechanical characteristics is known, providing insight on structural, mechanical, and electromechanical properties down to ~10 nm length scales.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Bulk and surface diffusion of neodymium in alpha-uranium: Ab initio calculations and kinetic Monte Carlo simulations

A fundamental understanding of lanthanide transport in metallic fuels is critical for high fidelity modeling of the fuel-cladding chemical interaction (FCCI) phenomenon, which can lead to the formation of brittle intermetallic compounds and premature failure of the cladding. Here we report a combined ab initio density functional theory (DFT) and kinetic Monte Carlo (KMC) study of the bulk diffusivity of Nd in α-U, fully taking into account the effect of radiation enhanced diffusion. In this work, the vacancy mechanism is considered to be the dominant mechanism for the bulk diffusion of Nd since a Nd interstitial is found to be intrinsically unstable in α-U. The surface diffusivity of a Nd adatom on α-U (001) surface has been further predicted using KMC simulations parameterized by DFT calculations. The present study suggests that Nd transport via the surface diffusion mechanism can be many orders of magnitude faster than bulk diffusion. Furthermore, the results from the present lower length scale study can be used to inform mesoscale phase-field simulations to determine the effective diffusion coefficient of Nd through α-U with a porous microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Spectroscopic Insight on Neodymium Solvation in Lithium Borohydride-Supported Electrolyte

Borohydride-based electrolytes have recently emerged as promising media for the electrodeposition of electropositive metals, including rare earth (RE) elements. While the presence of supporting alkali metal cations and RE counteranions provides essential electrochemical conductivity for achieving fast metal electrodeposition, interactions between the host ligand and solvated neodymium (Nd) complexes remain unclear. This study provides insights into the coordination structure of a concentrated and directly solvated Nd salt in a lithium borohydride-supported electrolyte. Our spectroscopic results indicate that the RE coordination environment is significantly influenced by the solvation mechanism, which can vary between metathesis and complexation pathways, primarily dictated by stoichiometric factors. Under dilute conditions, nearly complete metathesis of anions leads to a high coordination number for the host ligand (borohydride), consistent with the previously reported solvated Nd speciation in chlorine-free electrolytes. In contrast, concentrated dissolution of the Nd salt in the supported electrolyte is dominated by a complexation pathway featuring a Li-ion-paired complex with a low coordination number of the host ligand. Density functional theory (DFT) calculations indicated that the observed blue shift in the borohydride vibration was the result of an increase in electron density drawn into the terminal B–H interbond region from the hydride as the coordination changed from Li to Nd. In conjunction with DFT results, vibrational analyses allowed correlation of the experimental shifts associated with changes in Nd ligation and coordination spheres, further consolidating the prevalence of highly chloride-coordinated species under concentrated conditions. In conclusion, the outcomes of this work illuminate the distinctive and heterogeneous coordination structures that the electroactive RE species can adopt at high concentrations in lithium borohydride-supported electrolytes, as a key step to comprehend the reported metal electrodeposition performance in these media.

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

Equation of state for Hf, Ta, W, Re, Os, Ir, Pt, and Au to multi-terapascal pressures from density-functional theory

We present the zero-temperature equation of state (pressure dependence of compression) and phase stability predictions for the 5d-transition metals obtained from all-electron density-functional theory (DFT) calculations. The results compare favorably with experiments but extend beyond current experimental capabilities to 10 TPa. Our study reveals phase changes that are explained from the calculated electronic structure. The cubic face-centered and body-centered structures (fcc and bcc), together with two-, three-, and four-layered hexagonal structures, play major roles under compression. The results’ dependence on the electron exchange and correlation in the DFT approach is investigated, and it is shown that the impact of the choice, while significant at lower pressures, diminishes in the terapascal regime. We further illustrate that the normal parabolic trends in atomic volume and bulk modulus with atomic number, due to the occupation of bonding and anti-bonding 5d states, break down at TPa pressures, suggesting drastically different chemical bonding at these extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Towards Quantum-Chemical Level Calculations of SARS-CoV-2 Spike Protein Variants of Concern by First Principles Density Functional Theory

The spike protein (S-protein) is a crucial part of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), with its many domains responsible for binding, fusion, and host cell entry. In this review we use the density functional theory (DFT) calculations to analyze the atomic-scale interactions and investigate the consequences of mutations in S-protein domains. We specifically describe the key amino acids and functions of each domain, which are essential for structural stability as well as recognition and fusion processes with the host cell; in addition, we speculate on how mutations affect these properties. Such unprecedented large-scale ab initio calculations, with up to 5000 atoms in the system, are based on the novel concept of amino acid–amino acid-bond pair unit (AABPU) that allows for an alternative description of proteins, providing valuable information on partial charge, interatomic bonding and hydrogen bond (HB) formation. In general, our results show that the S-protein mutations for different variants foster an increased positive partial charge, alter the interatomic interactions, and disrupt the HB networks. We conclude by outlining a roadmap for future computational research of biomolecular virus-related systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Atomic Orbital Energy Matching vs Overlap in Actinide-Ligand Dative Bonding

Prebonding conditions and dative bond formation in actinide(IV)hexachloride complexes with U, Np, and Pu are studied by density functional theory (DFT) calculations, focusing on the interplay of atomic orbital (AO) overlap and AO energy matching in the formation of molecular orbitals (MOs) and the subsequent identification of dative bonds. The extent of donation is tracked via population analysis, MO localization, and bond-order criteria. DFT Fock matrices are used to setup models in which varying numbers of valence AOs interact to form dative bonds. The results confirm, among other effects, a contribution to metal−ligand covalency from the actinide (An) 6p shell. Better matching of An(5f) and Cl(3p) energies occurs as the An(5f) level stabilizes with increasing An effective nuclear charge. A near-degeneracy occurs in the case of the α-spin orbitals of the Pu system, but it is inconsequential. Altogether, better AO energy matching for An(5f) and Cl(3p) along the series U, Np, Pu is counter-balanced by decreasing AO overlap and decreasing availability of 5f acceptor orbitals, leading to similar bond orders and extents of donation in the three systems.

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

Kinetics of low-temperature methane activation on IrO 2 (1 1 0): Role of local surface hydroxide species

The ability of the IrO 2 (1 1 0) surface to promote CH 4 activation at low temperatures (~150 K) suggests possibilities for developing IrO 2 -based catalysts to selectively convert light alkanes to high-value chemicals. In this study, we present experimental results and microkinetic simulations, based on density functional theory (DFT) calculations, of temperature programmed reaction spectra (TPRS) obtained for CH 4 σ-complexes adsorbed on IrO 2 (1 1 0), and focus on clarifying how surface OH groups modify the branching between CH 4 desorption, activation and subsequent reactions during TPRS. Our DFT results predict that surface OH groups strongly destabilize CH 4 σ-complexes on IrO 2 (1 1 0) in addition to deactivating surface O-atoms that are needed to achieve CH 4 activation at low temperature. We demonstrate that a microkinetic model that incorporates the influence of surface OH groups on the CH 4 binding and reactivity reproduces experimental TPRS results which show that adsorbed CH 4 σ-complexes preferentially dissociate at low CH 4 coverage, but that an increasing fraction of the adsorbed CH 4 desorbs at low temperature (~120 K) with increasing initial CH 4 or H coverage. The simulations reveal that low-temperature CH4 desorption during TPRS arises primarily from CH 4 σ-complexes that are kinetically-trapped between adjacent surface OH groups, and are thus destabilized and unable to access reactive O-atoms. In addition, when the effect of adjacent OH groups is incorporated we find that the PBE-D3 functional incorporating dispersion provides CH 4 binding energies that agree more closely with that obtained from the TPRS experiments. Overall, our results demonstrate that the local effect of adjacent OH groups must be incorporated into any microkinetic models to properly capture the selectivity between extensive and partial oxidation in alkane conversion on IrO 2 (1 1 0) under reaction conditions.

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

Prediction of organic homolytic bond dissociation enthalpies at near chemical accuracy with sub-second computational cost

Bond dissociation enthalpies (BDEs) of organic molecules play a fundamental role in determining chemical reactivity and selectivity. However, BDE computations at sufficiently high levels of quantum mechanical theory require substantial computing resources. In this paper, we develop a machine learning model capable of accurately predicting BDEs for organic molecules in a fraction of a second. We perform automated density functional theory (DFT) calculations at the M06-2X/def2-TZVP level of theory for 42,577 small organic molecules, resulting in 290,664 BDEs. A graph neural network trained on a subset of these results achieves a mean absolute error of 0.58 kcal mol-1 (vs DFT) for BDEs of unseen molecules. We further demonstrate the model on two applications: first, we rapidly and accurately predict major sites of hydrogen abstraction in the metabolism of drug-like molecules, and second, we determine the dominant molecular fragmentation pathways during soot formation.

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

Pb-Apatite Framework as a Generator of Novel Flat-Band CuO-Based Physics

Based upon density functional theory (DFT) calculations, we present the basic electronic structure of CuPb 9 (PO 4 ) 6 O (Cu-doped lead apatite, aka LK-99), in two scenarios: (1) where the structure is constrained to the P3 symmetry and (2) where no symmetry is imposed. At the DFT level, the former is predicted to be metallic while the latter is found to be a charge-transfer insulator. In both cases the filling of these states is nominally d 9 , consistent with the standard Cu 2+ valence state, and Cu with a local magnetic moment of order 0.7 μ B . In the metallic case we find these states to be unusually flat (∼0.2 eV dispersion), giving a very high density of electronic states (DOS) at the Fermi level that we argue can be a host for novel electronic physics. The flatness of the bands is the likely origin of symmetry-lowering gapping possibilities that would remove the spectral weight from E F . Motivated by some initial experimental observations of metallic or semiconducting behavior, we propose that disorder (likely structural) is responsible for closing the gap. Here, we consider a variety of possibilities that could possibly close the charge-transfer gap but limit consideration to kinds of disorder that preserve electron count. Of the possible kinds we considered (spin disorder, O populating vacancy sites, and Cu on less energetically favorable Pb sites), the local Cu moment, and consequently the charge-transfer gap, remains robust. We conclude that disorder responsible for metallic behavior entails some kind of doping where the electron count changes. Further, we claim that the emergence of the flat bands should be due to weak wave function overlap between the orbitals on Cu and O sites, owing to the directional character of the constituent orbitals. Therefore, finding an appropriate host structure for minimizing hybridization between Cu and O while allowing them to still weakly interact should be a promising route for generating flat bands at E F which can lead to interesting electronic phenomena, regardless of whether LK-99 is a superconductor.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Synthesis, Electronic Structure, and Reactivity of a Planar Four‐Coordinate, Cobalt–Imido Complex

Abstract A four‐coordinate cobalt–imido complex, ( t Bu mPNP)Co=NMes ( t Bu mPNP=modified PNP pincer ligand) has been synthesized from addition of 2,4,6‐trimethylphenylazide (Mes–N 3 ) to the corresponding dinitrogen complex. The solid‐state structure determined by X‐ray diffraction established a rare, idealized planar geometry with a Co=N bond distance of 1.716(2) Å. Magnetic measurements revealed an S =1 ground state with CAS‐SCF calculations supporting radical character on the imide nitrogen. Thermolysis of the cobalt–imido compound induced selective insertion of the imido group into a Co−P bond and yielded a three‐coordinate cobalt complex with a distorted T‐shaped geometry. Transition state analysis conducted with DFT calculations established the thermodynamic stability of the P–N coupled product and provided insight into the exclusive selectivity.

Park, Yoonsu↗
DOE OSTI2025

Radiation Effects on Uranyl Tetrachloro Coordination Compounds: Impact of Lattice Water

Nuclear materials, such as uranium-bearing solids, are exposed to high levels of ionizing radiation throughout the nuclear fuel cycle; thus, it is important to develop a molecular-level understanding of how these materials behave and degrade in the presence of gamma (γ) irradiation. In the current study, three U(VI) tetrachloride complexes, M 2 [UO 2 Cl 4 ]·xH 2 O (where M = K + , Rb + , or Cs + and x = 0 or 2), and their respective chloride salts were exposed to 1–50 kGy of γ radiation using a 60 Co source. Irradiated materials were evaluated by using electron paramagnetic resonance (EPR) and Raman spectroscopy and were further explored by using density functional theory (DFT) methods. EPR spectra of the irradiated materials suggest the formation of a Cl-based radical for both the alkali salts and the uranyl tetrachloride compounds, and DFT calculations provide evidence that the Cl 2 –• radical is formed within these materials. The presence of water in the K + and Rb + compounds leads to additional spectroscopic signatures that could be traced back to water radiolysis and the formation of peroxide and superoxide species. DFT results support the formation of HO 2 • in the lattice and potentially the formation of a [UO 2 Cl 3 (O 2 )] 3– species, highlighting the impact of water within the hydrated material to alter U(VI) speciation by radiolysis.

36 MATERIALS SCIENCE↗
DOE OSTI2020

w19_OMEC Scientific Highlight [Slides]

Scientific Achievement: Revealed the effect of graphite substrate on H binding / hydrogen evolution reaction (HER) activity of organic molecular electrocatalysts (OMECs); and demonstrated the effect of molecular modification on calculated HER activity descriptor. Significance and Impact: Understanding of HER processes is important for discovering low-cost OMECs for hydrogen production and energy application; design of OMEC structures remains challenging due to the unexplored redox properties of organic molecules; and density functional theory (DFT) has been proposed to provide needed insights into HER activity for improved materials design. Research Details: DFT calculations on the binding energy of HER intermediates are used for predicting reduction potentials and and calculated H adsorption energy as activity descriptor; and LANL Institutional Computing (IC) resources are essential for successful project execution given the size of structures involved in modeling molecule-support interaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗