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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Understanding 2p core-level excitons of late transition metals by analysis of mixed-valence copper in a metal–organic framework

The L 2,3 -edge X-ray absorption spectra of late transition metals such as Cu, Ag, and Au exhibit absorption onsets lower in energy for higher oxidation states, which is at odds with the measured spectra of earlier transition metals. Time-dependent density functional theory calculations for Cu 2+ /Cu + reveal a larger 2p core-exciton binding energy for Cu 2+ , overshadowing shifts in single-particle excitation energies with respect to Cu + . Finally, we explore this phenomenon in a Cu + metal–organic framework with ~12% Cu 2+ defects and find that corrections with self-consistent excited-state total energy differences provide accurate XAS peak alignment.

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

Improved Kelbg Potentials for Z > 1 and Application to Carbon Plasmas

In this work, we present a general form for the electron‐ion diffractive potential derived from the quantum pair density matrix and fit to the improved Kelbg potential for atomic numbers up to $Z = 54$. We apply classical molecular dynamics using the improved Kelbg potential for carbon with various forms of the Pauli potential to compute internal energies and pressures for hot, dense plasma conditions. Our results are compared to an equation of state model based on path integral Monte Carlo and density functional theory simulations to examine the extent to which the improved Kelbg potential reproduces the internal energy and pressure of carbon plasmas. The regions of validity for carbon agree generally with those derived previously for hydrogen once pressure ionization effects are incorporated. Based on our carbon results and previously published hydrogen studies, we discuss the general applicability and limitations of these potentials for equation of state studies in warm dense matter and high energy density plasmas.

general physics↗
DOE OSTI2022

Predicting the Cetane Number, Sooting Tendency, and Energy Density of Terpene Fuel Additives

Abstract Discovering renewable fuels and fuel additives is paramount in reducing carbon emissions from internal combustion engines. Terpenes, a group of compounds that can be synthesized from plant matter and microorganisms, have gained significant interest in recent years as promising candidates for fuels/additives. Terpenes are a diverse class of compounds that contain rings and methyl branches, resulting in high energy densities and optimal cold weather behavior. Their variation in bond order, carbon chains, and functional groups lead to varying degrees of soot formation and performance in existing engines. The present work leverages predictive models, namely artificial neural networks, to predict the cetane number (CN), sooting tendency (quantified with yield sooting index, YSI), and energy density (quantified with lower heating value, LHV) of terpenes and hydrogenated terpenes whose sooting propensities were previously determined through experimental means. Predicted sooting propensities of these terpenes are compared with experimental values, and predicted cetane numbers and energy densities are used to comment on the compounds’ ability to act as fuels/additives. Expected prediction errors for CN, YSI, and LHV, defined by blind test set median absolute error, are within 5.56 cetane units, 3.63 yield sooting index units, and 0.77 MJ/kg respectively. Additionally, the present work investigates a variety of correlation/dependence metrics for property-property relationships, furthering our understanding of how combustion-relevant properties are related.

09 BIOMASS FUELS↗
DOE OSTI2022

Size-Dependent Energy of Ni Nanoparticles on Graphene Films on Ni(111) and Adhesion Energetics by Adsorption Calorimetry

The use of carbon supports for late transition-metal nanoparticle catalysts has grown substantially in recent years due to efforts to develop electrocatalysts for clean energy applications and catalysts for new aqueous-phase biomass-related conversions and due to the evolution of new carbon materials with unique properties (e.g., graphene, carbon nanotubes, and so forth). However, much less is known about the bonding energetics of catalytic metal nanoparticles on carbon supports in comparison with oxide supports, which are more common for thermal catalysis. Here, we report the growth morphology and heats of adsorption of Ni vapor deposited onto graphene/Ni(111) at 300 K and 100 K using metal vapor single-crystal adsorption calorimetry and He + low-energy ion scattering (LEIS). These results provide the Ni chemical potential versus particle size, and the Ni/graphene adhesion energy. LEIS intensities suggest that Ni grows as flat-topped face-centered cubic islands with a nearly constant thickness of ~1.5 nm when deposited at 300 K. At 100 K, Ni grows as smaller nanoparticles, well modeled as hemispherical hexagonal close-packed nanoparticles with a density of ~2 × 10 16 particles/m 2 . The Ni chemical potential as a function of average particle diameter in the 0.5 to 4 nm range at 100 K was determined from the heats of Ni gas adsorption. Further, by fitting the measured chemical potential as a function of diameter, we determined an adhesion energy of 3.6 J/m 2 for large Ni particles on graphene/Ni(111). This adhesion energy is in good agreement with previous scanning tunneling microscopy and density functional theory investigations of Ni/graphene/Ni(111).

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

