SEARCH · Search NASA

Results for “reference”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

DOE OSTI2025

High‐Speed Slot‐Die Coating with Donor‐Priority Rapid Aggregation Kinetics for Improved Morphology and Efficiency in Ecofriendly Organic Solar Cells

Abstract Solution‐processable organic solar cells (OSCs) represent a promising renewable photovoltaic technology with significant potential for eco‐compatible production. While high power conversion efficiencies (PCEs) have been achieved in OSCs, scaling this technology for high‐throughput manufacturing remains challenging. Key reason lies in the lack of efficient control strategies for the complex and long‐duration morphology evolution during high‐speed coating process with ecofriendly solvents. Here, a donor‐priority rapid aggregation process (DP‐RAP) scheme is proposed to solve this issue by adjusting the aggregation kinetics of donor and acceptor components. DP‐RAP enables blends with a nanoscale fiber network structure and favorable crystallinity, which contributes to balanced carrier transport and reduced recombination losses. As a result, the PCE is improved from 14.3% (reference) to 17.4% (DP‐RAP) for ultra‐high speed coated PM6:BTP‐eC9 devices in atmosphere, which is one of the highest values for non‐halogenated solvent‐processed solar cells at coating speeds of 500 mm s −1 . Moreover, the DP‐RAP based devices remain a stable PCE of approximately 17.4% across a broad range of coating speeds (20–500 mm s −1 ), illustrating its tolerance to the varied manufacturing conditions. This work highlights a promising avenue for the high‐speed, ecofriendly production of efficient OSCs, pushing the boundaries of practical manufacturing in renewable energy technologies.

Chemistry↗

DOE OSTI2024

Regulating the Third Metal to Design and Engineer Multilayered NiFeM (M: Co, Mn, and Cu) Nanofoam Anode Catalysts for Anion-Exchange Membrane Water Electrolyzers

In this study, alkaline anion-exchange membrane water electrolyzers (AEMWEs) for green hydrogen production have received intensive attention due to their feasibility of using earth-abundant platinum group metal (PGM)-free catalysts. Herein, the third metal is incorporated into NiFe-based catalysts to regulate their electronic structures and morphologies, aiming to achieve sufficient oxygen evolution reaction (OER) activity and performance in AEMWEs. The ternary NiFeM (M: Cu, Co, or Mn) catalysts are featured with multiple layered structures and nanofoam network morphologies, consisting of highly OER-active amorphous Ni-rich oxide shells and electrically conductive metallic alloy cores. The physical and electronic perturbations to the NiFe induced by a third element lead to a fine-tuning of the redox ability of the metal sites at the reaction centers, which breaks the scaling relationship between OH* and O* intermediates at the reaction centers. Thus, the unique structural configuration and electronic regulation simultaneously benefit catalytic activity and performance improvements. These NiFeM nanofoam catalysts demonstrated promising anode performance in actual AEMWEs, comparable to the IrO 2 reference, especially at high current densities. Notably, using various electrolytes (e.g., KOH solution or pure water) for AEMWEs exhibited a different performance trend among studied NiFeM catalysts, likely due to dynamic changes of catalysts under various OER environments.

25 ENERGY STORAGE↗

DOE OSTI2025

Molecular Interlocking Multidimensional Modulations of Cathode‐Electrolyte Interface for Constructing High Energy Density Quasi‐Solid‐State Batteries

Gel polymers are regarded as a promising candidate electrolyte for lithium-metal quasi-solid-state batteries, primarily due to their high ionic conductivity and solid-liquid synergistic properties. However, challenges such as interfacial side reactions, limitations in Li + transport caused by interfacial issues, and leaching of transition metals from the cathode have yet to be effectively solved. Herein, a novel gel electrolyte modulation strategy based on electrostatic filler assembly is proposed to address the issues of ineffective capacity utilization and inadequate cycling stability of high-energy-density cathode materials in solid-state lithium-ion batteries. It constructs a 3D interpenetrating charge-bridge network that effectively tackles the phase-separation challenge between fillers and electrolytes at the molecular level. Meanwhile, the molecular interlocking structure effectively inhibits the electrolyte erosion. More critically, it optimizes and stabilizes the cathode-electrolyte interface film, which facilitates the conduction of Li + -ions through a size-sieving mechanism. Consequently, this strategy enables effective adaptation across diverse high-energy-density cathode materials with satisfactory capacity performance (170.4 mAh g −1 at 4.5 V/1 C for LiNi 0.6 Co 0.2 Mn 0.2 O 2 and 194.0 mAh g −1 at 4.3 V/1 C for LiNi 0.9 Co 0.08 Mn 0.02 O 2 ). In conclusion, this investigation offers a straightforward and effective reference for addressing the critical challenges of ionic transport and interface stabilization in the design of gel electrolytes.

