Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Density functional study of the variants of inter-Coulombic decay resonances in the photoionization of Cl@C 60

Inter-Coulombic decay (ICD) resonances in the photoionization of Cl@C 60 endofullerene molecule are calculated using a perturbative density functional theory (DFT) method. This is the first ICD study of an open shell atom in a fullerene cage. Here, three classes of resonances are probed: (i) Cl inner vacancies decaying through C 60 outer continua, (ii) C 60 inner vacancies decaying through Cl outer continua, and (iii) inner vacancies of either system decaying through the continua of Cl-C 60 hybrid levels, the hybrid Auger-ICD resonances. Comparisons with Ar@C 60 results reveal that the properties of hybrid Auger-ICD resonances are affected by the extent of level hybridization.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Insights on carbon dioxide adsorption in a flexible 1-D coordination polymer from in situ X-ray scattering and density functional theory

A combination of in situ small-angle X-ray scattering (SAXS) microstructure characterization over scales ranging from 1 nm to 10 µm and powder X-ray diffraction (XRD) structure characterization under various gas/pressure/temperature conditions with density functional theory (DFT) calculations provides new insights for the CO2 sorption behavior of a one-dimensional porous coordination polymer: catena-bis­(di­benzoyl­methanato)(4,4′-bi­pyridyl)­nickel(II), denoted NiDBM-Bpy. The NiDBM-Bpy chains are held together by van der Waals forces, but the structure of guest-free NiDBM-Bpy is unsolved due to a lack of suitable crystals and high-quality powder XRD patterns. Nevertheless, SAXS and powder XRD can follow microstructural and structural changes as a function of gas pressure, composition and temperature. Both mixed-gas flow and static supercritical CO2 regimes are explored experimentally. DFT calculations are used to model the structural variation associated with XRD changes and hysteresis in the sorption isotherms. XRD and DFT calculations suggest that an orthorhombic Fddd structure with two CO2 per Ni emerges following a transition from the structure with lower molar volume and symmetry that exists without CO2 present.

Allen, Andrew J.↗
DOE OSTI2021

Electronic structure of Alq 3 and Liq using soft X-ray spectroscopy and density functional theory calculation

We show the electronic structure of aluminum tris-8-hydroxyquinoline (Alq 3 ) and 8-hydroquinolatolithium (Liq) was investigated using a combination of X-ray emission spectroscopy (XES), X-ray absorption spectroscopy (XAS), X-ray photoelectron spectroscopy (XPS), and spectral simulation with density functional theory. The chemical bonding states of Alq 3 and Liq were analyzed using XPS core-level spectra. The band gaps of Alq 3 and Liq were measured by combining the XAS and XES results. Additionally, resonant and non-resonant XES were used to measure the density of states for O sites in the molecules.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Extending Orbital-Optimized Density Functional Theory to L-Edge XPS and Beyond: Spin–Orbit Coupling via Nonorthogonal Quasi-Degenerate Perturbation Theory

Quantum mechanical calculations of core electron binding energies (CEBEs) are relevant to interpreting X-ray photoelectron spectroscopy (XPS). Orbital-optimized density functional theory (OO-DFT) accurately predicts K-edge CEBEs but is challenged by the presence of significant spin–orbit coupling (SOC) at L- and higher edges involving inner-shell orbitals with nonzero angular momentum. Here, to extend OO-DFT to L-edges and higher, our method utilizes scalar-relativistic, spin-restricted open-shell OO-DFT to construct a minimal, quasi-degenerate basis of core-hole states corresponding to a chosen inner-shell (e.g., ionizing all six possible 2p spin orbitals). Nonorthogonal configuration interaction (NOCI) is then used to obtain the matrix elements of the full Hamiltonian including SOC in this quasi-degenerate model space of determinants. Using a screened 1-electron SOC operator parametrized with the Dirac-Coulomb-Breit (DCB) Hamiltonian results in doublet splitting (DS) values for third row elements that are nearly in quantitative agreement with experiment. The resulting NOCI eigenvalues are shifted by the average of the (scalar) OO-DFT CEBEs to yield CEBEs (split by SOC) corrected for dynamic correlation. Comparing calculations on gas phase molecules with experimental results establishes that NO-QDPT with the SCAN functional (NO-QDPT/SCAN), using the DCB screened 1-electron SOC operator is accurate to about 0.2 eV for L-edge CEBEs of molecules containing third row atoms. However, this NO-QDPT approach becomes less accurate for fourth-row elements starting in the middle of the 3d transition metal series, with errors increasing as atomic number increases.

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

Hyperfine interactions for small systems including transition-metal elements using self-interaction corrected density-functional theory

