Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT calculation”

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 847 records · Page 47

Intricate short-range order in GeSn alloys revealed by atomistic simulations with highly accurate and efficient machine-learning potentials

GeSn alloys hold promise for silicon-compatible integrated applications in electronics, photonics, and topological quantum devices. However, understanding their intricate structures using density functional theory (DFT) calculations is hindered by spatiotemporal constraints. To overcome this limitation, we develop highly accurate and efficient machine-learning interatomic potentials based on a neuroevolution potential approach with farthest point sampling on a comprehensive DFT data set. The application of the developed machine-learning potential in large-scale atomistic simulations bridges the spatiotemporal gap between modeling and advanced characterizations, and facilitates the discovery of structural intricacies in GeSn alloys. Through extensive statistical sampling, we identify a type of short-range order (SRO) that is distinguished by both its structural signature and electronic band gap from the SRO structure previously predicted. Modeling based on a large simulation cell reveals the coexistence of nano SRO domains with various degrees of ordering, demonstrating a complex spatial heterogeneity of SRO structure. Finally, our study not only reinforces the significance of fine-level structural information in alloys, but it also constitutes an effective framework for exploring SRO in a broad range of complex alloys based on highly accurate and effective machine-learning potentials.

36 MATERIALS SCIENCE↗

Exploiting graphical processing units to enable quantum chemistry calculation of large solvated molecules with conductor-like polarizable continuum models

The conductor–like polarizable continuum model (C–PCM) with switching/Gaussian smooth discretization is a widely used implicit solvation model in quantum chemistry. We have previously implemented C–PCM solvation for Hartree–Fock (HF) and density functional theory (DFT) on graphical processing units (GPUs), enabling the quantum mechanical treatment of large solvated biomolecules. Here, we first propose a GPU–based algorithm for the PCM conjugate gradient linear solver that greatly improves the performance for very large molecules. The overhead for PCM–related evaluations now consumes less than 15% of the total runtime for DFT calculations on large molecules. Second, we demonstrate that our algorithms tailored for ground state C–PCM are transferable to excited state properties. Using a single GPU, our method evaluates the analytic gradient of the linear response PCM time–dependent density functional theory energy up to 80× faster than a conventional central processing unit (CPU)–based implementation. In addition, our C–PCM algorithms are transferable to other methods that require electrostatic potential (ESP) evaluations. For example, we achieve speed–ups of up to 130× for restricted ESP–based atomic charge evaluations, when compared to CPU–based codes. As a result, we also summarize and compare the different PCM cavity discretization schemes used in some popular quantum chemistry packages as a reference for both users and developers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unraveling the Origin of Structural Disorder in High Temperature Transition Al2O3: Structure of θ-Al2O3

The crystallography of transition Al2O3 has been extensively studied in the past due to the advantageous properties of the oxide in catalytic and a range of other technological applications. However, existing crystallographic models are insufficient to describe the structure of many important Al2O3 polymorphs due to their highly disordered nature. In this work, we investigate structure and disorder in high-temperature treated transition Al2O3, and provide a structural description for θ-Al2O3 by using a suite of complementary imaging, spectroscopy and quantum calculation techniques. Contrary to current understanding, our high-resolution imaging shows that θ-Al2O3 is a disordered composite phase of at least two different end members. By correlating imaging and spectroscopy results with DFT calculations, we propose a model that describes θ-Al2O3 as a disordered intergrowth of two crystallographic variants at the unit cell level. One variant is based on β-Ga2O3, and the other on a monoclinic phase that is closely-related to δ-Al2O3. The overall findings and interpretations afford new insight into the origin of poor crystallinity in transition Al2O3, and also provide new perspectives on structural complexity that can emerge from intergrowth of closely related structural polymorphs.

Environmental Molecular Sciences Laboratory↗

Development of a Multicenter Density Functional Tight Binding Model for Plutonium Surface Hydriding [Development of a Many-Body Density Functional Tight Binding Model for Plutonium Surface Hydriding]

