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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

NNL.Fe.qSNAP-ZBL.2024.1: A Fe Spectral Neighbor Analysis Potential for Radiation Damage Simulations

The NNL.Fe.qSNAP-ZBL.2024.1 machine-learned potential (MLP) has been generated to support the development of an elemental body-centered cubic (BCC) Fe athermal recombination corrected neutron damage model and simulations of primary recoil atom (PRA) cascades in BCC Fe. This MLP is a quadratic spectral neighbor analysis potential (qSNAP) hybridized with the universal Ziegler-Beirsack-Littmark (ZBL) potential at short-range and is named according to Naval Nuclear Laboratory MLP naming conventions (NNL.material-system.MLP-type.year.version). Training set calculations for Fe are presented along with the subsequent MLP fitting procedure. A key criterion of the fitting procedure is that ZBL describes the short-range interaction with minimal impact on the MLP. The MLP is compared to density functional theory (DFT) predicted properties relevant to radiation damage simulation, including threshold displacement energies, for validation. The NNL.Fe.qSNAP-ZBL.2024.1 potential is considered suitable for molecular dynamics (MD) simulations of radiation defects up to 800 K and PRA cascades in BCC Fe up to around 10 keV. The potential can additionally be used on a limited basis for recoils of 10–20 keV, within which range the emergence of structures outside the training set in cascade simulations may cause system instabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Examining the Impact of Local Constraint Violations on Energy Computations in DFT

ABSTRACT This work examines the impact of locally imposed constraints in Density Functional Theory (DFT). Using a metric referred to as the extent of violation index (EVI), we quantify how well exchange‐correlation functionals adhere to local constraints. Applying EVIs to a diverse set of molecules for GGA functionals reveals constraint violations, particularly for semi‐empirical functionals. We leverage EVIs to explore potential connections between these violations and errors in chemical properties. While no correlation is observed for atomization energies, a significant statistical correlation emerges between EVIs and total energies. Similarly, the analysis of reaction energies suggests weak positive correlations for specific constraints. However, definitive conclusions about error cancellation mechanisms cannot be made at this time. These observations revealed by EVIs may be useful for consideration when designing future generations of semilocal functionals.

