Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Energetic material”

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,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Oxidative Deboronation of Boronic Acids by Hydrogen Peroxide in Planta Generates Borate for Cross-Linking of Rhamnogalacturonan II

Vascular plants require boron to cross-link the rhamnogalacturonan-II (RG-II) domain of pectin to form functional cell walls. Boronic acids, which form reversible esters with cis-diols like borate, have been proposed to influence RG-II cross-linking, though the mechanism remains unclear. We used suspension-cultured rose cells adapted to grow without boron to investigate the effect of boronic acids on RG-II dimerization. When grown with phenylboronic acid (PBA) as the sole boron source, nearly all RG-II was crosslinked, whereas methylboronic acid (MBA) only partially restored cross-linking. In contrast, in vitro assays showed that homogeneous RG-II monomers did not dimerize with alkyl or aryl boronic acids unless supplemented with hydrogen peroxide (H2O2), which oxidatively converts boronic acids to boric acid. Real-time NMR spectroscopy and density functional theory calculations provided insight into the reaction mechanism and energetics of oxidation respectively. Together, our data show that exogenous boronic acids are a source of boric acid for plants, and that the deboronation reaction generates aryl or alkyl alcohol byproducts that can undergo further chemical modification in planta. The fate and potential roles of these byproducts in planta remain to be determined.

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

Influence of nanoscale interfaces on the dynamic deformation and spall failure of Cu–Fe alloy microstructures

Additively manufactured immiscible Cu–Fe alloys can exhibit a distribution of nanoscale interfaces due to the distribution of nanoscale clusters in the equilibrium and metastable phases. Molecular dynamics simulations investigate the role of such interfaces on the phase stability and transformation behavior during shock compression, as well as the mechanisms of damage nucleation during spall failure. The model multiphase Cu–Fe systems studied here comprise a distribution of Fe clusters in an FCC Cu matrix, as well as Cu clusters in a BCC Fe matrix. The length scales of the nanoscale clusters determine the energetics of the interfaces that can result in FCC/BCC (equilibrium) or FCC/FCC (metastable) interfaces in the Cu system matrix, and BCC/FCC (equilibrium) or BCC/BCC (metastable) interfaces in the Fe matrix. The MD simulations demonstrate that nanoscale metastable interface microstructures can induce plastic deformation and also stabilize phases or suppress phase transformations in metastable phase clusters. In contrast, equilibrium interfaces can influence phase transformation thresholds and serve as additional void nucleation sites during failure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Investigation of the distribution of laser-induced damage precursors in etched grating sidewalls fabricated using inductively coupled plasma reactive-ion etching

Inductively coupled plasma reactive-ion etching for grating fabrication is explored with emphasis on the type and distribution of process-induced defects that can act as precursors to laser-induced damage initiation. To provide flexibility in the ability to perform a variety of material characterization methods, grating-like features with large pitch were fabricated in silica monolayers. The results indicate that the etch process introduces surface fluorine and subsurface carbon-containing contamination, which reduces the laser-induced damage threshold. The local relative concentrations of the laser-damage precursors introduced during the etch step of the manufacturing process were investigated using angle-resolved damage testing, taking advantage of the shift of the location of the peak electric-field intensity incident on the grating sidewall. The results suggest an increased level of etch-related defects near the base of the sidewall, with directly fluorinated silicon identified as the primary defect.

