Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional calculation”

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 OSTI2016

Charge density wave order in 1D mirror twin boundaries of single-layer MoSe 2

Here, We provide direct evidence for the existence of isolated, one-dimensional charge density waves at mirror twin boundaries (MTBs) of single-layer semiconducting MoSe 2 . Such MTBs have been previously observed by transmission electron microscopy and have been predicted to be metallic in MoSe 2 and MoS 2 . Our low-temperature scanning tunnelling microscopy/spectroscopy measurements revealed a substantial bandgap of 100 meV opening at the Fermi energy in the otherwise metallic one-dimensional structures. We found a periodic modulation in the density of states along the MTB, with a wavelength of approximately three lattice constants. In addition to mapping the energy-dependent density of states, we determined the atomic structure and bonding of the MTB through simultaneous high-resolution non-contact atomic force microscopy. Density functional theory calculations based on the observed structure reproduced both the gap opening and the spatially resolved density of states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

General microscopic model of magnetoelastic coupling from first principles

Magnetoelastic coupling, i.e., the change of crystal lattice induced by a spin order, is not only scientifically interesting, but also technically important. In this work, we propose a general microscopic model from firstprinciples calculations to describe the magnetoelastic coupling and provide a way to construct the microscopic model from density functional theory calculations. Based on this model, we reveal that there exists a previously unexpected contribution to the electric polarization induced by the spin order in multiferroics due to the combined effects of magnetoelastic coupling and piezoelectric effect. Interestingly and surprisingly, we find that this latticedeformation contribution to the polarization is even larger than that from the pure electronic and ion-displacement contributions in BiFeO 3 . We conclude this model of magnetoelastic coupling can be generally applied to investigate the other magnetoelastic phenomena.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Double-Layer Kagome Metals Pt 3 Tl 2 and Pt 3 In 2

The connectivity and inherent frustration of the kagome lattice can produce interesting electronic structures and behaviors in compounds containing this structural motif. Here we report the properties of Pt 3 X 2 (X = In and Tl) that adopt a double-layer kagome net structure related to that of the topologically nontrivial high-temperature ferromagnet Fe 3 Sn 2 and the density wave hosting compound V 3 Sb 2 . We examined the structural and physical properties of single crystal Pt 3 Tl 2 and polycrystalline Pt 3 In 2 using X-ray and neutron diffraction, magnetic susceptibility, heat capacity, and electrical transport measurements, along with density functional theory calculations of the electronic structure. Our calculations show that Fermi levels lie in pseudogaps in the densities of states with several bands contributing to transport, and this is consistent with our Hall effect, magnetic susceptibility, and heat capacity measurements. Although electronic dispersions, characteristic of simple kagome nets with nearest-neighbor hopping, are not clearly seen, likely due to the extended nature of the Pt 5d states, we do observe moderately large and non-saturating magnetoresistance values and quantum oscillations in the magnetoresistance and magnetization associated with the kagome nets of Pt.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Embedded Cluster Calculations of Water Adsorption on the UO{sub 2} (111) Surface - 16503

The adsorption of water on the (111) surface of UO2 has been investigated with density functional theory calculations using the periodic electrostatic embedded cluster method. In particular, the adsorptions of a single water molecule, molecularly and dissociatively, as well as two water molecules, one molecularly and one dissociatively, have been studied. A strong dependence of the adsorption energy on basis set is found, and good agreement with energies calculated from previous theoretical studies using different methods is obtained. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Benchmarking Density Functional Theory Methods for Efficient Calculations of a Strongly Correlated Li 1– x Ni 1– y O 2−δ System