We detail the creation of a multicenter density functional tight binding (DFTB) model for hydrogen on δ-plutonium, using a framework of new Slater-Koster interaction parameters and a repulsive energy based on the Chebyshev Interaction Model for Efficient Simulation (ChIMES), where two- and three-center atomic interactions are represented by linear combinations of Chebyshev polynomials. Here, we find that our DFTB/ChIMES model yields a total electron density of states for bulk δ-Pu that compares well to that from Density Functional Theory, as well as to a grid of energy calculations representing approximate H 2 dissociation paths on the δ-Pu (100) surface. We then perform molecular dynamics simulations and minimum energy pathway calculations to determine the energetics of surface dissociation and subsurface diffusion on the (100) and (111) surfaces. Our approach allows for the efficient creation of multicenter repulsive energies with a relatively small investment in initial DFT calculations. Lastly, our efforts are particularly pertinent to studies that rely on quantum calculations for interpretation and validation, such as experimental determination of chemical reactivity both on surfaces and in condensed phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Characterizing Defects Inside Hexagonal Boron Nitride Using Random Telegraph Signals in van der Waals 2D Transistors

Single-crystal hexagonal boron nitride (hBN) is used extensively in many two-dimensional electronic and quantum devices, where defects significantly impact performance. Therefore, characterizing and engineering hBN defects are crucial for advancing these technologies. Here, we examine the capture and emission dynamics of defects in hBN by utilizing low-frequency noise (LFN) spectroscopy in hBN-encapsulated and graphene-contacted MoS 2 field-effect transistors (FETs). The low disorder of this heterostructure allows the detection of random telegraph signals (RTS) in large device dimensions of 100 μm 2 at cryogenic temperatures. Analysis of gate bias- and temperature-dependent LFN data indicate that RTS originates from a single trap species within hBN. By performing multi-space density functional theory (MS-DFT) calculations on a gated defective hBN/MoS 2 heterostructure model, we assign substitutional carbon atoms in boron sites as the atomistic origin of RTS. This study demonstrates the utility of LFN spectroscopy combined with MS-DFT analysis on a low-disorder all-vdW FET as a powerful means for characterizing the atomistic defects in single-crystal hBN.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

A critical examination of compound stability predictions from machine-learned formation energies

Machine learning has emerged as a novel tool for the efficient prediction of material properties, and claims have been made that machine-learned models for the formation energy of compounds can approach the accuracy of Density Functional Theory (DFT). The models tested in this work include five recently published compositional models, a baseline model using stoichiometry alone, and a structural model. By testing seven machine learning models for formation energy on stability predictions using the Materials Project database of DFT calculations for 85,014 unique chemical compositions, we show that while formation energies can indeed be predicted well, all compositional models perform poorly on predicting the stability of compounds, making them considerably less useful than DFT for the discovery and design of new solids. Most critically, in sparse chemical spaces where few stoichiometries have stable compounds, only the structural model is capable of efficiently detecting which materials are stable. The nonincremental improvement of structural models compared with compositional models is noteworthy and encourages the use of structural models for materials discovery, with the constraint that for any new composition, the ground-state structure is not known a priori. This work demonstrates that accurate predictions of formation energy do not imply accurate predictions of stability, emphasizing the importance of assessing model performance on stability predictions, for which we provide a set of publicly available tests.

36 MATERIALS SCIENCE↗

Theoretical investigations on the structure and properties of p-n-alkoxy benzoic acid based liquid crystals