Dent, R. [Univ. of Rochester, NY (United States). ↗
NASA NTRS1996

X-Ray Imaging Study

The space environment in which the Space Station Freedom and other space platforms will orbit is truly a hostile environment. For example, the currently estimated integral fluence for electrons above 1 Mev at 2000 nautical miles is above 2 x 1O(exp 10) electrons/sq cm/day and the proton integral fluence is above 1 x 10(exp 9) protons/sq cm/day. At the 200 - 400 nautical miles, which is more representative of the altitude which will provide the environment for the Space Station, each of these fluences will be proportionally less; however, the data indicates that the radiation environment will obviously have an effect on structural materials exposed to the environment for long durations. The effects of this combined environment is the issue which needs to be understood for the long term exposure of structures in space. At the same time, there will be substantial potential for collisions between the space platforms and space debris. The current NASA catalogue contains over 4500 objects floating in space which are not considered payloads. This debris can have significant effects on collision with orbiting spacecraft. In order to better understand the effect of these hostile phenomena on spacecraft, several types of studies are being performed to simulate at some level the effect of the environment. In particular the study of debris clouds produced by hypervelocity impact on the various surfaces anticipated on the Space Station is very important at this point in time. The need to assess the threat of such debris clouds on space structures is an on-going activity. The Space Debris Impact facility in Building 4612 provides a test facility to monitor the types of damage produced with hypervelocity impact. These facilities are used to simulate space environmental effects from energetic particles. Flash radiography or x-ray imaging has traditionally provided such information and as such has been an important tool for recording damage in situ with the event. The proper operation of the system can provide much useful information with respect to parametric analysis of the hypervelocity experiment. The following report outlines the procedures developed to optimize the operation of the x-ray imaging system and its operational characteristics.

OBrien, Susan K.↗
NASA NTRS2006

Cleaning Carbon Nanotubes by Use of Mild Oxygen Plasmas

Experiments have shown that it is feasible to use oxygen radicals (specifically, monatomic oxygen) from mild oxygen plasmas to remove organic contaminants and chemical fabrication residues from the surfaces of carbon nanotubes (CNTs) and metal/CNT interfaces. A capability for such cleaning is essential to the manufacture of reproducible CNT-based electronic devices. The use of oxygen radicals to clean surfaces of other materials is fairly well established. However, previously, cleaning of CNTs and of graphite by use of oxygen plasmas had not been attempted because both of these forms of carbon were known to be vulnerable to destruction by oxygen plasmas. The key to success of the present technique is, apparently, to ensure that the plasma is mild . that is to say, that the kinetic and internal energies of the oxygen radicals in the plasma are as low as possible. The plasma oxygen-radical source used in the experiments was a commercial one marketed for use in removing hydrocarbons and other organic contaminants from vacuum systems and from electron microscopes and other objects placed inside vacuum systems. In use, the source is installed in a vacuum system and air is leaked into the system at such a rate as to maintain a background pressure of .0.56 torr (.75 Pa). In the source, oxygen from the air is decomposed into monatomic oxygen by radio-frequency excitation of a resonance of the O2 molecule (N2 is not affected). Hence, what is produced is a mild (non-energetic) oxygen plasma. The oxygen radicals are transported along with the air molecules in the flow created by the vacuum pump. In the experiments, exposure to the oxygen plasma in this system was shown to remove organic contaminants and chemical fabrication residues from several specimens. Many high-magnification scanning electron microscope (SEM) images of CNTs were taken before and after exposure to the oxygen plasma. As in the example shown in the figure, none of these images showed evidence of degradation of CNT structures.

Petkov, Mihail↗
DOE OSTI2026

Unraveling interphase-driven failure pathways in LiMn0.6Fe0.4PO4/graphite pouch cells

LiMnxFe1−xPO4 (LMFP) is a promising high-voltage, thermally stable, and earth-abundant cathode material, yet its practical application is limited by interphase instability and Mn dissolution. In this work, we systematically evaluate LiMn0.6Fe0.4PO4/graphite pouch cells using three electrolyte formulations including control carbonate electrolyte, control + 2 wt% vinylene carbonate (VC), and control + 2 wt% VC + 1 wt% 1,3,2-dioxathiolane 2,2-dioxide (DTD), to establish how electrolyte composition governs interphase chemistry and long-term degradation. Electrochemical testing shows that both additives are preferentially reduced prior to ethylene carbonate (EC) during cell formation, generating robust cathode-electrolyte interphase (CEI) and solid-electrolyte interphase (SEI) layers that suppress gas evolution and raise the first-cycle coulombic efficiency to 89.3%. Additionally, the dual-additive electrolyte delivers the most stable performance, retaining over 85% capacity after 600 cycles while minimizing impedance growth under long-term cycling at C/3 and 40 °C. Soft X-ray absorption spectroscopy confirms that VC + DTD effectively suppresses electrolyte oxidation at the cathode surface, and micro-X-ray fluorescence shows substantially reduced Mn dissolution and deposition on the graphite anode. Density functional theory simulations further provided insights into the structural and energetic influences of alkoxide species on the cathode surface, proposing a Mn2+ extraction mechanism. The combined experimental and computational findings establish a mechanistic link between electrolyte composition and interphase evolution, highlighting the effectiveness of electrolyte engineering for extending the operational lifetime of LMFP-based lithium-ion batteries.

Chak, Chanmonirath Michael↗
DOE OSTI2022

Neutron vibrational spectroscopic evidence for short H∙∙∙H contacts in the RNiInH 1.4; 1.6 (R = Ce, La) metal hydride

Intermetallic metal hydrides are critical materials for hydrogen storage applications, however, metal hydrides with greater storage capacities are still needed. Within metal hydrides, the volumetric storage capacities are limited by the number of hydrogen-accommodating interstitial sites which can be simultaneously occupied given a minimum hydride nearest-neighbor distance of ~2.1 A, according to the Switendick-Westlake criterion. To date, violations of this criterion are rare. Perhaps the most well studied compounds violating this criterion are the RNiInHx compounds (R = Ce, La, Nd). Previous neutron diffraction studies on the deuterated species revealed the presence of Ni-D∙∙∙D-Ni-D∙∙∙D-Ni chains with anomalously close D∙∙∙D contacts of ~1.6 A. Yet there are no neutron vibrational spectroscopic investigations reported for these atypical hydrides. Here we use neutron vibrational spectroscopy (NVS) measurements to probe the hydrogen dynamics in LaNiInHx (x = 0.67, 1.6) and CeNiInH 1.4 . For x > 0.67, the presence of close H...H contacts yields two related features in the vibrational spectrum centered near ~90 meV corresponding to the oscillations of paired H atoms simultaneously occupying neighboring R 3 Ni tetrahedra. Notably, these features are energetically distinct from comparable vibrational motions for "unpaired" H atoms when x = 0.67. To compare, we also present powder neutron diffraction and NVS measurements for the newly characterized, chemically similar Sn compounds CeNiSnH, CeNiSnH 2 , and CeNiSnD 2 . These compounds also contain R 3 Ni tetrahedra, however, the H-occupied tetrahedra are well separated from each other with the closest H∙∙∙H distances exceeding 2.1 A, and the Switendick-Westlake criterion is not violated. Consequently, the spectral signature of the close H∙∙∙H contacts is absent in these hydrides.

08 HYDROGEN↗
DOE OSTI2022

Interfacial Unit-Dependent Catalytic Activity for CO Oxidation over Cerium Oxysulfate Cluster Assemblies

Atomically precise cerium oxo clusters offer a platform to investigate structure–property relationships that are much more complex in the ill-defined bulk material cerium dioxide. We investigated the activity of the MCe 70 torus family (M = Cd, Ce, Co, Cu, Fe, Ni, and Zn), a family of discrete oxysulfate-based Ce 70 rings linked by monomeric cation units, for CO oxidation. CuCe 70 emerged as the best performing MCe 70 catalyst among those tested, prompting our exploration of the role of the interfacial unit on catalytic activity. Temperature-programmed reduction (TPR) studies of the catalysts indicated a lower temperature reduction in CuCe 70 as compared to CeCe 70 . In situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) indicated that CuCe 70 exhibited a faster formation of Ce 3+ and contained CO bridging sites absent in CeCe 70 . Isothermal CO adsorption measurements demonstrated a greater uptake of CO by CuCe 70 as compared to CeCe 70 . The calculated energies for the formation of a single oxygen defect in the structure significantly decreased with the presence of Cu at the linkage site as opposed to Ce. Furthermore, this study revealed that atomic-level changes in the interfacial unit can change the reducibility, CO binding/uptake, and oxygen vacancy defect formation energetics in the MCe 70 family to thus tune their catalytic activity.

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