cathode-electrolyte interface↗

DOE OSTI2025

Breaking the Linear Scaling Relations for the Oxygen Reduction Reaction with a Dual‐Atom Catalyst Composed of a MnFe‐Porphyrrole Aerogel

Bimetallic catalysts offer enhanced catalytic performance through synergistic interactions between the two metals, allowing them to break the linear scaling relations and reach high electrocatalytic activity. This study presents bimetallic aerogel-based catalyst synthesized as a covalent, three-dimensional framework containing neighboring iron and manganese sites. The aerogel structure provides a high surface area and porosity, facilitating an ultra-high active site density and efficient mass transport. The MnFe porphyrrole's unique structure is obtained by alternately linking Mn-porphyrin and Fe-corrole complexes. It exhibited outstanding performance with an onset potential of 0.99 V RHE . Comparative studies with a free-base Fe porphyrrole catalyst (E onset 0.97 V RHE ) revealed that while Mn incorporation led to only a slight improvement in half-cell performance, it resulted in significantly enhanced performance in anion exchange membrane fuel cell. The MnFe catalyst achieved an OCV of 0.97 V and a peak power density of 0.27 W cm −2 , outperforming the free-base Fe counterpart. Using density functional theory calculations, we show that the higher ORR activity of MnFe-porphyrrole is due to charge transfer between Mn and Fe atoms, which is absent in the reference free-base Fe-porphyrrole. These findings underscore the advantages of bimetallic catalysts in improving ORR activity and fuel cell efficiency by leveraging synergistic effects.

Aerogel↗

DOE OSTI2025

Electronic Structure Distortions in Chromium Chelates Impair Redox Kinetics in Flow Batteries

Aminopolycarboxylate chelates are emerging as a promising class of electrolyte materials for aqueous redox flow batteries, offering tunable redox potentials, solubility, and pH stability through careful selection of ligands and transition metal ions. Despite their potential, the impact of molecular structure modifications on the electronic and electrochemical properties of these chelates remains underexplored. Here, in this study, we examine how introducing a hydroxyl group, often employed for its solubilizing properties, to the backbone of CrPDTA, a reference chelate material, significantly changes the thermodynamics and kinetics of the chelate's redox process. We correlate changes in molecular and electronic structures to different electrochemical responses resulting from the hydroxyl addition and show that the introduction of this functional group leads to a distortion in the octahedral coordination of chromium. Furthermore, increased anisotropic spin density and nonintegral oxidation state changes in the Cr metal center result in a larger barrier for electron transfer in CrPDTA‐OH. It is demonstrated that preserving a hexacoordinate chelate structure across a broad pH range is crucial for efficient flow battery application and it is emphasized that ligand modifications must avoid distorting the octahedral coordination of the transition metal.

25 ENERGY STORAGE↗

DOE OSTI2024

Metal Phosphide Nanoparticles Generated via a Molecular Precursor Route for Hydrotreatment of Methyl Laurate

Transition metal phosphide nanoparticles supported on silica were used as catalysts to investigate the hydrodeoxygenation of methyl laurate (used as a model compound for vegetable oils). Ni 2 P, Ni 1 Mo 1 P, and Ni 1.6 Mo 0.4 P were synthesized using a molecular precursor route. The nanoparticles were added to the silica support without any changes in their structure or particle size. As a reference and for comparison, MoP/SiO 2 was also prepared by the phosphite method. The prepared catalysts were characterized by X-ray diffraction, X-ray photoelectron spectroscopy, transmission electron microscopy, and chemical analysis. A higher conversion is reached with the MoP/SiO 2 monometallic catalyst, close to that with Ni 1 Mo 1 P/SiO 2 . However, this observation is misleading because the differences in dispersion obscure the result that the intrinsic activity (turnover frequency) of the bimetallic catalyst, prepared by a molecular precursor route, is significantly higher than that of the monometallic catalyst, prepared by the phosphite route, revealing a synergistic effect on the catalytic activity due to the formation of the bimetallic phosphide.

09 BIOMASS FUELS↗

DOE OSTI2025