Transition metal oxides (TMOs), such as LiNiO 2 , are promising candidates for energy storage and electronic devices due to their unique electronic properties, exceptional physical and chemical characteristics, and ability to adopt multiple oxidation states. However, accurately predicting their properties using mean-field density functional theory (DFT) is challenging due to the presence of strongly correlated d-electrons and the complex interplay between their structural, electronic, and magnetic responses. These challenges are further exacerbated by the need to model defects, surfaces, and interfaces, which require computationally efficient, large-scale simulations. To address these issues, we carry out a benchmark study on the Li 1–x NiO 2 system, evaluating the performance of several popular functionals. Our findings demonstrate that combining SCAN functional relaxation with single-step HSE calculations provides a practical and scalable computational strategy. This approach balances accuracy and efficiency, enabling high-throughput simulations of strongly correlated TMOs and improved predictive modeling capability of TMOs for practical applications.

25 ENERGY STORAGE↗
DOE OSTI2019

Lignin-KMC: A Toolkit for Simulating Lignin Biosynthesis

Lignin is a bountiful biopolymer of phenylpropanoid monomers that is critical for plant structure and function. Based on the abundance of lignin in the biosphere and interest in lignin valorization, a more comprehensive understanding of lignin biosynthesis is imperative. We introduce an open-source software toolkit, Lignin-KMC, that combines kinetic Monte Carlo and first-principles calculations of radical coupling events to model lignin biosynthesis in silico. Lignification is simulated using the Gillespie algorithm with rates derived from density functional theory calculations of individual fragment couplings. Using this approach, we confirm experimental findings regarding the impact of lignification conditions on the polymer structure such as (1) the positive correlation between sinapyl alcohol fraction and depolymerization yield and (2) the primarily benzodioxane linked structure of C-lignin. Furthermore, we identify the in planta monolignol supply rate as a possible control mechanism for lignin biosynthesis based on evolutionary stresses. These examples not only highlight the robustness of our modeling framework but also motivate future studies of new lignin types, unexplored monolignol chemistries, and lignin structure predictions, all with an overarching aim of developing a more comprehensive molecular understanding of native lignin, which, in turn, can accelerate the biological and chemistry communities interested in this important biopolymer.

09 BIOMASS FUELS↗
DOE OSTI2025

Defluorination Mechanisms and Real-Time Dynamics of Per- and Polyfluoroalkyl Substances on Electrified Surfaces

Per- and polyfluoroalkyl substances (PFAS) are persistent environmental contaminants found in groundwater sources and a wide variety of consumer products. In recent years, electrochemical approaches for the degradation of these harmful contaminants have garnered a significant amount of attention due to their efficiency and chemical-free modular nature. However, these electrochemical processes occur in open, highly non-equilibrium systems, and a detailed understanding of PFAS degradation mechanisms in these promising technologies is still in its infancy. To shed mechanistic insight into these complex processes, we present the first constant-electrode potential (CEP) quantum calculations of PFAS degradation on electrified surfaces. These advanced CEP calculations provide new mechanistic details about the intricate electronic processes that occur during PFAS degradation in the presence of an electrochemical bias, which cannot be gleaned from conventional density functional theory calculations. We complement our CEP calculations with large-scale ab initio molecular dynamics simulations in the presence of an electrochemical bias to provide time scales for PFAS degradation on electrified surfaces. Taken together, our CEP-based quantum calculations provide critical reaction mechanisms for PFAS degradation in open electrochemical systems, which can be used to prescreen candidate material surfaces and optimal electrochemical conditions for remediating PFAS and other environmental contaminants.

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

Defect levels and self-compensation in iodine-doped CdTe single crystals

This study systematically documents defect levels in n-type iodine-doped cadmium telluride (CdTe:I) crystals as measured by thermoelectric effect spectroscopy, Hall-effect measurements, and photoluminescence and calculated by density functional theory. The primary donor, I Te , was identified with an activation energy of ∼0.05 eV (measured)/∼0.13 eV (calculated). Deep acceptor states, V Cd and I Te –V Cd complex, exhibited activation energies of ∼0.1 eV (measured)/∼0.19 eV (calculated) and ∼0.12 eV (measured)/∼0.24 eV (calculated), respectively. Self-compensation and Fermi level pinning were observed at ∼0.7 eV (measured)/∼1.08 eV (calculated), contributing to high resistivity in as-grown samples. Post-growth Cd annealing effectively removes compensation centers, releases iodine donors, and significantly reduces resistivity, resulting in an ideal n-type solar cell material.