The interactions between the electronic magnetic moment and the nuclear spin moment, i.e., magnetic hyperfine (HF) interactions, play an important role in understanding electronic properties of magnetic systems and in realizing platforms for quantum information science applications. We investigate the HF interactions for atomic systems and small molecules, including Ti or Mn, by using Fermi–Löwdin orbital (FLO) based self-interaction corrected (SIC) density-functional theory. We calculate the Fermi contact (FC) and spin-dipole terms for the systems within the local density approximation (LDA) in the FLO-SIC method and compare them with the corresponding values without SIC within the LDA and generalized-gradient approximation (GGA), as well as experimental data. For the moderately heavy atomic systems (atomic number Z ≤ 25), we find that the mean absolute error of the FLO-SIC FC term is about 27 MHz (percentage error is 6.4%), while that of the LDA and GGA results is almost double that. Therefore, in this case, the FLO-SIC results are in better agreement with the experimental data. For the non-transition-metal molecules, the FLO-SIC FC term has the mean absolute error of 68 MHz, which is comparable to both the LDA and GGA results without SIC. For the seven transition-metal-based molecules, the FLO-SIC mean absolute error is 59 MHz, whereas the corresponding LDA and GGA errors are 101 and 82 MHz, respectively. Therefore, for the transition-metal-based molecules, the FLO-SIC FC term agrees better with experiment than the LDA and GGA results. We observe that the FC term from the FLO-SIC calculation is not necessarily larger than that from the LDA or GGA for all the considered systems due to the core spin polarization, in contrast to the expectation that SIC would increase the spin density near atomic nuclei, leading to larger FC terms.

Chemistry↗
DOE OSTI2021

Density Functional Tight-Binding Simulations Reveal the Presence of Surface Defects on the Quartz (101)–Water Interface

Understanding the structure and reactivity of quartz–water interfaces is critical for numerous applications in the geological, environmental, and biological sciences. However, disagreements on the atomic-level structure of the interfaces between experiments and simulations are hampering our ability to predict the surface reactivity. Here, we used density functional tight-binding (DFTB)-based molecular dynamics to simulate a series of quartz (101) surfaces having different types and densities of surface defects in water and compared them with the structures determined by X-ray reflectivity measurements. The DFTB simulations are able to reproduce previous classical and quantum mechanical predictions of the pristine quartz (101)–water interface that disagree with experimental observations. To remedy this situation, a set of defective quartz surfaces having various surface silicon (Si) vacancies were built as indicated by recent experimental studies. We found that the rotation of surface [SiO 4 ] tetrahedra near Si vacancies can lead to outward displacements of Si atoms similar to those observed in the experiments. The presence of additional surface Si vacancies caused inward relaxations of terminal oxygens through the formation of hydrogen bonds. Finally, the overall results indicate that the quartz (101)–water interface may include a mixture of geminal (≡Si–(OH) 2 )- and vicinal (≡Si–OH)-type silanol groups together with the presence of surface Si vacancies.

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

Survey on stochastic distribution systems: A full probability density function control theory with potential applications

Complex systems seen either in general engineering practice or economics are subjected to ever increased uncertainties that are mostly represented as random variables or parameters, and the characteristics of random variables are represented by their probability density functions (PDFs). Controlling their PDFs means to shape their stochastic distributions and in general it would provide a full treatment for system analysis and operational control and optimization. This leads to the development of stochastic distribution control (SDC) systems theory in the past decades, where the original aim of the controller design is to realize a shape control of the distributions of certain random variables in their PDFs sense for some engineering processes. Indeed, once the PDFs of these random variables or parameters are used to describe their distribution characters, the control task is to obtain control signals so that the output PDFs of stochastic systems are made to follow their target PDFs. The subject of SDC was initially originated for non-Gaussian stochastic control systems design but has found a wide spectrum of applications in general systems in terms of data-driven modeling, analysis, signal processing (filtering), data mining via multivariable statistics, decision-making (optimization) for systems subjected to uncertainties and even in economics. In this context, SDC constitutes an effective primer tool for complex system analysis, control and operational optimizations. In this review paper, a detailed survey of the developments on the research of SDC systems will be made together with their wide spectrum applications and future perspectives.

42 ENGINEERING↗
DOE OSTI2025

Density Functional Tight-Binding Enables Tractable Studies of Quantum Plasmonics

