Search NASA⌕ Search

SEARCH · Search NASA

Results for “optimized energy density functionals”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Accurate and balanced anisotropic Gaussian type orbital basis sets for atoms in strong magnetic fields

In high magnetic field calculations, anisotropic Gaussian type orbital (AGTO) basis functions are capable of reconciling the competing demands of the spherically symmetric Coulombic interaction and cylindrical magnetic (B field) confinement. However, the best available a priori procedure for composing highly accurate AGTO sets for atoms in a strong B field [W. Zhu et al., Phys. Rev. A 90, 022504 (2014)] yields very large basis sets. Their size is problematical for use in any calculation with unfavorable computational cost scaling. Here we provide an alternative constructive procedure. It is based upon analysis of the underlying physics of atoms in B fields that allow identification of several principles for the construction of AGTO basis sets. Aided by numerical optimization and parameter fitting, followed by fine tuning of fitting parameters, we devise formulae for generating accurate AGTO basis sets in an arbitrary B field. For the hydrogen iso-electronic sequence, a set depends on B field strength, nuclear charge, and orbital quantum numbers. For multi-electron systems, the basis set formulae also include adjustment to account for orbital occupations. Tests of the new basis sets for atoms H through C (1 ≤ Z ≤ 6) and ions Li+, Be+, and B+, in a wide B field range (0 ≤ B ≤ 2000 a.u.), show an accuracy better than a few μhartree for single-electron systems and a few hundredths to a few mHs for multi-electron atoms. The relative errors are similar for different atoms and ions in a large B field range, from a few to a couple of tens of millionths, thereby confirming rather uniform accuracy across the nuclear charge Z and B field strength values. Residual basis set errors are two to three orders of magnitude smaller than the electronic correlation energies in multi-electron atoms, a signal of the usefulness of the new AGTO basis sets in correlated wavefunction or density functional calculations for atomic and molecular systems in an external strong B field

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Gaseous cyclodextrin- closo -dodecaborate complexes χCD·B 12 X 12 2– (χ = α, β, and γ; X = F, Cl, Br, and I): electronic structures and intramolecular interactions

A fundamental understanding of cyclodextrin- closo -dodecaborate inclusion complexes is of great interest in supramolecular chemistry. Herein, we report a systematic investigation on the electronic structures and intramolecular interactions of perhalogenated closo -dodecaborate dianions B 12 X 12 2– (X = F, Cl, Br and I) binding to α-, β-, and γ-cyclodextrins (CDs) in the gas phase using combined negative ion photoelectron spectroscopy (NIPES) and density functional theory (DFT) calculations. The vertical detachment energy (VDE) of each complex and electronic stabilization of each dianion due to the CD binding (ΔVDE, relative to the corresponding isolated B 12 X 12 2– ) are determined from the experiments along α-, β- and γ-CD in the form of VDE (ΔVDE): 4.00 (2.10), 4.33 (2.43), and 4.30 (2.40) eV in X = F; 4.09 (1.14), 4.64 (1.69), and 4.69 (1.74) eV in X = Cl; 4.11 (0.91), 4.58 (1.38), and 4.70 (1.50) eV in X = Br; and 3.54 (0.74), 3.88 (1.08), and 4.05 (1.25) eV in X = I, respectively. All complexes have significantly higher VDEs than the corresponding isolated dodecaborate dianions with ΔVDE spanning from 0.74 eV at (α, I) to 2.43 eV at (β, F), sensitive to both host CD size and guest substituent X. DFT-optimized complex structures indicate that all B 12 X 12 2– prefer binding to the wide openings of CDs with the insertion depth and binding motif strongly dependent on the CD size and halogen X. Dodecaborate anions with heavy halogens, i.e. , X = Cl, Br, and I, are found outside of α-CD, while B 12 F 12 2– is completely wrapped by γ-CD. Partial embedment of B 12 X 12 2– into CDs is observed for the other complexes via multipronged B–X … H–O/C interlocking patterns. The simulated spectra based on the density of states agree well with those of the experiments and the calculated VDEs well reproduce the experimental trends. Molecular orbital analyses suggest that the spectral features at low binding energies originated from electrons detached from the dodecaborate dianion, while those at higher binding energies are derived from electron detachment from CDs. Energy decomposition analyses reveal that the electrostatic interaction plays a dominating role in contributing to the host–guest interactions for the X = F series partially due to the formation of a O/C–H … X–B hydrogen bonding network, and the dispersion forces gradually become important with the increase of halogen size.

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