Cd annealing↗
DOE OSTI2017

Formation of solid thorium monoxide at near-ambient conditions as observed by neutron reflectometry and interpreted by screened hybrid functional calculations

Oxidation of a ~1000 Å sputter-deposited thorium thin film at 150 °C in 100 ppm of flowing oxygen in argon produces the long-sought solid form of thorium monoxide. Changes in the scattering length density (SLD) distribution in the film over the 700-min experiment measured by in-situ, dynamic neutron reflectometry (NR) shows the densities, compositions and thickness of the various thorium oxides layers formed. Screened, hybrid density-functional theory calculations of potential thorium oxides aid interpretation, providing atomic-level picture and energetics for understanding oxygen migration. NR provided evidence of the formation of substoichiometric thorium oxide, ThO y (y < 1) at the interface between the unreacted thorium metal and its dioxide overcoat which grows inward, consuming the thorium at a rate of 2.1 Å/min while y increases until reaching 1:1 oxygen-to-thorium. Its presence indicates that kinetically-favored solid-phase ThO can be preferentially generated as a majority phase under the thermodynamically-favored ThO 2 top layer at conditions close to ambient.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Simulated field-modulated x-ray absorption in titania

Here, we present a method to compute the x-ray absorption near-edge structure (XANES) spectra of solid-state transition metal oxides using real-time time-dependent density functional theory, including spin–orbit coupling effects. This was performed on bulk-mimicking anatase titania (TiO 2 ) clusters, which allows for the use of hybrid functionals and atom-centered all electron basis sets. Furthermore, this method was employed to calculate the shifts in the XANES spectra of the Ti L-edge in the presence of applied electric fields to understand how external fields can modify the electronic structure, and how this can be probed using x-ray absorption spectroscopy. Specifically, the onset of t 2g peaks in the Ti L-edge was observed to red shift and the e g peaks were observed to blue shift with increasing fields, attributed to changes in the hybridization of the conduction band (3d) orbitals.

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

Antimony stable isotope fractionation during adsorption onto birnessite: A molecular perspective from X-ray absorption spectroscopy and density functional theory

Sorption of antimony (Sb) onto birnessite significantly influences the fate of Sb in oceanic and terrestrial environments and fractionates Sb isotopes. Nevertheless, little is known about Sb isotopic fractionation during its adsorption on birnessite. Here, in this study, we show the value of Δ 123 Sb adsorbed-aqueous increases from −0.398 to −0.332 ‰ in 1 h and then decreases and stabilizes at −0.384 ‰ in 72 h. The enrichment of the light Sb isotope is predominantly due to the distortion of the octahedral symmetry. X-ray absorption spectroscopy results indicate Sb first forms a double-corner-sharing complex on birnessite and then transforms to a double-edge-sharing complex during adsorption. The optimized bond distances for double-corner-sharing (3.37 Å) and double-edge-sharing (2.90 Å) complexes calculated using density functional theory (DFT) fits well with the structure (3.41 and 3.00 Å) revealed by X-ray absorption spectroscopy, respectively. The fractionation derived from reduced partition function ratios calculated using DFT aligns well with the experimental results. Therefore, the variation in Sb isotopic fractionation during adsorption is attributed to the evolving structure of Sb complexes on birnessite. Our results demonstrate the isotopic fractionation of Sb during adsorption on birnessite and provide a molecular-scale understanding of Sb behavior, contributing to the correct reconstruction of the Sb isotope composition of ancient seawater using ferromanganese crusts and nodules, and efforts to trace Sb migration in epigenetic mining environments.

Adsorption↗
DOE OSTI2021

P–V–T Equation of State of Iridium Up to 80 GPa and 3100 K