Interfacial Reactivity and Speciation Emerging from Na-Montmorillonite Interactions with Water and Formic Acid at 200 °C: Insights from Reactive Molecular Dynamics Simulations, Infrared Spectroscopy, and X-ray Scattering Measurements

Reactive organic fluid - mineral interactions at elevated temperatures contribute to the evolution of planetary matter. One of the less studied but important transformations in this regard involves the reactions of formic acid with naturally occurring clays such as sodium montmorillonite. To advance a mechanistic understanding of these interactions, we use ReaxFF reactive molecular dynamics simulations in conjunction with infrared (IR) spectroscopy and X-ray scattering experiments to investigate the speciation behavior of water-formic acid mixtures on sodium montmorillonite interfaces at 473 K and 1 atm. Using a newly developed reactive forcefield, we show that the experimental IR spectra of unreacted and reacted mixture can be accurately reproduced by ReaxFF/MD. We further benchmark the simulation predictions of sodium carbonate and bicarbonate formation in the clay interlayers using Small and Wide-Angle X-ray Scattering measurements. Subsequently, leveraging the benchmarked forcefield, we interrogate the pathway of speciation reactions with emphasis on carbonate, formate, and hydroxide groups elucidating the energetics, transition states, intermediates, and preferred products. Further, we also delineate the differences in reactivities and catalytic effects of clay edges, facets, and interlayers owing to their local chemical environments, which have far reaching consequences in their speciation behavior. The experimental and simulation approaches described in this study and the transferable forcefields can be applied translationally to advance the science of clay-fluid interactions for several applications including subsurface fluid storage and recovery and clay-pollutant dynamics.

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

