Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DENSITY FUNCTIONAL METHOD”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Spatial locality of electronic correlations in LiFeAs

In this work, we address the question of the degree of spatial nonlocality of the self-energy in the iron-based superconductors, a subject which is receiving considerable attention. Using LiFeAs as a prototypical example, we extract the self-energy from angular-resolved photoemission spectroscopy data. We use two distinct electronic structure references: density functional theory in the local density approximation and linearized quasiparticle self-consistent GW (LQSGW). We find that with the LQSGW reference, spatially local dynamical correlations provide a consistent description of the experimental data, and account for some surprising aspects of the data such as the substantial out-of-plane dispersion of the electron Fermi surface having dominant xz/yz character. Hence, correlations effects can be separated into static nonlocal contributions well described by LQSGW and dynamical local contributions. Hall effect and resistivity data are shown to be consistent with this description.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

The behavior of methane–water mixtures under elevated pressures from simulations using many-body potentials

Non-polarizable empirical potentials have been proven to be incapable of capturing the mixing of methane–water mixtures at elevated pressures. Although density functional theory-based ab initio simulations may circumvent this discrepancy, they are limited in terms of the relevant time and length scales associated with mixing phenomena. Here, we show that the many-body MB-nrg potential, designed to reproduce methane–water interactions with coupled cluster accuracy, successfully captures this phenomenon up to 3 GPa and 500 K with varying methane concentrations. Two-phase simulations and long time scales that are required to fully capture the mixing, affordable due to the speed and accuracy of the MBX software, are assessed. Constructing the methane–water equation of state across the phase diagram shows that the stable mixtures are denser than the sum of their parts at a given pressure and temperature. As a result, we find that many-body polarization plays a central role, enhancing the induced dipole moments of methane by 0.20 D during mixing under pressure. Overall, the mixed system adopts a denser state, which involves a significant enthalpic driving force as elucidated by a systematic many-body energy decomposition analysis.

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

First principles study of the electronic structure of the Ni 2 MnIn / InAs and Ti 2 MnIn / InSb interfaces

We present a first principles study of the electronic and magnetic properties of epitaxial interfaces between the Heusler compounds, Ti 2 MnIn and Ni 2 MnIn, and the III-V semiconductors, InSb and InAs, respectively. We use density functional theory (DFT) with a machine-learned Hubbard U correction determined by Bayesian optimization. Here, we evaluate these interfaces for prospective applications in Majorana-based quantum computing and spintronics. In both interfaces, states from the Heusler penetrate into the gap of the semiconductor, decaying within a few atomic layers. The magnetic interactions at the interface are weak and local in space and energy. Magnetic moments of less than 0.1 µB are induced in the two atomic layers closest to the interface. The induced spin polarization around the Fermi level of the semiconductor decays within a few atomic layers. The decisive factor for the induced spin polarization around the Fermi level of the semiconductor is the spin polarization around the Fermi level in the Heusler, rather than the overall magnetic moment. As a result, the ferrimagnetic narrow-gap semiconductor Ti 2 MnIn induces a more significant spin polarization in the InSb than the ferromagnetic metal Ni 2 MnIn induces in the InAs. This is explained by the position of the transition metal d states in the Heusler with respect to the Fermi level. Based on our results, these interfaces are unlikely to be useful for Majorana devices but could be of interest for spintronics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Semi-supervised Machine Learning Enables the Robust Detection of Multireference Character at Low Cost

Multireference (MR) diagnostics are common tools for identifying strongly correlated electronic structure that makes single-reference (SR) methods (e.g., density functional theory or DFT) insufficient for accurate property prediction. However, MR diagnostics typically require computationally demanding correlated wave function theory (WFT) calculations, and diagnostics often disagree or fail to predict MR effects on properties. To overcome these challenges, we introduce a semi-supervised machine learning (ML) approach with virtual adversarial training (VAT) of an MR classifier using 15 WFT and DFT MR diagnostics as inputs. In semi-supervised learning, only the most extreme SR or MR points are labeled, and the remaining point labels are learned. The resulting VAT model outperforms the alternatives, as quantified by the distinct property distributions of SR- and MR-classified molecules. Additionally, to reduce the cost of generating inputs to the VAT model, we leverage the VAT model’s robustness to noisy inputs by replacing WFT MR diagnostics with regression predictions in an MR decision engine workflow that preserves excellent performance. We demonstrate the transferability of our approach to larger molecules and those with distinct chemical composition from the training set. This MR decision engine demonstrates promise as a low-cost, high-accuracy approach to the automatic detection of strong correlation for predictive high-throughput screening.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Relationship between atomic and electronic structure in Ln-bearing oxides