Modified bornite materials with high electrochemical performance for sodium and lithium storage

Mineral-based functional materials have attracted enormous attention due to the advantages of high production, good material consistency, environmental friendliness, and low cost. However, the efficient utilization of natural minerals for high-performance energy storage is rarely reported. Here, we demonstrate a mineral-derived anode material Cu 5 FeS 4 coated with carbon for both sodium-ion and lithium-ion batteries. The findings show that ball milling bornite in the presence of polyvinylpyrrolidone followed by optimal-temperature carbonization results in S, N-codoped Cu 5 FeS 4 particles with the desired surface properties, enhanced electronic conductivity, and abundant active sites that deliver high reversible capacity and maintain superior long-term cycling performance at a high current density of 2.0 A g –1 . This work sheds light on potential energy storage use cases for naturally abundant mineral-based materials.

25 ENERGY STORAGE↗
DOE OSTI2018

Computational Study of Microhydration in Sulfonated Diels–Alder Poly(phenylene) Polymers

In this paper, the nature of microhydration in sulfonated Diels–Alder poly(phenylene) (SDAPP) polymer membranes is explored using ab initio and density functional theory (DFT) electronic structure calculations. The impact of the aromatic poly(phenylene) structure, including cooperative effects between multiple spatially adjacent sulfonic groups, on the hydration environment is addressed using a series of DFT B3LYP/6-311**-optimized structures for different SDAPP·nH 2 O clusters. In addition, larger SDAPP polymer fragments, along with selected hydrophilic domain structures extracted from molecular dynamic (MD) simulations, are also evaluated using ONIOM HF/PM6 semiempirical calculations. The SDAPP clusters reveal that spontaneous proton dissociation occurs at low levels of hydration to form sulfonic-acid-associated H 3 O + contact ion pairs (CIPs), which then evolve into solvated CIPs at higher hydration levels. For multiple sulfonic acid groups located on the poly(phenylene) side chains, the hydration energies are a function of the relative acid location and backbone configuration. Variations in the phenylene backbone torsional angles allow remote sulfonic acids to adopt an optimal separation to produce an extended hydrogen bonded network of waters between the SDAPP acids groups. Finally, these calculations provide a baseline to help describe the proton transport and hydration behavior of SDAPP membranes.

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

Fluorine-Tuned Carbon-Based Nickel Single-Atom Catalysts for Scalable and Highly Efficient CO 2 Electrocatalytic Reduction

Electrocatalytic CO 2 reduction is garnering significant interest due to its potential applications in mitigating CO 2 and producing fuel. However, the scaling up of related catalysis is still hindered by several challenges, including the cost of the catalytic materials, low selectivity, small current densities to maintain desirable selectivity. In this study, Fluorine (F) atoms were introduced into an N-doped carbon-supported single nickel (Ni) atom catalyst via facile polymer-assisted pyrolysis. This method not only maintains the high atom utilization efficiency of Ni in a cost-effective and sustainable manner but also effectively manipulates the electronic structure of the active Ni-N 4 site through F doping. The catalyst has also been further optimized by controlling the F states, including convalent and semi-ionic states, by adjusting the fluorine sources involved. Consequently, this catalyst with unique structure exhibited comparable electrocatalytic performance for CO 2 -to-CO conversion, achieving a Faradaic efficiency (FE) of over 99% across a wide potential range and an exceptional CO evolution rate of 9.5 x 10 4 h -1 at -1.16 V vs reversible hydrogen electrode (RHE). It also delivered a practical current of 400 mA cm -2 while maintaining more than 95% CO FE. Experimental analysis combined with density functional theory (DFT) calculations have also shown that F-doping modifies the electron configuration at the central Ni-N 4 sites. In conclusion, this modification lowers the energy barrier for CO 2 activation, thereby facilitating the production of the crucial *COOH intermediate.

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

Doping-induced antiferromagnetic bicollinear insulating state and superconducting temperature of monolayer FeSe systems