Routine investigations of plasmonic phenomena at the quantum level present a formidable computational challenge due to the large system sizes and ultrafast time scales involved. This Feature Article highlights the use of density functional tight-binding (DFTB), particularly its real-time time-dependent formulation (RT-TDDFTB), as a tractable approach to study plasmonic nanostructures from a quantum mechanical purview. We begin by outlining the theoretical framework and limitations of DFTB, emphasizing its efficiency in modeling systems with thousands of atoms over picosecond time scales. Applications of RT-TDDFTB are then explored in the context of optical absorption, nonlinear harmonic generation, and plasmon-mediated photocatalysis. We demonstrate how DFTB can reconcile classical and quantum descriptions of plasmonic behavior, capturing key phenomena such as size-dependent plasmon shifts and plasmon coupling in nanoparticle assemblies. Lastly, we showcase DFTB’s ability to model hot carrier generation and reaction dynamics in plasmon-driven H2 dissociation, underscoring its potential to model photocatalytic processes. Collectively, these studies establish DFTB as a powerful, yet computationally efficient tool to probe the emergent physics of materials at the limits of space and time.

Electrical energy↗
DOE OSTI2023

Towards Quantum-Chemical Level Calculations of SARS-CoV-2 Spike Protein Variants of Concern by First Principles Density Functional Theory

The spike protein (S-protein) is a crucial part of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), with its many domains responsible for binding, fusion, and host cell entry. In this review we use the density functional theory (DFT) calculations to analyze the atomic-scale interactions and investigate the consequences of mutations in S-protein domains. We specifically describe the key amino acids and functions of each domain, which are essential for structural stability as well as recognition and fusion processes with the host cell; in addition, we speculate on how mutations affect these properties. Such unprecedented large-scale ab initio calculations, with up to 5000 atoms in the system, are based on the novel concept of amino acid–amino acid-bond pair unit (AABPU) that allows for an alternative description of proteins, providing valuable information on partial charge, interatomic bonding and hydrogen bond (HB) formation. In general, our results show that the S-protein mutations for different variants foster an increased positive partial charge, alter the interatomic interactions, and disrupt the HB networks. We conclude by outlining a roadmap for future computational research of biomolecular virus-related systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Chemical bond effects in classical site density functional theory of inhomogeneous molecular liquids

Intra-molecular interactions or chemical bonds represent one of the main distinguishing characteristic of molecular fluids. Development of accurate and practical methods to treat these effects is one of the long standing problems in classical site density functional theory (SDFT). One particular instance when these issues become particularly severe is the case of classical interactions potentials with auxiliary sites or dummy atoms. In this situation current SDFT implementations, such as three-dimensional reference interaction site model (3D-RISM), lead to nonphysical results. We re-examine this issue in this work using our recent reformulation of SDFT. We put forward a simple practical solution to this problem, and illustrate its utility for the case of spherical solutes in diatomic liquids.

Chuev, Gennady N.↗
DOE OSTI2022

Regularized inversion of aerosol hygroscopic growth factor probability density function: application to humidity-controlled fast integrated mobility spectrometer measurements

Abstract. Aerosol hygroscopic growth plays an important role in atmospheric particle chemistry and the effects of aerosol on radiation and hence climate. The hygroscopic growth is often characterized by a growth factor probability density function (GF-PDF), where the growth factor is defined as the ratio of the particle size at a specified relative humidity to its dry size. Parametric, least-squares methods are the most widely used algorithms for inverting the GF-PDF from measurements of the humidified tandem differential mobility analyzer (HTDMA) and have been recently applied to the GF-PDF inversion from measurements of the humidity-controlled fast integrated mobility spectrometer (HFIMS). However, these least-squares methods suffer from noise amplification due to the lack of regularization in solving the ill-posed problem, resulting in significant fluctuations in the retrieved GF-PDF and even occasional failures of convergence. In this study, we introduce nonparametric, regularized methods to invert the aerosol GF-PDF and apply them to HFIMS measurements. Based on the HFIMS kernel function, the forward convolution is transformed into a matrix-based form, which facilitates the application of the nonparametric inversion methods with regularizations, including Tikhonov regularization and Twomey's iterative regularization. Inversions of the GF-PDF using the nonparameteric methods with regularization are demonstrated using HFIMS measurements simulated from representative GF-PDFs of ambient aerosols. The characteristics of reconstructed GF-PDFs resulting from different inversion methods, including previously developed least-squares methods, are quantitatively compared. The result shows that Twomey's method generally outperforms other inversion methods. The capabilities of Twomey's method in reconstructing the pre-defined GF-PDFs and recovering the mode parameters are validated.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Interlayer spacing in pillared and grafted MCM-22 type silicas: density functional theory analysis versus experiment

