Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

KSSOLV 2.0: An efficient MATLAB toolbox for solving the Kohn-Sham equations with plane-wave basis set

KSSOLV (Kohn-Sham Solver) is a MATLAB toolbox for performing Kohn-Sham density functional theory (DFT) calculations with a plane-wave basis set. KSSOLV 2.0 preserves the design features of the original KSSOLV software to allow users and developers to easily set up a problem and perform ground-state calculations as well as to prototype and test new algorithms. Furthermore, it includes new functionalities such as new iterative diagonalization algorithms, k-point sampling for electron band structures, geometry optimization and advanced algorithms for performing DFT calculations with local, semi-local, and hybrid exchange-correlation functionals. It can be used to study the electronic structures of both molecules and solids. So, we describe these new capabilities in this work through a few use cases. We also demonstrate the numerical accuracy and computational efficiency of KSSOLV on a variety of examples.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2021

Element-specific investigations of ultrafast dynamics in photoexcited Cu2ZnSnS4 nanoparticles in solution

Ultrafast, light-induced dynamics in copper–zinc–tin–sulfide (CZTS) photovoltaic nanoparticles are investigated through a combination of optical and x-ray transient absorption spectroscopy. Laser-pump, x-ray-probe spectroscopy on a colloidal CZTS nanoparticle ink yields element-specificity, which reveals a rapid photo-induced shift of electron density away from Cu-sites, affecting the molecular orbital occupation and structure of CZTS. We observe the formation of a stable charge-separated and thermally excited structure, which persists for nanoseconds and involves an increased charge density at the Zn sites. Combined with density functional theory calculations, the results provide new insight into the structural and electronic dynamics of CZTS absorbers for solar cells.

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

Establishing the Maximum Collectivity in Highly Deformed N = Z Nuclei

The lifetimes of the first excited 2 + states in the N = Z nuclei 80 Zr, 78 Y, and 76 Sr have been measured using the γ-ray line shape method following population via nucleon-knockout reactions from intermediate-energy rare-isotope beams. The extracted reduced electromagnetic transition strengths yield new information on where the collectivity is maximized and provide evidence for a significant, and as yet unexplained, odd-odd vs even-even staggering in the observed values. As a result, the experimental results are analyzed in the context of state-of-the-art nuclear density-functional model calculations.

59 ≤ A ≤ 89↗
DOE OSTI2023

Why Dissolving Salt in Water Decreases Its Dielectric Permittivity

The dielectric permittivity of salt water decreases on dissolving more salt. For nearly a century, this phenomenon has been explained by invoking saturation in the dielectric response of the solvent water molecules. Herein, we employ an advanced deep neural network (DNN), built using data from density functional theory, to study the dielectric permittivity of sodium chloride solutions. Notably, the decrease in the dielectric permittivity as a function of concentration, computed using the DNN approach, agrees well with experiments. Detailed analysis of the computations reveals that the dominant effect, caused by the intrusion of ionic hydration shells into the solvent hydrogen-bond network, is the disruption of dipolar correlations among water molecules. Accordingly, the observed decrease in the dielectric permittivity is mostly due to increasing suppression of the collective response of solvent waters.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Solving the electronic structure problem for over 100000 atoms in real space

Using a real-space high-order finite-difference approach, we investigate the electronic structure of large spherical silicon nanoclusters. Within Kohn-Sham density functional theory and using pseudopotentials, we report the self-consistent field convergence of a system with over 100000 atoms: a Si 107,641 ⁢H 9,084 nanocluster with a diameter of 16 nm. Our approach uses Chebyshev-filtered subspace iteration to speed up the convergence of the eigenspace, and blockwise Hilbert space-filling curves to speed up sparse matrix-vector multiplications, all of which are implemented in the parsec code. For the largest system, we utilized 2048 nodes (114 688 cores) on the Frontera machine in the Texas Advanced Computing Center. Our quantitative analysis of the electronic structure shows how it gradually approaches its bulk counterpart as a function of nanocluster size. The band gap is enlarged due to quantum confinement in nanoclusters, but decreases as the system size increases, as expected. In conclusion, our work serves as a proof of concept for the capacity of the real-space approach in efficiently parallelizing very large calculations using high-performance computer platforms, which can straightforwardly be replicated in other systems with more than 10 5 atoms.

0-dimensional systems↗
DOE OSTI2023

First-principles molten salt phase diagrams through thermodynamic integration

