Engineering Papers⌕ Search

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.

At least 451 records · Page 25

A Novel NASICON‐Type Na 4 MnCr(PO 4 ) 3 Demonstrating the Energy Density Record of Phosphate Cathodes for Sodium‐Ion Batteries

Abstract Sodium‐ion batteries (SIBs) have attracted incremental attention as a promising candidate for grid‐scale energy‐storage applications. To meet practical requirements, searching for new cathode materials with high energy density is of great importance. Herein, a novel Na superionic conductor (NASICON)‐type Na 4 MnCr(PO 4 ) 3 is developed as a high‐energy cathode for SIBs. The Na 4 MnCr(PO 4 ) 3 nanoparticles homogeneously embedded in a carbon matrix can present an extraordinary reversible capacity of 160.5 mA h g −1 with three‐electron reaction at ≈3.53 V during the Na + extraction/insertion process, realizing an unprecedentedly high energy density of 566.5 Wh kg −1 in the phosphate cathodes for SIBs. It is intriguing to reveal the underlying mechanism of the unique Mn 2+ /Mn 3+ , Mn 3+ /Mn 4+ , and Cr 3+ /Cr 4+ redox couples via X‐ray absorption near‐edge structure spectroscopy. The whole electrochemical reaction undergoes highly reversible single‐phase and biphasic transitions with a moderate volume change of 7.7% through in situ X‐ray diffraction and ex situ high‐energy synchrotron X‐ray diffraction. Combining density functional theory (DFT) calculations with the galvanostatic intermittent titration technique, the superior performance is ascribed to the low ionic‐migration energy barrier and desirable Na‐ion diffusion kinetics. The present work can offer a new insight into the design of multielectron‐reaction cathode materials for SIBs.

Zhang, Jian↗

Self-Assembly of Metallo-Supramolecules with Dissymmetrical Ligands and Characterization by Scanning Tunneling Microscopy

Asymmetrical and dissymmetrical structures are widespread and play a critical role in nature and life systems. In the field of metallo-supramolecular assemblies, it is still in its infancy for constructing artificial architectures using dissymmetrical building blocks. Herein, we report the self-assembly of supramolecular systems based on two dissymmetrical double-layered ligands. With the aid of ultra-high-vacuum, low-temperature scanning tunneling microscopy (UHV-LT-STM), we were able to investigate four isomeric structures corresponding to four types of binding modes of ligand LA with two major conformations complexes A. The distribution of isomers measured by STM and total binding energy of each isomer obtained by density functional theory (DFT) calculations suggested that the most abundant isomer could be the most stable one with highest total binding energy. Lastly, through shortening the linker between inner and outer layers and the length of arms, the arrangement of dissymmetrical ligand LB could be controlled within one binding mode corresponding to the single conformation for complexes B.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Order-disorder versus displacive transitions in Jahn-Teller active layered materials

Large anharmonic vibrations often play a crucial role in dynamically stabilizing crystalline phases whose structures are unstable at low temperature. Although the average structure of such phases can be measured through diffraction experiments, their local structure remains a challenge to characterize and understand. Dynamically stabilized phases are often classified as order/disorder or displacive based on the qualitative nature of their local structure. A robust understanding of how chemistry determines this distinction in behavior, however, is lacking. This article presents a parametric study of an anharmonic vibrational model that describes the transition from a cooperative to a noncooperative Jahn-Teller distortion in layered oxides—a class of materials widely used in Li-ion and Na-ion batteries. The results illustrate how the shape of the energy landscape determines the extent to which the high-temperature phase has order/disorder vs displacive character. In this work, we find that the nature of the high-temperature phase is determined by a competition between the strength of the elastic coupling between Jahn-Teller distortions at nearby sites and the energy scale driving the Jahn-Teller distortion. A comparison of the model to energy landscapes calculated with density-functional theory suggests that the high-temperature phases of the Jahn-Teller active layered compounds LiNiO 2 , NaNiO 2 , LiMnO 2 , and NaMnO 2 exhibit order/disorder character.

36 MATERIALS SCIENCE↗

DFT modelling of explicit solid–solid interfaces in batteries: methods and challenges

Density Functional Theory (DFT) calculations of electrode material properties in high energy density storage devices like lithium batteries have been standard practice for decades. By contrast, DFT modelling of explicit interfaces in batteries arguably lacks universally adopted methodology and needs further conceptual development. In this paper, we focus on solid–solid interfaces, which are ubiquitous not just in all-solid state batteries; liquid-electrolyte-based batteries often rely on thin, solid passivating films on electrode surfaces to function. We use metal anode calculations to illustrate that explicit interface models are critical for elucidating contact potentials, electric fields at interfaces, and kinetic stability with respect to parasitic reactions. The examples emphasize three key challenges: (1) the “dirty” nature of most battery electrode surfaces; (2) voltage calibration and control; and (3) the fact that interfacial structures are governed by kinetics, not thermodynamics. To meet these challenges, developing new computational techniques and importing insights from other electrochemical disciplines will be beneficial.

