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
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Kinetic Monte Carlo Simulation of Oxygen Diffusion in Ytterbium Disilicate

Silicon-based ceramic components for next-generation jet turbine engines offer potential weight savings, as well as higher operating temperatures, both of which lead to increased efficiency and lower fuel costs. Silicon carbide (SiC), in particular, offers low density, good strength at high temperatures, and good oxidation resistance in dry air. However, reaction of SiC with high-temperature water vapor, as found in the hot section of jet turbine engines in operation, can cause rapid surface recession, which limits the lifetime of such components. Environmental Barrier Coatings (EBCs) are therefore needed if long component lifetime is to be achieved. Rare earth silicates such as Yb2Si2O7 and Yb2SiO5 have been proposed for such applications; in an effort to better understand diffusion in such materials, we have performed kinetic Monte Carlo (kMC) simulations of oxygen diffusion in Ytterbium disilicate, Yb2- Si2O7. The diffusive process is assumed to take place via the thermally activated hopping of oxygen atoms among oxygen vacancy sites or among interstitial sites. Migration barrier energies are computed using density functional theory (DFT).

Computer Simulation↗
DOE OSTI2024

Atomic structure of different surface terminations of polycrystalline ZnPd

The intermetallic compound ZnPd has been found to have desirable characteristics as a catalyst for the steam reforming of methanol. The understanding of the surface structure of ZnPd is important to optimize its catalytic behavior. However, due to the lack of bulk single-crystal samples and the complexity of characterizing surface properties in the available polycrystalline samples using common experimental techniques, all previous surface science studies of this compound have been performed on surface alloy samples formed through thin-film deposition. In this study, we present findings on the chemical and atomic structure of the surfaces of bulk polycrystalline ZnPd studied by a variety of complementary experimental techniques, including scanning tunneling microscopy (STM), x-ray photoelectron spectroscopy (XPS), low energy electron microscopy (LEEM), photoemission electron microscopy (PEEM), and microspot low-energy electron diffraction ( μ -LEED). These experimental techniques, combined with density functional theory (DFT)-based thermodynamic calculations of surface free energy and detachment kinetics at the step edges, confirm that surfaces terminated by atomic layers composed of both Zn and Pd atoms are more stable than those terminated by only Zn or Pd layers. DFT calculations also demonstrate that the primary contribution to the tunneling current arises from Pd atoms, in agreement with the STM results. The formation of intermetallics at surfaces may contribute to the superior catalyst properties of ZnPd over Zn or Pd elemental counterparts. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2022

Xenon mobility in γ-uranium and uranium–molybdenum alloys

Diffusion in bcc uranium and U–Mo alloys is of great interest because fission gas and other fission products impact the performance of nuclear fuels. We investigate the mobility of xenon and molybdenum in bcc uranium (⁠$γ$-U) and metallic U–Mo alloys by calculating the migration energies of xenon and molybdenum for various local compositions using density functional theory. We also calculate the solute–vacancy binding energies of different solutes to vacancies in bcc uranium. We find that the solute–vacancy binding energy in bcc uranium is significantly higher than it is in other bcc metals (e.g., Fe and W). We also find that the migration energy of molybdenum is substantially higher than the migration energy of xenon, indicating that xenon is much more mobile than molybdenum in bcc uranium. The presence of molybdenum in the nearest-neighbor shell around a xenon atom typically increases the migration energy of xenon, which indicates a reduction of xenon mobility in U–Mo alloys compared to pure bcc uranium.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Polarization consistent dielectric screening in polarizable continuum model calculations of solvation energies