Experimental and Theoretical Studies of Ultrafast Vibrational Energy Transfer Dynamics in Energetic Materials

Energy transfer through anharmonically-coupled vibrations influences the earliest chemical steps in shockwave-induced detonation in energetic materials. A mechanistic description of vibrational energy transfer is therefore necessary to develop predictive models of energetic material behavior. We performed transient broadband infrared spectroscopy on hundreds of femtoseconds to hundreds of picosecond timescales as well as density functional theory and molecular dynamics simulations to investigate the evolution of vibrational energy distribution in thin film samples of pentaerythritol tetranitrate (PETN) , 1,3,5 - trinitroperhydro - 1,3,5 - triazine (RDX) , and 2,4,6 - triamino 1,3,5 - trinitrobenzene (TATB). Experimental results show dynamics on multiple timescales, providing strong evidence for coupled vibrations in these systems, as well as material-dependent evolution on tens to hundreds of picosecond timescales. Theoretical results also reveal pathways and distinct timescales for energy transfer through coupled vibrations in the three investigated materials, providing further insight into the mechanistic underpinnings of energy transfer dynamics in energetic material sensitivity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Reversing the Irreversible: Thermodynamic Stabilization of LiAlH 4 Nanoconfined Within a Nitrogen-Doped Carbon Host

A general problem when designing functional nanomaterials for energy storage is the lack of control over the stability and reactivity of metastable phases. Using the high-capacity hydrogen storage candidate LiAlH4 as an exemplar, we demonstrate an alternative approach to the thermodynamic stabilization of metastable metal hydrides by coordination to nitrogen binding sites within the nanopores of N-doped CMK-3 carbon (NCMK-3). The resulting LiAlH 4 @NCMK-3 material releases H 2 at temperatures as low as 126 °C with full decomposition below 240 °C, bypassing the usual Li 3 AlH 6 intermediate observed in bulk. Moreover, >80% of LiAlH 4 can be regenerated under 100 MPa H 2 , a feat previously thought to be impossible. Nitrogen sites are critical to these improvements, as no reversibility is observed with undoped CMK-3. Density functional theory predicts a drastically reduced Al–H bond dissociation energy and supports the observed change in the reaction pathway. Finally, the calculations also provide a rationale for the solid-state reversibility, which derives from the combined effects of nanoconfinement, Li adatom formation, and charge redistribution between the metal hydride and the host.

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

Aqueous solvation of the chloride ion revisited with density functional theory: impact of correlation and exchange approximations

The specificity of aqueous halide solvation is fundamental to a wide range of bulk and interfacial phenomena spanning from biology to materials science. Halide polarizability is thought to drive the ion specificity, and if so, it is essential to have an accurate description of the electronic properties of halide ions in water. To this end, the solvation of the chloride anion, Cl - has been reinvestigated with state-of-the-art density functional theory. Specifically, the PBE-D3, PBE0-D3, and SCAN functionals have been employed to probe the impact of correlation and exchange approximations. Anticipating the findings, adding exact exchange improves the electronic structure, but simultaneously significantly reduces the Cl - polarizability, resulting in an over-structured Cl–O radial distribution function (RDF) and longer water H-bond lifetimes to Cl - . SCAN does not yield as much improvement in the energetics of Cl - relative to bulk water, but does result in a smaller reduction of the polarizability and thus a less structured Cl–O RDF, which agrees better with experiment. Special consideration is therefore warranted in assessing the impact of exchange on the energy, charge density, and the charge density response when designing and testing hybrid functionals for aqueous halide solvation.

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

Uptake of Pb and the Formation of Mixed (Ba,Pb)SO 4 Monolayers on Barite During Cyclic Exposure to Lead-Containing Sulfuric Acid