25 ENERGY STORAGE↗

Facet-Dependent Photodegradation of Methylene Blue by Hematite Nanoplates in Visible Light

The expression of specific crystal facets in different nanostructures is known to play a vital role in determining the sensitivity toward photodegradation of organics, which can generally be ascribed to differences in surface structure and energy. Herein, we report the synthesis of hematite nanoplates with controlled relative exposure of basal (001) and edge (012) facets, enabling us to establish direct correlation between surface structure and the photocatalytic degradation efficiency of methylene blue (MB) in the presence of hydrogen peroxide. MB adsorption experiments showed that the capacity on (001) is about three times larger than on (012), and this is supported by density functional theory calculations suggesting adsorption energies about 1.5 times higher. However, this alone cannot explain the observed MB photodegradation rate, which is around 14.5 times faster on (001) than (012). The enhancement we attribute to a higher availability of photoelectron accepting surface Fe3+ sites on the (001) facet. This facilitates more efficient iron valence cycling and the heterogeneous photo-Fenton reaction yielding MB-oxidizing hydroxyl radicals at the surface. Our findings help establish a rational basis for design and optimization of hematite nanostructures as photocatalysts for environmental remediation.

Zong, Meirong↗

Hydroxyl radical‐driven transformations of bisphenol A and 2,4‐dinitroanisole: Experimental and computational analysis

Abstract This study used experimental and computational analysis to investigate the advanced oxidation of bisphenol A (BPA) and 2,4‐dinitroanisole (DNAN). The pseudo first‐order reaction rate constants depended on the molar peroxide ratio and were between 0.13 and 0.28 min −1 for BPA and between 0.018 and 0.032 min −1 for DNAN. The kinetic differences appear to be due in part to the energy requirements for oxidation, which depended on the reaction mechanism but were typically lower for BPA than they were for DNAN. Density functional theory (DFT) was used to develop transformation pathways that included experimentally‐detected byproducts. The most energetically favored pathway for BPA oxidation begins with the formation of hydroxylated derivatives, while for DNAN, the most energetically favorable degradation pathway begins with the substitution of the methoxy group. Overall, these findings demonstrate the power of combining experimental and computational tools to reveal transformation mechanisms during water treatment. Practitioner Points Advanced oxidation transformations for two emerging water pollutants, bisphenol A and dinitroanisole, was investigated. The observed reaction kinetics depended on molar peroxide ratio in a manner that is in keeping with previous findings. Density functional theory‐based analysis revealed reaction energy requirements and degradation pathways.

Engineering↗

Quantum Chemistry-Informed Active Learning to Accelerate the Design and Discovery of Sustainable Energy Storage Materials

Here we employed Density Functional Theory (DFT) to compute oxidation potentials of 1,400 homobenzylic ether molecules to search for the ideal sustainable redoxmer design. The generated data were used to construct an active learning model based on Bayesian optimization (BO) that targets candidates with desired oxidation potentials utilizing only a minimal number of DFT calculations. The active learning model demonstrated not only significant efficiency improvement over the random selection approach but also robust capability in identifying desired candidates in an untested set of 112,000 homobenzylic ether molecules. Our findings highlight the efficacy of quantum chemistry-informed active learning to accelerate the discovery of materials with desired properties from a vast chemical space.

25 ENERGY STORAGE↗

Enhanced Nanostructure Dynamics on Au(111) with Adsorbed Sulfur due to Au–S Complex Formation

Chemisorbed species can enhance the fluxional dynamics of nanostructured metal surfaces which has implications for applications such as catalysis. Scanning tunneling microscopy studies at room temperature reveal that the presence of adsorbed sulfur (S) greatly enhances the decay rate of 2D Au islands in the vicinity of extended step edges on Au(111). This enhancement is already significant at S coverages, θS, of a few hundredths of a monolayer (ML), and is most pronounced for 0.1–0.3 ML where the decay rate is increased by a factor of around 30. For θ S close to saturation at about 0.6 ML, sulfur induces pitting and reconstruction of the entire surface, and Au islands are stabilized. Enhanced coarsening at lower θS is attributed to the formation and diffusion across terraces of Au–S complexes, particularly AuS 2 and Au4S 4 , with some lesser contribution from Au3S4. This picture is supported by density functional theory analysis of complex formation energies and diffusion barriers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Activation of N 2 on Manganese Nitride-Supported Ni 3 and Fe 3 Clusters and Relevance to Ammonia Formation

