Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional calculations”

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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Single-Step Selective Oxidation of Methane by Iron-Oxo Species in the Metal–Organic Framework MFU-4l

The direct and selective conversion of methane to methanol can be considered a holy grail for catalysis research. In this work, we study a metal-organic framework known as MFU-4l, modified by design to include highly reactive iron-oxo species for the catalytic C-H bond activation of methane. We investigate the oxidation of methane and the further potential oxidation of the product methanol using N 2 O as an oxygen source and map the potential energy landscape of these reactions using density functional theory calculations. We show that the highest energy barrier encountered during the methane oxidation process is not the C-H bond breaking, but the activation of the iron center by N 2 O. Furthermore, the potential energy landscape for the C-H bond activation exhibits a large, high-energy plateau region instead of a sharp transition state, thus differing from the traditional radical rebound mechanism. Furthermore, this insight offers interesting potential routes to enhance the catalytic activity of the catalyst, to hinder unwanted deactivation pathways, and to reduce the activity towards the over-oxidation of the product.

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

Effect of the Chemical States of Copper on Methanol Decomposition and Oxidation

Here, the decomposition and oxidation reactions of CH 3 OH over metallic Cu(100) and Cu 2 O-covered Cu(100) surfaces are studied using a combination of in-situ ambient-pressure X-ray photoelectron spectroscopy, Auger electron spectroscopy, and density functional theory calculations. We identify the sequential chemical transformation pathways from bond cleavage to the formation of intermediates and final products under operational conditions. Accumulative surface adsorption of CH 3 O species on metallic Cu(100) impedes the decomposition of CH 3 OH. Co-dosing on metallic Cu(100) with low pressures of 1·10 -4 Torr CH 3 OH + 1·10 -4 Torr O 2 results in partial oxidation of CH 3 OH, where the chemisorbed O ads reduces surface sites available for CH 3 O adsorption, decreasing the surface activity for CH 3 OH decomposition. In contrast, the Cu 2 O overlayer formed under the elevated pressures of 0.33 Torr CH 3 OH + 0.66 Torr O 2 promotes the total oxidation of CH 3 OH into the final products of CO 2 and H 2 O, arising from the active reaction 2 between lattice O within Cu 2 O and intermediates of CH 3 O, CH 2 O, HCOO, and CO. Despite the more favorable O-H bond scission, C-O bond scission also occurs to result in surface accumulation of CH x on metallic Cu(100), blocking active sites for decomposition reactions of CH 3 OH and CH 3 O. By comparison, the CH x species on the Cu 2 O-covered Cu(100) undergo oxidation into CO 2 and H 2 O with lattice O in the Cu 2 O overlayer, thereby freeing active sites for the total oxidation of CH 3 OH. These results highlight the distinct roles of metallic Cu and Cu 2 O in the pathways of CH 3 OH decomposition and oxidation reactions, offering practical insights for the design of Cu-based catalysts with tailored reactivity and selectivity.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Influence of Alkali Metal Cations on the Oxygen Reduction Activity of Pt 5 Y and Pt 5 Gd Alloys

Electrolyte species can significantly influence the electrocatalytic performance. In this work, we investigate the impact of alkali metal cations on the oxygen reduction reaction (ORR) on active Pt 5 Gd and Pt 5 Y polycrystalline electrodes. Due to the strain effects, Pt alloys exhibit a higher kinetic current density of ORR than pure Pt electrodes in acidic media. In alkaline solutions, the kinetic current density of ORR for Pt alloys decreases linearly with the decreasing hydration energy in the order of Li + > Na + > K + > Rb + > Cs + , whereas Pt shows the opposite trend. To gain further insights into these experimental results, we conduct complementary density functional theory calculations considering the effects of both electrode surface strain and electrolyte chemistry. The computational results reveal that the different trends in the ORR activity in alkaline media can be explained by the change in the adsorption energy of reaction intermediates with applied surface strain in the presence of alkali metal cations. Our findings provide important insights into the effects of the electrolyte and the strain conditions on the electrocatalytic performance and thus offer valuable guidelines for optimizing Pt-based electrocatalysts.

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