The present study focused on the structural conformations, alkoxy chain lengths and mesogenic properties of two mole of alkoxy benzoic acid(nOBA) and one mole of suberic acid (SA) hydrogen bonded (nOBASA) complexes (n=8 to 10) by density functional theory (DFT) calculations and the Fourier Transform Infrared (FT-IR) spectrum. The intermolecular hydrogen bond formation was confirmed by the optimized geometric bond lengths and bond angles obtained by computation. Using the natural bond orbital (NBO) analysis, the stability of the molecule arising from hyper conjugative interactions and charge delocalization has been analyzed. Results obtained shows that the charge in electron density (ED) in σ*and π* antibonding orbital and second order delocalization energies E(2) authorizes the occurrence of intermolecular charge transfer. The molecular electrostatic potential (MEP) surface map is plotted over the optimized geometry of the molecule to obtain the chemical reactivity of the molecule. From the local charge distributions, the mesomorphic behavior and the nematic phase stabilities for each of the molecule have been predicted. Finally the calculated result is applied to simulated infrared spectra of 8OBASA mesogens which shows good agreement with the observed spectra. The comparison of the theoretical results obtained with the experimental ones shows the reliability of this DFT method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

High-pressure characterization of Ag 3 AuTe 2 : Implications for strain-induced band tuning

Recent band structure calculations have suggested the potential for band tuning in the chiral semiconductor Ag 3 AuTe 2 to zero upon application of negative strain. In this study, we report on the synthesis of polycrystalline Ag 3 AuTe 2 and investigate its transport and optical properties and mechanical compressibility. Transport measurements reveal the semiconducting behavior of Ag 3 AuTe 2 with high resistivity and an activation energy E a of 0.2 eV. The optical bandgap determined by diffuse reflectance measurements is about three times wider than the experimental E a ⁠. Despite the difference, both experimental gaps fall within the range of predicted bandgaps by our first-principles density functional theory (DFT) calculations employing the Perdew–Burke–Ernzerhof and modified Becke–Johnson methods. Furthermore, our DFT simulations predict a progressive narrowing of the bandgap under compressive strain, with a full closure expected at a strain of –4% relative to the lattice parameter. To evaluate the feasibility of gap tunability at such substantial strain, the high-pressure behavior of Ag 3 AuTe 2 was investigated by in situ high-pressure x-ray diffraction up to 47 GPa. Mechanical compression beyond 4% resulted in a pressure-induced structural transformation, indicating the possibility of substantial gap modulation under extreme compression conditions.

36 MATERIALS SCIENCE↗

Elaborating the phases and mechanical properties of multiphase alloy: Experimental two-dimensional mapping combined with theoretical calculations

Highlights: • Mechanical modulus and hardness distribution mapping is obtained for the first time by NanoBlitz 3D. • Phases mapping and exact element ratio are characterized for the nano multielement carbides by EPMA. • Carbides types are confirmed by acid etching and XRD refinement. • Young's modulus of multiple carbides is verified by DFT calculations and Voigt-Reuss-Hill model. - Abstract: The relationship between types of phases, microstructure and mechanical properties is always a critical challenge for the research and application of multiphase alloys. We report for the first time that NanoBlitz 3D method implemented in iNano is used to investigate the distribution of micro mechanical modulus and hardness for the entire complex multiphase alloy Fe-6.0W-5.0Mo-4.0Cr-2.0V-0.83C by constructing the two-dimensional mapping. In order to build the correlation between the microstructure and properties, the carbides in this tempering state steel are confirmed to be Fe{sub 3}(W, Mo){sub 3}C, Cr{sub 15.58}Fe{sub 7.42}C{sub 6} and V{sub 8}C{sub 7} by acid etching and XRD refinement. The phases mapping is obtained by electron probe micro-analysis. Exact element ratio is characterized for the nano-scale multielement carbides precipitates. The Young's modulus and hardness of the precipitated carbides are up to 330 GPa and 22 GPa, respectively. The distribution range of Young's modulus and hardness is during 150 to 250 GPa and 2.5 to 7.0 GPa for the matrix phases, respectively. The experimental Young's modulus and hardness results are verified by first-principles calculations combined with Voigt-Reuss-Hill approximation. In fact, this work provides a successful demonstration to three-dimensionally characterize the micro mechanical properties for all the ceramics and alloys.

36 MATERIALS SCIENCE↗

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↗

Rhodium Single–Atom Catalysts on Titania for Reverse Water Gas Shift Reaction Explored by First Principles Mechanistic Analysis and Compared to Nanoclusters