A polarization consistent framework, where dielectric screening is affected consistently in polarizable continuum model (PCM) calculations, is employed for the study of solvation energies. The computational framework combines a screened range-separated-hybrid functional (SRSH) with PCM calculations, SRSH-PCM, where dielectric screening is imposed in both PCM self-consistent reaction field (SCRF) iterations and the electronic structure Hamiltonian. We begin by demonstrating the impact of modifying the Hamiltonian to include such dielectric screening in SCRF iterations by considering the solutions of electrostatically embedded Hartree–Fock (HF) exact exchange equations. Long-range screened HF-PCM calculations are shown to capture properly the linear dependence of gap energy of frontier orbitals on the inverse of the dielectric constant, whereas unscreened HF-PCM orbital energies are fallaciously semi-constant with respect to the dielectric constant and, therefore, inconsistent with the ionization energy gaps. Similar trends affect density functional theory (DFT) calculations that aim to achieve predictive quality. Importantly, the dielectric screened calculations are shown to significantly affect DFT- and HF PCM-based solvation energies, where screened solvation energies are smaller compared to the unscreened values. Importantly, SRSH-PCM, therefore, appears to reduce the tendency of DFT-PCM to overestimate solvation energies, where we find the effect to increase with the dielectric constant and the polarity of the molecular solute, trends that enhance the quality of DFT-PCM calculations of solvation energy. Understanding the relationship of dielectric screening in the Hamiltonian and DFT-PCM calculations can ultimately benefit on-going efforts for the design of predictive and parameter free descriptions of solvation energies.

Chemistry↗
DOE OSTI2021

Energy barriers for Dy and H penetrating graphene on 6 H -SiC(0001) and freestanding bilayer graphene from first-principles calculations

Currently, intercalation of foreign guest atoms into two-dimensional (2D) layered van der Waals materials is an active research area motivated in part by the development of next-generation energy-storage technologies and optoelectronic devices. One such extensively studied 2D material is the graphene-on-SiC system. To realize and control the desired intercalated structures, it is fundamentally important to understand the kinetic process of intercalation. For the intercalation of a guest atom into graphene layers on SiC substrate, a critical kinetic parameter is the energy barrier of a guest atom penetrating the perfect graphene top layer into the gallery under it. However, accurate theoretical calculations for such penetration barriers are unavailable in literature. From our first-principles density functional theory calculations, we obtain the global energy barriers of 3.47 and 1.80 eV for single Dy and H atoms penetrating the graphene top layer on a graphene buffer layer supported by a Si-terminated 6H-SiC(0001) substrate, respectively. Furthermore, for comparison as well as for examining the lateral strain effects, we also obtain the global barriers of 5.05 and 1.50 eV for single Dy and H atoms penetrating freestanding bilayer graphene with a tensile strain of about 8.8% to match our model for supported graphene, as well as the global barriers of 7.21 and 4.18 eV for penetrating unstrained freestanding bilayer graphene, respectively. From corresponding minimum energy paths with multiple energy minima and saddle points, we can also obtain various local energy barriers and the global backward barrier from the graphene gallery back to the top surface.

2D materials↗
DOE OSTI2020

Energy density and spectrum of single-cycle and sub-cycle electromagnetic pulses

Based on the exact solution of the Maxwell equations in the form of a collapsing spherical vector wave specified by an arbitrary function of time, we have calculated the maximum energy density that can be achieved when focusing extremely short pulses of various shapes. It is shown that our earlier formula expressing the maximum energy density in terms of the spectrum parameters for Gaussian quasi-monochromatic pulses is approximately (with an accuracy of ±40 %) valid for the main types of extremely short pulses. (extremely short laser pulses)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Predicting Excitation Energies of Twisted Intramolecular Charge-Transfer States with the Time-Dependent Density Functional Theory: Comparison with Experimental Measurements in the Gas Phase and Solvents Ranging from Hexanes to Acetonitrile