Pillaring of synthetic layered crystalline silicates and aluminosilicates provides a strategy to enhance their adsorption and separation performance, and can facilitate the understanding of such behavior in more complex natural clays. In this study, we perform the first-principles density functional theory calculations for the pillaring of the pure silica polymorph of an MCM-22 type molecular sieve. Starting with a precursor material MCM-22P with fully hydroxylated layers, a pillaring agent, (EtO) 3 SiR, can react with hydroxyl groups (–OH) on adjacent internal surfaces, 2(–OH) + (EtO) 3 SiR + H 2 O → (–O) 2 SiOHR + 3EtOH, to form a pillar bridging these surfaces, or with a single hydroxyl, –OH + (EtO) 3 SiR + 2H 2 O → (–O)Si(OH) 2 R + 3EtOH, grafting to one surface. For computational efficiency, we replace the experimental organic ligand, R, by a methyl group. We find that the interlayer spacing in MCM-22 is reduced by 2.66 Å relative to weakly bound layers in the precursor MCM-22P. Including (–O) 2 SiR bridges for 50% (100%) of the hydroxyl sites in MCM-22P increases the interlayer spacing relative to MCM-22 by 2.52 Å (2.46 Å). For comparison, we also analyze the system where all –OH groups in MCM-22P are replaced by non-bridging grafted (–O)Si(OH) 2 R which results in a smaller interlayer spacing expansion of 2.17 Å relative to MCM-22. Our results for the interlayer spacing in the pillared materials are compatible with experimental observations for a similar MCM-22 type material with low Al content (Si : Al = 51 : 1) of an expansion relative to MCM-22 of roughly 2.8 Å and 2.5 Å from our x-ray diffraction and scanning transmission electron microscopy analyses, respectively. The latter analysis reveals significant variation in individual layer spacings.

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

Atomic-Level Structure of Mesoporous Hexagonal Boron Nitride Determined by High-Resolution Solid-State Multinuclear Magnetic Resonance Spectroscopy and Density Functional Theory Calculations

Mesoporous hexagonal boron nitride (p-BN) has received significant attention over the last decade as a promising candidate for water cleaning/pollutant removal and hydrogen storage applications. In this work, high-resolution solid-state NMR spectroscopy and plane-wave density-functional theory (DFT) calculations are used to obtain an atomic-level description of p-BN. 1 H– 15 N or 1 H– 14 N heteronuclear (HETCOR) correlation experiments recorded with either conventional NMR at room temperature or dynamic nuclear polarization surface-enhanced spectroscopy (DNP-SENS) at ca. 100 K reveal NB 2 H, NBH 2 , NBH 3 + species residing on the edges of BN sheets. Ultra-high field 35.2 T 11 B NMR spectroscopy was used to resolve 11 B NMR signals from BN 3 , BN 2 O x (OH) 1–x (x = 0–1), BNO x (OH) 2–x (x = 0–2), BO x (OH) 3–x (x = 0–3), and BO x (OH) 4–x – (x = 0–4). Importantly, 2D 11 B dipolar double-quantum–single-quantum homonuclear correlation spectra reveal that many pore/defect sites are composed of boron oxide/hydroxide clusters connected to the BN framework through BN 2 O units. 1D and 2D 11 B{ 15 N} HETCOR NMR experiments, in addition to plane-wave DFT calculations of nine different structural models, further confirm the assignment of all NMR signals. The detailed structure determination of the pore and edge/defect sites within p-BN should further enable the rational design and development of next-generation p-BN-based materials. In addition, the techniques outlined here should be applicable to determine structure within other porous and/or boron-based materials.

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

A density functional theory based comparative study of hybrid photoemissions from Cl@C 60 , Br@C 60 and I@C 60

Photoionization from atom-C 60 hybrid levels in halogen endofullerene molecules, Cl@C 60 , Br@C 60 and I@C 60 , are calculated using a linear response density functional method. Both the ordinary electron-configuration where the open shell halogen is at the center of C 60 and the stable configuration after the atom receives an electron from C 60 to form a closed shell anion are considered. Similar ground state hybridization is found for all three systems while, in general, a slight weakening of the effect is noticed after the electron transfer. At lower photon energies, cross sections of the outer hybrid levels attain identical shapes from enhancements driven by the C 60 plasmon resonances, while the higher energy emissions remain distinguishable from the differences in atomic responses. These results further show near insensitivity to the choice of a configuration. The inner hybrid cross sections in general exhibit similar overall structures, although differ in details between molecules. Furthermore, for these states the results significantly differ before and after the electron transfer – a feature that can be useful to experimentally determine the real configuration of the molecules via photoelectron spectroscopy.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Revealing symmetry-broken superconducting configurations by density functional theory