ReaxFF Reactive Force Field for Exploring Electronically Switchable Polarization in Zn 1– x Mg x O Ferroelectric Semiconductors

Cation misfit in traditional ferroelectric crystals offers a new material platform that can drive electronic components toward structural miniaturization and high-density integration, enabling deviation from von-Neumann architectures. Here, we explore ferroelectricity in Zn 1–x Mg x O, a nontraditional ferroelectric material with tunable properties. Using data from density-functional theory calculations, we have developed a ReaxFF reactive force field to explore the ferroelectric properties of Zn 1–x Mg x O and reveal the hysteresis behavior. We discover that ferroelectric switching can be observed at a critical thickness of 10 nm with a residual polarization of ~100 μC/cm 2 . Our analysis indicates that an increase in Mg-substitution correlates with a decrease in the coercive field. We also observe a strong temperature dependence of the coercive field in Zn 1–x Mg x O, with values decreasing as the temperature increases. Additionally, we find that the distribution of Mg atoms significantly impacts the coercive field, with a clustered distribution leading to a substantial increase. In particular, a decrease in coercive field values is observed when Mg atoms are randomly distributed, compared to uniform distribution. Furthermore, leveraging tunable hysteresis behavior offered by varied percentages and distribution of Mg-substitution provides valuable insights into the design of next-generation functional devices and will inspire further investigations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Defect Properties, Anion Ordering, and Photochromic Mechanism in Yttrium Oxyhydride

Yttrium oxyhydride (YHO) undergoes a reversible photochromic transition when exposed to ultraviolet light. However, the mechanism for this transformation is not fully understood, and the structure and precise chemical composition of YHO remain under debate. Here, we use first-principles density functional theory calculations with a hybrid functional to study the structure, chemical stability, and point defect properties of YHO. As experiments have shown, we find that YHO prefers a cubic structure, with H and O anions present in equal concentrations and located on tetrahedral sites. Stoichiometric and ordered YHO is chemically stable, but it has a wide band gap of 5.01 eV, considerably larger than that measured in experiments (2.4–3.8 eV). On the other hand, Y4H10O has a smaller band gap of 2.97 eV and also has a region of chemical stability; thus, the actual material may include some fraction of this H-rich structure. The defect chemistry of YHO is dominated by anionic antisite species (H O and O H ), with hydrogen interstitials (H i ) and vacancies (V H ) also present in reasonably high concentrations. We show that antisite disorder lowers the band gap relative to the perfectly ordered structure, bringing the magnitude of the gap into closer agreement with experiment. Based on our calculations of defect migration and the positions of defect states relative to the band edges, we link the onset of photochromic behavior to the reaction H O – → V O 0 + H i – , which follows photoexcitation of a H O + defect. H i – can subsequently migrate away and be trapped by additional H O + defects, contributing to the persistence of the reaction, while the resultant oxygen vacancy, V O 0 , introduces an occupied defect state that leads to optical absorption at visible wavelengths. Our results can explain reported discrepancies between experimental and computational results for YHO, and they allow us to propose specific atomic-scale processes that can lead to photochromism. In conclusion, understanding these mechanisms is key for unlocking YHO’s application in devices ranging from smart windows and optoelectronics to electrochemical synapses for neural networks.

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

Structural Evolution of Mixed-Addenda Keggin Polyoxometalate Anions with Atom-by-Atom Substitution