In the present study, the high-pressure high-temperature equation of the state of iridium has been determined through a combination of in situ synchrotron X-ray diffraction experiments using laser-heating diamond-anvil cells (up to 48 GPa and 3100 K) and density-functional theory calculations (up to 80 GPa and 3000 K). The melting temperature of iridium at 40 GPa was also determined experimentally as being 4260 (200) K. The results obtained with the two different methods are fully consistent and agree with previous thermal expansion studies performed at ambient pressure. The resulting thermal equation of state can be described using a third-order Birch–Murnaghan formalism with a Berman thermal-expansion model. The present equation of the state of iridium can be used as a reliable primary pressure standard for static experiments up to 80 GPa and 3100 K. A comparison with gold, copper, platinum, niobium, rhenium, tantalum, and osmium is also presented. On top of that, the radial-distribution function of liquid iridium has been determined from experiments and calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Selective Hydrogenation of Furfural Acetone over a Cu Catalyst: Combined Theoretical and Experimental Study

The selective hydrogenation of biomass-derived hydroxymethyl furfural (HMF)-acetone-HMF (HAH) presents an alternative route to producing highly functionalized polyesters and polyurethanes. While HAH undergoes furan ring hydrogenation over Pd, Ru, and Ni catalysts, furan ring hydrogenation is not observed over Cu. Herein, we combine reaction kinetics experiments and density functional theory calculations to elucidate the selective hydrogenation behavior of HAH over Cu catalysts. We identified furfural acetone (FAc) as a suitable surrogate for modeling HAH hydrogenation over Cu surfaces and performed reaction kinetics experiments between temperatures of 313–393 K and a H 2 partial pressure of 55 bar. Similar to the behavior of HAH, hydrogenation of FAc follows a consecutive two-step hydrogenation pathway over Cu and does not undergo furan ring hydrogenation. The apparent activation energy barriers for hydrogenation of the aliphatic double bond (0.58 eV) and carbonyl (0.43 eV) of FAc measured in a continuous flow reactor setup are consistent with those reported in prior studies for batch HAH hydrogenation. Reaction orders with respect to each reactant, including H 2 and FAc, were determined to be nearly one. Density functional theory (DFT; GGA-PBE-D3) calculations on Cu(111) showed that the hydrogenation of the aliphatic double bond of FAc is more facile than the hydrogenation of the furan ring, which displays weak interactions with the Cu surface. We determined an apparent activation energy barrier for FAc hydrogenation (0.58 eV) that agreed with our DFT predictions (highest barrier for FAc hydrogenation of 0.57 eV). Finally, our DFT calculations further show that weak interactions between the furan ring and Cu surface are responsible for the selective hydrogenation behavior.

Biomass↗
DOE OSTI2014

Local versus global electronic properties of chalcopyrite alloys: X-ray absorption spectroscopy and ab initio calculations

Element-specific unoccupied electronic states of Cu(In, Ga)S{sub 2} were studied as a function of the In/Ga ratio by combining X-ray absorption spectroscopy with density functional theory calculations. The S absorption edge shifts with changing In/Ga ratio as expected from the variation of the band gap. In contrast, the cation edge positions are largely independent of composition despite the changing band gap. This unexpected behavior is well reproduced by our calculations and originates from the dependence of the electronic states on the local atomic environment. The changing band gap arises from a changing spatial average of these localized states with changing alloy composition.

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

Machine learning-enhanced design of lead-free halide perovskite materials using density functional theory

The investigation of emerging non-toxic perovskite materials has been undertaken to advance the fabrication of environmentally sustainable lead-free perovskite solar cells. This study introduces a machine learning methodology aimed at predicting innovative halide perovskite materials that hold promise for use in photovoltaic applications. The seven newly predicted materials are as follows: CsMnCl 4 , Rb 3 Mn 2 Cl 9 , Rb 4 MnCl 6 , Rb 3 MnCl 5 , RbMn 2 Cl 7 , RbMn 4 Cl 9 , and CsIn 2 Cl 7 . The predicted compounds are first screened using a machine learning approach, and their validity is subsequently verified through density functional theory calculations. CsMnCl 4 is notable among them, displaying a bandgap of 1.37 eV, falling within the Shockley-Queisser limit, making it suitable for photovoltaic applications. Here, through the integration of machine learning and density functional theory, this study presents a methodology that is more effective and thorough for the discovery and design of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Ab initio density functional theory study of non-polar (101{sup ¯}0), (112{sup ¯}0) and semipolar (202{sup ¯}1) GaN surfaces