Defect modeling in semiconductors: the role of first principles simulations and machine learning

Abstract Point defects in semiconductors dictate their electronic and optical properties. Vacancies, interstitials, substitutional defects, and defect complexes can form in the semiconductor lattice and significantly impact its performance in applications such as solar absorption, light emission, electronics, and catalysis. Understanding the nature and energetics of point defects is essential for the design and optimization of next-generation semiconductor technologies. Here, we provide a comprehensive overview of the current state of research on point defects in semiconductors, focusing on the application of density functional theory (DFT) and machine learning (ML) in accelerating the prediction and understanding of defect properties. DFT has been instrumental in accurately calculating defect formation energies, charge transition levels, and other defect-related properties such as carrier recombination rates and lifetimes, and ion migration barriers. ML techniques, particularly neural networks, have emerged as powerful tools for enabling rapid prediction of defect properties at DFT-accuracy in order to overcome the expense of using large supercells and advanced functionals. We begin this article with a discussion of different types of point defects and complexes, their impact on semiconductor properties, and the experimental and DFT approaches typically used for their characterization. Through multiple case studies, we explore how DFT has been successfully applied to understand defect behavior across a variety of semiconductors, and how ML approaches integrated with DFT can efficiently predict defect properties and facilitate the discovery of new materials with tailored defect behavior. Overall, the advent of ‘DFT+ML’ promises to drive advancements in semiconductor technology, catalysis, and renewable energy applications, paving the way for the development of high-performance semiconductors which are defect-tolerant or have desirable dopability.

Rahman, Md Habibur (ORCID:000000027705984X)↗
DOE OSTI2026

A cell-centered AMR-ALE framework for 3D multi-material hydrodynamics. Part I: Lagrangian and indirect Euler AMR algorithms