Polyoxometalates (POMs) are molecular metal oxides with distinctive electronic properties that make them promising materials for applications in energy, sensors, and memory devices. One of the most promising methods of tuning the stability, photochromic, redox, and electron-spin properties of POMs is through the substitution of the metal “addenda” atoms that, along with oxygen, constitute their cage-like structures. Because traditional synthesis methods typically produce a distribution of POMs, the isolation and characterization of multimetallic POMs with predetermined stoichiometry remains challenging. The presence of multiple energetically accessible isomers further complicates the experimental characterization and theoretical modeling of multimetallic POMs. Herein, we leverage the distinguishing mass-selection capabilities of ion soft landing to prepare stoichiometrically selected Keggin PMo x W 12-x O 40 3- (x = 0 – 6, 8, 10, and 12) POMs on self-assembled monolayer surfaces free of the solvent molecules and counterions that often confound characterization of complex species at interfaces. The structures of the supported POMs are characterized with atom-by-atom precision using in situ infrared (IR) reflection absorption spectroscopy complemented by detailed density functional theory calculations. Our joint experimental and theoretical results reveal an almost linear shift in the positions of the IR bands towards lower wavenumbers with an increase in the number of lighter molybdenum atoms compared to heavier W atoms in PMoxW 12-x O 40 3- . The theoretical calculations also indicate that numerous isomeric structures may be populated at the experimental conditions and, consequently, contribute to the overall IR spectra. In conclusion, our findings indicate that in addition to the number of substituted addenda atoms and the presence of multiple isomeric structures, interactions with the surface play an important role in determining the IR spectra and structure of supported bimetallic POMs.

