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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

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↗

DOE OSTI2024

Chemical Environment and Structural Variations in High Entropy Oxide Thin Film Probed with Electron Microscopy

For this work, we employ analytical transmission electron microscopy (TEM) to correlate the structural and chemical environment variations within a stacked epitaxial thin film of the high entropy oxide (HEO) Mg 0.2 Co 0.2 Ni 0.2 Cu 0.2 Zn 0.2 O (J14), with two layers grown at different substrate temperatures (500 and 200 °C) using pulsed laser deposition (PLD). Electron diffraction and atomically resolved STEM imaging reveal the difference in out-of-plane lattice parameters in the stacked thin film, which is further quantified on a larger scale using four-dimensional STEM (4D-STEM). In the layer deposited at a lower temperature, electron energy loss spectroscopy (EELS) mapping indicates drastic changes in the oxidation states and bonding environment for Co ions, and energy-dispersive X-ray spectroscopy (EDX) mapping detects more significant cation deficiency. Ab initio density functional theory (DFT) calculations validate that vacancies on the cation sublattice of J14 result in significant electronic and structural changes. The experimental and computational analyses indicate that low temperatures during film growth result in cation deficiency, an altered chemical environment, and reduced lattice parameters while maintaining a single phase. Our results demonstrate that the complex correlation of configurational entropy, kinetics, and thermodynamics can be utilized for accessing a range of metastable configurations in HEO materials without altering cation proportions, enabling further engineering of functional properties of HEO materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

On analysis and stochastic modeling of the particle kinetic energy equation in particle-laden isotropic turbulent flows

In this study, we analyze three-dimensional particle-laden, isotropic turbulence to develop an understanding of inertial particle dynamics from a kinetic energy perspective. Data trends implying inhomogeneous sampling of the flow by particles are identified and used to support a proposed particle behavior: particles appear to accumulate in regions of low flow kinetic energy over time because they lose kinetic energy and slow down in such regions, ultimately causing them to spend more time there. To elucidate this behavior, we derive a particle kinetic energy equation from the particle momentum equation, which incorporates inertial effects through the Schiller–Naumann drag correlation. Upon extracting fundamental physics from this equation, hypotheses regarding the role of the Stokes number in the temporal change of particle kinetic energy and the previously proposed particle behavior are evaluated using simulation data considering three Stokes numbers. Finally, a Fokker–Planck equation is used to derive the steady-state probability density function of the particle kinetic energy. The model fits the simulation data well and provides a tool for further investigation into understanding preferential concentration, as well as a reduced order model for predicting particle kinetic energy in turbulent flows.

42 ENGINEERING↗

DOE OSTI2025

Dissociation of gas-phase anisole induced by low-energy electron interactions: understanding patterns of aromatic bond cleavage

Abstract In view of elucidating the fragmentation patterns of aromatic systems induced by low-energy electron interactions, dissociative electron attachment (DEA) to gas-phase anisole was performed. Anionic fragments resulting from this DEA process were detected by a quadrupole mass spectrometer, and ion yields of those fragments as a function of incident electron energy were rendered. Our study showed the formation of CH 3 − , HCC − , and OCH 3 − fragments, suggesting that various dissociation channels proceed out of DEA to anisole. We employed density functional theory to compute thermodynamic threshold energies for each potential dissociation channel. Those theoretical calculations supported the prediction that the CH 3 − and OCH 3 − fragments form via mechanisms of single-bond cleavage; the HCC − fragments may form through two-, three-, or four-body dissociation channels that entail hydrogen transfers and the cleavage of multiple aromatic bonds. The experimental resonance energies that form the CH 3 − , HCC − , and OCH 3 − fragments were 6.0 eV, 5.8 and 9.7 eV, and 9.8 eV, respectively. Given the classification of anisole as a monosubstituted aromatic species, our results explain generalizable patterns of electron-mediated dissociation in aromatic systems.

Finley, Jacob (ORCID:0009000606963352)↗

DOE OSTI2021

Ab Initio modeling of Near-Edge EELS spectra for chemisorbed molecules

Abstract Electron energy loss spectroscopy (EELS) has recently been applied to probe chemisorbed molecules on metal nanostructures, but a fundamental understanding of the correlation between these spectra and the electronic structures of the adsorbates has been limited. We report here on the insights afforded by time-dependent density functional theory to decipher the energy loss near edge structure (ELNES) of EELS spectra associated with chemisorption. These first-principles calculations simulate the ELNES-EELS spectra for chemisorbed CO on various facets of Au and Pt. Computational predictions of key signatures such as the ‘red shift’ and reductions in the peak intensity of the 2 π * and 6 σ * peaks, as compared to free CO in the gas phase, are validated in comparison to experimentally collected EELS spectra. These signatures are revealed to arise from changes in the electronic structure in terms of unoccupied density of states associated with the chemisorption process. They are consistent with a Blyholder model that incorporates donation and back-donation of electrons. They are also characteristic of the chemisorption process, such as the choice of metal, site of adsorption and the coverage and distribution of adsorbates. Our simulations thus provide guidelines for the use of ELNES-EELS to characterize the atomic structure and adsorption property of nanostructured surfaces and facilitate the development of advanced nanomaterials for catalytic applications.

Materials Science↗

DOE OSTI2021

Quantum-classical hybrid algorithm for the simulation of all-electron correlation

While chemical systems containing hundreds to thousands of electrons remain beyond the reach of quantum devices, hybrid quantum-classical algorithms present a promising pathway toward a quantum advantage. Hybrid algorithms treat the exponentially scaling part of the calculation-the static correlation-on the quantum computer and the non-exponentially scaling part-the dynamic correlation-on the classical computer. While a variety of algorithms have been proposed, the dependence of many methods on the total wave function limits the development of easy-to-use classical post-processing implementations. Here, we present a novel combination of quantum and classical algorithms, which computes the all-electron energy of a strongly correlated molecular system on the classical computer from the 2-electron reduced density matrix (2-RDM) evaluated on the quantum device. Significantly, we circumvent the wave function in the all-electron calculations by using density matrix methods that only require input of the statically correlated 2-RDM. Although the algorithm is completely general, we test it with two classical density matrix methods, the anti-Hermitian contracted Schrödinger equation (ACSE) and multiconfiguration pair-density functional theories, using the recently developed quantum ACSE method for simulating the statically correlated 2-RDM. Furthermore, we obtain experimental accuracy for the relative energies of all three benzyne isomers and thereby demonstrate the ability of the developed algorithm to achieve chemically relevant and accurate results on noisy intermediate-scale quantum devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