Many applications of physics and engineering involve wide ranges of time and spatial scales. The numerical simulation of localized small scales such as shock waves and material interfaces requires a large number of computational cells in these regions. For these applications, Lagrangian and Arbitrary-Lagrangian-Eulerian (ALE) related methods are engaging since the moving mesh feature naturally brings mesh cells on shock discontinuities and material interfaces are carefully captured. In addition, Adaptive-Mesh-Refinement (AMR) strategies aim to optimize computational resources by concentrating finer mesh cells only in areas of interest while using coarser cells elsewhere. A key but challenging AMR requirement consists in efficiently distributing the computational effort to achieve high accuracy without the prohibitive computational costs associated with uniformly fine grids. Here, in this document, the coupling of the p4est AMR library with a cell-centered Lagrangian scheme is presented with the goal to perform reliable 3D Lagrangian-AMR and indirect Euler-AMR multi-material simulations. In particular, it is shown that starting from a 3D indirect ALE code, the memory management and load balancing requirements can be delegated to an external library (here the p4est library) to unlock ALE-AMR capabilities. First, we present a strategy to transcribe the octant-based connectivity of the 3D AMR framework with that of an unstructured mesh of polygonal cells used in Lagrangian hydrodynamics. Then, we show how refinement and coarsening operations must be adapted to the particular Lagrangian framework to ensure the conservation of volume during those steps. Finally, several numerical test cases are presented that demonstrate the capabilities of the Lagrangian-AMR and indirect Euler-AMR algorithms.

3D cell-centered Lagrangian numerical scheme↗
DOE OSTI2025

Magnetic Phase Diagram of 2D Magnet Fe 5– x GeTe 2 Revealed by Cryogenic Lorentz 4D-STEM

A two-dimensional van der Waals ferromagnet Fe 5–x GeTe 2 (F5GT) exhibits room-temperature magnetic transition and magnetic anisotropy. Studies have identified diverse magnetic states, including stripe domains, skyrmionic type-I and topologically trivial type-II bubbles, depending on the magnetic and thermal history. Yet, the underlying micromagnetic energetics driving these states remains unclear. Here, we establish the magnetic phase diagram of bulk F5GT using cryogenic Lorentz four-dimensional scanning transmission electron microscopy. We reveal that stripe domains spontaneously form upon zero-field cooling, while metastable magnetic bubbles are stabilized by external fields. Transitions from type-I to type-II bubbles are driven by an oblique external field. Micromagnetic simulations confirm that these transitions arise from the interplay of uniaxial anisotropy, dipolar interactions, and external fields, without requiring significant Dzyaloshinskii–Moriya Interaction (DMI). In conclusion, these findings clarify the micromagnetic origin of spin textures of F5GT and establish it as a tunable platform for nanoscale topological magnetism.

2D Magnets↗
DOE OSTI2022

First-Principles Elucidation of Initial Dehydrogenation Pathways in Mg(BH 4 ) 2

Complex borohydrides such as Mg(BH 4 ) 2 offer one of highest capacities to chemically store hydrogen for onboard applications; however, it suffers greatly from kinetic constraints that prevent realization of full capacity and reversibility. Understanding these kinetic limitations solely from experiments is extremely challenging due to the unusual complexity of various competing elemental reaction steps involved during the de/rehydrogenation reaction. This work aims to map out the energetics associated with initial dehydrogenation of Mg(BH 4 ) 2 from first-principles simulations and to identify the preferred reaction pathways. Our calculations suggest the rate-limiting step during BH 4 – –B 3 H 8 – conversion is the formation of the B 2 H 7 – intermediate. We further emphasize and clarify that the B 3 H 8 – and H – intermediates, formed during initial Mg(BH 4 ) 2 decomposition, appear as molecular species that are embedded in the Mg–BH 4 –Mg matrix as evidenced in the nuclear magnetic resonance measurements and not as bulk MgH 2 and Mg(B 3 H 8 ) 2 as previously assumed in theoretical predictions of the thermodynamics.