Barite (BaSO 4 ) is a common additive in lead-acid batteries, where it acts as a nucleating agent to promote the reversible formation and dissolution of PbSO 4 during battery cycling. However, little is known about the molecular-scale mechanisms the nucleation and cyclic evolution of PbSO 4 over a battery’s lifetime. In this study, we explore the responses of a barite (001) surface to cycles of high and low lead concentrations in the presence of 100 mM sulfuric acid using in-situ atomic force microscopy and high-resolution X-ray reflectivity. Here we find that PbSO 4 epitaxial films readily nucleate on the barite surface, even from solutions that are undersaturated relative to bulk PbSO 4 . Despite this, barite (001) proves to be an ineffective nucleator of bulk PbSO 4 , as multilayer growth is suppressed even in highly supersaturated solutions. Instead, we find evidence that Pb 2+ ions can directly exchange with Ba 2+ to create mixed (Ba,Pb)SO 4 surfaces. These chemically mixed surfaces do not host PbSO4 monolayers as readily as pristine barite, and the original reactivity is not regained until a fresh surface is reestablished by aggressive etching. Our results can be partly explained by traditional models of Stranski-Krastanov (S-K) growth, in which monolayer films are stabilized by a reduction in surface energy, but multilayer growth is inhibited by epitaxial strain. Complementary density functional theory calculations confirmed the basic energetic-terms of traditional S-K models, but also showed evidence form more complex, thickness-dependent, energetics than would be predicted from the S-K models. The experimental results can be better understood by extending the S-K model to consider the formation of mixed surfaces and films, which can have reduced strain and interfacial energies relative to pure films, while also being stabilized by entropy of mixing. These insights into non-stoichiometric heteroepitaxy will enable better predictions of how barite affects PbSO 4 nucleation in battery environments.

25 ENERGY STORAGE↗
DOE OSTI2024

An experimental and computational view of the photoionization of diol–water clusters

In the interstellar medium, diols and other prebiotic molecules adsorb onto icy mantles surrounding dust grains. Water in the ice may affect the reactivity and photoionization of these diols. Ethylene glycol (EG), 1,2-propylene glycol, and 1,3-propylene glycol clusters with water clusters were used as a proxy to study these interactions. The diol–water clusters were generated in a continuous supersonic molecular beam, photoionized by synchrotron-based vacuum ultraviolet light from the Advanced Light Source, and subsequently detected by reflectron time-of-flight mass spectrometry. The appearance energies for the detected clusters were determined from the mass spectra, collected at increasing photon energy. Clusters of both diol fragments and unfragmented diols with water were detected. The lowest energy geometry optimized conformers for the observed EG–water clusters and EG fragment–water clusters have been visualized using density functional theory (DFT), providing insight into hydrogen bonding networks and how these affect fragmentation and appearance energy. As the number of water molecules clustered around EG fragments (m/z 31 and 32) increased, the appearance energy for the cluster decreased, indicating a stabilization by water. This trend was supported by DFT calculations. Fragment clusters from 1,2-propylene glycol exhibited a similar trend, but with a smaller energy decrease, and no trend was observed from 1,3-propylene glycol. In conclusion, we discuss and suggest that the reactivity and photoionization of diols in the presence of water depend on the size of the diol, the location of the hydroxyl group, and the number of waters clustered around the diol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

The initiation and growth of delaminations induced by matrix microcracks in laminated composites

A recent variational mechanics analysis of microcracking damage in cross-ply laminates of the form /(S)/90n/s, where (S) is any orthotropic sublaminate much stiffer than /90n/, has been extended to account for the presence of delaminations emanating from the tips of microcracks in the /90 2n/T sublaminate. The new two-dimensional stress analysis is used to calculate the total strain energy, effective modulus, and longitudinal thermal expansion coefficient for a laminate having microcracks and delaminations. These results are used to calculate the energy release rate for the initiation and growth of a delamination induced by a matrix microcrack. At low crack densities, /(S)/90n/s laminates are expected to fail by microcracking and to show little or no delamination. At some critical crack density, which is a function of laminate structure and material properties, the energy release rate for delamination exceeds that for microcracking and delamination is predicted to dominate over microcracking. A quasi-three-dimensional model is used to predict the propagation of arbitrarily shaped delamination fronts. All predictions agree with experimental observations.

Nairn, J. A.↗
DOE OSTI2021

The role of surfaces in flexoelectricity

We develop the relationship between the strain derivative of the mean-inner potential and surface contributions to flexoelectricity, identifying the true surface-specific component of the flexoelectric response of finite samples. Density functional theory calculations on a range of experimentally observed, low energy SrTiO 3 , MgO, and Si surfaces demonstrate that the mean-inner potential and its contributions to flexoelectricity are sensitive to small differences in surface structure, chemistry, and adsorbates. We also introduce a method to estimate mean-inner potential contributions to flexoelectricity using electron scattering factors and use this approximation to predict total flexoelectric responses for a variety of insulators. Strategies to experimentally disentangle bulk and surface flexoelectric terms are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

An efficient protocol for excited states of large biochromophores