Precise prediction of phase diagrams in molecular dynamics simulations is challenging due to the simultaneous need for long time and large length scales and accurate interatomic potentials. Here, we show that thermodynamic integration from low-cost force fields to neural network potentials trained using density-functional theory (DFT) enables rapid first-principles prediction of the solid–liquid phase boundary in the model salt NaCl. We use this technique to compare the accuracy of several DFT exchange–correlation functionals for predicting the NaCl phase boundary and find that the inclusion of dispersion interactions is critical to obtain good agreement with experiment. Importantly, our approach introduces a method to predict solid–liquid phase boundaries for any material at an ab initio level of accuracy, with the majority of the computational cost at the level of classical potentials.

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

Theoretical Calculation Guided Design of Single-Atom Catalysts toward Fast Kinetic and Long-Life Li–S Batteries

Lithium–sulfur (Li–S) batteries are enticing next-generation energy storage technologies due to their high theoretical energy density, environmental friendliness, and low cost. Yet, low conductivity of sulfur species, dissolution of polysulfides, poor conversion from sulfur reduction, and lithium sulfide (Li 2 S) oxidation reactions during discharge–charge processes hinder their practical applications. Herein, under the guidance of density functional theory calculations, we have successfully synthesized large-scale single atom vanadium catalysts seeded on graphene to achieve high sulfur content (80 wt % sulfur), fast kinetic (a capacity of 645 mAh g –1 at 3 C rate), and long-life Li–S batteries. Both forward (sulfur reduction) and reverse reactions (Li 2 S oxidation) are significantly improved by the single atom catalysts. This finding is confirmed by experimental results and consistent with theoretical calculations. The ability of single metal atoms to effectively trap the dissolved lithium polysulfides (LiPSs) and catalytically convert the LiPSs/Li 2 S during cycling significantly improved sulfur utilization, rate capability, and cycling life. Our work demonstrates an efficient design pathway for single atom catalysts and offers solutions for the development of high energy/power density Li–S batteries.

25 ENERGY STORAGE↗
DOE OSTI2024

Violation of the Cauchy–Born rule in multi-principal element alloys

Multi-principal element alloys are a novel class of materials that are formed by combining multiple elements in high concentrations and show exceptional properties compared to conventional alloys. These alloys have high configurational entropy due to inherent atomic disorder. The Cauchy–Born rule is a popular homogenization method for linking atomistic to continuum length scales. In this Letter, we use ab initio density functional theory calculations to report that the Cauchy–Born rule, which holds in conventional alloys in the absence of defects, is not valid in multi-principal element alloys. Further, the violation of the Cauchy–Born rule in these alloys is attributed to the presence of atomic disorder due to which the deformations are inhomogeneous. Our results also show the deviations in deformation are related to the magnitude of the stretch and shear by power laws.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Global description of $β$ – decay with the axially deformed Skyrme finite-amplitude method: Extension to odd-mass and odd-odd nuclei

Here we use the finite-amplitude method (FAM), an efficient implementation of the quasiparticle random phase approximation, to compute $β$-decay rates with Skyrme energy-density functionals for 3983 nuclei, essentially all the medium-mass and heavy isotopes on the neutron-rich side of stability. We employ an extension of the FAM that treats odd-mass and odd-odd nuclear ground states in the equal filling approximation. Our rates are in reasonable agreement both with experimental data where available and with rates from other global calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Tunable magnetic order in low-symmetry SeO 3 ligand linked TM 3 (SeO 3 ) 3 H 2 O ( TM =Mn , Co, and Ni) compounds

The ability to design crystal lattices with specific sublattice geometries has long been sought after especially in the context of harnessing magnetic frustration to elicit emergent physics. One approach which has seen some success recently is to design magnetic sublattices chemically through the use of nonmagnetic linker ligands as scaffolding for the magnetic ions. In this work we report the magnetic properties of one such family of materials, the transition metal (TM) selenite hydrates with chemical formula TM 3 (SeO 3 ) 3 H 2 O. These materials link highly distorted TMO 6 octahedra via nonmagnetic [SeO 3 ] 2+ linkers. Studying members with TM = Mn, Co, and Ni we use magnetic susceptibility and neutron powder diffraction to identify antiferromagnetic order in all three compounds with moderate frustration indexes. A comparison of the magnetic structures suggests that while the overall structure of the SeO 3 scaffolding remains unchanged the TM effects changes in the magnetic properties through tuning the internal bonding parameters and the single-ion physics. Changing from Mn to Co is found to be particularly consequential manifesting changes in both the ordered-moment direction and in the direction of the ordering vector. Field-dependent measurements of the susceptibility and heat capacity reveal metamagnetic transitions in both Mn 3 (SeO 3 ) 3 H 2 O and Co 3 (SeO 3 ) 3 H 2 O indicating nearby magnetic ground states accessible under relatively small applied fields. Density functional theory calculations broadly confirm these results, showing both a sensitivity of the magnetic structure to the TM and its local environment. Although no spin liquid behavior is achieved, these results suggest the fruitfulness of such synthesis philosophies and encourage future work to engender higher frustration in these materials via doping, field, pressure, or larger linker ligands.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Prediction of hydration energies of adsorbates at Pt(111) and liquid water interfaces using machine learning