Electronically excited states characterized by intramolecular charge transfer play an essential role in many biological processes and optical devices. The ability to make quantitative ab initio predictions of the relative energetics involved is a challenging yet desirable goal, especially for large molecules in solution. In this work, we present a data set of 61 experimental measurements of absorption and emission processes, both in the gas phase and in solvents representing a broad range of polarities, which involve intramolecular charge transfer mediated by a nonzero, "twisted" dihedral angle between one or more donor and acceptor subunits. Among a variety of density functionals investigated within the framework of linear-response theory, the "optimally tuned" LRC-ωPBE functional, which utilizes a system-specific yet nonempirical procedure to specify the range-separation parameter, emerges as the preferred choice. For the entire set of excitation energies, involving changes in dipole moment ranging from 4 to >20 Debye, the mean signed and absolute errors are 0.02 and 0.18 eV, respectively (compared, e.g., to -0.30 and 0.30 for PBE0, 0.44 and 0.47 for LRC-ωPBEh, 0.83 and 0.83 for ωB97X-V). We analyze the performance of polarizable continuum solvation models available in Q-Chem that partition the solvent response into fast and slow time scales, and clear trends emerge when measurements corresponding to the four small 4-(dimethylamino)benzonitrile (DMABN)-like molecules and a charged species are excluded. We make the case that the large errors found only for small molecules in the gas phase and weak solvents cannot be expected to improve via the optimal tuning procedure, which enforces a condition that is exact only in the well-separated donor-acceptor limit, and present empirical evidence implicating the outsized importance for small donor-acceptor systems of relaxation effects that cannot be accounted for by the linear-response time-dependent density functional theory within the adiabatic approximation. Finally, we demonstrate the utility of the optimally tuned density functional approach by targeting the charge-transfer states of a large biomimetic model system for light-harvesting structures in Photosystem II.

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

The Structural Characterisation and DFT-Aided Interpretation of Vibrational Spectra for Cyclo(l-Cys-d-Cys) Cyclic Dipeptide in a Solid State

Cyclic dipeptides with two intramolecular peptide bonds forming a six-membered 2,5-diketopiperazine ring are gaining significant attention due to their biological and chemical properties. Small changes in the local geometry of such molecules (from cis to trans) can lead to significant structural differences. This work presents the results of a study of cyclo(l-Cys-d-Cys), a dipeptide comprising two cysteine molecules in opposite chiral configurations, with the functional groups situated at both sides of the diketopiperazine ring. X-ray diffraction (XRD) experiment revealed that the molecule crystallises in the $P$-1 space group, which includes the centre of inversion. The IR and Raman vibrational spectra of the molecule were acquired and interpreted in terms of the potential energy distribution (PED) according to the results of density functional theory (DFT) calculations. The DFT-assisted analysis of energy frameworks for the hydrogen bond network within molecular crystals was performed to support the interpretation of X-ray structural data. The optimisation of the computational model based on three-molecule geometry sections from the crystallographic structure, selected to appropriately reflect the intermolecular interactions responsible for the formation of 1D molecular tapes in cyclo(l-Cys-d-Cys) crystal, allowed for better correspondence between theoretical and experimental vibrational spectra. This work can be considered the first complete structural characterisation of cyclo(l-Cys-d-Cys), complemented via vibrational spectroscopy results with full band assignment aided with the use of the DFT method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Toward alcohol synthesis from CO hydrogenation on Cu(111)-supported MoS 2 – predictions from DFT+KMC

In the quest for cheap and efficient catalysts for alcohol synthesis from syngas, a material of interest is single-layer MoS2 owing to its low cost, abundancy, and flexible structure. Because of the inertness of its basal plane, however, it is essential to find ways that make it catalytically active. Herein, by means of density functional theory based calculations of reaction pathways and activation energy barriers and accompanying kinetic Monte Carlo simulations, we show that while S vacancy row structures activate the MoS 2 basal plane, further enhancement of chemical activity and selectivity can be achieved by interfacing the MoS 2 layer with a metallic support. When defect-laden MoS2 is grown on Cu(111), there is not only an increase in the active region (surface area of active sites) but also charge transfer from Cu to MoS 2 , resulting in a shift of the Fermi level such that the frontier states (d orbitals of the exposed Mo atoms) appear close to it, making the MoS 2 /Cu(111) system ready for catalytic activity. Finally, our calculated thermodynamics of reaction pathways lead to the conclusion that the Cu(111) substrate promotes both methanol and ethanol as the products, while kinetic Monte Carlo simulations suggest a high selectivity toward the formation of ethanol.

2D materials↗
DOE OSTI2023

Band Energy Dependence of Defect Formation in the Topological Semimetal Cd3As2