First-principles calculations of ionic-liquid-gated (Li,Fe)OHFeSe suggest a metal-insulator transition at a nominal Li/Fe ratio of 75/25 in the (Li,Fe)OH layer. While doping increases Fermi energy, the formation of an antiferromagnetic bicollinear phase is the key for the transition. It is expected that this insulating phase also exists in other FeSe systems upon heavy electron doping, and its presence can hinder the increase of superconducting temperature. Here these results offer clues on how to optimize superconductivity amid its interplay with magnetic properties in FeSe systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Adaptive Data-Driven Deep-Learning Surrogate Model for Frontal Polymerization in Dicyclopentadiene

Frontal polymerization (FP) is a self-sustaining curing process that enables rapid and energy-efficient manufacturing of thermoset polymers and composites. Computational methods conventionally used to simulate the FP process are time-consuming, and repeating simulations are required for sensitivity analysis, uncertainty quantification, or optimization of the manufacturing process. Here, in this work, we develop an adaptive surrogate deep-learning model for FP of dicyclopentadiene (DCPD), which predicts the evolution of temperature and degree of cure orders of magnitude faster than the finite-element method (FEM). The adaptive algorithm provides a strategy to select training samples efficiently and save computational costs by reducing the redundancy of FEM-based training samples. The adaptive algorithm calculates the residual error of the FP governing equations using automatic differentiation of the deep neural network. A probability density function expressed in terms of the residual error is used to select training samples from the Sobol sequence space. The temperature and degree of cure evolution of each training sample are obtained by a 2D FEM simulation. The adaptive method is more efficient and has a better prediction accuracy than the random sampling method. With the well-trained surrogate neural network, the FP characteristics (front speed, shape, and temperature) can be extracted quickly from the predicted temperature and degree-of-cure fields.

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

To address surface reaction network complexity using scaling relations machine learning and DFT calculations

Surface reaction networks involving hydrocarbons exhibit enormous complexity with thousands of species and reactions for all but the very simplest of chemistries. We present a framework for optimization under uncertainty for heterogeneous catalysis reaction networks using surrogate models that are trained on the fly. The surrogate model is constructed by teaching a Gaussian process adsorption energies based on group additivity fingerprints, combined with transition-state scaling relations and a simple classifier for determining the rate-limiting step. The surrogate model is iteratively used to predict the most important reaction step to be calculated explicitly with computationally demanding electronic structure theory. Applying these methods to the reaction of syngas on rhodium(111), we identify the most likely reaction mechanism. Lastly, propagating uncertainty throughout this process yields the likelihood that the final mechanism is complete given measurements on only a subset of the entire network and uncertainty in the underlying density functional theory calculations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Power Generation from a Radiative Thermal Source Using a Large-Area Infrared Rectenna

Electrical power generation from a moderate-temperature thermal source by means of direct conversion of infrared radiation is important and highly desirable for energy harvesting from waste heat and micropower applications. Here, we demonstrate direct rectified power generation from an unbiased large-area nanoantenna-coupled tunnel diode rectifier called a rectenna. Using a vacuum radiometric measurement technique with irradiation from a temperature-stabilized thermal source, a generated power density of 8 nW / cm 2 is observed at a source temperature of 450°C for the unbiased rectenna across an optimized load resistance. The optimized load resistance for the peak power generation for each temperature coincides with the tunnel diode resistance at zero bias and corresponds to the impedance matching condition for a rectifying antenna. Current-voltage measurements of a thermally illuminated large-area rectenna show current zero crossing shifts into the second quadrant indicating rectification. Photon-assisted tunneling in the unbiased rectenna is modeled as the mechanism for the large short-circuit photocurrents observed where the photon energy serves as an effective bias across the tunnel junction. Furthermore, the measured current and voltage across the load resistor as a function of the thermal source temperature represents direct current electrical power generation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

HfNb_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-niobium (Hf-Nb). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Nb. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect BCC lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,029 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Alloying and Defect Control within Chalcogenide Perovskites for Optimized Photovoltaic Application