Efficient energy transport in photosynthetic antenna is a long-standing source of inspiration for artificial light harvesting materials. However, characterizing the excited states of the constituent chromophores poses a considerable challenge to mainstream quantum chemical and semiempirical excited state methods due to their size and complexity and the accuracy required to describe small but functionally important changes in their properties. In this paper, we explore an alternative approach to calculating the excited states of large biochromophores, exemplified by a specific method for calculating the Q y transition of bacteriochlorophyll a, which we name Chl-xTB. Using a diagonally dominant approximation to the Casida equation and a bespoke parameterization scheme, Chl-xTB can match time-dependent density functional theory’s accuracy and semiempirical speed for calculating the potential energy surfaces and absorption spectra of chlorophylls. We demonstrate that Chl-xTB (and other prospective realizations of our protocol) can be integrated into multiscale models, including concurrent excitonic and point-charge embedding frameworks, enabling the analysis of biochromophore networks in a native environment. We exploit this capability to probe the low-frequency spectral densities of excitonic energies and interchromophore interactions in the light harvesting antenna protein LH2 (light harvesting complex 2). The impact of low-frequency protein motion on interchromophore coupling and exciton transport has routinely been ignored due to the prohibitive costs of including it in simulations. Our results provide a more rigorous basis for continued use of this approximation by demonstrating that exciton transition energies are unaffected by low-frequency vibrational coupling to exciton interaction energies.

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

Efficient computational design of two-dimensional van der Waals heterostructures: Band alignment, lattice mismatch, and machine learning

Here, we develop a computational database, website applications (web-apps), and machine-learning (ML) models to accelerate the design and discovery of two-dimensional (2D) heterostructures. Using density functional theory (DFT) based lattice parameters and electronic band energies for 674 nonmetallic exfoliable 2D materials, we generate 226 779 possible bilayer heterostructures. We classify these heterostructures into type-I, -II, and -III systems according to Anderson’s rule, which is based on the relative band alignments of the noninteracting monolayers. We find that type II is the most common and type III the least common heterostructure type. We subsequently analyze the chemical trends for each heterostructure type in terms of the Periodic Table of constituent elements. The band alignment data can also be used for identifying photocatalysts and high-work-function 2D metals for contacts. We validate our results by comparing them to experimental data as well as hybrid-functional predictions. Additionally, we carry out DFT calculations of a few selected systems (MoS 2 /WSe 2 , MoS 2 /h-BN, and MoSe 2 /CrI 3 ), to compare the band-alignment description with the predictions from Anderson’s rule. We develop web-apps to enable users to virtually create combinations of 2D materials and predict their properties. Additionally, we use ML tools to predict band-alignment information for 2D materials. The web-apps, tools, and associated data will be distributed through the JARVIS-HETEROSTRUCTURE website. Our analysis, results, and the developed web-apps can be applied to the screening and design applications, such as finding alternative photocatalysts, photodetectors, and high-work-function (WF) 2D-metal contacts.

2-dimensional systems↗
NASA NTRS1988

Pressure-strain-rate events in homogeneous turbulent shear flow