Cadmium Arsenide (Cd3As2) is a prototypical Dirac semimetal that manifests topological properties in a 3D bulk material. In defect-free Cd3As2, the Fermi level lies at a minimum in the density of states at the Dirac point, but experimentally it forms with excess electron carriers and an elevated EF, thereby masking the topological features. To computationally study the self-doping of Cd3As2, we combine density functional theory (DFT) calculations for defect formation energies with quasi-particle self-consistent GW (QSGW) electronic structure calculations. We demonstrate an innate dependence of the point defect formation energies on carrier concentrations and use the QSGW calculated density of states to extrapolate formation energies to arbitrary electron concentrations. This approach allows the quantitative modeling of thermodynamic defect equilibria in topological semimetals and is used to predict how Cd3As2 growth conditions affect the position of EF relative to the Dirac point.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗
DOE OSTI2023

An Extremely Stable and Soluble NH 2 -Substituted Anthraquinone Electrolyte for Aqueous Redox Flow Batteries

Aqueous redox flow batteries require long-term stable redox molecules for electrical energy storage. Anthraquinones, especially ether bond-decorated ones, experience two dominant decomposition pathways, including side-chain loss and anthrone formation. Here, with the aid of DFT calculations, we designed an anthraquinone (3-NH 2 -2-2PEAQ) bearing an ether substituent and a neighboring NH 2 group, which suppresses both of these decomposition mechanisms and exhibits a high solubility of 1.1 M. When paired with ferrocyanide in a full cell, this anthraquinone, at a concentration of 1.0 M, can operate at greater than 1 V with an extremely low capacity fade rate of 0.01%/day and a Coulombic efficiency above 99.8% while cycling for over a month. The synergistic effects of the ether and parent amino substituents are extremely sensitive to the precise substitution pattern of the anthraquinone. This study demonstrates the effectiveness of the judicious use of DFT-based prediction in rapidly identifying promising candidate electrolytes from a large chemical space.

25 ENERGY STORAGE↗
DOE OSTI2020

Perspectives on van der Waals Density Functionals: The Case of TiS 2

The van der Waals interaction is of foundational importance for a wide variety of physical systems. In particular, van der Waals forces lie at the heart of potential device technologies that may be realized from the functional organization of layered two-dimensional (2D) nanomaterials. Furthermore, for intermediate to large-scale applications modeling, van der Waals density functionals have become the de facto choice for first-principles calculations. In particular, the vdW-DF family of functionals have provided a systematic approach to this theoretically challenging problem. While much progress has been made, there remains room for improvement in the microscopic description of vdW forces from these density functionals. In this work, we compute benchmark results for the binding energy and the electronic density response to binding in TiS2 via accurate diffusion quantum Monte Carlo calculations. We compare these benchmark data to results obtained from local, semilocal, and van der Waals functionals. In particular, we gauge the quality of the original vdW-DF/vdW-DF2 functionals, as well as updated variants such as vdW-DF-C09, vdW-DF-optB88, vdW-DF-optB86b, and vdW-DF2-B86R. We find a close relationship between the accuracy of predicted interlayer separation distances and binding energies for TiS 2 , with the vdW-DF-optB88 functional performing very well in terms of both quantities. In general, the more recently developed functionals are systematic improvements over older ones. However, when considering the response of the electron density to binding, we find that local-density approximation (LDA) and PBEsol generally outperform the vdW-DF functionals in describing the interlayer charge accumulation with vdW-DF-C09 variants performing the best overall.

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

Localization in the SCAN meta-generalized gradient approximation functional leading to broken symmetry ground states for graphene and benzene

Density functional theory calculations play a central role in understanding chemical and solid-state systems. Progress depends on density functionals that accurately reproduce both energies, for thermochemistry, and properly describe ground states and other properties that are of interest. The Cr dimer, benzene and graphene are especially important benchmark systems for quantum chemistry and condensed matter physics. The Strongly Constrained and Appropriately Normed (SCAN) functional, which is an advanced meta-generalized gradient approximation functional that significantly improves molecular energies is shown to perform poorly for the Cr dimer. This is connected with its poor performance for itinerant solid-state magnets and is a consequence of over localization of electrons, thus illustrating an analogy between the Cr dimer and itinerant magnets. The Cr dimer is a notoriously difficult system for density functionals. However, we additionally find that SCAN predicts an incorrect symmetry broken ground state for 2D graphene and for the benzene molecule, which is surprising considering that ground states of these are known to be well described even by the simplest local density approximation. We show that SCAN overly favors localized spin polarized states, which is a serious deficiency of this approach. Thus, the challenge of finding density functionals that accurately treat both localized and delocalized electronic systems remains.

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