Metastable states of matter are of great interest as they offer the promise of novel functionality. They are often a natural consequence of exposure to nonequilibrium environments. In oxides, metastability can take the form of new polymorphs, chemical disorder, and even amorphization. While significant attention has been given to the impact those changes have on the atomic properties of the material, the corresponding changes in the electronic structure have received less attention. Here, using density functional theory, we consider how the electronic structure varies with potential metastable structures in two classes of lanthanide-bearing oxides—pyrochlores and interlanthanide sesquioxides. We find that the changes depend strongly on both the crystal structure and crystal chemistry of the compound with, for example, disordering and amorphization either increasing or decreasing the bandgap depending on the chemistry. For the 𝐴 2 ⁢𝐵 2 ⁢O 7 pyrochlores, we find different dependencies of the bandgap on the 𝐴 = 𝐿⁢𝑛 cations as the 𝐵 cation is changed, which we relate to the nature of the density of states at the conduction band minimum for different 𝐵 chemistries. Our calculations are validated by electron energy loss spectroscopy measurements for two pyrochlore compounds in which amorphization does reduce the bandgap, consistent with our calculations on these two compounds. In conclusion, our results highlight the relationship between atomic and electronic structure and how radiation can be used to modify and potentially control the electronic properties of oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Twisted Diindeno‐Fused Dibenzo[ a , h ]anthracene Derivatives and their Dianions

Abstract We report a facile synthesis of diindeno‐fused dibenzo[a,h]anthracene derivatives ( DIDBA‐2Cl , DIDBA‐2Ph , and DIDBA‐2H) with different degrees of non‐planarity using three substituents (chloro, phenyl, and hydrogen) of various sizes. The planarization of their cores, as evidenced by the decreased end‐to‐end torsional angles, was confirmed by X‐ray crystallography. Their enhanced energy gaps with twisting were investigated by a combination of spectroscopic and electrochemical methods with density functional theory, which showed a transition from singlet open‐shell to closed‐shell configuration. Moreover, their doubly reduced states, DIDBA‐2Ph 2− and DIDBA‐2H 2− , were achieved by chemical reduction. The structures of dianions were identified by X‐ray crystallographic analysis, which elucidated that the electron charging further distorted the backbones. The electronic structure of the dianions was demonstrated by experimental and theoretical approaches, suggesting decreased energy gaps with larger non‐planarity, different from the neutral species.

Gu, Yanwei↗
DOE OSTI2023

Twisted Diindeno‐Fused Dibenzo[ a , h ]anthracene Derivatives and their Dianions

Abstract We report a facile synthesis of diindeno‐fused dibenzo[a,h]anthracene derivatives ( DIDBA‐2Cl , DIDBA‐2Ph , and DIDBA‐2H) with different degrees of non‐planarity using three substituents (chloro, phenyl, and hydrogen) of various sizes. The planarization of their cores, as evidenced by the decreased end‐to‐end torsional angles, was confirmed by X‐ray crystallography. Their enhanced energy gaps with twisting were investigated by a combination of spectroscopic and electrochemical methods with density functional theory, which showed a transition from singlet open‐shell to closed‐shell configuration. Moreover, their doubly reduced states, DIDBA‐2Ph 2− and DIDBA‐2H 2− , were achieved by chemical reduction. The structures of dianions were identified by X‐ray crystallographic analysis, which elucidated that the electron charging further distorted the backbones. The electronic structure of the dianions was demonstrated by experimental and theoretical approaches, suggesting decreased energy gaps with larger non‐planarity, different from the neutral species.

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

Structural, Electronic, and Photophysical Insights into a Few Atom Copper-Sulfur Cluster in the Solid and Solution States

Coinage-metal chalcogenide clusters are widely studied for their attractive photoluminescence properties. Copper chalcogenides are especially promising, but are often confined to solid-state investigations due to their limited solution stability and the difficulty of synthesizing stable, well-defined clusters. Here, we investigate copper–sulfur clusters incorporating a small number of Cu atoms to elucidate fundamental atomic interactions, ground- and excited-state characteristics, and photophysical behavior in both solid and solution. We have synthesized the Cu6(4,6-dimethyl-2-mercaptopyrimidine)6 cluster in both neutral and charged states, Cu6 and Cu6-2H2+, respectively, by selective ligand protonation. The molecular structures are determined using single-crystal X-ray diffraction, while Cu K-edge X-ray absorption spectroscopy is used to probe Cu electronic structure differences arising from the ligand modification. Steady-state and pump-probe optical spectroscopy is used to investigate photophysical properties, interpreted using density functional theory methods. Both clusters exhibit good stability in the solid state and in solution and show characteristic near-infrared emission with microsecond lifetimes. Overall, the Cu6S6 clusters display favorable charge–transfer characteristics and show potential for further use in driving photochemical transformations.