Assessing the Performance of Exchange‐Correlation Density Functionals in Describing the Iron‐Catalyzed Ammonia Synthesis System

Density functional theory (DFT) has been widely employed for elucidating mechanistic aspects of heterogeneous catalysis. However, the accuracy of DFT calculations relies heavily on selecting an exchange-correlation (XC) functional that correctly describes the electronic structure of materials involved in the reactions. This study assesses the accuracy of several XC density functionals for modeling the iron-catalyzed ammonia synthesis system. In the assessment of functional accuracy, experimental references are compared to DFT-calculated values for the formation energy of gas-phase ammonia and nitrogen, bulk Fe/Fe-nitride (γ′-Fe 4 N) lattice constants and cohesive/formation energies, and nitrogen and ammonia binding energies on Fe(100), Fe(111), Fe(110), and γ′-Fe 4 N(111). It is observed that the experimental value for each of these descriptors is accurately modeled by at least one functional. RPBE alone provides reliable estimates for both the lattice constant and cohesive energy of Fe and γ′-Fe 4 N. Temperature-programmed desorption experiments led to estimates for N and NH 3 adsorption across several Fe-based facets that are best captured by RPBE. These results highlight the importance of choosing an appropriate XC functional that accurately describes experimental systems and offer insights into effectively modeling the interactions between nitrogen and ammonia on Fe-based surfaces.

adsorption↗

DOE OSTI2024

Large language model evaluation for high–performance computing software development

We apply AI-assisted large language model (LLM) capabilities of GPT-3 targeting high-performance computing (HPC) kernels for (i) code generation, and (ii) auto-parallelization of serial code in C ++, Fortran, Python and Julia. Our scope includes the following fundamental numerical kernels: AXPY, GEMV, GEMM, SpMV, Jacobi Stencil, and CG, and language/programming models: (1) C++ (e.g., OpenMP [including offload], OpenACC, Kokkos, SyCL, CUDA, and HIP), (2) Fortran (e.g., OpenMP [including offload] and OpenACC), (3) Python (e.g., numpy, Numba, cuPy, and pyCUDA), and (4) Julia (e.g., Threads, CUDA.jl, AMDGPU.jl, and KernelAbstractions.jl). Kernel implementations are generated using GitHub Copilot capabilities powered by the GPT-based OpenAI Codex available in Visual Studio Code given simple + + prompt variants. To quantify and compare the generated results, we propose a proficiency metric around the initial 10 suggestions given for each prompt. For auto-parallelization, we use ChatGPT interactively giving simple prompts as in a dialogue with another human including simple “prompt engineering” follow ups. Results suggest that correct outputs for C++ correlate with the adoption and maturity of programming models. For example, OpenMP and CUDA score really high, whereas HIP is still lacking. We found that prompts from either a targeted language such as Fortran or the more general-purpose Python can benefit from adding language keywords, while Julia prompts perform acceptably well for its Threads and CUDA.jl programming models. Finally, we expect to provide an initial quantifiable point of reference for code generation in each programming model using a state-of-the-art LLM. Overall, understanding the convergence of LLMs, AI, and HPC is crucial due to its rapidly evolving nature and how it is redefining human-computer interactions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Gas–Phase Preparation of the 14π Hückel Polycyclic Aromatic Anthracene and Phenanthrene Isomers (C 14 H 10 ) via the Propargyl Addition–BenzAnnulation (PABA) Mechanism

Polycyclic aromatic hydrocarbons (PAHs) imply the missing link between resonantly stabilized free radicals and carbonaceous nanoparticles, commonly referred to as soot particles in combustion systems and interstellar grains in deep space. Whereas gas phase formation pathways to the simplest PAH – naphthalene (C 10 H 8 ) – are beginning to emerge, reaction pathways leading to the synthesis of the 14π Hückel aromatic PAHs anthracene and phenanthrene (C 14 H 10 ) are still incomplete. Here, by utilizing a chemical microreactor in conjunction with vacuum ultraviolet (VUV) photoionization (PI) of the products followed by detection of the ions in a reflectron time-of-flight mass spectrometer (ReTOF-MS), the reaction between the 1'- and 2'-methylnaphthyl radicals (C 11 H 9 •) with the propargyl radical (C 3 H 3 •) accesses anthracene (C 14 H 10 ) and phenanthrene (C 14 H 10 ) via the Propargyl Addition–BenzAnnulation (PABA) mechanism in conjunction with a hydrogen assisted isomerization. Furthermore, the preferential formation of the thermodynamically less stable anthracene isomer compared to phenanthrene suggests a kinetic, rather than a thermodynamics control of the reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Converting Second‐Order Saddle Points to Transition States: New Principles for the Design of 4π Photoswitches