The atomic structures of non-polar GaN(101{sup ¯}0), (112{sup ¯}0) and semipolar GaN(202{sup ¯}1), (202{sup ¯}1{sup ¯}) surfaces were studied using ab initio calculations within density functional theory. The bulk-like truncated (1 × 1) structure with buckled Ga-N or Ga-Ga dimers was found stable on the non-polar GaN(101{sup ¯}0) surface in agreement with previous works. Ga-N heterodimers were found energetically stable on the GaN(112{sup ¯}0)-(1 × 1) surface. The formation of vacancies and substitution site defects was found unfavorable for non-polar GaN surfaces. Semipolar GaN(202{sup ¯}1)-(1 × 1) surface unit cells consist of non-polar (101{sup ¯}0) and semipolar (101{sup ¯}1) nano-facets. The (101{sup ¯}1) nano-facets consist of two-fold coordinated atoms, which form N-N dimers within a (2 × 1) surface unit cell on a GaN(202{sup ¯}1) surface. Dimers are not formed on the GaN(202{sup ¯}1{sup ¯}) surface. The stability of the surfaces with single (101{sup ¯}0) or (101{sup ¯}1) nano-facets was analyzed. A single non-polar (101{sup ¯}0)-(1 × 1) nano-facet was found stable on the GaN(202{sup ¯}1) surface, but unstable on the GaN(202{sup ¯}1{sup ¯}) surface. A single (101{sup ¯}1) nano-facet was found unstable. Semipolar GaN surfaces with (202{sup ¯}1) and (202{sup ¯}1{sup ¯}) polarity can be stabilized with a Ga overlayer at Ga-rich experimental conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Theoretical estimate of the sensitivity of the CUORE detector to solar axions

In this paper we calculate the potential sensitivity of the CUORE detector to axions produced in the Sun through the Primakoff process and detected by the inverse coherent Bragg-Primakoff process. The conversion rate is calculated using density functional theory for the electron density and realistic expectations for the energy resolution and background of CUORE. Monte Carlo calculations for 5 y × 741 kg = 3705 kg y of exposure are analyzed using time correlation of individual events with the theoretical time-dependent counting rate and lead to an expected limit on the axion-photon coupling g{sub aγγ}<3.83 × 10{sup −10} GeV{sup −1} for axion masses less than 100 eV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Determination of the individual atomic site contribution to the electronic structure of 3,4,9,10-perylene-tetracarboxylic-dianhydride (PTCDA)

We have studied the element and orbital-specific electronic structure of thin films of 3,4,9,10-perylene-tetracarboxylic-dianhydride (PTCDA) using a combination of synchrotron radiation-exited resonant x-ray emission spectroscopy, x-ray absorption spectroscopy, x-ray photoelectron spectroscopy, as well as density functional theory calculations. Resonant and non-resonant x-ray emission spectroscopies were used to measure the C and O 2p partial densities of state in PTCDA. Furthermore, resonant x-ray emission at the C and O K-edges is shown to be able to measure the partial densities of states associated with individual atomic sites. The flat molecular orientation of PTCDA on various substrates is explained in terms of the carbonyl O atom acting as a hydrogen-bond acceptor leading to multiple in-plane intermolecular C=O···H−C hydrogen bonding between carbonyl groups and the perylene core of the neighboring PTCDA molecules. We support this conclusion by comparison of our calculations to measurements of the electronic structure using element-, site-, and orbital-selective C and O K-edge resonant x-ray emission spectroscopy, and photoemission spectroscopy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