Copper-sulfur clusters↗
DOE OSTI2021

What's Left for a Computational Chemist To Do in the Age of Machine Learning?

Machine learning (ML) has become a central focus of the computational chemistry community. In this paper, I will first discuss my personal history in the field. Then I will provide a broader view of how this resurgence in ML interest echoes and advances upon earlier efforts. Although numerous changes have brought about this latest wave, one of the most significant is the increased accuracy and efficiency of low-cost methods (e.g., density functional theory or DFT) that have made it possible to generate large data sets for ML models. ML has also been used to bypass, guide, or improve DFT. The field of computational chemistry thus finds itself at a crossroads as ML both augments and supersedes traditional efforts. I will present what I believe the role of the computational chemist will be in this evolving landscape, with specific focus on my experience in the development of autonomous workflows in computational materials discovery for open-shell transition-metal chemistry.

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

“Self-trapping” in solar cell hybrid inorganic-organic perovskite absorbers

In the simplest picture, a “self-trapped” polaron forms when an excess electron or hole deforms a crystal lattice, creating a potential well with bound states. Properties of self-trapped polarons in methylammonium lead iodide perovskite (MAPbI 3 ), which is widely used as solar cell absorber, are of great interest, and are a subject of ongoing investigations and debates concerning the existence of large polarons with the co-presence of metastable self-trapping. Herein, we employ a self-interaction-free density functional theory method to investigate the stability of small polarons in tetragonal MAPbI 3 phase. The electron small polaron is found to be unstable, while the hole small polaron is found to be metastable at realistic operation temperatures of solar cells. Further, the hole polaron is found to have a hole band close to the conduction band, which in conjunction with its metastability suggests that small polarons will have an appreciable effect on charge-carrier recombinations in MAPbI 3 . Further, in this study, we posit that the existence of the metastable polarons in addition to the large polarons may explain the experimentally observed non-monotonic temperature dependence of bimolecular charge-carrier recombination rate in tetragonal MAPbI 3 phase.

14 SOLAR ENERGY↗
DOE OSTI2025

Effects of iron carbide crystal phases and dopants on the conversions of CO 2 into ethylene

The density functional theory method was used to investigate the conversions of CO 2 to ethylene formation on two common iron carbide surfaces: Fe 3 C(0 1 0) and Fe 5 C 2 (1 1 1). Based on the structure relaxation of reaction intermediates and the elementary reaction transition states. We deduced the most competitive reaction pathways for ethylene production. The main CO 2 -to-ethylene routes and the competition of side products, CO and CH 4 , are discussed. Our analyses showed that CO 2 conversion is surface structure sensitive, whereas CH 4 and C2+ hydrocarbon formations depend on the reactivity of native C atoms in the carbides. To modify the intrinsic catalyst performance, mixing dopants in Fe catalysts is an effective strategy. Furthermore, we demonstrate that doping Zn and Zr can alter the local electronic structure and enhance CO 2 adsorption on the catalyst surface.

CO2 hydrogenation↗
DOE OSTI2021

Modeling of aqueous species interaction energies prior to nucleation in cement-based gel systems

Arguably the most ubiquitous construction material in modern civilization, concrete is enabling the development of megacities around the globe together with increasing living standards in developing nations. However, the exact formation mechanisms of the strength-giving calcium-rich gels remain a topic of debate. Using density functional modeling, we simulate the fundamental solution-based building blocks of cement hydrates (calcium ions and silicate and aluminate monomers) and their propensity to form pair-wise complexes with bonding environments characteristic of those found in calcium-silicate-hydrate, calcium-alumino-silicate-hydrate and sodium-containing calcium-alumino-silicate-hydrate gels, as assessed from Gibbs free energies of chemical reactions. By accurately simulating the high pH pore solution chemistry in Portland cements and related systems, along with discrete solvation of the species, we hypothesize potential early age formation routes of the gels and discuss limitations and future work associated with this approach.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Intrinsic reactivity and dissolution characteristics of tetracalcium aluminoferrite