Dual-site models were constructed to represent manganese nitride (Mn 4 N)-supported Ni 3 and Fe 3 clusters for NH 3 synthesis. Density functional theory calculations produced an energy barrier of approximately 0.55 eV for N–N bond activation at the interfacial nitrogen vacancy sites (N v ); also, the hydrogenation and removal of interfacial N is promoted by earth-abundant Ni and Fe metals. Steady-state microkinetic modeling revealed that the turnover frequencies of NH 3 production follow an order of Fe 3 @Mn 4 N ≈ Ni 3 @Mn 4 N > Mn 4 N > Fe >> Ni. Moreover, we present clear evidence that, before NH 3 formation, NH migrates from N v onto the metallic sites. Using N binding energy (BE N ) and the transition-state energy of N 2 activation (E TS ) as descriptors, we concluded that the beneficial effects owing to interfacial N v sites are the most pronounced when BE N is either too strong or too weak while E TS is high; otherwise, excessive N v sites may hinder catalyst performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unmasking Static Correlation Error in Hybrid Kohn–Sham Density Functional Theory

To uncover the way Hartree–Fock exchange brings static correlation error into hybrid Kohn–Sham density functional theory, we compare the potential energy curves of four diatomic molecules, namely, H 2 , F 2 , HF, and NaF, using both restricted and unrestricted Kohn–Sham theory. Here, we find that increasing the percentage of Hartree–Fock exchange significantly increases the static correlation error because Hartree–Fock exchange replaces the ability of local exchange in density functionals to account for the localization effects associated with static correlation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Non-Electricity Based Renewable Fuels: Theory and Computation for Solar Thermochemical Hydrogen

Dominated by photovoltaics and wind, current renewable energy sources generate mostly electricity, but 80% of the global final energy consumption occurs in form of fuels. Therefore, direct solar fuel generation would be a major breakthrough for the energy transition. Solar thermochemical hydrogen (STCH) is one of the very few potential routes towards scalable renewable fuels, but currently suffers from lack of an oxide working material that could optimally perform energy conversion within the thermodynamic boundary conditions. Theory and computation can contribute in two distinct ways, through materials search and discovery, but also by providing detailed mechanistic models for specific systems so to advance our understanding of possible design strategies. To enable high-throughput materials screening, we developed a defect graph neural network (dGNN) machine learning approach,[1] which accelerates the prediction of defect formation energies by replacing the tedious density functional theory (DFT) supercell calculations for all possible defect sites. This approach enables high-throughput database screening of oxides, which was integrated with thermodynamic modeling to extract the reduction entropies as additional selection criterion for STCH. Once potential candidate materials are identified, detailed models can guide materials design by predicting performance characteristics. One challenge is to quantitatively predict thermochemical equilibria at high concentrations when the redox active defects start to interact with each other, thereby impeding the formation of additional defects. Introducing a model for the free energy of defect interaction, parametrized on the basis of DFT data, we simulated the complete STCH redox cycle for (Sr,Ce)MnO3 alloys, achieving near-quantitative agreement with experimental data.[2] The analysis of these simulations reveals how defect interactions diminish the reduction entropy and H2 yield, suggesting to include these interactions in design considerations. Finally, we revisit the popular van't Hoff method for analyzing reduction enthalpies and entropies. This method is not ideal, as it involves a temperature-dependent convolution of gas-phase and solid-state entropies, causing uncertainties in the same order of magnitude as the physical quantities of interest. To avoid this problem, we suggest a simple alternative approach which can be applied to experimental and simulated data alike.

first-principles calculations↗

Nonextensive hydrodynamics of boost-invariant plasmas

We use quasiparticle anisotropic hydrodynamics to study the non-conformal and non-extensive dynamics of a system undergoing boost-invariant Bjorken expansion. To introduce nonextensivity, we use an underlying Tsallis distribution with a time-dependent nonextensivity parameter q. By taking moments of the quasiparticle Boltzmann equation in the relaxation-time approximation, we obtain dynamical equations which allow us to determine the time evolution of all microscopic parameters including q. We compare numerical solutions for bulk observables obtained using the nonextensive evolution with results obtained using quasiparticle anisotropic hydrodynamics with a Boltzmann distribution function (q → 1). We show that the evolution of the temperature, pressure ratio, and scaled energy density, are quite insensitive to which distribution function is assumed. However, we find significant differences in the early-time evolution of the bulk pressure which are observed for even small deviations from the Boltzmann distribution function. Finally, we discuss the existence of non-conformal hydrodynamic attractors for the longitudinal and transverse pressures, the bulk and shear viscous corrections, and the nonextensivity parameter q.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Computational methods based on density functional theory for reactions and processes involving electronic spin (Final Technical Report)