Through density functional theory calculations, we show that the alloy perovskite system BaZr 1-x Ti x S 3 (x < 0.25) is a promising candidate for producing high power conversion efficiency (PCE) solar cells with ultrathin absorber layers. To maximize the minority carrier lifetime, which is important for achieving high PCE, the defect calculations show that BaZr 1-x Ti x S 3 films should be synthesized under moderate (i.e., near stoichiometric) growth conditions to minimize the formation of deep-level defects. The perovskite BaZrS 3 is also found to exhibit ambipolar self-doping properties, indicating the ability to form homo p–n junctions. However, our theoretical calculations and experimental solid-state reaction efforts indicate that the doped perovskite BaZr 1-x Ti x S 3 (x > 0) may not be stable under thermal equilibrium growth conditions. Calculations of decomposition energies suggest that introducing compressive strain may be a plausible approach to stabilize BaZr 1-x Ti x S 3 thin films.

14 SOLAR ENERGY↗
DOE OSTI2020

Accelerated Modeling of Lithium Diffusion in Solid State Electrolytes using Artificial Neural Networks

Abstract Previous efforts to understand structure‐function relationships in high ionic conductivity materials for solid state batteries have predominantly relied on density functional theory (DFT‐) based ab initio molecular dynamics (MD). Such simulations, however, are computationally demanding and cannot be reasonably applied to large systems containing more than a hundred atoms. Here, an artificial neural network (ANN) is trained to accelerate the calculation of high accuracy atomic forces and energies used during such MD simulations. After carefully training a robust ANN for four and five element systems, nearly identical lithium ion diffusivities are obtained for Li 10 GeP 2 S 12 (LGPS) when benchmarking the ANN‐MD results with DFT‐MD. Applying the ANN‐MD approach, the effect of chlorine doping on the lithium diffusivity is calculated in an LGPS‐like structure and it is found that a dopant concentration of 1.3% maximizes ionic conductivity. The optimal concentration balances the competing consequences of effective atomic radii and dielectric constants on lithium diffusion and agrees with the experimental composition. Performing simulations at the resolution necessary to model experimentally relevant and optimal concentrations would be infeasible with traditional DFT‐MD. Systems that require a large number of simulated atoms can be studied more efficiently while maintaining high accuracy with the proposed ANN‐MD framework.

Rao, Karun K.↗
DOE OSTI2018

Overcoming ammonia synthesis scaling relations with plasma-enabled catalysis

Correlations between the energies of elementary steps limit the rates of thermally catalysed reactions at surfaces. Here, we show how these limitations can be circumvented in ammonia synthesis by coupling catalysts to a non-thermal plasma. We postulate that plasma-induced vibrational excitations in N 2 decrease dissociation barriers without influencing subsequent reaction steps. We develop a density-functional-theory-based microkinetic model to incorporate this effect, and parameterize the model using N 2 vibrational excitations observed in a dielectric-barrier-discharge plasma. We predict plasma enhancement to be particularly great on metals that bind nitrogen too weakly to be active thermally. Ammonia synthesis rates observed in a dielectric-barrier-discharge plasma reactor are consistent with predicted enhancements and predicted changes in the optimal metal catalyst. The results provide guidance for optimizing catalysts for application with plasmas.

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

Mechanism for spontaneous oxygen and hydrogen evolution reactions on CoO nanoparticles

Overall photocatalytic water splitting with a high efficiency of ~5% has recently been observed for CoO nanoparticle suspensions in the absence of an applied bias or co-catalyst. Although experimental measurements indicate that the overall photocatalytic water splitting is caused by optimal band edge alignments with respect to the redox potentials of water, the mechanism by which H 2 and O 2 simultaneously evolve on these nanoparticles is unknown. Here, we used first-principles density functional theory (DFT) calculations to elucidate the mechanisms for the charge separation and H 2 and O 2 evolution on CoO nanoparticles under illumination in aqueous solution. We demonstrated that electrons are driven to the CoO(100) facet and holes are driven to the hydroxylated CoO(111) facet (OH*–CoO(111)) as a result of the built-in potential arising from the difference in the band edge positions on the two facets. Furthermore, based on a set of criteria, depending on if the photoexcited electrons and holes have sufficient energy to overcome the kinetic barrier along the H 2 and O 2 evolution reaction pathways, respectively, on the relevant surface facet, we show that H 2 evolution preferentially occurs on the CoO(100) facet, while O 2 evolution occurs on the OH*–CoO(111) surface. Our understanding of the overall water splitting mechanism on CoO nanoparticles provides a general explanation for the experimentally observed overall water splitting phenomena on a variety of self-standing photocatalysts, including γ-Ga 2 O 3 , Cu 2 O, and KTaO 3 , without an external driving potential or co-catalyst. In addition, we provide a new strategy for designing novel photocatalysts with high efficiency by controlling their surface configurations and morphologies.

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