08 HYDROGEN↗
DOE OSTI2022

Data-driven design of novel halide perovskite alloys

The great tunability of the properties of halide perovskites presents new opportunities for optoelectronic applications as well as significant challenges associated with exploring combinatorial chemical spaces. Here, in this work, we develop a framework powered by high-throughput computations and machine learning for the design and prediction of mixed cation halide perovskite alloys. In a chemical space of ABX 3 perovskites with a selected set of options for A, B, and X species, pseudo-cubic structures with B-site mixing are simulated using density functional theory (DFT) and several properties are computed, including stability, lattice constant, band gap, vacancy formation energy, refractive index, and optical absorption spectrum, using both semi-local and hybrid functionals. Neural networks (NN) are used to train predictive models for every property using tabulated elemental properties of A, B, and X site species as descriptors. Starting from a DFT dataset of 229 points, we use the trained NN models to predict the structural, energetic, electronic, and optical properties of an extensive dataset of 17955 compounds, and perform high-throughput screening in terms of stability, band gap, and defect tolerance, to obtain 392 promising compounds that are ranked as potential absorbers according to their photovoltaic figure of merit. Compositional trends in the screened set of attractive mixed cation halide perovskites are revealed and additional computations are performed on selected compounds. The data-driven design framework developed here is promising for designing novel mixed compositions and can be extended to a wider perovskite chemical space in terms of A, B, and X species, different kinds of mixing at the A, B, or X sites, non-cubic phases, and other properties of interest.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Defects in epitaxial Ru(0001) on Al2O3(0001): Dislocations, stacking faults, and deformation twins

In this study, defects in epitaxial Ru(0001) films on c-plane sapphire, with nominal thicknesses of 10–80 nm, deposited at 350 °C and step-annealed to 950 °C, were characterized using transmission electron microscopy. The variation of Ru and sapphire lattice parameters with temperature is such that the misfit strain for the observed 30° rotated-honeycomb epitaxial relationship is essentially constant with temperature at 1.5%, resulting in a biaxial stress of 10.0 GPa and an energy density of 150 MJ m –3 in unrelaxed films. Stress relaxation occurs by the formation of defects. For the 20–80 nm thick films, the defects are a- and c-type dislocations and stacking faults, argued to be of I 2 type. In addition, the films show the surprising presence of $\{11\bar{2}1\}1/3\langle11\bar{2}\bar{6}\rangle$ deformation twins. The 10 nm-thick films were found to be defect free. The critical thickness for misfit strain relaxation via the formation of threading and misfit dislocations is computed as 7±2 nm, depending on the choice of the dislocation core radius. Energetic analysis of twin formation, using both the infinite-matrix and the finite-matrix (Mori–Tanaka) approaches, provides values of the twin aspect ratios, assumed to be ellipsoidal, and shows that the latter but not the former approach can qualitatively explain the formation of the observed twins. In addition to providing the maximum strain relief compared to other potential twin types, $\{11\bar{2}1\}1/3\langle11\bar{2}\bar{6}\rangle$ twins do not require lattice shuffles and have a boundary that is a special boundary, namely, a 35° tilt boundary with a-type dislocations every other {0002} plane, that may also favor their formation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Evidence for anti-synergism between ion-assisted etching and in-plasma photoassisted etching of silicon in a high-density chlorine plasma