Khanna, Vaibhav [Department of Chemistry Universit↗
DOE OSTI2024

Interpretable, extensible linear and symbolic regression models for charge density prediction using a hierarchy of many-body correlation descriptors

Here, density functional theory (DFT) is routinely used to make electronic structure predictions for high-throughput screening of materials and molecules for technologically relevant areas, like the identification of better catalysts, electronic materials, and drug discovery. However, the DFT formalism is limited by (a) its poor (quadratic-to-quartic) scaling, and (b) the need to perform repeated eigenvalue computations of the electronic Hamiltonian as part of its self-consistent field (SCF) iteration procedure to obtain the converged ground state electron density, ρ (r). Approaches that directly predict ρ (r) of a structure with high accuracy can accelerate conventional SCF calculations and can also be used in linearly scaling methods such as orbital-free DFT. To this end, we present a procedure to predict the ground state electron density of molecular and periodic three-dimensional systems directly from the atomic structure with a particular emphasis on physical interpretability. In our framework, ρ (r) is modeled using many-body correlation descriptors that accurately capture the effects of local atomic arrangements in the neighborhood of a grid point. Our use of a linear regression scheme to fit to charge density data enables transparent analysis of the relative contributions of various types of local atomic correlations. By systematically including increasingly complex correlations, our model is shown to accurately predict ρ (r) for a variety of chemically and electronically diverse systems — amorphous Ge, Al(001) slab, crystalline Ga 2 O 3 , molecular benzene, and polyethylene. We then demonstrate a symbolic regression-based protocol to construct easily computable, interpretable features from lower-order correlations that significantly improves our electron density predictions with effectively no increase in the computational cost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Rhenium-Sulfido and -Dithiolato Corroles: Reflections on Chalcophilicity

The high-temperature (~180 °C) reaction between free-base meso-triarylcorroles and Re 2 (CO) 10 , followed by exposure to PCl 3 and thiols (or elemental sulfur), affords rhenium-sulfido (ReS) corroles in 67-76% yields. The use of shorter reaction times, lower temperatures (~130 °C), and a dithiol (e.g., ethane-1,2-dithiol) also allows the isolation of rhenium-dithiolato corroles, presumptive intermediates on the path to ReS corroles. The ReS corroles exhibit high thermal stability and two reversible oxidations and reductions in their cyclic voltammograms, with redox potentials nearly identical to those observed for analogous ReO corroles. The electrochemical HOMO-LUMO gaps of the complexes, at 2.2 eV, are consistent with ligand-centered oxidation and reduction. The UV-vis spectra of ReS corroles, on the other hand, differ significantly from those of their ReO counterparts. Scalar-relativistic DFT calculations suggest that this difference reflects low-energy LUMO+2 and LUMO+3 levels, consisting of Re-S π-antibonding interactions; the ReO corroles, in contrast, exhibit a larger LUMO+1/LUMO+2 gap, as expected for a relatively classical Gouterman-type metalloporphyrin analogue. The high stability of ReS corroles is consistent with geochemists' view of rhenium as a moderately chalcophilic element (i.e., one that partitions into sulfide melts) as well as with a recent quantitative analysis of thiophilicity, which indicates that rhenium's oxophilicity and thiophilicity are essentially evenly balanced.

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

Assessing the Sensitivity of Pourbaix Diagrams to Computational Protocols: Electrochemical Stability of Ni Oxides as a Case Study

Pourbaix diagrams stand as a useful tool in assessing and visualizing materials’ electrochemical stability and are widely used for electrocatalyst design. However, their reliability hinges on the accuracy of the chemical potentials of involved phases, which may bear uncertainties and can be significantly impacted by decision-making steps in the computational protocol. Here, this study introduces a robust sensitivity analysis framework, exemplified through a detailed examination of the computational Pourbaix diagram of Ni, the oxides of which are used as high-activity and cost-friendly catalysts for many electrochemical reactions. Quantities of interest derived from the Pourbaix diagram include the appearance and stability domain of the catalytically active Ni oxide phases along with the onset electrochemical potentials of phase transitions. These metrics can guide the design of operational conditions for Ni oxide electrocatalysts. We find that the employed DFT exchange-correlation functional has the most significant influence on the computed Pourbaix diagram. Uncertainties on crystal structures, along with their related ab initio energetics, are also found to affect the size of the phase stability domain. Higher-order coupling among input parameters is found to play a crucial role in influencing the appearance and distribution of Ni phases in the diagram. Our findings suggest a need to consider variations and uncertainties associated with the computational procedures on predicted Pourbaix diagrams for materials design.

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

Colloidal Synthesis of Palladium Nanocluster‐Decorated Cs 3 Sb 2 Cl 9 Perovskite Heterostructural Nanorods for Enhanced CO 2 Photoreduction

Developing efficient and sustainable photocatalysts for CO 2 reduction remains a significant challenge, particularly with environmentally benign materials. Here, in this study, we report the first one-step synthesis of metal–lead-free perovskite heterostructural nanocrystals by decorating Cs 3 Sb 2 Cl 9 perovskite nanorods with size-controlled Pd nanoclusters via a one-step hot-injection method. The resulting Pd-Cs 3 Sb 2 Cl 9 heteronanorods (HNRs) exhibit strong interfacial electronic coupling, enhanced charge separation, and excellent colloidal stability. Transient absorption spectroscopy and DFT calculations reveal a built-in electric field that drives directional electron transfer from the perovskite host to the Pd domains. Under UV irradiation, the Pd-Cs 3 Sb 2 Cl 9 HNRs demonstrate excellent CO 2 photoreduction activity with high CH 4 selectivity, achieving a record apparent quantum yield (AQY) of 2.62% among halide perovskite nanocrystal-based systems with a large electronic yield of 689.3 ± 12.2 µmol·g cat −1 . In situ spectroscopic monitoring and Gibbs free energy analysis further unveil a Pd-facilitated reaction pathway involving stabilization of key intermediates. This work introduces a new class of lead-free perovskite-based heterostructures through a facile one-step synthesis strategy and offers a new design principle for next-generation photocatalysts for solar fuel production.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Magnetic Properties Tuning via Broad Range Site Deficiency in Square Net Material UCu x Bi 2

HfCuSi 2 -type pnictogen compounds have recently been shown to be a versatile platform for designing materials with topologically nontrivial band structures. However, these phases require strict control over the electron count to tune the Fermi level, which can only be achieved in compositions with A 2+ M 2+ Pn 2 and A 3+ M + Pn 2 (A = lanthanides, M = transition metals, Pn = pnictogens P–Bi) charge distribution. While such lanthanide compounds have been thoroughly studied as candidate magnetic topological materials, their heavy element analogs with uranium and bismuth remain largely underexplored. In this report, we present the synthesis of UCu x Bi 2 single crystals and study their magnetic properties. Detailed structural analysis revealed that flux-grown crystals always form as a site-deficient UCu x Bi 2 composition, where x varies between 0.20 and 0.64. Magnetic property measurements revealed a dependence of the magnetic coupling on the Cu site deficiency, linearly changing the Néel temperature from 51 K for UCu 0.60 Bi 2 to 118 K for UCu 0.30 Bi 2 . Moreover, higher Cu concentration promotes a metamagnetic transition in highly magnetically anisotropic UCu 0.60 Bi 2 single crystals. We show that DFT calculations can successfully model site deficiency in the UCu x Sb 2 and UCu x Bi 2 systems. This work paves the way toward using the site deficiency to tune the Fermi level in more ubiquitous A 3+ M 2+ x Pn 2 phases, which previously have not been considered topological candidate materials due to unfavorable electron count.

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

Field-Coupled Water Splitting with Metal-Free Donor–Acceptor Covalent Organic-Framework Junctions

Advancing metal-free electrocatalysts for hydrogen and oxygen evolution reactions (HER/OER) across acidic and alkaline media requires coordinated control of intermediate binding thermodynamics, interfacial charge delivery, and near-electrode transport dynamics. Here, we design amide-linked benzene–triazine covalent organic frameworks (BTA/TzTA-Hz COFs) and integrate them with carbon nanotubes (CNTs) to form COF–CNT junctions that establish a built-in interfacial electric field. Density functional theory (DFT) and electrostatic potential maps indicate complementary active motifs, with benzene-proximal fragments associated with HER and triazine-proximal motifs associated with OER. CNT integration shifts the contact-potential difference by ≈0.20 V, while operando electrochemical impedance spectroscopy suggests partially separable high-frequency junction-charging and lower-frequency Faradaic/transport responses. A 300 mT static magnetic field lowers the HER and OER overpotentials by tens of millivolts. Under anodic bias, the effective interfacial charging capacitance increases, and Mott–Schottky analysis shows an apparent ∼0.15 V flat-band shift with an essentially unchanged slope. Together, these observations are consistent with field-perturbed interfacial charging and altered bias partitioning. Field-dependent impedance and bubble imaging are consistent with magnetohydrodynamic convection that promotes bubble detachment and near-electrode mass transport for both half-reactions, and they reveal an OER-specific high-frequency perturbation under anodic bias. Under field, the heterostructure reaches an OER onset overpotential of ∼261 mV and requires an overpotential of 366 mV at 10 mA cm –2 in alkaline electrolyte. These results illustrate how reticular-framework chemistry, junction engineering, and both built-in and applied fields can program reactivity through interfacial electrostatics and near-electrode transport in organic-framework electrocatalysts.

Garcia-Enriquez, Lissette [Univ. of Texas at El Pa↗
DOE OSTI2025

Temperature-Dependent Mechanical Properties of Ni-Based Concentrated Alloys: Insights from First-Principles Calculations

The present work focuses on predicting temperature-dependent mechanical properties of Ni-based concentrated alloys Ni 18 Cr 10 Co 10 Fe 6 M 4 (abbreviated by X 44 M 4 , with M = Al, V, Mn, Fe, Nb, Mo, and W) using density functional theory (DFT). These predictions are based on shear (plastic) and elastic deformations, utilizing the special quasirandom structure (SQS), the phonon-based quasiharmonic approach (QHA), and the quasistatic approach. The resulting properties include coefficient of thermal expansion via QHA, ideal shear strength (τ IS ), and stable and unstable stacking fault energies (γ SF and γ US ) through pure alias shear deformation, and elastic constants (c ij ), bulk modulus (B 0 ), and shear modules (G 0 ) via elastic deformation. Notably, predicting accurate γ SF is challenging due to uncertainties that can exceed the γ SF values. τ IS and γ US exhibit a strong linear relationship, enabling the accurate prediction of γ US based on the precisely determined τ IS . All mechanical properties of X 44 M 4 decrease with increasing temperature, except for some γ SF cases such as X 44 M 4 with M = V, Mn, Fe, Mo, and W. Among the X 44 M 4 alloys, X 44 Nb 4 exhibits the lowest τ IS , γ US , and G 0 values, and the highest B 0 /G 0 ratio, while X 44 Mn 4 has the lowest B 0 and B 0 /G 0 ratio. We found that volume is a crucial descriptor for understanding and modeling mechanical properties (except B0 and maybe also γ SF ) affected by alloying elements and temperature. Ni-based dilute alloys (e.g., Ni 11 M 1 and Ni 31 M 1 ) and concentrated alloys (e.g., X 44 M 4 ) show similar trends in mechanical properties influenced by alloying elements and temperature, simplifying the analysis and design of Ni-based alloys.

Elastic properties↗
NASA NTRS2019

Analysis of At-Altitude LTE Power Spectra for C2 Communications for UAS Traffic Management

The National Aeronautics and Space Administration’s (NASA) Unmanned Aircraft Systems Traffic Management (UTM) project works to develop tools and technologies essential for safely enabling civilian low-altitude small Unmanned Aerial Systems (sUAS, also known as drones) operations. This paper presents results of work completed in the paper [1] presented at the 2018 ICNS conference where proposed approaches were explored for evaluating and analyzing sUAS Command and Control (C2) links based on commercial cellular networks. This paper focuses on the UTM Project’s Technology Capability Level 3 (TCL-3) test results which address the communications portion identified within the same paper. A software defined radio (SDR) was flown as a sUAS payload to capture received signal spectrum in Long Term Evolution (LTE) frequency bands of interest. The purpose was to measure the RF environment at UTM altitudes to characterize the interference potential. The SDR payload was flown at various stationary altitudes where the LTE over-the-air complex (I/Q) samples were captured by the SDR and later post-processed. The SDR received inputs through an omnidirectional antenna. The complex samples captured were an aggregate of transmissions received from all line-of-sight (LOS) towers within the geographic area for the specific radio frequency bandwidth the SDR is programmed to capture. Using this approach, the complex samples captured do not distinguish between the various eNodeB's (Long Term Evolution (LTE) transmitting towers). The complex samples were post processed via a Discrete Fourier Transform (DFT) algorithm to view the captured spectrum along with the power levels across the captured LTE bandwidth. This SDR payload process of capturing complex samples was done at two different regions within the US: 1) NASA's Ames Research Center (ARC) in Moffett Field, CA, and 2) Griffiss Airfield in Rome, NY. The data capture at the ARC site was done at two physical locations within the Ames campus where many stationary altitude captures where done as high as 800 ft. above ground level (AGL). The data captured at the Griffiss Airport (also known as the NY Corridor Site) were acquired at one location with three specific stationary altitude levels – {Ground Level (GL), 300 ft., and 400 ft.}. The LTE spectrum power levels were captured for two LTE carriers, AT&T and Verizon, at both sites where their respective spectra and power levels were measured and compared at various altitudes. The overall results show that there is an increase in LTE spectrum power levels at higher altitudes for drones. A detailed analysis of this data and conclusions drawn from the results are presented in this paper.

Kerczewski, Robert J.↗
DOE OSTI2025

Single-Phase Spinel NiCo 2 O 4 as Highly Active and Stable Electrocatalysts for Urea Oxidation Reaction in Urea Electrolysis

Exploring and designing a stable and active catalyst for the urea electro-oxidation reaction (UOR, CO(NH 2 ) 2 + 6OH – → CO 2 + N 2 + 5H 2 O + 6e – ) is crucial for the long-term sustainability of ecological systems and clean energy production. We found that spinel NiCo 2 O 4 is a stable and active electrocatalyst for UOR at a relatively low anodic potential without triggering the competing oxygen evolution reaction (OER). A urea electrolysis cell (CO(NH 2 )2 + H 2 O → CO 2 + N 2 + 2H 2 ) utilizing a spinel NiCo 2 O 4 anode and a commercial Pt cathode was further characterized through galvanostatic polarization tests, demonstrating excellent structural stability at various current densities. Post-mortem analysis of long-term urea electrolysis measurements suggested that NiCo 2 O 4 electrocatalysts maintained a stable spinel structure. However, redistribution of Ni 3+ to Ni 2+ valence on the catalyst surface was observed, in contrast to the intact Co valence, indicating that (i) Ni sites are active toward urea adsorption and sequential electro-oxidation; (ii) while urea oxidation proceeds primarily through the direct electro-oxidation mechanism, chemical reactions between the Ni 3+ site and urea occur during long-term electrochemical UOR operation. Density functional theory (DFT) simulations were used to calculate the adsorption energies of urea molecules on NiO, Co 3 O 4 , and NiCo 2 O 4 , revealing the importance of regulating the configuration of adsorbed urea molecules on the NiCo 2 O 4 surface.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Investigation of Americium-Containing Phosphates, Silicates, Borates, Molybdates, and Fluorides Synthesized via High-Temperature Flux Crystal Growth

The crystal chemistry of americium-containing extended structures was investigated, and several classes of americium-containing solid-state oxide materials were obtained in single-crystal form via high-temperature flux crystal growth. This enabled the structural characterization of rare examples of ternary, quaternary, and penternary americium-containing silicates K 3 Am (Si 2 O 7 ) and Cs 6 Am 2 Si 21 O 48 , phosphates Na 3 Am (PO 4 ) 2 and K 3 Am (PO 4 ) 2 , borates Ba 3 Am 2 (BO 3 ) 4 and AmBO 3 , borate halides Ca 5 Am(BO 3 ) 4 Cl, molybdates Li 0.5 Am 0.5 MoO 4 , and fluorides CsAm 2 F 7 . Using these crystallographic data, the ionic radii of Am 3+ with coordination numbers of six (0.975 Å), seven (1.052 Å), and nine (1.162 Å) were established. A maximum entropy method (MEM) analysis was performed on the single-crystal X-ray diffraction data that were collected for K 3 Nd(PO 4 ) 2 /K 3 Am(PO 4 ) 2 , K 3 NdSi 2 O 7 /K 3 AmSi 2 O 7 , and NdBO 3 /AmBO 3 , to qualitatively compare the ionicities of the Nd–O and Am–O bonds. In conclusion, Raman spectroscopy data were collected on single crystals of K 3 Am(PO 4 ) 2 and compared to the calculated Raman spectrum of K 3 Am(PO 4 ) 2 obtained from DFT calculations.

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

Ethylene Polymerizations Catalyzed by Fluorinated “Sandwich” Diimine-Nickel and Palladium Complexes

Nickel and palladium complexes bearing “sandwich” diimine ligands with perfluorinated aryl caps have been synthesized, characterized, and explored in ethylene polymerization reactions. The X-ray crystallographic analysis of the precatalysts 16 and 6b show differences from their non-fluorinated analogues 17 and 19, with the perfluorinated aryl caps centered precisely over the nickel and palladium centers which results in higher buried volumes of the metal centers relative to the non-fluorinated analogs. The sandwich diimine palladium complexes 5a and 5b containing perfluorinated aryl caps polymerize ethylene in a controlled fashion with activities that are substantially increased compared with their non-fluorinated analogues. Migratory insertion rates in relevant methyl ethylene complexes agree with activities exhibited in bulk polymerization experiments. DFT studies suggest that facility of ethylene rotation from its preferred orientation perpendicular to the Pd-alkyl bond into a parallel in-plane conformation contributes to the higher polymerization activity for 5b relative to 18a. For these palladium systems polymer molecular weights can be controlled via hydrogen addition (hydrogenolysis) which is unusual for late transition metal-catalyzed olefin polymerizations with no catalyst deactivation occurring. Sandwich diimine nickel complexes 6a and 6b with perfluorinated aryl caps show ethylene polymerization activities that are about half that of classical tetraisopropyl-substituted catalyst 2 but again are more active than the analogous nonfluorinated sandwich complexes. Ethylene polymerizations exhibit living behavior and branched ultra-high molecular weight polyethylenes (UHMWPEs) with very low molecular weight distributions (less than 1.1) are obtained. The activated nickel catalysts are stable in the absence of monomer and show good long-term stability at 25 °C.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Computational Analysis of the Energetic Stability of High-Entropy Structures of a Prototypical Lanthanide-Based Metal–Organic Framework

High-entropy materials are characterized by their complex compositions, typically comprising five or more elements in near-equiatomic proportions. Applying this concept to metal ions in metal−organic frameworks (MOFs) has paved the way for exploring a new class of high-entropy MOFs. While the compositional strategy of high-entropy materials leverages configurational entropy to aid thermodynamic stability, it also poses significant analytical challenges due to the vast compositional landscape and diverse phases that these materials can adopt. We present a computational study of several complexities associated with selecting potential high-entropy versions of a prototype lanthanidebased MOF. We compute the energetics of metal mixing of these heterometallic MOFs using density functional theory (DFT) and machine learning interatomic potential (MLIP) methods. The use of MLIP methods allows a systematic exploration of the convex hull of thermodynamically stable MOF structures containing up to 5 distinct metals.

Chemical structure↗
DOE OSTI2025

Transient kinetic insights into selective propene oxidation over industrial bismuth molybdate catalysts

Selective oxidation of propene to acrolein over industrial multicomponent bismuth molybdate (BMO) catalysts significantly depends on reaction conditions that include operating parameters and catalyst state. Here, this work investigates selective propene oxidation in the intrinsic kinetic paradigm of Temporal Analysis of Products (TAP) reactor by systematically varying catalyst redox state, temperature, and oxygen-to-propene feed ratio. A 93 % propene conversion with an acrolein yield of 80 % is achieved at elevated temperatures (450 °C) on oxidized catalysts under oxygen-rich conditions (O 2 :C 3 H 6 = 10). However, these conditions diminish acrolein-to-CO 2 selectivity due to enhanced total oxidation to CO 2 . In contrast, a reduced catalyst state, moderate temperature (350 °C), and lower oxygen feed (O 2 :C 3 H 6 = 1) nearly doubles the acrolein to-CO 2 selectivity, albeit at a lower acrolein yield (33 %). Transient kinetic studies together with a kinetic model that simplify the major product formation pathways in lumped non-elementary forms reveal that the availability of surface oxygen species plays a pivotal role in governing reaction pathways. Additionally, density functional theory (DFT) calculations on pure BMO catalysts inform the role of surface redox states on propene and oxygen activation barriers. Lattice oxygen at acrolein-selective sites drives both acrolein and CO 2 formation, while adsorbed oxygen at activation sites favors unselective CO 2 generation. This work establishes a critical relationship between transient product selectivity and surface oxygen availability, which is strongly influenced by catalyst redox state, feed ratio, and reaction temperature. These insights underscore the importance of dynamic reactor operation strategies and offer a foundation for designing next-generation propene oxidation processes with tunable acrolein selectivity.

42 - ENGINEERING↗
DOE OSTI2024

Revisiting Artifacts of Kohn–Sham Density Functionals for Biosimulation

We revisit the problem of unphysical charge density delocalization/fractionalization induced by the self-interaction error of common approximate Kohn–Sham (KS) density functional theory functionals on simulation of small to medium-sized proteins in a vacuum. Aside from producing unphysical electron densities and total energies, the vanishing of the HOMO–LUMO gap associated with the unphysical charge delocalization leads to an unphysical low-energy spectrum and catastrophic failure of most popular solvers for the KS self-consistent field (SCF) problem. We apply a robust quasi-Newton SCF solver to obtain solutions for some of these difficult cases. The anatomy of the charge delocalization is revealed by the natural deformation orbitals obtained from the density matrix difference between the Hartree–Fock and KS solutions; the charge delocalization not only can occur between charged fragments (such as in zwitterionic polypeptides) but also involves neutral fragments. The vanishing-gap phenomenon and troublesome SCF convergence are both attributed to the unphysical KS Fock operator eigenspectra of molecular fragments (e.g., amino acids or their side chains). Analysis of amino acid pairs suggests that the unphysical charge delocalization can be partially ameliorated by the use of some range-separated hybrid functionals but not by semilocal or standard hybrid functionals. Last, we demonstrate that solutions without the unphysical charge delocalization can be located even for semilocal KS functionals highly prone to such defects, but such solutions have non-Aufbau character and are unstable with respect to mixing of the non-overlapping “frontier” orbitals. Caution should be exercised when unexpectedly small (or vanishing) HOMO–LUMO gaps and atypical SCF convergence patterns (e.g., oscillatory) are observed in KS DFT simulations in any context (bio or otherwise).

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

Quantitative Determination of Electronic Effects in Free Radicals through Open-Shell Hammett Substituent Constants

Hammett substituent constants (σ), which quantify the electronic effects of functional groups, are widely used for predicting the properties of organic compounds and investigating reaction mechanisms. While these values have been obtained for a wide range of closed-shell substituents, measurements of analogous values for open-shell substituents are rare due to challenges associated with their short lifetimes. Here, in this report, we developed a combined experimental and computational approach for quantifying the electronic properties of open-shell substituents based on changes in nitrile vibrational frequencies (ν(C≡N)). By coupling pulse radiolysis and time-resolved infrared spectroscopy (PR-TRIR), we measured ν(C≡N) IR bands of 30 para - and meta -substituted benzonitriles bearing C-, N-, and S-centered radicals. A linear scaling relationship was obtained between these experimental values and values obtained from DFT calculations. Using these computed values, two different Hammett constants, σ m , and σ p + , were determined for a series of C-, N-, O-, S-, Si-, and B-centered radicals. The differences between σm and σp+ values enable the separate evaluation of inductive and resonance effects in these open-shell substituents. The results suggest that there are three classes of radicals: one is electron withdrawing (σ m , σ p + > 0), one is electron donating (σm, σp+ < 0), and one is inductively withdrawing but resonance donating (σ m > 0, σ p + < 0). Our study represents a general approach to the analysis of the electronic properties of open-shell species and has potential applications in a wide range of molecular processes involving free radical intermediates.

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

Hydrated Metal and Metal-Nitrate Complexes in Water: Full Lanthanide(III) Series plus Miscellaneous Metal Ions

This is a dataset of hydrated metal complexes and metal–nitrate hydrated complexes intended for public use, reproducibility, and downstream structural analysis. A key feature is coverage across the full lanthanide(III) series (La–Lu), enabling systematic comparisons of coordination motifs and bonding trends across the entire lanthanide sequence. In addition to the lanthanides, the dataset also includes other metal ions such as UO2(VI), Fe(II), and Fe(III). The dataset provides optimized geometries for hydrated and nitrate-containing hydrated complexes, together with representative ab initio molecular dynamics (AIMD) trajectories saved in standard XYZ formats. The accompanying NWChem input decks enable reproduction of the reported calculations and provide a starting point for extending the simulations to related coordination environments. Computationally, DFT calculations employ the B3LYP functional with DFT-D3BJ dispersion corrections and a COSMO continuum solvent model (dielectric constant 78.4) to represent solvation beyond the explicitly treated first hydration shell. AIMD simulations are performed with the NWChem qmd module at 298 K, integrating nuclear motion with the velocity-Verlet algorithm and controlling temperature using a Nosé–Hoover thermostat. Trajectories are approximately 4.8 ps in length and are used primarily to assess short-time stability of candidate coordination motifs, including (for lanthanides) differences between 8- versus 9-water coordination and comparisons between nitrate-bound and nitrate-free hydrated complexes.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