Prabhakaran, Venkateshkumar [Pacific Northwest Nat↗
DOE OSTI2024

Assessment of Long-Term Degradation of Adsorbents for Direct Air Capture by Ozonolysis

Porous adsorbents are a promising class of materials for the direct air capture of CO 2 (DAC). Practical implementation of adsorption-based DAC requires adsorbents that can be used for thousands of adsorption–desorption cycles without significant degradation. We examined the potential degradation of adsorbents by a mechanism that appears to have not been considered previously, namely, ozonolysis by trace levels of ozone from ambient air. We focused on amine-appended metal–organic frameworks, specifically amine-functionalized Mg 2 (dobpdc), as a representative DAC adsorbent. Estimates based on the number of amine sites in these adsorbents and the ozone concentration in air suggest that degradation by ozone may be relevant over thousands of adsorption–desorption cycles if reactions with adsorbed ozone are fast. We used density functional theory calculations to estimate reaction rates for amine groups and carbon–carbon double bonds in amine-functionalized Mg 2 (dobpdc).

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

Soft X-ray Emission Spectroscopy of Phosphorus Compounds for Energy Conversion and Storage

Phosphorus (P) is a ubiquitous component of materials for energy conversion and storage. Despite the effectiveness of soft X-ray emission spectroscopy (XES) to characterize these complex materials, XES studies of P-containing compounds are still largely missing. In this investigation, the electronic structure and XES signatures from the P L 2,3 -edge are probed for nine solid-state P-containing compounds characterized by various oxidation states and chemical environments: GaP (3−) , InP (3−) , red-P (0) , H 3 P (3+) O 3 , Na 2 H 2 P 2 (4+) O 6 , H 3 P (5+) O 4 , KH 2 P (5+) O 4 , Na 2 HP (5+) O 4 , and InP (5+) O 4 . The XES spectrum of each material exhibits distinct spectral fingerprints, which provide a robust basis for the speciation of complex P-containing samples. Density functional theory calculations shed light on the electronic structure of selected materials and particularly on the effect of phosphorus-oxygen hybridization. This study provides a comprehensive understanding of P L 2,3 XES spectral fingerprints across a wide range of oxidation states and chemical environments. These data sets and insights pave the way for dedicated future work, such as in situ and operando investigations of P-containing compounds in energy-related applications, where speciation of P compounds is needed.

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

Long-Lived Charge-Transfer Excitons in a Graphene–PTCDI–TiOPc Trilayer Heterostructure

Excitation transfer across the interfaces between graphene, perylenetetracarboxylic diimide (PTCDI), and titanyl phthalocyanine (TiOPc) was studied by using transient absorption and photoluminescence spectroscopy. Both photoluminescence quenching and transient absorption measurements confirm the presence of a type-II interface between PTCDI and TiOPc. While the graphene/PTCDI interface is expected to exhibit type-I behavior, transient absorption measurements indicate that only electrons transfer from PTCDI to graphene, with no evidence of hole transfer. Density functional theory calculations reveal significant ground-state electron transfer from graphene to PTCDI, resulting in band bending that prevents excited holes from transferring from PTCDI to graphene. This feature is exploited in a trilayer heterostructure of graphene/PTCDI/TiOPc, where the spatial separation of photoexcited electrons and holes in graphene and TiOPc, respectively, leads to the formation of long-lived photoexcitations with a lifetime of approximately 500 ps. Furthermore, spatially resolved transient absorption measurements reveal the immobile nature of these excitations, confirming that they are charge-transfer excitons rather than free electrons and holes. Furthermore, these results provide valuable insights into the complex interlayer photoexcitation transfer properties and demonstrate precise control over the layer population and the recombination lifetime of photocarriers in such hybrid heterostructures.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Modeling Framework to Predict Melting Dynamics at Microstructural Defects in TNT-HMX High Explosive Composites

Many high explosive (HE) formulations are composite materials whose microstructure is understood to impact functional characteristics. Interfaces are known to mediate the formation of hot spots that control their safety and initiation. Here, to study such processes at molecular scales, we developed all-atom force fields (FFs) for Octol, a prototypical HE formulation comprised of TNT (2,4,6-trinitrotoluene) and HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). We extended a FF for TNT and recasted it in a form that can be readily combined with a well-established FF for HMX. The resulting FF was extensively validated against experimental results and density functional theory calculations. We applied the new combined TNT-HMX FF to predict and rank surface and interface energies, which indicate that there is an energetic driver for coarsening of microstructural grains in TNT-HMX composites. Finally, we assess the impact of several microstructural environments on the dynamic melting of TNT crystal under ultrafast thermal loading. We find that both free surfaces and planar material interfaces are effective nucleation points for TNT melting. However, MD simulations show that TNT crystal is prone to superheating by at least 50 K on subnanosecond time scales and that the degree of superheating is inversely correlated with surface and interface energy. The modeling framework presented here will enable future studies on hot spot formation processes in accident scenarios that are governed by strong coupling between microstructural interfaces, material mechanics, momentum and energy transport, phase transitions, and chemistry.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Determining the Conformation of Supported Complexes Using an 17 O TEDOR-like NMR Experiment

Dynamic nuclear polarization surface-enhanced nuclear magnetic resonance (NMR) spectroscopy has enabled the determination of the three-dimensional configuration of surface sites, in particular supported metal complexes of relevance to single-site heterogeneous catalysis. These approaches have chiefly leveraged the application of NMR double-resonance experiments that either reveal the complex conformation via point-to-point intramolecular distances between spin-labeled atoms or the complex-surface orientation via distances between the spins and the surface plane. Either method typically requires expensive isotope labeling and each reports on different structural features. The application of an experiment that simultaneously reveals both types of distances with chemical resolution would be ideal. Here, in this article, we describe an 17 O{ 1 H} pseudo-3D correlation experiment that achieves this goal. Specifically, Si–O–Si and Si–O–M oxygens are well-resolved by 17 O NMR; therefore, distances can be simultaneously measured radially, between Si– 17 O–M and the 1 H’s of the ligands, and vertically to the Si– 17 O–Si linkages of the silica support. We demonstrate the experiment using supported yttrium and zirconium complexes. Good agreement is obtained when comparing the experimental results to theoretical predictions from density functional theory calculations, highlighting the reliability of this relatively simple experiment.

alkyls↗
DOE OSTI2025