A detailed study of the intercomponent energy transfer processes by the pressure-strain-rate in homogeneous turbulent shear flow is presented. Probability density functions (pdf's) and contour plots of the rapid and slow pressure-strain-rate show that the energy transfer processes are extremely peaky, with high-magnitude events dominating low-magnitude fluctuations, as reflected by very high flatness factors of the pressure-strain-rate. A concept of the energy transfer class was applied to investigate details of the direction as well as magnitude of the energy transfer processes. In incompressible flow, six disjoint energy transfer classes exist. Examination of contours in instantaneous fields, pdf's and weighted pdf's of the pressure-strain-rate indicates that in the low magnitude regions all six classes play an important role, but in the high magnitude regions four classes of transfer processes, dominate. The contribution to the average slow pressure-strain-rate from the high magnitude fluctuations is only 50 percent or less. The relative significance of high and low magnitude transfer events is discussed.

Brasseur, James G.↗
DOE OSTI2021

Ground State Properties of the Wide Band Gap Semiconductor Beryllium Sulfide (BeS)

We report the results from self-consistent calculations of electronic, transport, and bulk properties of beryllium sulfide (BeS) in the zinc-blende phase, and employed an ab-initio local density approximation (LDA) potential and the linear combination of atomic orbitals (LCAO). We obtained the ground state properties of zb-BeS with the Bagayoko, Zhao, and Williams (BZW) computational method, as enhanced by Ekuma and Franklin (BZW-EF). Our findings include the electronic energy bands, the total (DOS) and partial (pDOS) densities of states, electron and hole effective masses, the equilibrium lattice constant, and the bulk modulus. The calculated band structure clearly shows that zb-BeS has an indirect energy band gap of 5.436 eV, from Γ to a point between Γ and X, for an experimental lattice constant of 4.863 Å. This is in excellent agreement with the experiment, unlike the findings of more than 15 previous density functional theory (DFT) calculations that did not perform the generalized minimization of the energy functional, required by the second DFT theorem, which is inherent to the implementation of our BZW-EF method.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Reliable p K a Prediction through Efficient Incorporation of Anharmonicity within the Nuclear–Electronic Orbital Framework

Accurate pK a prediction is critical for understanding chemical reactivity and molecular properties across a wide range of applications. Computational approaches usually invoke a harmonic treatment of the vibrational modes for zero-point energies, as well as thermal and entropic contributions. Herein, we present a general protocol for relative pK a prediction that incorporates the significant anharmonic effects using nuclear–electronic orbital (NEO) theory. This protocol is validated against experimental data for a range of molecules in acetonitrile, including protonated nitrogen bases, nitrophenols, anilines, and diamines, as well as cobalt electrocatalysts. For simple acids, the NEO approach offers only a slight improvement over conventional density functional theory with the standard harmonic vibrational treatment, whereas for hydrogen-bonded acids, the NEO approach offers more significantly improved performance at a comparable computational cost. This accessible methodology provides a practical route for accurate pKa prediction in challenging systems and is extendable to related thermodynamic properties such as hydricities and proton-coupled redox potentials.

Density functional theory↗
DOE OSTI2025

Atomic ionization: sd energy imbalance and Perdew–Zunger self-interaction correction energy penalty in 3d atoms

To accurately describe the energetics of transition metal systems, density functional approximations (DFAs) must provide a balanced description of s- and d- electrons. One measure of this is the sd transfer error, which has previously been defined as E ( 3 d n − 1 4 s 1 ) − E ( 3 d n − 2 4 s 2 ) . Theoretical concerns have been raised about this definition due to its evaluation of excited-state energies using ground-state DFAs. A more serious concern appears to be strong correlation in the 4s 2 configuration. Here, we define a ground-state measure of the sd energy imbalance, based on the errors of s- and d-electron second ionization energies of the 3d atoms, that effectively circumvents the aforementioned problems. We find an improved performance as we move from the local spin density approximation (LSDA) to the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) to the regularized and restored Strongly Constrained and Appropriately Normed (r 2 SCAN) meta-GGA for first-row transition metal atoms. However, we find large (∼2 eV) ground-state sd energy imbalances when applying a Perdew–Zunger 1981 self-interaction correction. This is attributed to an “energy penalty” associated with the noded 3d orbitals. A local scaling of the self-interaction correction to LSDA results in a balance of s- and d-errors.

Science & Technology - Other Topics↗
DOE OSTI2025

Nanoscale and Element-Specific Lattice Temperature Measurements Using Core-Loss Electron Energy-Loss Spectroscopy

Measuring nanoscale local temperatures, particularly in vertically integrated and multicomponent systems, remains challenging. Spectroscopic techniques like X-ray absorption and core-loss electron energy-loss spectroscopy (EELS) are sensitive to lattice temperature, but understanding thermal effects is nontrivial. This work explores the potential for nanoscale and element-specific core-loss thermometry by comparing the Si L2,3 edge’s temperature- dependent redshift against plasmon energy expansion thermometry (PEET) in a scanning TEM. Using density functional theory (DFT), time-dependent DFT, and the Bethe−Salpeter equation, we ab initio model both the Si L 2,3 and plasmon redshift. We find that the core-loss redshift occurs due to bandgap reduction from electron−phonon renormalization. Our results indicate that despite lower core-loss signal intensity compared to plasmon features, core-loss thermometry has key advantages and can be more accurate through standard spectral denoising. Specifically, we show that the Varshni equation easily interprets the core-loss redshift for semiconductors, which avoids plasmon spectral convolution for PEET in complex junctions and interfaces. We also find that core-loss thermometry is more accurate than PEET at modeling thermal lattice expansion in semiconductors, unless the specimen’s temperature-dependent dielectric properties are fully characterized. Furthermore, core-loss thermometry has the potential to measure nanoscale heating in multicomponent materials and stacked interfaces with elemental specificity at length scales smaller than the plasmon’s wave function.

Bethe–Salpeter equation↗