This award supports one post-doctoral researcher for 1.5 years. Publications that acknowledge this grant: Refs. 1–14. Refs. 3,10,12 assess current methodology for the evaluation of magnetic exchange couplings in transition metal complexes. In particular, Ref. 10 validates the use of an approximate (non-iterative) Green’s function approach for the calculation of magnetic exchange couplings and will be the foundation for Thrust 2 in this proposal. Refs. 3 and 12 focus on widely used density functional approaches based on the standard energy differences methodology for the particular case of oxo-bridged Fe(III) complexes. Refs. 2,4–7,11 apply current methodologies to problems of practical interest in molecular magnetism. Ref. 13 presents a methodology to explicitly simulate the dynamics of open quantum systems within density functional theory (DFT) calculations based on the Liouville-von Neumann equation of motion for quantum systems driven out-of-equilibrium. Ref. 8 uses non-collinear spin DFT to explain the mechanical behavior of magnetic mono-atomic Pt wires produced in break-junction experiments in the presence of a magnetic field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Elucidation of Bottom-up Growth of CaCO3 Involving Prenucleation Clusters from Structure Predictions and Decomposition of Globally Optimized (CaCO3)n Nanoclusters

Low-energy minima structures for (CaCO3)n, n = 28, are predicted using bottom-up genetic algorithms in conjunction with density functional theory electronic structure calculations, in comparison with the frozen and relaxed top-down clusters generated by cuts from the calcite, vaterite, and aragonite crystal structures. Similarities in structural motifs for the bottom-up and relaxed top-down are revealed using a fragment recognition technique. Fragment energy decomposition analysis shows that the bottom-up and relaxed top-down clusters belong to two classes of amorphous clusters with distinct intracluster energy distributions, despite their structural similarity. The bottom-up clusters with >20 formula units are surface stabilized with negative surface energy densities. In contrast, the top-down clusters are interior stabilized with positive surface energy densities. We prove that the sign of the surface energy density determines whether the nucleation reaction energy as a function of nuclear size has a maximum or a minimum. The surface-stabilized bottom-up clusters are proposed to be a type of prenucleation cluster at the minimum of the nucleation reaction energy. A mechanism for mineralization of CaCO3 involving prenucleation clusters and nonclassical growth pathway is proposed on the basis of our theoretical findings, which is consistent with previous titration experiments.

calcium carbonate, prenucleation cluster, biominer↗

First-Principles Modeling of the Repassivation of Corrosion Resistant Alloys: Part I. O and Cl Adsorption Energy

To understand the mechanism of repassivation, a key process that controls the performance of corrosion resistant alloys (CRAs), it is useful to obtain some fundamental parameters that describe the atomistic processes. Oxygen and chlorine adsorption energies were systematically calculated with density functional theory for Ni-22Cr alloys with seven different solute elements: Co, Cu, Fe, Mn, Mo, Ru and W. All elements other than Cu enhance both O and Cl adsorption to alloy surfaces. Cr has the strongest effect among these elements. Surface configurations that give stronger O adsorption also tend to give stronger Cl adsorption. While this makes it difficult to differentiate the effect of these solute elements on the oxidation/repassivation tendency of alloys, these fundamental properties provide the foundation for higher scale models that could be used to study the oxidation process during corrosion and potentially inform the design of CRAs. Bader charge analysis reveals a strong relationship between the adatom adsorption energy and the Bader charges of both O and Cl adatom. In conclusion, the work function of alloy surfaces is increased by O while decreased by Cl, suggesting that Cl adsorption can have an intrinsic negative impact on corrosion resistance by decreasing the barrier to electron transfer.

36 MATERIALS SCIENCE↗

Defect Stability in CdTe Based on Formation Energies and Migration Barriers

Native point defects are thought to play a key role in CdTe, either as compensation centers in intentionally doped material, as a source of conductivity in nominally undoped material, or as electron–hole recombination centers. Here, the discussion of their concentration and impact has often centered only on formation energies and transition levels. Using hybrid density functional calculations, including the effects of spin–orbit coupling (SOC), we discuss the stability of native point defects in CdTe based on their formation energies and migration barriers. We show that although Cd interstitials are the lowest energy donor defects, they are unstable at room temperature due to a low migration barrier. They are important for maintaining charge neutrality during growth or annealing at high temperatures, but once the material is brought to room temperature, they are not frozen in as often assumed and are expected to anneal out, leaving the other more stable defects to determine the conductivity. Taking this into account in the solution of the charge neutrality equation, we are able to predict the conductivity type and carrier concentrations that are in good agreement with experimental observations.

14 SOLAR ENERGY↗