The thermocatalytic reduction of CO2 by H2 often proceeds via two competing reaction mechanisms—the reverse water gas shift reaction (rWGSR, CO2 + H2 ⇌ CO + H2O) and methanation (CO2 + 4H2 ⇌ CH4 + 2H2O). Atomically dispersed Rh1 catalysts on TiO2 show high selectivity toward the rWGSR compared with larger Rh nanoclusters, but the origin of this size-dependent selectivity has not been fully explained. Here we report density functional theory (DFT) calculations and microkinetic simulations that clarify the Rh1 active sites and rWGSR pathway on anatase TiO2(101), as well as the high rWGSR selectivity of Rh1 compared with supported Rhx (x = 2–8 atoms) nanoclusters. DFT-computed formation energies, vibrational frequency analysis, and microkinetic modeling suggest three plausible active sites: Rh1 on titania (Rh1/TiO2(101)), Rh1 with a nearby hydroxyl group (Rh1OH/TiO2(101)), and Rh1 near an oxygen vacancy at a three-fold coordinated site (Rh1 near O3cvac). Predicted turnover frequencies and apparent activation barriers for Rh1 indicate a faster reaction involving CO2 dissociation assisted by a support oxygen vacancy via Rh1 near O3cvac, as well as slower reactions involving Rh1OH/TiO2(101) or Rh1/TiO2(101) through a COOH intermediate. Furthermore, these Rh1 sites are selective toward CO rather than CH4 because of the weak adsorption of CO, large barrier for C-O bond dissociation, and the lack of nearby metal sites for H2 dissociation, in contrast to Rhx nanoclusters, including Rh2 dimers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of Oxygen Vacancy on the Magnetic Properties of Ni-Doped SnO{sub 2} Nanoparticles

We studied the ferromagnetism of Ni-doped SnO{sub 2} with and without oxygen vacancy (V{sub O}) by experiments and calculations. Sn{sub 0.96}Ni{sub 0.04}O{sub 2} nanocrystalline powder prepared by the sol-gel method exhibits room-temperature ferromagnetism, and the vacuum annealing reduces its saturation magnetizations. The density functional theory (DFT) calculations show that the doped Ni atom can introduce a local magnetic moment of about 2.0 μ{sub B} in SnO{sub 2} (110) surface, which is ascribed to the spin-polarization of Ni 3d and O 2p electrons. The type of preferential magnetic coupling of Ni-doped SnO{sub 2} (110) surface depends upon the distributions of Ni dopants, and ferromagnetic coupling is more favorable in most cases. Introduction of V{sub O} obviously weakens the ferromagnetic interaction of Ni-doped SnO{sub 2} (110) surface, which can give a reasonable explanation for our experimental observations.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Tight-binding approximations to time-dependent density functional theory — A fast approach for the calculation of electronically excited states

We propose a new method of calculating electronically excited states that combines a density functional theory based ground state calculation with a linear response treatment that employs approximations used in the time-dependent density functional based tight binding (TD-DFTB) approach. The new method termed time-dependent density functional theory TD-DFT+TB does not rely on the DFTB parametrization and is therefore applicable to systems involving all combinations of elements. We show that the new method yields UV/Vis absorption spectra that are in excellent agreement with computationally much more expensive TD-DFT calculations. Errors in vertical excitation energies are reduced by a factor of two compared to TD-DFTB.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Mechanistic Insights into the Structure-Dependent Selectivity of Catalytic Furfural Conversion on Platinum Catalysts

