Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT analysis”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Hydride Migration within RhH 2 Ag 19 Superatom: A Combined Neutron Diffraction and DFT Analysis

Abstract An investigation combining neutron diffraction and DFT allows determining the most likely hydride migration pathway within the icosahedral metal framework of [RhH 2 Ag 19 {S 2 P(O n Pr) 2 } 12 ] (RhH 2 Ag 19 ). Starting from the experimentally derived solid‐state structures, a computational analysis is able to reveal an energetically favorable migration pathway with a maximum energy barrier of 4.2 kcal mol −1 . The two hydrides migrate simultaneously within the Rh@Ag 12 icosahedral core, traversing several positional isomers. This study expands the understanding of hydride dynamics in nanoclusters and provides critical insights into the structural flexibility of the superatom framework. These findings have significant implications for hydrogen storage, catalysis, and the design of advanced hydride‐containing materials.

Chemistry↗
DOE OSTI2025

Mechanistic Studies of Ligand Substitution, Linkage Isomerism, and Insertion Reactions in Electron Rich Pd(II) Complexes of a Zwitterionic Diimine Ligand

We have prepared cationic palladium complexes possessing a new zwitterionic ligand bis-N,N’–1-(2,4,6-triphenylpyridyl) oxalamide [(N ^ N)Pd(Me)(L)] + [BArF] - , (BArF=3,5-(CF 3 ) 2 C 6 H 3 , L=NCMe, CO). The structure of [(N ^ N)Pd(Me)(CO)] + [BArF] - was determined by X-ray diffraction analysis. Energy Decomposition Analysis (EDA) indi-cates this N ^ N zwitterionic ligand is more electron-donating relative to bidentate diimine ligands. Low temperature NMR analysis shows the existence of linkage isomers with the N ^ N isomer the most stable. Structures were assigned using NMR and DFT analysis. Barriers to interconversion of isomers are ΔG ‡ = 10-12 kcal/mol. Kinetics of acetoni-trile displacement from [(N ^ N)Pd(Me)(NCCH 3 )] + [BArF] - by CD 3 CN, ethylene and t Bu 3 P were measured and mechanisms of exchange determined. The ethylene complex, [(N ^ N)Pd(Me)(C 2 H 4 )] + was generated at -45 °C, and the barrier of migrato-ry insertion determined at 0 °C (ΔG ‡ = 23.4 kcal/mol) and compared to related diimine complexes. The methyl carbonyl complex undergoes migratory insertion in the presence of CO at -70 to -55 °C (ΔG ‡ = ca. 15.7 kcal/mol) to yield the acyl carbonyl complex. Furthermore, the neutral bis-trimethylsilylmethyl complex [(N ^ N)Pd(CH 2 SiMe 3 ) 2 was prepared and characterized by X-ray diffraction analysis. It displays dynamic behavior at very low temperatures in the NMR spectrum (-90 °C, ΔG ‡ =7.9 kcal/mol) which, supported by DFT analysis, is ascribed to rotation of the bulky -CH 2 SiMe 3 groups.

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

Surface Phase Stability of Fe 2 O 3 (001) in Hydrogen Reducing Environments: A DFT and XPS Analysis

Here, this study combines density functional theory (DFT) and ab initio thermodynamics calculations with X-ray photoelectron spectroscopy (XPS) investigations to identify the reduction properties of the Fe 2 O 3 (001) surface with implications for corrosion resistance, hydrogen transport, and energy safety. Ab initio thermodynamics modeling predicts fully hydroxylated surface stability across a broad range of pressures (1 × 10 –23 to 1 × 10 5 mbar) and temperatures below 700 K, consistent with previous experimental studies. Above 800 K, exposures to 1 × 10 –4 mbar H 2 , 1 × 10 –4 mbar O 2 , or 1 × 10 –4 mbar H 2 + 1 × 10 –4 mbar O 2 each yield unique XPS signals indicating a loss of −OH coverage, aligning with DFT predictions. Insight into the mechanism of reduction as a function of H 2 exposure is provided, as well as conditions that promote further reduction toward Fe 3 O 4 . Theoretical and experimental investigations indicate the ability to maintain the Fe 2 O 3 protective layer of iron oxides that have been exposed to H 2 environments by including trace amounts of aqueous O2.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Fe‐Triazolate Metal–Organic Frameworks as Water Oxidation Catalysts with Dual Photoanode Functionality

Artificial photosynthesis is an emerging technology that achieves renewable fuels, such as hydrogen, from sunlight. Its realization depends on finding highly active and stable catalysts of water splitting and photoactive materials for light absorption. To be scalable, these should contain only abundant elements. Here, for the first time, Fe-triazolate (Fe(ta) 2 ) and its metal substituted derivatives (Fe-Metal(ta) 2 ) Metal-organic frameworks (MOFs) are characterized as new dual-function materials for photo-absorption and water oxidation catalysis in acidic media. The materials were studied by a range of structural, spectroscopic, and computational density functional theory (DFT) techniques. Fe(ta) 2 and Fe-Mn(ta) 2 were found to be highly active and stable in chemical and photochemical water oxidation, and in addition function as photoanodes, with photo-electrocatalytic currents (∼2.00 x 10 −3 Acm −2 at + 1.4 V vs. Ag/AgCl) at pH = 1. The possibility of a unique catalytic mechanism where O─O bond formation is possible from the coupling of two adjacent Fe IV = O fragments was demonstrated by DFT analysis. Thus, Fe-triazolate MOF has been established as a new, stable, scalable, versatile, and efficient platform for sustainable energy conversion in the realm of artificial photosynthesis.

Artificial photosynthesis↗
DOE OSTI2025

Atomic-resolution imaging as a mechanistic tool for studying single-site heterogeneous catalysis

Heterogeneous catalysts dominate the chemical industry but typically feature diverse, incompletely defined active sites. Thus, describing structure-activity relationships, unlike homogeneous catalysts, remains challenging. In contrast, molecularly defined single-site heterogeneous catalysts (SSHCs), using appropriate tools, are poised to address these challenges and provide new avenues for catalysis research and development. The present study explores eco-friendly H 2 production mediated by discrete MoO 2 sites supported on carbon nanohorns (CNHs) and active for alcohol dehydrogenation. While informative, detailed ensemble EXAFS/XANES/XPS, kinetic measurements, and DFT analysis alone cannot provide a full molecular picture of the reaction pathway. Here, using single-molecule atomic-resolution time-resolved electron microscopy (SMART-EM), we propose the identification of four key catalytic intermediates anchored to CNHs and uncover a new reaction pathway involving alkoxide/hemiacetal equilibration and acetal oligomerization. Furthermore, these intermediates are inferred through a combination of theory and SMART-EM, showcasing the potential of SMART-EM as a complementary tool for exploring mechanistic hypotheses in catalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The presence of Ru metal modifies the behavior of deuterium bound to pyridinic nitrogen in doped graphene-like materials

We used a combination of x-ray photoelectron spectroscopy (XPS) and density functional theory (DFT) to investigate how the presence of a ruthenium (Ru) metal substrate modifies the N-D bond strength of the pyridinic nitrogen species in nitrogen (N)-doped graphene (Gr) on Ru(0001) compared to metal-free highly-oriented pyrolytic graphite (HOPG) substrate. The N-dopants were introduced through low-energy N+/N2+ irradiation. N is embedded in the carbon materials in two predominant configurations: as graphitic N (GN: N substituted in the hexagonal C lattice) and pyridinic N (PN: substitutional N adjacent to a C vacancy). XPS data showed that upon atomic deuterium (D) exposure at 220 K, the PN peak shifted to a higher binding energy by +1.2 eV for HOPG and +1.0 eV for N-doped graphene on Ru(0001), while the GN peak remained unchanged, indicating that the D atoms bound solely to pyridinic N. 85% of PN sites on HOPG can be saturated with D atoms, whereas only ~30% of pyridinic N sites are able to bind D atoms in N-doped Gr/Ru(0001). Our DFT analysis shows that this difference is due to the coordination of PN to Ru atoms, which necessitates bond cleavage of the N-Ru interaction prior to D atom adsorption. D begins desorbing from N-HOPG at ~573 K and is fully desorbed by ~973 K, whereas desorption from N-doped graphene on Ru(0001) begins at ~290 K, with complete desorption observed at ~700 K, indicating that the Ru metal weakens the N-D bond strength. We also studied a high-surface-area, layered, and porous N-doped carbon material, internally labeled NC900, which was synthesized by pyrolysis of glucose and graphitic carbon nitride (g-C3N4) at 1173 K. D exposure caused a +1.1 eV PN peak shift in NC900, with no change in the GN peak. Moreover, the NC900 exhibited the same desorption behavior as HOPG, demonstrating that well-defined model systems can effectively capture the behavior of more complex N-doped carbon materials.

Alupothe Gedara, Buddhika S.↗
DOE OSTI2026

Cosolvent-tuned interactions in ionic liquids: A vibrational and quantum-chemical study of ethylene glycol ratio effects

Ionic liquids (ILs) are attractive media for CO 2 capture but remain limited by viscosity and cost. Blending ILs with ethylene glycol (EG) is a practical route to mitigate these constraints, yet the molecular origins of cosolvent effects and their dependence on composition are not well resolved. We combine Fourier-transform infrared (FT-IR) spectroscopy with quantum-chemical (DFT) analysis to elucidate how the IL:EG molar ratio modulates intermolecular interactions and electronic structure. Computed vibrational frequencies enable mode assignment and deconvolution of overlapping bands, revealing systematic, ratio-dependent shifts and broadenings in (i) EG O–H stretching, (ii) cation and EG C–H stretchings (imidazolium C2–H, C4–H, C5–H, methyl and ethyl groups, -CH2 of EG), (iii) anion signature modes (e.g., CN motifs), and (iv) EG C–O and C–C stretchings, consistent with the redistribution of hydrogen-bonding networks. Molecular electrostatic potential (MESP) maps quantify attenuation of extreme potential regions with increasing EG, indicating progressive screening of cation–anion electrostatic interactions. Quantum Theory of Atoms in Molecules (QTAIM) identifies emergent bond critical points between EG and the IL ions, while Reduced Density Gradient–Noncovalent Interaction (RDG–NCI) analysis differentiates strong directional hydrogen bonds from dispersive contacts across compositions. Together, these results show that EG fraction controls a switch from predominantly ion–ion to mixed ion–EG coordination, altering local polarity and polarizability that underlie the observed FT-IR trends. The framework provides composition–structure–spectrum relationships that can guide rational selection of IL:EG ratios to balance favorable molecular interactions with practical performance targets in scalable CO 2 capture systems.

DAC↗
DOE OSTI2025

Impact of SO 2 on NiFe Nanoparticle Exsolution and Dissolution from LaFe 0.9 Ni 0.1 O 3 Perovskite Oxides

Ni-doped LaFeO 3 perovskite oxide is a promising cathode material for solid oxide electrolysis cells (SOECs) designed for CO 2 /H 2 O coelectrolysis. Here, the performance of LaFe 0.9 Ni 0.1 O 3 is being investigated under real-world conditions that include exposure to acid gases, such as SO 2 , relevant to SOEC operation. Experiments show that LaFe 0.9 Ni 0.1 O 3 exsolves NiFe nanoparticles, along with the formation of surface SO 4 2– and SO 3 2– after being exposed to 200 ppm of SO 2 . This suggests that the ionic diffusion of Ni 3+ and Fe 3+ between the bulk and the surface remains unaffected throughout the exsolution–dissolution–exsolution cycle. Thermochemical water splitting has been employed as a probe reaction to evaluate the catalytic properties of the exsolved NiFe nanoparticles. These nanoparticles demonstrated improved hydrogen production compared to bare perovskite oxide substrates. However, after exposure to SO 2 , the formation of Fe-rich NiFe nanoparticles led to poor thermocatalytic performance and rapid deactivation of the perovskite at elevated temperatures. Density functional theory (DFT) analysis was utilized to validate the experimental findings, indicating a significantly negative reaction energy for water splitting over exsolved Fe, as well as stronger binding of SO 2 to Fe than to Ni. Computational analysis further suggests that the presence of surface sulfate promotes the formation of Fe-rich NiFe nanoparticles, aligning with the experimental results. Overall, this study clarifies how SO 2 affects the structure of SOEC perovskite oxide candidate materials. Future engineering efforts should focus on enhancing nanoparticle exsolution and sulfur resistance, which is crucial for improving the hydrogen production capacity of La-based perovskite oxides for electro- and thermocatalytic water splitting in real environments containing acid gases.

Najimu, Musa [Univ. of Southern California, Los An↗
DOE OSTI2025

Defect Diffusion Graph Neural Networks for Materials Discovery in High-Temperature Energy Applications

Here, the migration of crystallographic defects dictates material properties and performance for a plethora of technological applications. Density functional theory (DFT)-based nudged elastic band (NEB) calculations are a powerful computational technique for predicting defect migration activation energy barriers, yet they become prohibitively expensive for high-throughput screening of defect diffusivities. Without introducing hand-crafted (i.e., chemistry- or structure-specific) descriptors, we propose a generalized deep learning approach to train surrogate models for NEB energies of vacancy migration by hybridizing graph neural networks with transformer encoders and simply using pristine host structures as input. With sufficient training data, computationally efficient and simultaneous inference of vacancy defect thermodynamics and migration activation energies can be obtained to compute temperature-dependent vacancy diffusivities and to down-select candidates for more thorough DFT analysis or experiments. Thus, as we specifically demonstrate for potential water-splitting materials, candidates with desired defect thermodynamics, kinetics, and host stability properties can be more rapidly targeted from open-source databases of experimentally validated or hypothetical materials.

14 SOLAR ENERGY↗
DOE OSTI2024

Spectroscopic Signature of Metal-hydroxo and Peroxo Species in K-edge X-ray Absorption Spectra

Metal-dioxygen species are important intermediates formed during dioxygen activations by metalloenzymes in various biological processes, by catalysts in fuel cells, and prior to O 2 evolution by photosystem II. In this work, we focus on manganese-porphyrin complexes using tetramesitylporphyrin ligand (TMP) to explore changes in Mn K-edge X-ray absorption spectroscopy (XAS) associated with the formation of Mn-hydroxide and Mn–O 2 peroxide species. With limited spectroscopic characterization of these compounds, Mn K β X-ray emission spectroscopy (XES), XAS, density functional theory (DFT), and time-dependent DFT (TD-DFT) analysis will enhance our understanding of their complex electronic structure. We show that the shape of the pre-edge in the K-edge Mn X-ray absorption near-edge structure (XANES) can serve as a spectroscopic signature of the Mn III -peroxo formation and thus can be used to track the presence of the side-on peroxide as an intermediate in time-resolved or in situ experiments. In conclusion, our results will help to further summarize the spectroscopic fingerprints for peroxo and hydroxo species, addressing the challenge of identifying the reactive metal species in catalytic reactions.

Extended X-ray absorption fine structure↗
DOE OSTI2025

Tetracene Functionalized Si(111) Achieves Enhanced Solar-to-Chemical Energy Conversion via Molecular Acceptor States

The properties of semiconductor|liquid interfaces play a critical role in determining the efficiency of solar-to-hydrogen (STH) conversion. Here, we investigate how molecular functionalization of Si(111) and Si(111)|TiO 2 surfaces impacts photoelectrochemical (PEC) hydrogen production efficiency. We find that functionalization of ∼3% of the atop sites of Si(111) with either 9-anthracene (Anth) or 5-tetracene (Tet), with the remaining sites passivated by methyl groups, provides substrates with high electronic quality and low surface oxide densities, as determined by X-ray photoelectron spectroscopy (XPS) measurements. Surface photovoltage (SPV) spectroscopy shows that surfaces modified with Anth or Tet exhibit an increased photovoltage, with Tet-functionalized surfaces yielding an additional 192 meV relative to methyl-terminated Si(111), indicating improved charge separation for Si-Tet. Further improvement in onset potential was achieved by replacing a nitrogen-containing TiO 2 atomic layer deposition (ALD) precursor (TDMAT) with a precursor lacking nitrogen (TTIP), which eliminates the parasitic defect band in the TiO 2 overlayer (p-Si(111)-Tet|TTIP-TiO 2 |Pt: V OC = +0.283 ± 0.041 V vs RHE). Density functional theory (DFT) analysis demonstrates that compared with Anth-modified Si(111), the Tet-modified surface exhibits more hybridized Si(111)-Tet states closer to the silicon band edges. Mercury contact current–voltage (I–V, dark) measurements quantified the relative interfacial density of states of Si-Tet, Si-Anth and Si-Me surfaces─revealing that the interfacial state density was highest for Si-Tet. This suggests that such hybridized interfaces serve to capture better photoexcited charge, which enables facile electron transfer to molecular acceptors in solution. Altogether, the data indicate that beneficial hybrid molecular LUMO surface states interacting with the Si conduction band edge results in improved hydrogen evolution (HER) performance for p-Si devices.

Group theory↗
NASA NTRS1994

Blocking reduction of Landsat Thematic Mapper JPEG browse images using optimal PSNR estimated spectra adaptive postfiltering

Two representative sample images of Band 4 of the Landsat Thematic Mapper are compressed with the JPEG algorithm at 8:1, 16:1 and 24:1 Compression Ratios for experimental browsing purposes. We then apply the Optimal PSNR Estimated Spectra Adaptive Postfiltering (ESAP) algorithm to reduce the DCT blocking distortion. ESAP reduces the blocking distortion while preserving most of the image's edge information by adaptively postfiltering the decoded image using the block's spectral information already obtainable from each block's DCT coefficients. The algorithm iteratively applied a one dimensional log-sigmoid weighting function to the separable interpolated local block estimated spectra of the decoded image until it converges to the optimal PSNR with respect to the original using a 2-D steepest ascent search. Convergence is obtained in a few iterations for integer parameters. The optimal logsig parameters are transmitted to the decoder as a negligible byte of overhead data. A unique maxima is guaranteed due to the 2-D asymptotic exponential overshoot shape of the surface generated by the algorithm. ESAP is based on a DFT analysis of the DCT basis functions. It is implemented with pixel-by-pixel spatially adaptive separable FIR postfilters. PSNR objective improvements between 0.4 to 0.8 dB are shown together with their corresponding optimal PSNR adaptive postfiltered images.

Linares, Irving↗
DOE OSTI2026

Mo 4 FeGa 17.25– x Ge x : Complementary Point Substitutions, Buffering Frameworks, and Merging of the 18-n and Octet Bonding Schemes

We present the discovery of Mo 4 FeGa 17.25–x Ge x (x ∼ 2.1(4), based on determination of Ge content), a complex gallide featuring a variety of point substitution phenomena. Its crystal structure is derived from the Ti 2 Ni type, in which a diamond network of face-sharing octahedra is interpenetrated by a second diamond network of vertex-sharing stella quadrangula. However, while in the ideal Ti 2 Ni type these two frameworks run uninterrupted through the crystal, Mo 4 FeGa 17.25–x Ge x shows three variations. First, the inner tetrahedron of every other stella quadrangula is replaced with a main group atom, creating tetrahedra reminiscent of the Zintl phase NaTl. Next, a selection of Ga 6 octahedra are filled with Fe atoms, in a manner analogous to the stuffed AuCu 3 -type phases. Finally, the refined crystal structure shows that in each unit cell a single Ga2 atom in the octahedral network is substituted with a dumbbell of Ga/Ge atoms. Electronic structure calculations on ordered models of Mo 4 FeGa 17.25–x Ge x reveal a narrow band near the Fermi energy, which is explained in terms of the 18-n and octet bonding schemes using reversed approximation Molecular Orbital analysis. A DFT-chemical pressure analysis connects the tetrahedron/atom and atom/dumbbell substitutions to the relief of atomic packing tensions and highlights soft atomic motions within the octahedral framework and driving forces for atom/dumbbell substitution. The combination of soft vibrational modes, disorder, and a narrow band gap could make this phase of interest for potential thermoelectric properties.

Chemical structure↗
DOE OSTI2025

Hydrogen Bond Benchmark: Focal‐Point Analysis and Assessment of DFT Functionals

We performed a hierarchical, convergent ab initio benchmark study and systematically analyzed the performance of density functional approximations for describing hydrogen bonds in small neutral, cationic, and anionic complexes, as well as in larger systems involving amide, urea, deltamide, and squaramide moieties. Focal point analyses (FPA), extrapolating to the ab initio limit, were carried out using correlated wave function methods up to CCSDT(Q) for the small complexes and CCSD(T) for the larger systems, together with correlation-consistent Gaussian basis sets up to the complete basis set limit. Optimized geometries and vibrational frequencies were obtained at the CCSD(T) level. The resulting FPA hydrogen-bond energies converge within a few tenths of a kcal mol −1 . These reference data were used to evaluate 60 density functionals (including 12 dispersion-corrected), spanning the local-density approximation (LDA), generalized gradient approximations (GGAs), meta-GGAs, hybrids, meta-hybrids, double-hybrids, and range-separated hybrids. Overall, the meta-hybrid M06-2X provides the best performance for both hydrogen bond energies and geometries, while the dispersion-corrected GGAs BLYP-D3(BJ) and BLYP-D4 also yield accurate hydrogen-bond data and can serve as cost-effective options for studying large and complex systems.

coupled cluster theory↗
DOE OSTI2025

Controllable oxygen vacancy defect engineering of BiVO 4 porous structures for room temperature NH 3 detection

Controlling structural features of sensing material while judiciously introducing vacancy defect states for revamping the electronic properties of the sample to obtain its superior gas sensing performance, is quite rare. Herein, we report for the first time, the room temperature (RT) ammonia (NH 3 ) detection of peanut-like porous bismuth vanadate (BiVO 4 ) with a stable monoclinic phase. The developed BiVO 4 possess abundant oxygen vacancies and porosity by virtue of calcinations (400–800 ℃). BiVO 4 calcined at 400 ℃ exhibits high selectivity towards NH 3 with a maximum response of 1421 @ 270 ppm at RT, which is 2.5 fold enhanced compared to without calcined BiVO 4 (response of 547 @ 270 ppm NH 3 ). The oxygen defective BiVO 4 appeared highly durable and stable even under high humid conditions (∼60 %). Besides, the porosity of BiVO 4 not only enhances the specific surface area (19.8 m 2 /g) but also results in fast diffusion of NH 3 molecules, leading to a reduction in the decay time (34 s for 90 ppm NH 3 ). The density functional theory (DFT) uncovers that the oxygen vacancy formation in BiVO 4 augments the NH 3 sensing capabilities by enhancing the adsorption energy of NH 3 . This work provides insight into the sensing mechanism of increased response caused by defect engineering and porosity, which will be favourable for fabricating high-performance NH 3 sensors at RT.

DFT analysis↗
DOE OSTI2025

Computational insights into hydrogen adsorption energies on medium-entropy oxides

High entropy oxides (HEOs) have emerged as promising catalysts for several important chemical transformations including alkane activation. Hydrogen adsorption energy (HAE) has been used as a key descriptor for many reactions including methane C–H activation and hydrogen evolution reactions. Hence, understanding the relationship between HAEs and the surface chemistry of HEO surfaces could lay the foundation for meaningful correlations among methane C–H activation, HAE, and the complex, local environment of HEO surfaces. Here, we used a medium-entropy oxide as a prototypical system – Mg 0.25 Ni 0.25 Cu 0.25 Zn 0.25 O with a rock-salt structure – to interrogate these relationships. We sampled 2000 different surfaces of its (100) plane and calculated the HAEs at randomly chosen surface O sites using density functional theory (DFT). Our analysis of the 2000 data points reveals that the HAEs at the surface O sites are significantly influenced by the local environment around the adsorption sites, particularly the nature of the metal atom directly below the surface O site where H adsorbs. After comparing several popular graph-neural-network-based machine learning models, we found that the DimeNet++ model performed best achieving satisfactory accuracy in predicting HAEs for both Mg 0.25 Ni 0.25 Cu 0.25 Zn 0.25 O and slightly varied compositions. Our work underscores the promise of such models and the need for further refinement to address the complexity of HEOs.

Song, Haohong [Vanderbilt Univ., Nashville, TN (Un↗
DOE OSTI2025

Origin of Stabilization of Ligand-Centered Mixed Valence Ruthenium Azopyridine Complexes: DFT Insights for Neuromorphic Applications

Redox-driven conductance changes are critical processes in molecular- and coordination-complex-based memristive thin films and devices that are envisioned for neuromorphic technologies, but fundamental mechanisms of conductance switching are not fully understood. Here, we explore charge disproportionation (CD) processes in [Ru II L 2 ](PF 6 ) 2 molecular systems that intrinsically involve interfragment charge transfer (IFCT). Using a combination of ab initio molecular dynamics simulation (AIMD), time-dependent density functional theory (TD-DFT), and density functional theory (DFT) calculations, we investigate the electron transfer mechanisms and the roles of temperature and cell volumetric expansion in facilitating the counterion movements and electronic transitions required for low-cost IFCT and charge redistribution. A detailed analysis of the density of states and TD-DFT calculations highlights that unpaired electrons play a crucial role in low-energy transitions, with the azo (N=N) groups of the ligand serving as the primary sites for electronic transport between molecular fragments, further stabilizing the asymmetric state. Localization of added electrons on azo ligands occurs with negligible change at the Ru centers, supported by atomic volume expansions up to +4.74 bohr 3 , and goes along with a progressive reduction of the HOMO−LUMO gap across redox states, suggesting enhanced conductivity. The TD-DFT analysis reveals a dominant IFCT excitation at 2082.76 nm in the doubly reduced (22) state, while a stabilization energy of 1.20 eV of the asymmetric (13) state relative to the symmetric (22) state is predicted by constrained DFT. Periodic DFT and AIMD simulations emulating a molecular film show that the stabilization of the asymmetric state, relative to a symmetric one, translates in net charge separation values (order of ∼0.33 e) that are strongly linked to increased counterion mobility (average counterion displacements exceeding 0.7 Å per atom during CD events) and the involvement of azo groups in electron redistribution. These findings, which align with previously reported experimental and computational data, provide key insights into the IFCT mechanisms and electronic transport facilitated by azo groups, with important implications for redox-driven memristive and neuromorphic technologies.

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

A Robust Methodology to Elucidate Kinetics of Room Temperature Electrochemical Propane Adsorption on Platinum

Electrocatalytic activation of alkanes can further decarbonize chemical manufacturing by leveraging affordable renewable electricity and readily available shale gas reserves in the United States. Earlier works have identified the unique role of Pt in adsorbing and activating alkanes, like propane, at room temperature in acidic, aqueous electrolytes, revealing spontaneous formation of deeply dehydrogenated propane-derived surface species with an intact C 3 - backbone. Although an adsorption mechanism was hypothesized, it has not been explicitly investigated to date, preventing the quantification of kinetic rate parameters. A robust methodology to investigate and benchmark propane adsorption kinetics on Pt is critical for the rational design of electrocatalysts that exhibit higher selectivity toward desired partially oxidized products. Herein, we analyze an oxidative current transience that appears during the adsorption of propane on Pt in aqueous electrochemical conditions and develop a methodology that elucidates the adsorption mechanism and enables quantification of rate parameters such as order dependences and apparent activation barriers. This method yields an expected first-order dependence with respect to propane concentration at low coverage and reveals a second-order dependence with respect to the concentration of surface active sites. Additionally, the apparent activation barrier for propane adsorption was calculated using an Arrhenius analysis of the current transience under temperature control. The experimentally measured activation barrier of 35 kJ mol –1 is in excellent agreement with the theoretical barrier calculated by density functional theory (DFT). The kinetic analysis was extended, via the use of transition state theory, to extract entropy and enthalpy of activation, yielding consistent results with the proposed two-step adsorption mechanism and DFT calculations. These results demonstrate reliable quantification of kinetic parameters for electrocatalytic activation of C–H bonds in alkanes that can be employed for rational catalyst development for a versatile range of electrocatalytic conditions.

alkane activation↗