Separation of Linear and Cyclic Siloxanes in Pure Silica Zeolites

We present a computational assessment of pure-silica zeolites for separating linear and cyclic siloxanes. We developed a force field (FF) for pure silica zeolites, when combined with our previously developed FF for siloxanes using standard Lorentz–Berthelot combining rules, shows good agreement with dispersion-corrected density functional theory calculations. We used molecular dynamics simulations to investigate diffusion of siloxanes in pure silica zeolites and identified a pure silica zeolite with the structure code FAU that enables kinetic separation of linear and cyclic siloxanes. FAU allows the diffusion of linear siloxanes (L2–L6) while excluding cyclic siloxanes. D4 siloxane does not diffuse in any of the investigated zeolites, eliminating the potential of pure silica zeolites to achieve equilibrium-based separations of linear and cyclic siloxanes.

adsorption↗
DOE OSTI2026

High Pressure X-ray Diffraction and Equation of State of Hydrazine

Synchrotron X-ray diffraction has been used to investigate the structure and equation of state (EOS) of hydrazine (N 2 H 4 ) up to 54.3 GPa at 298 K. The diffraction patterns could be fit to a monoclinic unit-cell structure and put strong constraints on previously reported phase transitions documented by vibrational spectroscopy over this pressure range. Pressure–volume ( P–V ) data were fit using a Vinet EOS, yielding parameters: V 0 = 45.2 Å 3 /molecule (fixed), K 0 = 11.8(7) GPa, and K 0 ′ = 6.5(2). Previously measured high-pressure vibrational frequency shifts were used to estimate the vibrational free energy and model P–V–T isotherms from 0 to 1200 K. The results of the P–V–T isotherms are compared to existing shock Hugoniot data on hydrazine and 298 K isotherms for assemblages of possible decomposition products. This comparison suggests dissociation at high density under shock loading. Good correspondence was found between the static lattice EOS as calculated by the model and the previously reported EOS as calculated by density functional theory. Finally, these results resolve existing uncertainties about the EOS and crystal symmetry of hydrazine at high pressure and provide valuable baseline information on this important energetic material.

diffraction↗
DOE OSTI2020

Ca Cobaltites as Potential Cathode Materials for Rechargeable Ca-Ion Batteries: Theory and Experiment

Rechargeable Ca-based batteries can potentially achieve higher energy capacity than Li-ion batteries. The development of Ca-ion batteries, however, remains in its infancy, especially due to the challenge of finding cathode materials with high reversible capacity. In this work, we investigate properties of Ca cobaltites as Ca-ion intercalation cathodes, by means of density functional theory calculations validated by synthesis and electrochemical measurements. Computationally, we accessed thermodynamic, diffusion properties, energy capacity, and the voltage profiles for four Ca cobaltite compounds of different stoichiometry and Co oxidation states: Ca 3 Co 2 O 6 , CaCo 2 O 4 , Ca 2 Co 2 O 5 , and [Ca 2 CoO 3 ][CoO 2 ] 1.62 . We found good stability and relatively low migration barriers of some Ca cobaltites during cycling with the layered CaCo 2 O 4 having the lowest Ca migration barrier of 0.7 eV and the highest theoretical capacity. To validate our calculations, we synthesized Ca 3 Co 2 O 6 , CaCo 2 O 4 , and [Ca 2 CoO 3 ][CoO 2 ] 1.62 with the Pechini method and, subsequently, to test the electrochemical extraction of Ca.

25 ENERGY STORAGE↗
DOE OSTI2020

Chemisorption-Driven Roughening of Hydrothermally Grown KTa1–xNbxO3 Nanoparticles