Carboxylate binding prefers two cations to one

Almost all studies of specific ion binding by carboxylates (–COO - ) have considered only a single cation, but clustering of ions and ligands is a common phenomenon. We apply density functional theory to investigate how variations in the number of acetate ligands in binding to two monovalent cations affects ion binding preferences. We study a series of monovalent (Li + , Na + , K + , Cs + ) ions relevant to experimental work on many topics, including ion channels, battery storage, water purification and solar cells. We find that the preferred optimal structure has 3 acetates except for Cs + , which has 2 acetates. The optimal coordination of the cation by the carboxylate O atoms is 4 for both Na + and K + , and 3 for Li + and Cs + . There is a 4-fold coordination minimum just a few kcal mol -1 higher than the optimal 3-fold structure for Li+. For two cations, multiple minima occur in the vicinity of the lowest free energy state. Here we find that, for Li, Na and K, the preferred optimal structure with two cations is favored over a mixture of single cation complexes, providing a basis for understanding ionic cluster formation that is relevant for engineering proteins and other materials for rapid, selective ion transport.

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

Validating first-principles molecular dynamics calculations of oxide/water interfaces with x-ray reflectivity data

Metal oxide/water interfaces play a crucial role in many electrochemical and photocatalytic processes, such as photoelectrochemical water splitting, the creation of fuel from sunlight, and electrochemical CO 2 reduction. First-principles electronic structure calculations can reveal unique insights into these processes, such as the role of the alignment of the oxide electronic energy levels with those of liquid water. An essential prerequisite for the success of such calculations is the ability to predict accurate structural models of these interfaces, which in turn requires careful experimental validation. Here we report a general, quantitative validation protocol for first-principles molecular dynamics simulations of oxide/aqueous interfaces. The approach makes direct comparisons of interfacial x-ray reflectivity (XR) signals from experimental measurements and those obtained from ab initio simulations with semilocal and van der Waals functionals. The protocol is demonstrated here for the case of the Al 2 O 3 (001)/water interface, one of the simplest oxide/water interfaces. We discuss the technical requirements needed for validation, including the choice of the density functional, the simulation cell size, and the optimal choice of the thermodynamic ensemble. Our results establish a general paradigm for the validation of structural models and interactions at solid/water interfaces derived from first-principles simulations. Furthermore, while there is qualitative agreement between the simulated structures and the experimental best-fit structure, direct comparisons of simulated and measured XR intensities show quantitative discrepancies that derive from both bulk regions (i.e., alumina and water) as well as the interfacial region, highlighting the need for accurate density functionals to properly describe interfacial interactions. Our results show that XR data are sensitive not only to the atomic structure (i.e., the atom locations) but also to the electron-density distributions in both the substrate and at the interface.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Statistical characteristics of cloud variability. Part 2: Implication for parameterizations of microphysical and radiative transfer processes in climate models

The effects of subgrid cloud variability on grid-average microphysical rates and radiative fluxes are examined by use of long-term retrieval products at the Tropical West Pacific, Southern Great Plains, and North Slope of Alaska sites of the Department of Energy's Atmospheric Radiation Measurement program. Four commonly used distribution functions, the truncated Gaussian, Gamma, lognormal, and Weibull distributions, are constrained to have the same mean and standard deviation as observed cloud liquid water content. The probability density functions are then used to upscale relevant physical processes to obtain grid-average process rates. It is found that the truncated Gaussian representation results in up to 30% mean bias in autoconversion rate, whereas the mean bias for the lognormal representation is about 10%. The Gamma and Weibull distribution function performs the best for the grid-average autoconversion rate with the mean relative bias less than 5%. For radiative fluxes, the lognormal and truncated Gaussian representations perform better than the Gamma and Weibull representations. The results show that the optimal choice of subgrid cloud distribution function depends on the nonlinearity of the process of interest, and thus, there is no single distribution function that works best for all parameterizations. Examination of the scale (window size) dependence of the mean bias indicates that the bias in grid-average process rates monotonically increases with increasing window sizes, suggesting the increasing importance of subgrid variability with increasing grid sizes.

54 ENVIRONMENTAL SCIENCES↗