Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional theory (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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Minimum conditions for accurate modeling of urea production via co-electrolysis

Co-electrolysis of carbon oxides and nitrogen oxides promise to simultaneously help restore the balance of the C and N cycles while producing valuable chemicals such as urea. However, co-electrolysis processes are still largely inefficient and numerous knowledge voids persist. Here, we provide a solid thermodynamic basis for modelling urea production via co-electrolysis. First, we determine the energetics of aqueous urea produced under electrochemical conditions based on experimental data, which enables an accurate assessment of equilibrium potentials and overpotentials. Next, we use density functional theory (DFT) calculations to model various co-electrolysis reactions producing urea. The calculated reaction free energies deviate significantly from experimental values for well-known GGA, meta-GGA and hybrid functionals. These deviations stem from errors in the DFT-calculated energies of molecular reactants and products. In particular, the error for urea is approximately -0.25 ± 0.10 eV. Finally, we show that all these errors introduce large inconsistencies in the calculated free-energy diagrams of urea production via co-electrolysis, such that gas-phase corrections are strongly advised.

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

Half‐sandwich ruthenium complex with a very low overpotential and excellent activity for water oxidation under acidic conditions

Abstract Molecular catalysts are acknowledged for the ability to design reaction sites within well‐defined structures to achieve high catalytic activities. However, in many cases, molecular catalysts undergo structural changes to some other form(s), which are finally the real catalysts. Here, we report two half‐sandwich ruthenium complexes; [Ru([9]aneN3)(bpy)Cl] + 3 and [Ru([9]aneN3)(pic)OH 2 ] 2+ 4 , based on the macrocyclic 1,4,7‐triazacyclononane ([9]aneN3) ligand for water oxidation (WO). The two complexes have similar core but different ancillary ligands, which greatly affected their stability as well as activity for WO. Complex 3 retained high stability and excellent activity (turnover number [TON] = 1250) in chemical WO and first‐order reaction kinetics with respect to [Ce IV ] with a calculated rate constant ( k cat ) of 34.59 s −1 . Further, the complex demonstrated very low overpotential of ~210 mV in electrochemical WO. At an overpotential of only 400 mV, turnover frequency (TOF) of complex 3 was electrochemically estimated to be 131.2 s −1 . In contrast, complex 4 underwent picolinate ligand dissociation, as a deactivation pathway, to form the tri‐aqua derivative. Density functional theory (DFT) calculations are used to explain the dissociation mechanism of picolinate ligand in complex 4 , which happens through a stepwise dissociation mechanism.

Younus, Hussein A.↗
DOE OSTI2022

Energy and mechanical properties predictions in Fe-Ni binary system by ab initio calculations

Here, the fcc random solid solutions in Fe-Ni binary system were investigated using ab initio calculations. Density functional theory (DFT) was applied along with supercells generated by the Special Quasi-Random Structure (SQS) method to predict the enthalpy of mixing and magnetic and elastic properties. The results were compared with relevant experimental and computational data produced by other groups. The effect of SQS supercell size and relaxation procedure on the predicted properties of the alloys were also studied. Our calculations show that 16-atoms SQS model provides adequate approach to the mixing energy of the alloy. However, calculation of mechanical properties requires larger supercell. Excellent agreement between experimental and calculated elastic properties was found in fcc region of Fe-Ni system.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Spatial Control of Substitutional Dopants in Hexagonal Monolayer WS 2 : The Effect of Edge Termination

The ability to control the density and spatial distribution of substitutional dopants in semiconductors is crucial for achieving desired physicochemical properties. Substitutional doping with adjustable doping levels has been previously demonstrated in 2D transition metal dichalcogenides (TMDs); however, the spatial control of dopant distribution remains an open field. In this work, edge termination is demonstrated as an important characteristic of 2D TMD monocrystals that affects the distribution of substitutional dopants. Particularly, in chemical vapor deposition (CVD)-grown monolayer WS 2 , it is found that a higher density of transition metal dopants is always incorporated in sulfur-terminated domains when compared to tungsten-terminated domains. Here two representative examples demonstrate this spatial distribution control, including hexagonal iron- and vanadium-doped WS 2 monolayers. Density functional theory (DFT) calculations are further performed, indicating that the edge-dependent dopant distribution is due to a strong binding of tungsten atoms at tungsten-zigzag edges, resulting in the formation of open sites at sulfur-zigzag edges that enable preferential dopant incorporation. Based on these results, it is envisioned that edge termination in crystalline TMD monolayers can be utilized as a novel and effective knob for engineering the spatial distribution of substitutional dopants, leading to in-plane hetero-/multi-junctions that display fascinating electronic, optoelectronic, and magnetic properties.

2D transition metal dichalcogenides↗
DOE OSTI2023

The Intermetallic Reactivity Database: Compiling Chemical Pressure and Electronic Metrics toward Materials Design and Discovery

Here, the advent of high-throughput Density Functional Theory (DFT) calculations has supported the creation of large databases containing the quantitative output necessary for constructing theoretical phase diagrams and predicting physical properties. In this Article, we present a complementary resource, the Intermetallic Reactivity Database (IRD), focused on the chemical bonding features of solid-state structures and indicators of potential structural transformations. Each IRD entry augments common features, such as band structures and density of states (DOS) distributions, with chemically motivated information including DFT-Chemical Pressure (CP) schemes and visualizable representations of the atomic charges. Together, these data types enable the rationalization and prediction of potential structural phenomena encountered in intermetallic chemistry, as we illustrate with four examples: the origins of the Y 2 Ni 2 Mg structure in terms of CP features of its parent structures, the anticipation of intergrowth phases from the net atomic CPs collected in Al-containing binary phases, the correlation between trends in the CP schemes of CaCu 5 -type phases and experimentally observed structural variations, and finally, the development of theoretical methodology with the testing of a streamlined method generating DFT-CP schemes. Altogether, these examples highlight how the IRD supports the creation of models of structural chemistry that extend beyond the bounds of its entries.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multiphase tin equation of state using density functional theory

In this work, we perform density functional theory (DFT) calculations of five solid phases and the liquid phase of tin. The calculations include cold curves of the five solid phases, phonon calculations in the quasiharmonic approximation over a range of volumes for each solid phase, and DFT-based molecular dynamics (DFT-MD) simulations of the liquid phase, including those of the melt curve using the Z method. Using the DFT results, we construct a tabular multiphase sesame equation of state for tin, referred to as sesame 2162. Comparisons to experimental data are made and show a high level of agreement in isobaric data, isothermal data, shock data, and phase boundary measurements, including measurements of the melt curve. The 2162 EOS will be useful for hydrodynamics simulations and has been designed with an eye toward hydrodynamics simulations that incorporate materials strength models and allow for modeling of the kinetics of phase transitions.

3-dimensional systems↗
DOE OSTI2024

Exploring the Computational Aspects of Propylene Oligomerization Catalysis Using M′ 2 M Type Trimetallic MOF Nodes

Metal-organic frameworks have emerged as promising materials in the field of catalysis. They offer an optimal ground for screening catalysts and tailoring their catalytic properties. Here, in this work, via density functional theory (DFT) calculations, we investigated the catalytic activity of the trimetallic MOF nodes, M' 2 M for propylene oligomerization, by varying the active metal M from Sc to Cu with M' being Fe, aiming to grasp the impact of altering the active atoms on the catalyst's activity. Additionally, we examined how substituting the spectator atom, M', with other transition metals, i.e., from Sc to Cu, affects these energy barriers, keeping Ni as the active metal. We proposed several cases with lower or comparable energy barriers to the experimentally reported Fe 2 Ni trimetallic MOF node. In addition, we found a correlative relationship between spin-density from natural population analysis and energy barriers in the realm of C-C bond formation, whereby an elevation in spin-density is found to be inversely proportional to the magnitude of the energy barriers. Moreover, we calculated the energy barriers for C-C coupling and beta-hydride elimination using multireference NEVPT2 calculations on top of the CASSCF wave function to validate the rate-determining step that is predicted by DFT.

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

Towards extreme fast charging of 4.6V LiCoO 2 via mitigating high-voltage kinetic hindrance

High-voltage LiCoO 2 (LCO) is an attractive cathode for ultra-high energy density lithium-ion batteries (LIBs) in the 3C markets. However, the sluggish lithium-ion diffusion at high voltage significantly hampers its rate capability. Herein, combining experiments with density functional theory (DFT) calculations, we demonstrate that the kinetic limitations can be mitigated by a facial Mg 2+ +Gd 3+ co-doping method. The as-prepared LCO shows significantly enhanced Li-ion diffusion mobility at high voltage, making more homogenous Li-ion de/intercalation at a high-rate charge/discharge process. The homogeneity enables the structural stability of LCO at a high-rate current density, inhibiting stress accumulation and irreversible phase transition. When used in combination with a Li metal anode, the doped LCO shows an extreme fast charging (XFC) capability, with a superior high capacity of 193.1 mAhg -1 even at the current density of 20C and high-rate capacity retention of 91.3% after 100 cycles at 5C. In conclusion, this work provides a new insight to prepare XFC high-voltage LCO cathode materials.

25 ENERGY STORAGE↗
DOE OSTI2022

Machine learning for impurity charge-state transition levels in semiconductors from elemental properties using multi-fidelity datasets

Quantifying charge-state transition energy levels of impurities in semiconductors is critical to understanding and engineering their optoelectronic properties for applications ranging from solar photovoltaics to infrared lasers. While these transition levels can be measured and calculated accurately, such efforts are time-consuming and more rapid prediction methods would be beneficial. Here, we significantly reduce the time typically required to predict impurity transition levels using multi-fidelity datasets and a machine learning approach employing features based on elemental properties and impurity positions. We use transition levels obtained from low-fidelity (i.e., local-density approximation or generalized gradient approximation) density functional theory (DFT) calculations, corrected using a recently proposed modified band alignment scheme, which well-approximates transition levels from high-fidelity DFT (i.e., hybrid HSE06). Further, the model fit to the large multi-fidelity database shows improved accuracy compared to the models trained on the more limited high-fidelity values. Crucially, in our approach, when using the multi-fidelity data, high-fidelity values are not required for model training, significantly reducing the computational cost required for training the model. Our machine learning model of transition levels has a root mean squared (mean absolute) error of 0.36 (0.27) eV vs high-fidelity hybrid functional values when averaged over 14 semiconductor systems from the II–VI and III–V families. As a guide for use on other systems, we assessed the model on simulated data to show the expected accuracy level as a function of bandgap for new materials of interest. Finally, we use the model to predict a complete space of impurity charge-state transition levels in all zinc blende III–V and II–VI systems.

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

Thermal fluctuations and carrier localization induced by dynamic disorder in MAPbI 3 described by first-principles based tight-binding model

Halide perovskites are strongly influenced by large amplitude anharmonic lattice fluctuations at room temperature. We develop a tight binding model for dynamically disordered MAPbI$_3$ based on density functional theory (DFT) calculations to calculate electronic structure for finite temperature crystal structures at the length scale of thermal disorder and carrier localization. The model predicts individual Hamiltonian matrix elements and band structures with high accuracy, owing to the inclusion of additional matrix elements and descriptors for non-Coulombic interactions. We apply this model to electronic structure at length and time scales inaccessible to first principles methods, finding an increase in band gap, carrier mass, and the sub-picosecond fluctuations in these quantities with increasing temperature as well as the onset of carrier localization in large supercells induced by thermal disorder at 300 K. Finally, we identify the length scale $L^*= 5$ nm as the onset of localization in the electronic structure, associated with associated with decreasing band edge fluctuations, increasing carrier mass, and Rashba splitting approaching zero.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Does Chirality Influence the Stability of Amino Acid –Cu Complexes in the Salt-Induced Peptide Formation Reaction? Insights from Density Functional Theory Calculations

The polymerization of amino acids into peptides and ultimately proteins is critical to evolution of life. Ribosomally-synthesized proteins are homochiral, composed entirely of L-amino acids. Abiotic syntheses of amino acids, however, produce racemic mixtures(equal amounts of L and D enantiomers), meaning that the pool of monomers available for abiotic protein synthesis would be racemic or nearly so. There are no prebiotically plausible mechanisms known to yield homochiral peptides from racemic pools of amino acids. While some mechanisms capable of producing or amplifying enantiomeric excesses have been observed, these excesses are small or not applicable to all biologically relevant amino acids. As a result, a prebiotically-feasible mechanism to create enantiopure pools of amino acids has not been identified. However, it is not clear that these enantio-enriched or enantiopure pools are necessary to produce homochiral peptides. Steric and electronic effects, or a combination of both, may influence how different enantiomers interact with each other, potentially leading to homochiral peptides from racemic mixtures. Here, we use density functional theory (DFT) calculations to assess the stability of homochiral and heterochiral reactive complexes produced during polymerization via the salt-induced peptide formation (SIPF) reaction. Experimental studies of the SIPF reaction have shown that it enables the formation of peptides under diverse environmental conditions, making it a plausible pathway to peptide formation on early Earth. In the SIPF reaction, NaCl acts as a condensation reagent while a divalent metal cation, primarily Cu2+, facilitates polymerization by complexing and activating amino acids forming a monochlorocuprate complex. To assess if hetero-or homochiral monochlorocuprate complexes were energetically favored, we compared the stability of LL, DD, and LD enantiomers of Cu2+–(alanine)2and Cu2+–(valine) 2 complexes. Gaussian 09 density DFT energy-minimization calculations were made for each Cu –amino acid complex in the cis and trans configuration. Models were energy minimized using the BH and HLYP/6-31++G(d,p) level of theory and free energies were compared to determine the most stable configuration.

A C Fox↗
DOE OSTI2021

Electrochemical CO 2 Reduction Reaction over Cu Nanoparticles with Tunable Activity and Selectivity Mediated by Functional Groups in Polymeric Binder

Electrochemical carbon dioxide reduction reaction (CO 2 RR) using copper (Cu)-based catalysts has received significant attention mainly because Cu is an element capable of producing hydrocarbons and oxygenates. One possible way to control the CO 2 RR performance at the electrode interface is by modifying catalysts with specific functional groups of different polymeric binders, which are necessary components in the process of electrode fabrication. However, the modification effect of the key functional groups on the CO 2 RR activity and selectivity is poorly understood over Cu-based catalysts. In this work, the role of functional groups (e.g., -COOH and -CF 2 groups) in hydrophilic and hydrophobic polymeric binders on the CO 2 RR of Cu-based catalysts is investigated using a combination of electrochemical measurements, in-situ characterization and density functional theory (DFT) calculations. DFT results reveal that functional groups influence the binding energies of key intermediates involved in both CO 2 RR and the competing hydrogen evolution reaction, consistent with experimental observation of binder-dependent product distributions among formic acid, CO, CH 4 , and H 2 . This study provides a fundamental understanding that the selection of desired polymeric binders is a useful strategy for tuning the CO 2 RR activity and selectivity.

25 ENERGY STORAGE↗
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 OSTI2021

Crystal Structure of Colloidally Prepared Metastable Ag 2 Se Nanocrystals

Structural polymorphism is known for many bulk materials; however, on the nanoscale metastable polymorphs tend to form more readily than in the bulk, and with more structural variety. One such metastable polymorph observed for colloidal Ag 2 Se nanocrystals has traditionally been referred to as the “tetragonal” phase. While there are reports on the chemistry and properties of this metastable polymorph, its crystal structure, and therefore electronic structure, has yet to be determined. We report that an anti-PbCl 2 -like structure type (space group P2 1 /n) more accurately describes the powder X-ray diffraction and X-ray total scattering patterns of colloidal Ag 2 Se nanocrystals prepared by several different methods. Density functional theory (DFT) calculations indicate that this anti-PbCl 2 -like Ag 2 Se polymorph is a dynamically stable, narrow-band-gap semiconductor. Here, the anti-PbCl 2 -like structure of Ag 2 Se is a low-lying metastable polymorph at 5–25 meV/atom above the ground state, depending on the exchange-correlation functional used.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Molecular chemisorption of N 2 on IrO 2 (110)

We investigated the adsorption of N 2 on stoichiometric and O-rich IrO 2 (110) surfaces using temperature programmed desorption (TPD) experiments and density functional theory (DFT) calculations. TPD shows that N 2 desorbs predominantly from the stoichiometric-IrO 2 (110) surface in a well-defined peak at 270 K for N 2 coverages below about 0.5 ML, and that a shoulder centered near 235 K develops in the N 2 TPD traces as the coverage approaches saturation, indicating that adsorbed N 2 molecules become destabilized at high N 2 coverage. Experiments of N 2 adsorption onto O-rich IrO 2 (110) surfaces provide evidence that N 2 adsorbs exclusively on the coordinatively-unsaturated Ir atoms (Ircus) of the surface, and that pre-adsorbed O-atoms (“on-top” oxygen) stabilize adsorbed N 2 molecules, causing the main N 2 TPD peak to shift toward higher temperature with increasing oxygen coverage. Consistent with prior results, our DFT calculations predict that an N 2 molecule preferentially adsorbs into an upright configuration on an Ircus atom of the IrO 2 (110) surface and achieves a binding energy of about 100 kJ/mol. The computed binding energy agrees well with our experimental estimate of ~90 kJ/mol for low N 2 coverages on stoichiometric IrO 2 (110). Lastly, the DFT calculations also quantitatively reproduce the observed stabilization of N 2 by co-adsorbed on-top O-atoms and predict the destabilization of N 2 on IrO 2 (110) as the N 2 adlayer becomes crowded at high coverage.

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

Computational Design of Miniproteins as SARS-CoV-2 Therapeutic Inhibitors

A rational therapeutic strategy is urgently needed for combating SARS-CoV-2 infection. Viral infection initiates when the SARS-CoV-2 receptor-binding domain (RBD) binds to the ACE2 receptor, and thus, inhibiting RBD is a promising therapeutic for blocking viral entry. In this study, the structure of lead antiviral candidate binder (LCB1), which has three alpha-helices (H1, H2, and H3), is used as a template to design and simulate several miniprotein RBD inhibitors. LCB1 undergoes two modifications: structural modification by truncation of the H3 to reduce its size, followed by single and double amino acid substitutions to enhance its binding with RBD. We use molecular dynamics (MD) simulations supported by ab initio density functional theory (DFT) calculations. Complete binding profiles of all miniproteins with RBD have been determined. The MD investigations reveal that the H3 truncation results in a small inhibitor with a -1.5 kcal/mol tighter binding to RBD than original LCB1, while the best miniprotein with higher binding affinity involves D17R or E11V + D17R mutation. DFT calculations provide atomic-scale details on the role of hydrogen bonding and partial charge distribution in stabilizing the minibinder:RBD complex. This study provides insights into general principles for designing potential therapeutics for SARS-CoV-2.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

High Thermoelectric Performance in Chalcopyrite Cu 1-x Ag x GaTe 2 –ZnTe: Nontrivial Band Structure and Dynamic Doping Effect

The understanding of thermoelectric properties of ternary I–III–VI 2 type (I = Cu, Ag; III = Ga, In; and VI = Te) chalcopyrites is less well developed. Although their thermal transport properties are relatively well studied, the relationship between the electronic band structure and charge transport properties of chalcopyrites has been rarely discussed. In this study, we reveal the unusual electronic band structure and the dynamic doping effect that could underpin the promising thermoelectric properties of Cu 1–x Ag x GaTe 2 compounds. Density functional theory (DFT) calculations and electronic transport measurements suggest that the Cu 1–x Ag x GaTe 2 compounds possess an unusual non-parabolic band structure, which is important for obtaining a high Seebeck coefficient. Moreover, a mid-gap impurity level was also observed in Cu 1–x Ag x GaTe 2 , which leads to a strong temperature-dependent carrier concentration and is able to regulate the carrier density at the optimized value for a wide temperature region and thus is beneficial to obtaining the high power factor and high average ZT of Cu 1–x Ag x GaTe 2 compounds. We also demonstrate a great improvement in the thermoelectric performance of Cu 1–x Ag x GaTe 2 by introducing Cu vacancies and ZnTe alloying. The Cu vacancies are effective in increasing the hole density and the electrical conductivity, while ZnTe alloying reduces the thermal conductivity. As a result, a maximum ZT of 1.43 at 850 K and a record-high average ZT of 0.81 for the Cu 0.68 Ag 0.3 GaTe 2 –0.5%ZnTe compound are achieved.

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

Simulations and Experiments Reveal Effect of Nanopores on Helium Diffusion in Quartz

The diffusion properties of noble gases in minerals are widely used to reconstruct the thermal histories of rocks. In this work, we combine density functional theory (DFT) calculations with laboratory experiments to investigate controls on helium diffusion in quartz. DFT calculations for perfect α-quartz predict substantially lower activation energies and frequency factors for helium diffusion than observed in laboratory experiments, especially in the [001] direction. These results imply that no helium could be retained in quartz at Earth surface temperatures, which conflicts with observations of partial cosmogenic 3 He retention over geologic time scales. Here, we implement a model of helium diffusion in α-quartz modulated by nanopore defects that disrupt energetically favorable diffusion pathways. In this model, we find that laboratory-determined diffusivities can be most closely reproduced when a helium atom encounters ~70 nanopore sites per million interstitial sites. The results of our model indicate that diffusion of helium in natural quartz, like other noble gases in other minerals, can be significantly modulated by extended defects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