SEARCH · Engineering Papers
Results for “energy 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.
Simplified tuning of long-range corrected density functionals for use in symmetry-adapted perturbation theory
Long considered a failure, second-order symmetry-adapted perturbation theory (SAPT) based on Kohn–Sham orbitals, or SAPT0(KS), can be resurrected for semiquantitative purposes using long-range corrected density functionals whose asymptotic behavior is adjusted separately for each monomer. Here, as in other contexts, correct asymptotic behavior can be enforced via “optimal tuning” based on the ionization energy theorem of density functional theory, but the tuning procedure is tedious, expensive for large systems, and comes with a troubling dependence on system size. Here, we show that essentially identical results are obtained using a fast, convenient, and automated tuning procedure based on the size of the exchange hole. In conjunction with “extended” (X)SAPT methods that improve the description of dispersion, this procedure achieves benchmark-quality interaction energies, along with the usual SAPT energy decomposition, without the hassle of system-specific tuning.
The Fermi–Löwdin self-interaction correction for ionization energies of organic molecules
(Semi)-local density functional approximations (DFAs) suffer from self-interaction error (SIE). When the first ionization energy (IE) is computed as the negative of the highest occupied orbital (HO) eigenvalue, DFAs notoriously underestimate them compared to quasi-particle calculations. Here, the inaccuracy for the HO is attributed to SIE inherent in DFAs. We assessed the IE based on Perdew-Zunger self-interaction correction on 14 small to moderate-sized organic molecules relevant in organic electronics and polymer donor materials. Though self-interaction corrected DFAs were found to significantly improve the IE relative to the uncorrected DFAs, they overestimate. However, when the self-interaction correction is interiorly scaled using a function of the iso-orbital indicator zs , only the regions where SIE is significant get a correction. We discuss these approaches and show how these methods significantly improve the description of the HO eigenvalue for the organic molecules.
Machine-learning-aided density functional theory calculations of stacking fault energies in steel
A combined large-scale first principles approach with machine learning and materials informatics is proposed to quickly sweep the chemistry-composition space of advanced high strength steels (AHSS). AHSS are composed of iron and key alloying elements such as aluminum and manganese. A systematic exploration of the distribution of aluminum and manganese atoms in iron is used to investigate low stacking fault energies configurations using first principles calculations. To overcome the computational cost of exploring the composition space, this process is sped up using an automated machine learning tool: DeepHyper. Here our results predict that it is energetically favorable for Al to stay away from a stacking fault, but Mn atoms do not affect the stacking fault energy and can stay in the vicinity of the fault. The distribution of Al and Mn atoms in systems containing stacking faults and the effects of their interactions on the equilibrium distribution are systematically analyzed.
Evaluating Material Design Principles for Calcium-Ion Mobility in Intercalation Cathodes
Multivalent-ion batteries offer an alternative to Li-based technologies, with the potential for greater sustainability, improved safety, and higher energy density, primarily due to their rechargeable system featuring a passivating metal anode. Although a system based on the Ca 2+ /Ca couple is particularly attractive given the low electrochemical plating potential of Ca 2+ , the remaining challenge for a viable rechargeable Ca battery is to identify Ca cathodes with fast ion transport. In this work, a high-throughput computational pipeline is adapted to (1) discover novel Ca cathodes in a largely unexplored space of empty intercalation hosts and (2) develop material design rules for Ca-ion mobility. One candidate from the screening, W 2 O 3 (PO 4 ) 2 , is confirmed to have a low Nudged Elastic Band (NEB) barrier of 168 meV within a one-dimensional (1D) ion percolation topology. This candidate is subsequently synthesized and electrochemically tested, achieving reversible Ca cycling with a capacity of 25 mA h/g. To further accelerate the screening for promising Ca intercalation electrodes, machine learning (ML) Random Forest (RF) and Extreme Gradient Boosting (XGB) classification models are created with local environment descriptors based on a large, structurally and chemically diverse dataset of minimum energy pathways, spanning over 5,000 density functional theory (DFT) site energy calculations. Accuracies of 92% are achieved, material design metrics are quantified, ML force-fields are leveraged in an accelerated iteration of the screening, and a total of 27 novel Ca cathode materials are highlighted for further investigation.
Quantum chemical calculations of lithium-ion battery electrolyte and interphase species
Abstract Lithium-ion batteries (LIBs) represent the state of the art in high-density energy storage. To further advance LIB technology, a fundamental understanding of the underlying chemical processes is required. In particular, the decomposition of electrolyte species and associated formation of the solid electrolyte interphase (SEI) is critical for LIB performance. However, SEI formation is poorly understood, in part due to insufficient exploration of the vast reactive space. The Lithium-Ion Battery Electrolyte (LIBE) dataset reported here aims to provide accurate first-principles data to improve the understanding of SEI species and associated reactions. The dataset was generated by fragmenting a set of principal molecules, including solvents, salts, and SEI products, and then selectively recombining a subset of the fragments. All candidate molecules were analyzed at the ω B97X-V/def2-TZVPPD/SMD level of theory at various charges and spin multiplicities. In total, LIBE contains structural, thermodynamic, and vibrational information on over 17,000 unique species. In addition to studies of reactivity in LIBs, this dataset may prove useful for machine learning of molecular and reaction properties.
Excitation energy of fission fragments within nuclear time-dependent density functional theory
The number and properties of the neutrons and photons emitted in nuclear fission are directly related to the excitation energy of the fission fragments when they are formed at scission. Though not observable experimentally because of the extremely short timescales, the excitation energy of fission fragments can be predicted by microscopic theory based on time-dependent density functional theory (TDDFT). Initial results on the value of the total kinetic energy of fission reactions were very promising, but could not probe all possible fragmentations. Here, in this work, we perform large-scale TDDFT calculations in 240 Pu enabled by the development of a new TDDFT solver. We obtain TDDFT trajectories covering nearly all possible fragmentations. We find that the total kinetic energy is close to experimental values only for the most likely fission while it is severely underestimated at both small and large asymmetries. This conclusion seems rather independent of the parametrization of the energy functional, both in its particle-hole and particle-particle channels.
OrbNet Denali: A machine learning potential for biological and organic chemistry with semi-empirical cost and DFT accuracy
We present OrbNet Denali, a machine learning model for an electronic structure that is designed as a drop-in replacement for groundstate density functional theory (DFT) energy calculations. The model is a message-passing graph neural network that uses symmetryadapted atomic orbital features from a low-cost quantum calculation to predict the energy of a molecule. OrbNet Denali is trained on a vast dataset of 2.3 × 10 6 DFT calculations on molecules and geometries. This dataset covers the most common elements in biochemistry and organic chemistry (H, Li, B, C, N, O, F, Na, Mg, Si, P, S, Cl, K, Ca, Br, and I) and charged molecules. OrbNet Denali is demonstrated on several well-established benchmark datasets, and we find that it provides accuracy that is on par with modern DFT methods while offering a speedup of up to three orders of magnitude. For the GMTKN55 benchmark set, OrbNet Denali achieves WTMAD-1 and WTMAD-2 scores of 7.19 and 9.84, on par with modern DFT functionals. For several GMTKN55 subsets, which contain chemical problems that are not present in the training set, OrbNet Denali produces a mean absolute error comparable to those of DFT methods. For the Hutchison conformer benchmark set, OrbNet Denali has a median correlation coefficient of R 2 = 0.90 compared to the reference DLPNO-CCSD(T) calculation and R 2 = 0.97 compared to the method used to generate the training data (ωB97X-D3/def2- TZVP), exceeding the performance of any other method with a similar cost. Similarly, the model reaches chemical accuracy for noncovalent interactions in the S66x10 dataset. For torsional profiles, OrbNet Denali reproduces the torsion profiles of ωB97X-D3/def2-TZVP with an average mean absolute error of 0.12 kcal/mol for the potential energy surfaces of the diverse fragments in the TorsionNet500 dataset.
Exchange–correlation entropy from the generalized thermal adiabatic connection
Warm dense matter is a highly energetic phase characterized by strong correlations, thermal effects, and quantum mechanical electrons. Thermal density functional theory is commonly used in simulations of this challenging phase, driving the development of temperature-dependent approximations to the exchange–correlation free energy. Approaches using the adiabatic connection formula are well known at zero temperature and have been recently leveraged at non-zero temperatures as well. In this work, a generalized thermal adiabatic connection (GTAC) formula is proposed, introducing a fictitious temperature parameter. This allows extraction of the exchange–correlation entropy S XC using simulated interaction strength scaling. Here, this procedure uses a Hellmann–Feynman approach to express the exchange–correlation entropy in terms of a temperature- and interaction strength-dependent exchange–correlation potential energy. In addition, analysis of S XC as a function of interaction strength suggests new forms for approximations, and GTAC itself offers a new framework for exploring both the exact and approximate interplay of temperature, density, and interaction strength across a wide range of conditions.
Minimum-Energy Conical Intersections by Compressed Multistate Pair-Density Functional Theory
Compressed multistate pair-density functional theory (CMS-PDFT) is a multistate version of multiconfiguration pair-density functional theory that can capture the correct topology of coupled potential energy surfaces (PESs) around conical intersections. Here, in this work, we develop interstate coupling vectors (ISCs) for CMS-PDFT in the OpenMolcas and PySCF/mrh electronic structure packages. Yet, the main focus of this work is using ISCs to calculate minimum-energy conical intersections (MECIs) by CMS-PDFT. This is performed using the projected constrained optimization method in OpenMolcas, which uses ISCs to restrain the iterations to the conical intersection seam. We optimize the S 1 /S 0 MECIs for ethylene, butadiene, and benzene and show that CMS-PDFT gives smooth PESs in the vicinities of the MECIs. Furthermore, the CMS-PDFT MECIs are in good agreement with the MECI calculated by the more expensive XMS-CASPT2 method.
Data-Driven Recommendation of Optimal Tuning Scheme for Range-Separated Hybrid Functionals in Solution-Phase UV/Vis Absorption Energy Prediction
Time-dependent density functional theory (TDDFT) combined with range-separated hybrid (RSH) functionals and a tuned range-separation parameter γ offers a computationally economical approach for high-throughput excited- state property predictions. The γ-tuning procedure in the gas phase is well established. However, no agreement on the best γ- tuning procedure has been made when considering the solvent effect with implicit solvent models like the polarizable continuum model (PCM). To answer that question, this study created a diverse dataset with 937 molecules with experimental solutionphase UV/vis absorption spectra. Three γ-tuning methods, the gasphase γ-tuning (GPγT), the partial vertical γ-tuning (PVγT), and the strict vertical γ-tuning (SVγT), were evaluated for the ωPBEh functional over the entire dataset. Additional benchmarks are done for the optimally tuned screened range-separated hybrid combined with the PCM approach (SRSH-PCM) and the solvation-mediated tuning procedure (sol-med-OT). Our findings revealed that the optimal γ-values obtained by the PVγT and the SVγT are significantly smaller than the GPγT. This trend holds consistently across all molecules in our dataset, and we explained the origin of this phenomenon. TDDFT calculations with PVγTand SVγT-tuned γ-values and default global Fock exchange fraction achieve superior performance compared to those using GPγTtuned or default γ and slightly outperform SRSH-PCM and sol-med-OT with similar or lesser computational cost. Furthermore, we found that the smaller γ-values from SVγT captured the expected 1/(εR) asymptotic behavior in the solution phase, resulting in accurate prediction of solution-phase CT excitations, consistent with the screened asymptote behavior encoded in SRSH-PCM. These results show that SVγT is the best scheme for high-throughput UV/vis absorption spectrum calculations using the ωPBEh functional from a data-driven perspective.
Efficient Pd on carbon catalyst for ammonium formate dehydrogenation: Effect of surface oxygen functional groups
Formate solution is one of a kind promising liquid organic hydrogen carrier but suffers from kinetic challenges. This study investigates the kinetics of the surfacefunctionalized Pd on carbon catalysts for formate dehydrogenation and the impact of O-functional groups. Here, the fraction of the distinguished O-functional groups was modulated by the different concentrated HNO 3 solutions treatment or by H 2 reduction. This study shows that the O-functional groups play critical roles in dispersing Pd nanoparticles and decreasing the activation energy for dehydrogenation. Density functional theory (DFT) calculation reveals that most O-functional groups enhance formate adsorption on the Pd active site. However, the existence of C=O groups consumes reducing agents and hinders the formation of metallic Pd. The electron transfer from Pd to oxygen functional groups is unfavorable to dehydrogenation. The as-prepared Pd5/re-ACA (reduced activated carbon washed by acid) exhibited significant activity with a higher turnover frequency of 13,511 h -1 than commercial Pd/C.
Thermal hybrid exchange-correlation density functional for improving the description of warm dense matter
Finite-temperature density functional theory (FT-DFT) has become a standard tool for first principles calculations of the properties of warm dense matter (WDM) relevant to high-energy density physics (HEDP) applications. In this work we introduce theoretical grounds of thermal hybrid exchange-correlation (XC) functionals within the generalized Mermin-Kohn-Sham scheme for an improved description of WDM. Building on the previously developed KSDT (Karasiev-Sjostrom-Dufty-Trickey) [Karasiev et al., Phys. Rev. Lett. 112, 076403 (2014)] local density approximation (LDA) and the KDT16 (Karasiev-Dufty-Trickey 2016) [Karasiev et al., Phys. Rev. Lett. 120, 076402 (2018)] generalized-gradient approximation (GGA) XC free-energy density functionals, we construct a new thermal hybrid XC functional, referred to here as KDT0. The KDT0 model at low temperature reduces to the popular ground-state PBE0 hybrid due to properties of the used KDT16 density functional approximation. Application to static calculations of electronic band gap and band structure at a wide range of temperatures for various systems of interest to HEDP show that KDT0 provides a significant improvement to the lower LDA and GGA rung XC functionals and to the ground-state PBE0 hybrid.
Prussian blue analogues as platform materials for understanding and developing oxygen evolution reaction electrocatalysts
Transition metal based materials containing Fe have drawn great attention as oxygen evolution reaction (OER) catalysts. The nature of the electrocatalytic active species remains under debate due to the ambiguous physicochemical properties of the catalyst materials, such as the oxidation states and crystal structures. Here, in order to address this issue, transition metal Prussian blue analogues (TM-PBA, Na(TM)(Fe)(CN) 6 , TM = V, Fe, Co, and Ni) with an isomorphous structure are investigated for OER catalysis. Our combined experimental measurements and density functional theory (DFT) calculations reveal that TM-PBAs exhibit volcano-like OER activity with Ni-PBA located near the top of the volcano. Such a volcano-like activity profile can be attributed to the distinctive binding energy difference between *O and *OH on different TM-PBAs surfaces. This research demonstrates that TM-PBAs can be used as platform materials for understanding structure-property-activity relationships in OER catalysts.
Effect of applied potential on metal surfaces: Surface energy, Wulff shape and charge distribution
Here we use grand canonical density functional theory to predict the surface energies, Wulff shapes, charge distributions and catalytically active sites of different metal surfaces under electrochemical conditions. We propose a method for computing surface energies from grand canonical density functional theory (GC-DFT) calculations of periodic slab models and use it to compute the surface energies of the facets of Pt, Cu, and Ag crystals to predict their Wulff shapes under electrochemical conditions. GC-DFT predicts that, for the pure metals studied, solvation only slightly affects the Wulff shape while applied potentials considerably affect the surface energies and corresponding Wulff shapes. We used Bader charge analysis of GC-DFT computed electron densities to investigate the effect of applied potential on the distribution of electron density over the atoms of the surfaces of Pt, Cu, Ag, and the 75–25 Ag-Pt and Au-Ni alloys. This analysis shows that, under an applied potential, the electron density is unevenly distributed over the surface atoms and that the charges of atoms more exposed to solvent are more sensitive to bias. Our results show that the most sensitive atom to bias can be used to identify the most favorable adsorption site and thus, the active sites of electrochemical reactions, which is computationally less demanding than calculating the adsorption energies on all possible adsorption sites.
New reactions of diazene and related species for modelling combustion of amine fuels
Potential energy surfaces for reactions involving N 2 H 2 isomers of diazene (diimide) have been explored using density functional theory, with energies based on coupled-cluster theory. A focus is on processes that create or consume these species, and isomerisation between the E (trans) and Z (cis) forms of HNNH. Furthermore, these include isomerisation and dissociation pathways for HNNH, addition of H atoms to form N 2 H 3 , abstraction by H atoms yielding short-lived NNH, and abstraction reactions of H with N 2 H 3 . Transition state and capture theories are applied for high-pressure-limiting behaviour, while low-pressure and falloff regions are characterised via the methods of Troe and coworkers. Rate constants and thermochemistry are provided to improve models of diamine chemistry, relevant to the combustion of NH 3 especially at high concentrations, high pressures or under reducing conditions. Results indicate that amine radical recombination mainly yields the E HNNH isomer, while H-abstraction from N 2 H 3 results in E HNNH and H 2 NN. However, at elevated temperature E → Z isomerisation becomes competitive, and Z HNNH, being more reactive, acts to enhance the diazene consumption rate.
Impact of Pendent Ammonium Groups on Solubility and Cycling Charge Carrier Performance in Nonaqueous Redox Flow Batteries
The synthesis, characterization, electrochemical performance, and theoretical modeling of two base-metal charge carrier complexes incorporating a pendent quaternary ammonium group, [Ni(bppn-Me3)][BF4], 3′, and [Fe(PyTRENMe)][OTf]3, 4’, are described. Both complexes were produced in high yield and fully characterized using NMR, IR, and UV–vis spectroscopies as well as elemental analysis and single-crystal X-ray crystallography. The solubility of 3′ in acetonitrile showed a 283% improvement over its neutral precursor, whereas the solubility of complex 4’ was effectively unchanged. Cyclic voltammetry indicates an ∼0.1 V positive shift for all waves, with some changes in reversibility depending on the wave. Bulk electrochemical cycling demonstrates that both 3′ and 4’ can utilize the second more negative wave to a degree, whereas 4’ ceases to have a reversible positive wave. Flow cell testing of 3′ and 4’ with Fc as the posolyte reveals little improvement to the cycling performance of 3′ compared with its parent complex, whereas 4’ exhibits reductions in capacity decay when cycling either negative wave. Postcycling CVs indicate that crossover is the likely source of capacity loss in complexes 3, 3′, and 4’ because there is little change in the CV trace. Density functional theory calculations indicate that the ammonium group lowers the HOMO energy in 3′ and 4’, which may impart stability to cycling negative waves while making positive waves less accessible. The incorporation of a positively charged species can improve solubility, stored electron density, and capacity decay depending on the complex, features critical to high energy density redox flow battery performance.
Fermi–Löwdin orbital self-interaction correction of adsorption energies on transition metal ions
Density functional theory (DFT)-based descriptions of the adsorption of small molecules on transition metal ions are prone to self-interaction errors. Here, we show that such errors lead to a large over-estimation of adsorption energies of small molecules on Cu + , Zn + , Zn 2+ , and Mn + in local spin density approximation (LSDA) and Perdew, Burke, Ernzerhof (PBE) generalized gradient approximation calculations compared to reference values computed using the coupled-cluster with single, doubles, and perturbative triple excitations method. These errors are significantly reduced by removing self-interaction using the Perdew–Zunger self-interaction correction (PZ-SIC) in the Fermi–Löwdin Orbital (FLO) SIC framework. In the case of FLO-PBE, typical errors are reduced to less than 0.1 eV. Furthermore, analysis of the results using DFT energies evaluated on self-interaction-corrected densities [DFT(@FLO)] indicates that the density-driven contributions to the FLO-DFT adsorption energy corrections are roughly the same size in DFT = LSDA and PBE, but the total corrections due to removing self-interaction are larger in LSDA.