A coherent theory for the superconductivity of both conventional and unconventional superconductors is currently lacking. Here we show that superconductivity arises from the formation of a symmetry-broken superconducting configuration (SCC) due to atomic perturbation of the normal conducting configuration (NCC). This electron–phonon interaction creates straight one-dimensional tunnels (SODTs) for charge density of electrons and/or holes as revealed by the calculations based on density functional theory (DFT). The SODTs act as resistance-free superhighways and are correlated to the Cooper pairs in the Bardeen–Cooper–Schrieffer (BCS) theory. The formation of SODTs implies that the electron–phonon interaction in the BCS theory can be represented by the difference in charge densities between SCC and NCC predicted by DFT. The present work highlights that in conventional superconductors, SODTs are embedded within the bulk materials and are easily destroyed by phonon vibrations, resulting in a low critical superconducting temperature (T C ). Conversely, in unconventional superconductors such as YBa 2 Cu 3 O 7 (YBCO 7 ), SODTs are protected by a layered pontoon structure with very weak bonding to the bulk materials, maintaining SODTs’ stability at higher temperatures and leading to a much higher T C . The present approach is validated for 14 conventional superconductors of 18 pure elements and MgB 2 examined in this work, including the presently predicted superconductivity in Cu, Ag, Au, Sb, Bi, and MgB 2 at 0 K and 0 GPa, and one unconventional superconductor of YBCO 7 . Our discovery indicates that DFT can be a practical tool for predicting superconductors, enabling a systematic search for new superconducting materials in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Density functional theory studies of transition metal carbides and nitrides as electrocatalysts

Transition metal carbides and nitrides are interesting non-precious materials that have been shown to replace or reduce the loading of precious metals for catalyzing several important electrochemical reactions. The purpose of this review is to summarize density functional theory (DFT) studies, describe reaction pathways, identify activity and selectivity descriptors, and present a future outlook in designing carbide and nitride catalysts for the hydrogen evolution reaction (HER), oxygen evolution reaction (OER), oxygen reduction reaction (ORR), nitrogen reduction reaction (N 2 RR), CO 2 reduction reaction (CO 2 RR) and alcohol oxidation reactions. Furthermore, this topic is of high interest to scientific communities working in the field of electrocatalysis and this review should provide theoretical guidance for the rational design of improved carbide and nitride electrocatalysts.

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

Hybrid density functional study of band gap engineering of SrTi O 3 photocatalyst via doping for water splitting

Perovskite SrTiO 3 (STO) is an attractive photocatalyst for solar water splitting but suffers from a limited photoresponse in the visible spectral range due to its wide band gap. By means of hybrid density functional theory calculations, we systematically study engineering its band gap via doping 4$\textit{d}$ and 5$\textit{d}$ transition-metals $\textit{M}$ (M = Zr, Nb, Mo, Tc, Ru, Rh, Pd, Hf, Ta, W, Re, Os, Ir, and Pt) and chalcogen elements $\textit{Y}$ ($\textit{Y}$ = S and Se).We find that transition-metal dopant M either has no effect on the STO band gap or introduces detrimental midgap states except for Pd and Pt that are able to reduce the STO band gap. In contrast, doping S and Se significantly reduces STO’s direct band gap, thus, leading to appreciable optical absorption transitions in the visible spectral range. In this work, our findings provide that Pd-, S-, and Se-doped STO are potential promising photocatalysts for water splitting under visible light irradiation thereby providing insightful theoretical guides for experiments to improve the photocatalytic activity of STO.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Time-dependent density-functional-theory calculations of the nonlocal electron stopping range for inertial confinement fusion applications

Nonlocal electron transport is important for understanding laser-target coupling for laser-direct-drive (LDD) inertial confinement fusion (ICF) simulations. Current models for the nonlocal electron mean free path in radiation-hydrodynamic codes are based on plasma-physics models developed decades ago; improvements are needed to accurately predict the electron conduction in LDD simulations of ICF target implosions. Here we utilized time-dependent density functional theory (TD-DFT) to calculate the electron stopping power (SP) in the so-called conduction-zone plasmas of polystyrene in a wide range of densities and temperatures relevant to LDD. Compared with the modified Lee-More model, the TD-DFT calculations indicated a lower SP and a higher stopping range for nonlocal electrons. We fit these electron SP calculations to obtain a global analytical model for the electron stopping range as a function of plasma conditions and the nonlocal electron kinetic energy. This model was implemented in the one-dimensional radiation-hydrodynamic code LILAC to perform simulations of LDD ICF implosions, which are further compared with simulations by the standard modified Lee-More model. In conclusion, results from these integrated simulations are discussed in terms of the implications of this TD-DFT-based mean-free-path model to ICF simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