Aqueous phase heterogeneous catalysis is important to various industrial processes, including biomass conversion, Fischer–Tropsch synthesis, and electrocatalysis. Accurate calculation of solvation thermodynamic properties is essential for modeling the performance of catalysts for these processes. Explicit solvation methods employing multiscale modeling, e.g., involving density functional theory and molecular dynamics have emerged for this purpose. Although accurate, these methods are computationally intensive. This study introduces machine learning (ML) models to predict solvation thermodynamics for adsorbates on a Pt(111) surface, aiming to enhance computational efficiency without compromising accuracy. In particular, ML models are developed using a combination of molecular descriptors and fingerprints and trained on previously published water–adsorbate interaction energies, energies of solvation, and free energies of solvation of adsorbates bound to Pt(111). These models achieve root mean square error values of 0.09 eV for interaction energies, 0.04 eV for energies of solvation, and 0.06 eV for free energies of solvation, demonstrating accuracy within the standard error of multiscale modeling. Feature importance analysis reveals that hydrogen bonding, van der Waals interactions, and solvent density, together with the properties of the adsorbate, are critical factors influencing solvation thermodynamics. Furthermore, these findings suggest that ML models can provide rapid and reliable predictions of solvation properties. This approach not only reduces computational costs but also offers insights into the solvation characteristics of adsorbates at Pt(111)–water interfaces.

Adsorption↗
DOE OSTI2021

Experimental and calculated infrared spectra of disubstituted naphthoquinones

In recent years there has been interest in incorporating substituted 1,4-naphthoquinones (NQs) into the A 1 binding site in photosystem I (PSI) photosynthetic protein complexes. This interest in part stems from the considerably altered bioenergetics of electron transfer that occur in PSI with such substitutions. Time resolved FTIR studies of PSI complexes with disubstituted NQs incorporated have and currently are being undertaken, and with this in mind it is worth considering FTIR absorption spectra of these disubstituted NQs in solution. Here we present FTIR absorbance spectra for 2-bromo-3-methyl-1,4-naphthoquinone (BrMeNQ), 2-chloromethyl-3-methyl-1,4-naphthoquinone (CMMeNQ) and 2-ethylthio-3-methyl-1,4-naphthoquinone (ETMeNQ) in tetrahydrofuran (THF). The FTIR spectra of these di-substituted naphthoquinones (NQs) were compared to FTIR spectra of 2-methyl-3-phytyl-1,4-naphthoquinone [phylloquinone (PhQ)], 2,3-dimethyl-1,4-naphthoquinone (DMNQ), and 2-methyl-1,4-naphthoquinone (2MNQ). To aid in the assignment of bands in the experimental spectra, density functional theory (DFT) based vibrational frequency calculations for all the substituted NQs in solution were undertaken. The calculated and experimental spectra agree well. By calculating normal mode potential energy distributions, unambiguous quantitative band assignments were made. The calculated and experimental spectra together make predictions about what may be observable in time resolved FTIR difference spectra obtained using PSI with the different NQs incorporated. As a result, time resolved FTIR difference spectra are presented that support these predictions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Atomically dispersed zeolite-supported rhodium complex: Selective and stable catalyst for acetylene semi-hydrogenation

In this study, supported rhodium catalysts are known to be unselective for semi-hydrogenation reactions. Here, by tuning the electronic structure of supported mononuclear rhodium sites determined by the metal nuclearity and the electron-donor properties of the support, we report that atomically dispersed HY zeolite-supported rhodium with reactive acetylene ligands affords a stable ethylene selectivity > 90 % for acetylene semi-hydrogenation at 373 K and atmospheric pressure, even when ethylene is present in a large excess over acetylene. Infrared and X-ray absorption spectra and measurements of rates of the catalytic reaction complemented with calculations at the level of density functional theory show how the catalyst performance depends on the electronic structure of the rhodium, influenced by the support as a ligand that is a weak electron donor.

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