Abstract Molecular photoswitches have demonstrated potential for storing solar energy at the molecular level, with power densities comparable to commercial batteries and hydroelectric energy storage. However, development of efficient photoswitches is hindered by limitations in cyclability and optical properties of existing materials. We here demonstrate that certain limitations in photoswitches based on electrocyclizations stem from the issue of controlling competition between Woodward‐Hoffmann allowed and forbidden pathways. Our approach moves beyond the traditional view of activation barriers and reveals that second‐order saddle points are crucial in dictating the competition between disrotatory and conrotatory pathways. These insights suggest new opportunities to manipulate the competition between these pathways through geometric constraints, fundamentally altering the connectivity of the potential energy surface. Our study also emphasizes the necessity of multi‐reference methods and the need to conduct higher‐dimensional explorations for competing pathways beyond photoswitch design.

Chemistry↗

DOE OSTI2024

Allometric relationships and trade‐offs in 11 common M editerranean‐climate grasses

Abstract Biomass allocation in plants is the foundation for understanding dynamics in ecosystem carbon balance, species competition, and plant–environment interactions. However, existing work on plant allometry has mainly focused on trees, with fewer studies having developed allometric equations for grasses. Grasses with different life histories can vary in their carbon investment by prioritizing the growth of specific organs to survive, outcompete co‐occurring plants, and ensure population persistence. Further, because grasses are important fuels for wildfire, the lack of grass allocation data adds uncertainty to process‐based models that relate plant physiology to wildfire dynamics. To fill this gap, we conducted a greenhouse experiment with 11 common California grasses varying in photosynthetic pathway and growth form. We measured plant sizes and harvested above‐ and belowground biomass throughout the life cycle of annual species, while for the establishment stage of perennial grasses to quantify allometric relationships for leaf, stem, and root biomass, as well as plant height and canopy area. We used basal diameter as a reference measure of plant size. Overall, basal diameter is the best predictor for leaf and stem biomass, height, and canopy area. Including height as another predictor can improve model accuracy in predicting leaf and stem biomass and canopy area. Fine root biomass is a function of leaf biomass alone. Species vary in their allometric relationships, with most variation occurring for plant height, canopy area, and stem biomass. We further explored potential trade‐offs in biomass allocation across species between leaf and fine root, leaf and stem, and allocation to reproduction. Consistent with our expectation, we found that fast‐growing plants allocated a greater fraction to reproduction. Additionally, plant height and specific leaf area negatively influenced the leaf‐to‐stem ratio. However, contrary to our hypothesis, there were no differences in root‐to‐leaf ratio between perennial and annual or C 4 and C 3 plants. Our study provides species‐specific and functional‐type‐specific allometry equations for both above‐ and belowground organs of 11 common California grass species, enabling nondestructive biomass assessment in California grasslands. These allometric relationships and trade‐offs in carbon allocation across species can improve ecosystem model predictions of grassland species interactions and environmental responses through differences in morphology.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Clarifying the trophic state concept to advance macroscale freshwater science and management

For over a century, ecologists have used the concept of trophic state (TS) to characterize an aquatic ecosystem's biological productivity. However, multiple TS classification schemes, each relying on a variety of measurable parameters as proxies for productivity, have emerged to meet use‐specific needs. Frequently, chlorophyll a, phosphorus, and Secchi depth are used to classify TS based on autotrophic production, whereas phosphorus, dissolved organic carbon, and true color are used to classify TS based on both autotrophic and heterotrophic production. Both classification approaches aim to characterize an ecosystem's function broadly, but with varying degrees of autotrophic and heterotrophic processes considered in those characterizations. Moreover, differing classification schemes can create inconsistent interpretations of ecosystem integrity. For example, the US Clean Water Act focuses exclusively on algal threats to water quality, framed in terms of eutrophication in response to nutrient loading. This usage lacks information about non‐algal threats to water quality, such as dystrophication in response to dissolved organic carbon loading. Consequently, the TS classification schemes used to identify eutrophication and dystrophication may refer to ecosystems similarly (e.g., oligotrophic and eutrophic), yet these categories are derived from different proxies. These inconsistencies in TS classification schemes may be compounded when interdisciplinary projects employ varied TS frameworks. Even with these shortcomings, TS can still be used to distill information on complex aquatic ecosystem function into a set of generalizable expectations. The usefulness of distilling complex information into a TS index is substantial such that usage inconsistencies should be explicitly addressed and resolved. To emphasize the consequences of diverging TS classification schemes, we present three case studies for which an improved understanding of the TS concept advances freshwater research, management efforts, and interdisciplinary collaboration. To increase clarity in TS, the aquatic sciences could benefit from including information about the proxy variables, ecosystem type, as well as the spatiotemporal domains used to classify TS. As the field of aquatic sciences expands and climatic irregularity increases, we highlight the importance of re‐evaluating fundamental concepts, such as TS, to ensure their compatibility with evolving science.