The effects of structure and size on the selectivity of catalytic furfural conversion over supported Pt catalysts in the presence of hydrogen have been studied using first principles density functional theory (DFT) calculations and microkinetic modeling. Four Pt model systems, i.e., periodic Pt(111), Pt(211) surfaces, as well as small nanoclusters (Pt13 and Pt55) are chosen to represent the terrace, step, and corner sites of Pt nanoparticles. Our DFT results show that the reaction routes for furfural hydrogenation and decarbonylation are strongly dependent on the type of reactive sites, which lead to the different selectivity. On the basis of the size-dependent site distribution rule, we correlate the site distributions as a function of the Pt particle size. Our microkinetic results indicate the critical particle size that controls the furfural selectivity is about 1.0 nm, which is in good agreement with the reported experimental value under reaction conditions. This work was supported by National Basic Research Program of China (973 Program) (2013CB733501) and the National Natural Science Foundation of China (NSFC-21306169, 21176221, 21136001, 21101137 and 91334103). This work was also partially supported by the US Department of Energy (DOE), the Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences & Biosciences. Pacific Northwest National Laboratory (PNNL) is a multiprogram national laboratory operated for DOE by Battelle. Computing time was granted by the grand challenge of computational catalysis of the William R. Wiley Environmental Molecular Sciences Laboratory (EMSL). EMSL is a national scientific user facility located at Pacific Northwest National Laboratory (PNNL) and sponsored by DOE’s Office of Biological and Environmental Research.

biomass conversion↗

Influence of La/W ratio on electrical conductivity of lanthanum tungstate with high La/W ratio

The proton-conducting properties of lanthanum tungstates (LWOs) with high La/W ratios were investigated using electrochemical measurements and quantum chemical calculations. Single phases of LWOs with high La/W ratios (6.3≤La/W≤6.7) were synthesized by high-temperature sintering at around 1700 °C. The electrical conductivity of LWO increased with increasing La/W ratio in the single-phase region. The LWO synthesized at the optimum sintering temperature and time, and with the optimum La/W ratio gave the maximum conductivity, i.e., 2.7×10{sup −3} S cm{sup −1} with La/W=6.7 at 500 °C. Density functional theory calculations, using the nudged elastic band method, were performed to investigate the proton diffusion barrier. The results suggest that the proton diffusion paths around La sites have the lowest proton diffusion barrier. These findings improve our understanding of LWO synthesis and the proton-conducting mechanism and provide a strategy for improving proton conduction in LWOs. - Graphical abstract: The LWOs with high La/W ratios were synthesized for the first time. The optimum La/W ratio gave the maximum conductivity with La/W=6.7 at 500 °C. The proton diffusion paths were also considered with density functional theory calculations. - Highlights: • The proton-conducting properties of lanthanum tungstates (LWOs) were investigated. • Single phase LWOs with high La/W ratios (6.3≤La/W≤6.7) were synthesized successfully. • LWOs with the high La/W ratios showed high proton conductivity. • The DFT calculation suggested the lowest proton diffusion barrier in the path around La sites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Formic Acid Dissociative Adsorption on NiO(111): Energetics and Structure of Adsorbed Formate

Here, the dissociative adsorption of carboxylic acids on oxide surfaces is important for understanding adsorbed carboxylates, which are important as intermediates in catalytic reactions, for the organo-functionalization of oxide surfaces, and in many other aspects of oxide surface chemistry. We present here the first direct experimental measurement of the heat of dissociative adsorption of any carboxylic acid on any single-crystal oxide surface. The enthalpy of the dissociative adsorption of formic acid, the simplest carboxylic acid, to produce adsorbed formate and hydrogen (as a surface hydroxyl) on a (2 × 2)-NiO(111) surface is measured by single crystal adsorption calorimetry. The differential heat of adsorption decreases with formic acid coverage from 202 to 99 kJ/mol at saturation (0.25 ML). The structure of the adsorbed products is clarified by density functional theory (DFT) calculations, which provide energies in reasonable agreement with the calorimetry. These calculations show that formic acid readily dissociates on both the oxygen and Ni terminations of the octapolar NiO(111) surfaces, donating its acid H to a surface lattice oxygen, while HCOO adsorbs preferentially with bridging-type geometry near the M-O 3 /O-M 3 sites. The calculated energetics at low coverages agrees well with experimental data, while larger differences are observed at high coverage (0.25 ML). The large decrease in experimental heat of adsorption with coverage can be brought into agreement with the DFT energies if we assume that both types of octapolar surface terminations (O- and Ni-) are present on the starting surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