The chemical and physical properties of nanoparticle surfaces have significant effects on their growth processes and the resulting morphology. Hydrothermally grown KNbO3, KTaO3, and KTa1–xNbxO3 were studied to examine the complex relationships between surface composition, phase, chemistry, and energetics and how these may be used to model and thereby control nanoparticle growth mechanisms. Two different composition-dependent growth modes were identified, where one type formed smooth surface facets, while the other resulted in roughened nanoparticle morphologies. Electron microscopy characterization, density functional theory calculations, and mathematical growth models were used to illuminate the role of surface properties and chemisorption on the nanoparticle growth morphology. Surface energy reduction by chemisorption can increase the rate of terrace nucleation, driving the roughening of these lower energy nanoparticle surfaces.

Ly, Tiffany↗
DOE OSTI2020

Acetonitrile Transition Metal Interfaces from First Principles

Acetonitrile is among the most commonly used nonaqueous solvents in catalysis and electrochemistry. Here, we study its interfaces with multiple facets of the metals Ag, Cu, Pt, and Rh using density functional theory calculations; the structures reported shed new light on experimental observations and underscore the importance of solvent–solvent interactions at high coverage. We investigate the relationship of potential of zero charge (PZC) to metal work function, reporting results in agreement with experimental measurements. We develop a model to explain the effects of solvent chemisorption and orientation on the PZC to within a mean absolute deviation of 0.08–0.12 V for all facets studied. Our electrostatic field dependent phase diagram agrees with spectroscopic observations and sheds new light on electrostatic field effects. This work provides new insight into experimental observations on acetonitrile metal interfaces and provides guidance for future studies of acetonitrile and other nonaqueous solvent interfaces with transition metals.

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

Observation of molecular hydrogen produced from bridging hydroxyls on anatase TiO2(101)

Anatase TiO2 is used extensively in a wide range of catalytic and photocatalytic processes and is a promising catalyst for hydrogen production. Here, we show that molecular hydrogen was produced from bridging hydroxyls, HOb, on the (101) surface of single crystal anatase (a-TiO2(101)). Electron bombardment at 30 K produced bridging oxygen vacancies in the surface. Deuterated bridging hydroxyls, DOb, were subsequently formed via dissociation of adsorbed D2O and confirmed by infrared reflection-absorption spectroscopy. During temperature programmed desorption (TPD) spectroscopy, D2 desorption was observed at 520 K. Density functional theory calculations show that both H2 and H2O production from HOb are endothermic at 0 K on a-TiO2(101), but H2 (H2O) desorption is entropically driven above 230 K (800 K). The calculated activation barrier for H2 desorption is 1.40 eV, which is similar to the desorption energy obtained from analysis of the D2 TPD spectra. The H2 desorption likely proceeds in two steps: H atom diffusion on the surface and then recombination. The calculations indicate that diffusion of the hydrogen atoms into the subsurface could also be an important process. We thank Drs. Z. Dohnalek and R. Rousseau for stimulating discussions. This work was supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, & Biosciences and performed in EMSL, a national scientific user facility sponsored by the Department of Energy's Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory (PNNL). PNNL is a multiprogram national laboratory operated for DOE by Battelle. We also acknowledge the WPI computational facilities for providing resources to perform the DFT studies.

Deskins, N A.↗
DOE OSTI2021

Pressure-Induced Enhancement of Thermoelectric Figure of Merit and Structural Phase Transition in TiNiSn

Half-Heusler thermoelectric materials are potential candidates for high thermoelectric efficiency. Here, we report high-pressure thermoelectric and structural property measurements, density functional theory calculations on the half-Heusler material TiNiSn, and an increase of 15% in the relative dimensionless figure of merit, ZT, around 3 GPa. Thermal and electrical properties were measured utilizing a specialized sample cell assembly designed for the Paris–Edinburgh large-volume press to a maximum pressure of 3.5 GPa. High-pressure structural measurements performed up to 50 GPa in a diamond-anvil cell indicated the emergence of a new high-pressure phase around 20 GPa. A first-principles structure search performed using an ab initio random structure search approach identified the high-pressure phase as an orthorhombic type, in good agreement with the experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