Dark-exciton energy splitting in monolayer W Se 2 : Insights from time-dependent density functional theory

Here, we present here a formalism based on time-dependent density functional theory to describe characteristics of both intra- and intervalley excitons in semiconductors, the latter of which had remained a challenge. Through the usage of an appropriate exchange-correlation kernel (nanoquanta), we trace the energy difference between the intra- and intervalley dark excitons in monolayer (1L) W⁢ Se 2 to the domination of the exchange part in the exchange-correlation energies of these states. Furthermore, our calculated transition contribution maps establish the momentum resolved weights of the electron-hole excitations in both bright and dark excitons thereby providing a comprehensive understanding of excitonic properties of 1L W⁢ Se 2 . We find that the states consist of hybridized excitations around the corresponding valleys which leads to brightening of the dark excitons, i.e., significantly decreasing their lifetime which is reflected in the PL spectrum. Using many-body perturbation theory, we calculate the phonon contribution to the energy band gap and the linewidths of the excited electrons, holes, and (bright) exciton to find that as the temperature increases the band gap significantly decreases, while the linewidths increase. Our work paves for describing the ultrafast charge dynamics of transition metal dichalcogenide within an ab initio framework.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Effect of heterostructure engineering on electronic structure and transport properties of two-dimensional halide perovskites

Organic-inorganic halide perovskite solar cells have attracted much attention due to their low-cost fabrication, flexibility, and high-power conversion efficiency. More recent efforts show that the reduction from three- to two-dimensions (2D) of organic–inorganic halide perovskites promises an exciting opportunity to tune their electronic properties. Here, we explore the effect of reduced dimensionality and heterostructure engineering on the intrinsic material properties, such as energy stability, bandgap and transport properties of 2D hybrid organic–inorganic halide perovskites using first-principles density functional theory. We show that the energy stability of engineered perovskite heterostructures is significantly enhanced. The heterostructures with improved stability also show excellent transport properties similar to their bulk counterparts. Furthermore, these layered chemistries demonstrate the advantage of a broad range of tunable bandgaps and high-absorption coefficient in the visible spectrum. The proposed 2D heterostructured material holds potential for nano-optoelectronic devices as well as for effective photovoltaics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Unusual aliovalent doping effects on oxygen non-stoichiometry in medium-entropy compositionally complex perovskite oxides

Aliovalent doping of perovskite oxides can tune the oxygen vacancy formation energy. This work discovers normal vs. abnormal aliovalent doping effects on redox behaviors in medium-entropy compositionally complex perovskite oxides (CCPOs) (La 1–x Sr x )(Mn 1/3 Fe 1/3 Ti 1/3 )O 3–δ (LS_MFT) vs. (La 1–x Sr x )(Mn 1/3 Fe 1/3 Cr 1/3 )O 3–δ (LS_MFC). In the LS_MFC series, the oxygen non-stoichiometry range Δδ (= δ red – δ ox ) linearly depends on the Sr molar ratio x, while the LS_MFT series shows a V-shape dependence of Δδ on x. This unusual observation is investigated and explained based on the analysis of the energy loss near edge structure (ELNES) in STEM electron energy loss spectroscopy, along with density functional theory (DFT) calculations. In LS_MFC, Cr–L 2,3 , Mn–L 2,3 and Fe–L 2,3 peaks have a similar linear shift to higher energy with increasing x, which indicates higher oxidation states of Cr, Mn, and Fe with lower oxygen vacancy formation energies. In LS_MFT, the V-shape of the Δδ vs. x curve is caused by the stable Ti 4+ state and a V-shape Mn/Fe valency dependence on x. Furthermore, this study suggests the possible existence of different (including unexpected) coupled aliovalent doping effects in CCPOs with multiple B-site redox active elements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