Etching of p-Si in 60 mTorr 10%Cl2/90%Ar Faraday-shielded inductively coupled high density plasmas was investigated under both ionassisted etching (IAE) and photoassisted etching (PAE) conditions. Real-time etching rates and after-etching Si surface chemical compositions were obtained by laser interferometry and vacuum-transfer x-ray photoelectron spectroscopy (XPS), respectively. Precisely controlled ion energy distributions (IEDs) were generated by applying pulsed negative DC bias on the conductive sample stage. Above a 36 eV threshold at a total flow rate of 250 SCCM, the IAE rate increased with the square root of the ion energy. In contrast to the DC bias, etching under RF bias did not exhibit a threshold ion energy because of the wide IED. XPS spectra revealed that the surface layer under PAE conditions had a significantly lower chlorine content, composed of only SiCl. Under IAE conditions, however, silicon dangling bonds (Si•), SiCl 2 , and SiCl 3 were found on the surface, in addition to SiCl, with a relative abundance of SiCl > SiCl 2 , and SiCl 3 . The absence of higher chlorides and Si• under PAE conditions suggested that vacuum ultraviolet photons and above threshold-energy ions interact with the surface very differently. By varying the duty cycle of the pulsed DC bias, it was found that the IAE rate scaled with the energetic ion dose, but only for low duty cycles. For higher duty cycles, the apparent IAE yield fell off with an increasing Cl coverage on the surface, as the duty cycle went up, which pointed to a negative synergy (antisynergism) between PAE and IAE as the explanation. This antisynergism was further supported by the observed decrease of the total etching rate with an increasing period of the pulsed DC bias. A plausible mechanism is that increasing the pulsing period causes more near-surface damage, creating more recombination centers that lead to a higher loss rate of electron-hole pairs through recombination, thereby reducing the PAE rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

First-principles DFT modeling of nitrobenzene adsorption on the Ag(111) surface at varying monolayer coverages

Here, the adsorption behavior of nitrobenzene on the Ag(111) surface as a function of coverage is investigated using density functional theory. Adsorption energies and optimized geometries are analyzed together with isolated intermolecular interaction calculations and electronic structure analysis, including Bader charge partitioning and projected density of states, to disentangle the roles of adsorbate–surface bonding and through-space adsorbate–adsorbate interactions. At low and intermediate coverages (θ = 1/9 and θ = 2/9), bidentate adsorption at top sites is favored due to strong adsorbate–surface interactions, with additional stabilization at θ = 2/9 arising from favorable intermolecular separations. At higher coverage (θ = 1/3), bidentate adsorption is destabilized by strong intermolecular repulsion, and monodentate adsorption becomes energetically preferred as rotational freedom allows more favorable intermolecular spacing. Charge density differences, Bader charge, and density of states analyses show that charge transfer from the Ag surface is localized primarily on the nitro group and increases with coverage and adsorption denticity, although this increase does not directly correlate with adsorption strength at high coverage due to competing intermolecular interactions. These results demonstrate that surface coverage can induce a transition in preferred adsorption denticity driven by intermolecular interactions, highlighting the importance of adsorbate packing in organic molecule adsorption at metal interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Density Functional Theory Study of Iron–Oxygen Divacancies in Magnetite (Fe 3 O 4 ) and Hematite (Fe 2 O 3 )

Density functional theory (DFT) calculations are employed to investigate the formation energies, charge redistribution, and binding energies of iron–oxygen divacancies in magnetite (Fe 3 O 4 ) and hematite (Fe 2 O 3 ). For magnetite, we focus on the low-temperature phase to explore variations with local environments. Building on previous DFT calculations of the variations in formation energies for oxygen vacancies with local charge and spin order in magnetite, we extend this analysis to include octahedral iron vacancies before analyzing the iron–oxygen divacancies. We also assessed the relative stability of iron–oxygen divacancies by comparing their formation energies with those of individual vacancies. Our findings reveal a significant energetic driving force for the formation of divacancy clusters, particularly in magnetite, where divacancies in the +1 charge state exhibit formation energies comparable to those of neutral iron vacancies under oxidizing conditions. In hematite, the results indicate a strong tendency for oxygen vacancies to bind to iron vacancies. These results highlight the significance of iron–oxygen vacancy complexes in the transport properties of iron oxides, with particular relevance to diffusion mechanisms under irradiation conditions.

36 MATERIALS SCIENCE↗