Tetracalcium aluminoferrite is one of the main phases of Portland cement, the dissolution characteristics of which however have not been fully revealed. This work adopted first-principles calculations and Mössbauer spectroscopy as well as X-ray diffraction analysis to clarify the intrinsic reactivity of different species of tetracalcium aluminoferrite crystals. Ca ions are the most reactive species of tetracalcium aluminoferrite and Al ions are slightly more reactive than Fe ions. Moreover, the tetrahedral Al and Fe ions are more reactive than the octahedral Al and Fe ions. These indicate dissolution characteristics of tetracalcium aluminoferrite crystals, namely, Ca ions dissolve first and the tetrahedral species dissolve faster than the octahedral species. After the initial dissolution, the dissolved ions nucleate and grow on the defective surface of C{sub 4}AF and the C{sub 4}AF residues are fully hydroxylated. The further hydration thereafter becomes the complex dissolution with the nucleation and growth of hydration products.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Machine learning assisted modeling of mixing timescale for LES/PDF of high-Karlovitz turbulent premixed combustion

Accurate modeling of mixing in the transported probability density function (PDF) method remains a great challenge, especially for turbulent premixed combustion under extreme conditions such as high Karlovitz number Ka. Recently, a power-law based mixing timescale model was developed for the large-eddy simulations (LES)/PDF modeling of high-Ka number turbulent premixed flames. It is found in this work that the power-law mixing timescale model is highly sensitive to the model parameters. It is thus critically needed to develop accurate calibration of these model parameters. The empirical specification of the model parameters developed in Zhang et. al. is found to be inadequate for accurate modeling of the mixing timescale. Here, machine learning is introduced as an attractive alternative in this work for the specification of the model parameters. A high-Ka number DNS jet flame is used as the training and validation of the machine learning models. The choices of the input parameters are discussed and compared for the machine learning models. The effect of differential molecular diffusion on mixing is examined by including the effect of the Lewis number in the training of the machine learning models. The performance of different machine learning algorithms is compared for the specification of the mixing model parameters. Overall, excellent performance of the machine learning models is observed for assisting the mixing modeling. The feasibility, interpretability, applicability, generality, and portability of using machine learning are discussed in general to provide a perspective on applying data-driven machine learning for turbulent combustion modeling studies.

42 ENGINEERING↗
DOE OSTI2022

High temperature stability and transport characteristics of hydrogen in alumina via multiscale computation

Here, the impact of hydrogen charge states on the stability and transport characteristics of hydrogen interstitials in alumina polymorphs is evaluated by multiscale computational methods including density functional theory (DFT), ab initio molecular dynamics (AIMD) and machine learned force fields. Thermodynamic calculations show that the protonic H i +1 interstitial is the most stable defect species for most values of the electronic bandgap in both and amorphous alumina (Al 2 O 3 ). Further, active learned Gaussian approximation potentials (GAP) were developed using AIMD data to study temperature dependent long time proton diffusion in alumina. Diffusivity calculations from GAP-MD simulations are found to be comparable with of the AIMD data, while being ~340 times faster and scalable to larger systems. Comparisons with diffusivity values for other interstitial charge states (H i 0 and H i -1 ) and published experimental literature indicate that H i +1 diffusion is the likely mechanism of hydrogen transport. A good agreement is obtained between H i +1 diffusivity calculated in α-Al 2 O 3 from DFT: 5.05 10 -3 exp(-0.81 eV/k B /T) cm 2 /s and reported experiment: 9.7X10 -4 exp(-0.83 eV/k B /T) cm 2 /s. Computationally and experimentally calculated energy barriers (0.81 and 0.83 eV respectively) only differ by 2.5%. Similarly, the pre-exponential diffusion coefficients only differ by 0.5 orders of magnitude. Moreover, the diffusivity of H i +1 in amorphous Al 2 O 3 in the 1000–2000 K range is calculated to be 2.53X10 -2 exp(-0.89 eV/k B /T), just one order of magnitude higher than the corresponding value in α-Al 2 O 3 . This suggests that local structural disorder does not significantly affect the energy landscape and diffusion behavior of H i +1 in Al 2 O 3 . Overall, these results show promise for the application of alumina polymorphs as hydrogen permeation barriers.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Analysis and prediction of reaction kinetics using the degree of rate control