classification↗

DOE OSTI2023

2,5‐Dimercapto‐1,3,4‐Thiadiazole (DMCT)‐Based Polymers for Rechargeable Metal–Sulfur Batteries

Organosulfur materials are a sustainable alternative to the present-day layered oxide cathodes in lithium-based batteries. One such organosulfur material that was intensely explored from the 1990s to early 2010s is 2,5-dimercapto-1,3,4-thiadiazole (DMCT). However, research interest declined as the electrode reactions with DMCT were assumed to be too sluggish to be practical. Armed with the advances in metal–sulfur batteries, we revisit DMCT-based materials in the form of poly[tetrathio-2,5-(1,3,4-thiadiazole)], referred to as pDMCT-S. With an appropriate choice of electrode design and electrolyte, pDMCT-S cathode paired with a Li-metal anode shows a capacity of 715 mA h g −1 and a Coulombic efficiency of 97.7% at a C/10 rate, thus quelling the concerns of sluggish reactions. Surprisingly, pDMCT-S shows significantly improved long-term cyclability compared to a sulfur cathode. Investigations into the origin of the stability reveals that the discharge product Li-DMCT in its mesomeric form can strongly bind to polysulfides, preventing their dissolution into the electrolyte and shuttling. This unique mechanism solves a critical problem faced by sulfur cathodes. Encouragingly, this mechanism results in a stable performance of pDMCT-S with Na-metal cells as well. In conclusion, this study opens the potential for exploring other organic materials that have inherent polysulfide sequestering capabilities, enabling long-life metal–sulfur batteries.

2,5-dimercapto-1,3,4-thiadiazole↗

DOE OSTI2024

Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations

One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

1D assumption↗

DOE OSTI2024

JAMES BUTTLE REVIEW: Interflow, subsurface stormflow and throughflow: A synthesis of field work and modelling

Interflow, throughflow and subsurface stormflow are interchangeable terms that refer to the lateral subsurface flow above a restricting layer of lower hydraulic conductivity that occurs during and following storm events. Interflow (used here) is a more dominant process in steeper catchments with high infiltration capacity soils overlying a more impermeable soil or geologic layer. Interflow as a runoff process was first recognised in the early 1900s, yet hydrologists still struggle to predict its occurrence, persistence, importance, interaction with other streamflow generation processes, and potential to connect to valleys and streams during and following storms. We review the history of interflow research and address some of the challenges in understanding its role in runoff production. We argue that characterising the controls on interflow initiation and occurrence relies on detailed field observations of subsurface properties, which exist only in limited experimental settings. This data shortcoming contributes to our inability to predict interflow or determine its contribution to streamflow more broadly. There remain many opportunities to advance our understanding of interflow that include both modelling and experimental or observational approaches in hydrology.

hillslope hydrology↗

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 OSTI2025

Projector‐Based Quantum Embedding Study of Iron Complexes

Projection‐based embedding theory (PBET) is used to calculate and assess the challenging spin‐crossover energies for a selection of small Fe‐containing systems by embedding the metal center into the frozen potential of the ligands. MP2, CCSD, and CCSD(T) are embedded in potentials from the SCAN and r 2 SCAN functionals and compared with the canonical values for the constituent methods and previously reported reference values. Considering the PBET calculations as a correction for the underlying DFT, the embedding calculations are able to provided improvement for most cases. In some cases, the PBET methods are able to compensate for limitations in the wave function methods and produce results similar to more rigorous calculations from the literature. For the systems with spin‐crossover energies near zero, the current methodology fails to provide consistent improvement. In conclusion, the isolated recalculation of the electronic structure around the metal center when embedded into a DFT treatment of the ligand field shows promise as a pragmatic and lower cost treatment compared to the canonical treatment of the whole system of the difficult class of spin‐crossover complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