Combining Computational Modeling with Reaction Kinetics Experiments for Elucidating the In Situ Nature of the Active Site in Catalysis

Microkinetic modeling based on density functional theory (DFT) derived energetics is important for addressing fundamental questions in catalysis. The quantitative fidelity of microkinetic models (MKMs), however, is often insufficient to conclusively infer the mechanistic details of a specific catalytic system. Here, this can be attributed to a number of factors such as an incorrect model of the active site for which DFT calculations are performed, deficiencies in the hypothesized reaction mechanism, inadequate consideration of the surface environment under reaction conditions, and intrinsic errors in the DFT exchange-correlation functional. Despite these limitations, we aim at developing a rigorous understanding of the reaction mechanism and of the nature of the active site for heterogeneous catalytic chemistries under reaction conditions. By achieving parity between experimental and modeling outcomes through robust parameter estimation and by ensuring coverage-consistency between DFT calculations and MKM predictions, it is possible to systematically refine the mechanistic model and, thereby, our understanding of the catalytic active site in situ.

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

Probing pH-Dependent Dehydration Dynamics of Mg and Ca Cations in Aqueous Solutions with Multi-Level Quantum Mechanics/Molecular Dynamics Simulations

The dehydration of aqueous calcium and magnesium cations is the most fundamental process controlling their reactivity in chemical and biological phenomena, such as the formation of ionic solids or passing through ion channels. It holds particular relevance in light of recent advancements in the development of carbon capture techniques that rely on mineralization for long-term carbon storage. Specifically, dehydration of Ca 2+ and Mg 2+ is a key step in proposed carbon capture processes aiming to exploit the relatively high concentration of dissolved carbon dioxide in seawater via the formation of carbonate minerals from solvated Ca 2+ and Mg 2+ cations for sequestration and storage. Nevertheless, atomic-scale understanding of the dehydration of aqueous Ca 2+ and Mg 2+ cations remains limited. Here, we utilize rare event sampling via density functional theory molecular dynamics and embedded wavefunction theory calculations to elucidate the dehydration dynamics of aqueous Ca 2+ and Mg 2+ . Emphasis is placed on the investigation of the effect pH has on the stability of the different coordination environments. Here, our results reveal significant differences in the dehydration dynamics of the two cations and provide insight into how they may be modulated by pH changes.

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

Iodine Vacancies do not Cause Nonradiative Recombination in Halide Perovskites

The iodine vacancy (V I ) has frequently been discussed as a strong nonradiative recombination center in halide perovskites. This proposition was mainly based on the presence of charge-state transition levels in the band gap, as found in early first-principles calculations. In this work, we perform accurate hybrid-density-functional calculations for V I in CsPbI 3 , CsSnI 3 , and CsGeI 3 and find that V I does not have any transition levels in the band gap in CsPbI 3 , in contrast to the results from calculations based on semilocal functionals. The iodine vacancy V I does introduce levels in the band gap in CsSnI 3 and CsGeI 3 , but our explicitly computed nonradiative capture coefficients demonstrate that V I has a negligible impact on nonradiative recombination. Our study corrects a misunderstanding of the role of V I in the iodide-based perovskites, and shifts the focus toward identifying and mitigating actual recombination centers in order to further improve the optoelectronic performance.

14 SOLAR ENERGY↗
NASA NTRS1993

Geometric model from microscopic theory for nuclear absorption

A parameter-free geometric model for nuclear absorption is derived herein from microscopic theory. The expression for the absorption cross section in the eikonal approximation, taken in integral form, is separated into a geometric contribution that is described by an energy-dependent effective radius and two surface terms that cancel in an asymptotic series expansion. For collisions of light nuclei, an expression for the effective radius is derived from harmonic oscillator nuclear density functions. A direct extension to heavy nuclei with Woods-Saxon densities is made by identifying the equivalent half-density radius for the harmonic oscillator functions. Coulomb corrections are incorporated, and a simplified geometric form of the Bradt-Peters type is obtained. Results spanning the energy range from 1 MeV/nucleon to 1 GeV/nucleon are presented. Good agreement with experimental results is obtained.

John, Sarah↗