“Degree of rate control” (DRC) analysis provides a quantitative approach for analysing the kinetics of multi-step reaction mechanisms that has been widely applied to both heterogeneous and homogeneous catalysis research, as well as electrocatalysis. The DRC of any given transition state or intermediate is defined as a partial derivative such that it approximately equals the fractional increase in net rate to the product of interest per differential decrease in its standard-state free energy for that species (÷RT), holding constant the standard-state free energies of all other transition states and intermediates. Even very complex mechanisms usually have only a few species with non-zero DRCs and are thus the species whose interactions with the catalyst most strongly affect the net rate. These key DRC values thus offer a simple and intuitive route to optimize catalyst materials, especially with the assistance of computational methods like density functional theory (DFT). These high-DRC species are also the species whose energetics must be most accurately measured or calculated to achieve an accurate kinetic model for any reaction mechanism. In simple cases with a single “rate-determining step”, the DRC for its transition state (TS) is + 1. Catalyst-bound intermediates, on the other hand, often have negative DRCs equal to a small integer times their fractional occupancy of catalyst sites. The apparent activation energy equals a weighted average of the standard-state enthalpies (relative to reactants) of all the species (intermediates, transition states and products) in the reaction mechanism, each weighted by its DRC (+RT). It has been shown that the apparent transfer coefficient in electrocatalysis, an inverted form of the Tafel slope, is a weighted average of the number of electrons transferred to generate each intermediate or product species in the mechanism, weighted again by the DRC. Quantitative analysis of kinetic isotope effects (KIEs, or the ratio of net rates for different reactant isotopes) in complex mechanisms has shown that the logarithm of the KIE equals the weighted average over all species in the mechanism of the difference between the two isotopes in their standard-state free energies (÷RT), again weighted by the DRC. The reaction orders with respect to fluid-phase concentrations of reactants, products and intermediates have also been proven to be directly related to DRCs. Thus, there are numerous experimental observables which equate to short linear combinations of DRCs, so that combinations of experimental measurements might provide access to DRC values. Since its invention for transition states in 1994 and its generalization to include intermediates in 2009, the DRC has thus far mainly been calculated for microkinetic models based either entirely on DFT or on DFT with the key energies (i.e., those for high-DRC species) being fine-tuned to match experiments. The relationships summarized in this work provide new opportunities for using experiments earlier in the development and optimization of microkinetic models that require input from computational catalysis.

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

Ultrafast X-ray Spectroscopy of Intersystem Crossing in Hexafluoroacetylacetone: Chromophore Photophysics and Spectral Changes in the Face of Electron-Withdrawing Groups

Intersystem crossings between singlet and triplet states represent a crucial relaxation pathway in photochemical processes. Herein, we probe the intersystem crossing in hexafluoro-acetylacetone with ultrafast X-ray transient absorption spectroscopy at the carbon K-edge. In this study, we observe the excited state dynamics following excitation with 266 nm UV light to the 1 ππ* (S 2 ) state with element and site-specificity using a broadband soft X-ray pulse produced by high harmonic generation. These results are compared to X-ray spectra computed from orbital optimized density functional theory methods. It is found that the electron-withdrawing fluorine atoms decongest the X-ray absorption spectrum by enhancing separation between features originating from different carbon atoms. This facilitates the elucidation of structural and electronic dynamics at the chromophore. The evolution of the core-to-valence resonances at the carbon K-edge reveals an ultrafast population transfer between the 1 nπ* (S 1 ) and 3 ππ* (T 1 ) states on a 1.6 ± 0.4 ps time scale, which is similar to the 1.5 ps time scale earlier observed for acetylacetone [ J. Am. Chem. Soc. 2017, 139, 16576—16583, DOI: 10.1021/jacs.7b07532]. It therefore appears that terminal fluorination has little influence on the intersystem crossing rate of the acetylacetone chromophore. In addition, the significant role of hydrogen-bond opened and twisted rotational isomers is elucidated in the excited state dynamics by comparison of the experimental transient X-ray spectra with theory.

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

Revisiting the Hole Size in Double Helical DNA with Localized Orbital Scaling Corrections

The extent of electronic wave function delocalization for the charge carrier (electron or hole) in double helical DNA plays an important role in determining the DNA charge transfer mechanism and kinetics. The size of the charge carrier’s wave function delocalization is regulated by the solvation induced localization and the quantum delocalization among the π stacked base pairs at any instant of time. Using a newly developed localized orbital scaling correction (LOSC) density functional theory method, we accurately characterized the quantum delocalization of the hole wave function in double helical B-DNA. This approach can be used to diagnose the extent of delocalization in fluctuating DNA structures. Our studies indicate that the hole state tends to delocalize among 4 guanine–cytosine (GC) base pairs and among 3 adenine–thymine (AT) base pairs when these adjacent bases fluctuate into degeneracy. The relatively small delocalization in AT base pairs is caused by the weaker π–π interaction. This extent of delocalization has significant implications for assessing the role of coherent, incoherent, or flickering coherent carrier transport in DNA.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