Search NASA⌕ Search

SEARCH · Search NASA

Results for “hybrid functional”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Reinforcement learning based hybrid bond-order coarse-grained interatomic potentials for exploring mesoscale aggregation in liquid–liquid mixtures

Exploring mesoscopic physical phenomena has always been a challenge for brute-force all-atom molecular dynamics simulations. Although recent advances in computing hardware have improved the accessible length scales, reaching mesoscopic timescales is still a significant bottleneck. Coarse-graining of all-atom models allows robust investigation of mesoscale physics with a reduced spatial and temporal resolution but preserves desired structural features of molecules, unlike continuum-based methods. Here, we present a hybrid bond-order coarse-grained forcefield (HyCG) for modeling mesoscale aggregation phenomena in liquid–liquid mixtures. The intuitive hybrid functional form of the potential offers interpretability to our model, unlike many machine learning based interatomic potentials. We parameterize the potential with the continuous action Monte Carlo Tree Search (cMCTS) algorithm, a reinforcement learning (RL) based global optimizing scheme, using training data from all-atom simulations. The resulting RL-HyCG correctly describes mesoscale critical fluctuations in binary liquid–liquid extraction systems. cMCTS, the RL algorithm, accurately captures the mean behavior of various geometrical properties of the molecule of interest, which were excluded from the training set. The developed potential model along with the RL-based training workflow could be applied to explore a variety of other mesoscale physical phenomena that are typically inaccessible to all-atom molecular dynamics simulations.

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

Characterization of Deformational Isomerization Potential and Interconversion Dynamics with Ultrafast X-ray Solution Scattering

Dimeric complexes composed of d 8 square planar metal centers and rigid bridging ligands provide model systems to understand the interplay between attractive dispersion forces and steric strain in order to assist the development of reliable methods to model metal dimer complexes more broadly. [Ir 2 (dimen) 4 ] 2+ (dimen = para-diisocyanomenthane) presents a unique case study for such phenomena, as distortions of the optimal structure of a ligand with limited conformational flexibility counteract the attractive dispersive forces from the metal and ligand to yield a complex with two ground state deformational isomers. Here, in this work, we use ultrafast X-ray solution scattering (XSS) and optical transient absorption spectroscopy (OTAS) to reveal the nature of the equilibrium distribution and the exchange rate between the deformational isomers. The two ground state isomers have spectrally distinct electronic excitations that enable the selective excitation of one isomer or the other using a femtosecond duration pulse of visible light. We then track the dynamics of the nonequilibrium depletion of the electronic ground state population—often termed the ground state hole—with ultrafast XSS and OTAS, revealing a restoration of the ground state equilibrium in 2.3 ps. This combined experimental and theoretical study provides a critical test of various density functional approximations in the description of bridged d 8 –d 8 metal complexes. The results show that density functional theory calculations can reproduce the primary experimental observations if dispersion interactions are added, and a hybrid functional, which includes exact exchange, is used.

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

Self-consistent implementation of locally scaled self-interaction-correction method

Recently proposed local self-interaction correction (LSIC) method is a one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange–correlation and Coulomb energy densities. The LSIC method is exact for the one-electron densities, also recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to the well-known Perdew–Zunger SIC (PZSIC) method as a special case. This article presents the self-consistent implementation of the LSIC method using the ratio of Weizsäcker and Kohn–Sham kinetic energy densities as an iso-orbital indicator. The atomic forces as well as the forces on the Fermi-Löwdin orbitals are also implemented for the LSIC energy functional. Results show that LSIC with the simplest local spin density functional predicts atomization energies of the AE6 dataset better than some of the most widely used generalized-gradient-approximation (GGA) functional and barrier heights of the BH6 database better than some of the most widely used hybrid functionals. The LSIC method [a mean absolute error (MAE) of 0.008 Å] predicts bond lengths of a small set of molecules better than the PZSIC-LSDA (MAE 0.042 Å) and LSDA (0.011 Å). Furthermore, this work shows that accurate results can be obtained from the simplest density functional by removing the self-interaction-errors using an appropriately designed SIC method.

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

Power Split Supercharging: A Mild Hybrid Approach to Boost Fuel Economy

This work investigates an innovative low-voltage (<60 V) hybrid device that enables engine boosting and downsizing in addition to mild hybrid functionalities such as regenerative braking, start-stop, and torque assist. A planetary gear set and a brake permit the power split supercharger (PSS) to share a 9 kW motor between supercharging the engine and direct torque supply to the crankshaft. In contrast, most e-boosting schemes use two separate motors for these two functionalities. This single motor structure restricts the PSS operation to only one of the supercharging or parallel hybrid modes; therefore, an optimized decision making strategy is necessary to select both the device mode and its power split ratio. An adaptive equivalent consumption minimization strategy (A-ECMS), which uses the battery state of charge (SoC) history to adjust the equivalence factor, is developed for energy management of the PSS. The A-ECMS effectiveness is compared against a dynamic programming (DP) solution with full drive cycle preview through hardware-in-the-loop experiments on an engine dynamometer testbed. The experiments show that the PSS with A-ECMS reduces vehicle fuel consumption by 18.4% over standard FTP75 cycle, compared to a baseline turbocharged engine, while global optimal DP solution decreases the fuel consumption by 22.8% compared to the baseline.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Optimizing accuracy and efficacy in data-driven materials discovery for the solar production of hydrogen

The production of hydrogen fuels, via water splitting, is of practical relevance for meeting global energy needs and mitigating the environmental consequences of fossil-fuel-based transportation. Water photoelectrolysis has been proposed as a viable approach for generating hydrogen, provided that stable and inexpensive photocatalysts with conversion efficiencies over 10% can be discovered, synthesized at scale, and successfully deployed. While a number of first-principles studies have focused on the data-driven discovery of photocatalysts, in the absence of systematic experimental validation, the success rate of these predictions may be limited. We address this problem by developing a screening procedure with co-validation between experiment and theory to expedite the synthesis, characterization, and testing of the computationally predicted, most desirable materials. Here, starting with 70 150 compounds in the Materials Project database, the proposed protocol yielded 71 candidate photocatalysts, 11 of which were synthesized as single-phase materials. Experiments confirmed hydrogen generation and favorable band alignment for 6 of the 11 compounds, with the most promising ones belonging to the families of alkali and alkaline-earth indates and orthoplumbates. This study shows the accuracy of a nonempirical, Hubbard-corrected density-functional theory method to predict band gaps and band offsets at a fraction of the computational cost of hybrid functionals, and outlines an effective strategy to identify photocatalysts for solar hydrogen generation.

08 HYDROGEN↗
DOE OSTI2024

Density functional theory of material design: fundamentals and applications—II

Abstract This is the second and the final part of the review on density functional theory (DFT), referred to as DFT-II. In the first review, DFT-I, we have discussed wavefunction-based methods, their complexity, and basics of density functional theory. In DFT-II, we focus on fundamentals of DFT and their implications for the betterment of the theory. We start our presentation with the exact DFT results followed by the concept of exchange-correlation (xc) or Fermi-Coulomb hole and its relationship with xc energy functional. We also provide the exact conditions for the xc-hole, xc-energy and xc-potential along with their physical interpretation. Next, we describe the extension of DFT for non-integer number of electrons, the piecewise linearity of total energy and discontinuity of chemical potential at integer particle numbers, and derivative discontinuity of the xc potential, which has consequences on fundamental gap of solids. After that, we present how one obtains more accurate xc energy functionals by going beyond the LDA. We discuss the gradient expansion approximation (GEA), generalized gradient approximation (GGA), and hybrid functional approaches to designing better xc energy functionals that give accurate total energies. However, these functionals fail to predict properties like the ionization potential and the band gap. Thus, we next describe different methods of modelling these potentials and results of their application for calculation of the band gaps of different solids to highlight accuracy of different xc potentials. Finally, we conclude with a glimpse on orbital-free density functional theory and the machine learning approach.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Efficient computational design of two-dimensional van der Waals heterostructures: Band alignment, lattice mismatch, and machine learning

Here, we develop a computational database, website applications (web-apps), and machine-learning (ML) models to accelerate the design and discovery of two-dimensional (2D) heterostructures. Using density functional theory (DFT) based lattice parameters and electronic band energies for 674 nonmetallic exfoliable 2D materials, we generate 226 779 possible bilayer heterostructures. We classify these heterostructures into type-I, -II, and -III systems according to Anderson’s rule, which is based on the relative band alignments of the noninteracting monolayers. We find that type II is the most common and type III the least common heterostructure type. We subsequently analyze the chemical trends for each heterostructure type in terms of the Periodic Table of constituent elements. The band alignment data can also be used for identifying photocatalysts and high-work-function 2D metals for contacts. We validate our results by comparing them to experimental data as well as hybrid-functional predictions. Additionally, we carry out DFT calculations of a few selected systems (MoS 2 /WSe 2 , MoS 2 /h-BN, and MoSe 2 /CrI 3 ), to compare the band-alignment description with the predictions from Anderson’s rule. We develop web-apps to enable users to virtually create combinations of 2D materials and predict their properties. Additionally, we use ML tools to predict band-alignment information for 2D materials. The web-apps, tools, and associated data will be distributed through the JARVIS-HETEROSTRUCTURE website. Our analysis, results, and the developed web-apps can be applied to the screening and design applications, such as finding alternative photocatalysts, photodetectors, and high-work-function (WF) 2D-metal contacts.

2-dimensional systems↗
NASA NTRS2018

The GMAO Hybrid Ensemble-Variational Atmospheric Data Assimilation System: Version 2.0

This document describes the implementation and usage of the Goddard Earth Observing System (GEOS) Hybrid Ensemble-Variational Atmospheric Data Assimilation System (Hybrid EVADAS). Its aim is to provide comprehensive guidance to users of GEOS ADAS interested in experimenting with its hybrid functionalities. The document is also aimed at providing a short summary of the state-of-science in this release of the hybrid system. As explained here, the ensemble data assimilation system (EnADAS) mechanism added to GEOS ADAS to enable hybrid data assimilation applications has been introduced to the pre-existing machinery of GEOS in the most non-intrusive possible way. Only very minor changes have been made to the original scripts controlling GEOS ADAS with the objective of facilitating its usage by both researchers and the GMAO's near-real-time Forward Processing applications. In a hybrid scenario two data assimilation systems run concurrently in a two-way feedback mode such that: the ensemble provides background ensemble perturbations required by the ADAS deterministic (typically high resolution) hybrid analysis; and the deterministic ADAS provides analysis information for recentering of the EnADAS analyses and information necessary to ensure that observation bias correction procedures are consistent between both the deterministic ADAS and the EnADAS. The nonintrusive approach to introducing hybrid capability to GEOS ADAS means, in particular, that previously existing features continue to be available. Thus, not only is this upgraded version of GEOS ADAS capable of supporting new applications such as Hybrid 3D-Var, 3D-EnVar, 4D-EnVar and Hybrid 4D-EnVar, it remains possible to use GEOS ADAS in its traditional 3D-Var mode which has been used in both MERRA and MERRA-2. Furthermore, as described in this document, GEOS ADAS also supports a configuration for exercising a purely ensemble-based assimilation strategy which can be fully decoupled from its variational component. We should point out that Release 1.0 of this document was made available to GMAO in mid-2013, when we introduced Hybrid 3D-Var capability to GEOS ADAS. This initial version of the documentation included a considerably different state-of-science introductory section but many of the same detailed description of the mechanisms of GEOS EnADAS. We are glad to report that a few of the desirable Future Works listed in Release 1.0 have now been added to the present version of GEOS EnADAS. These include the ability to exercise an Ensemble Prediction System that uses the ensemble analyses of GEOS EnADAS and (a very early, but functional version of) a tool to support Ensemble Forecast Sensitivity and Observation Impact applications.

Todling, Ricardo↗
DOE OSTI2022

Machine learning for impurity charge-state transition levels in semiconductors from elemental properties using multi-fidelity datasets

Quantifying charge-state transition energy levels of impurities in semiconductors is critical to understanding and engineering their optoelectronic properties for applications ranging from solar photovoltaics to infrared lasers. While these transition levels can be measured and calculated accurately, such efforts are time-consuming and more rapid prediction methods would be beneficial. Here, we significantly reduce the time typically required to predict impurity transition levels using multi-fidelity datasets and a machine learning approach employing features based on elemental properties and impurity positions. We use transition levels obtained from low-fidelity (i.e., local-density approximation or generalized gradient approximation) density functional theory (DFT) calculations, corrected using a recently proposed modified band alignment scheme, which well-approximates transition levels from high-fidelity DFT (i.e., hybrid HSE06). Further, the model fit to the large multi-fidelity database shows improved accuracy compared to the models trained on the more limited high-fidelity values. Crucially, in our approach, when using the multi-fidelity data, high-fidelity values are not required for model training, significantly reducing the computational cost required for training the model. Our machine learning model of transition levels has a root mean squared (mean absolute) error of 0.36 (0.27) eV vs high-fidelity hybrid functional values when averaged over 14 semiconductor systems from the II–VI and III–V families. As a guide for use on other systems, we assessed the model on simulated data to show the expected accuracy level as a function of bandgap for new materials of interest. Finally, we use the model to predict a complete space of impurity charge-state transition levels in all zinc blende III–V and II–VI systems.

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

Distinguishing erbium dopants in Y 2 O 3 by site symmetry: Ab initio theory of two spin-photon interfaces

Here, we present a first-principles study of the defect formation and electronic structure of erbium (Er)–doped yttria (Y 2 O 3 ). This is an emerging material for spin-photon interfaces in quantum information science due to the narrow-linewidth optical emission from Er dopants at standard telecommunication wavelengths and their potential for quantum memories and transducers. We calculate formation energies of neutral and negatively and positively charged Er dopants and find the charge-neutral configuration to be the most stable, consistent with experiment. Of the two substitutional sites of Er for Y, the C 2 (more relevant for quantum memories) and C 3i (more relevant for quantum transduction), we identify the former as possessing the lowest formation energy. The electronic properties are calculated using the Perdew-Burke-Ernzerhof functional along with the Hubbard U parameter and spin-orbit coupling, which yields a ~ 6 μ B orbital and a ~ 3 μ B spin magnetic moment, and 11 electrons in the Er 4 f shell, confirming the formation of charge-neutral Er 3+ . This standard density functional theory approach underestimates the band gap of the host and lacks a first-principles justification for U . To overcome these issues we performed screened hybrid functional calculations, including a negative U for the 4 f orbitals, with mixing (α) and screening (w) parameters. These produced robust electronic features with slight modifications in the band gap and the 4 f splittings depending on the choice of tuning parameters. We also computed the many-particle electronic excitation energies and compared them with experimental values from photoluminescence.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Data-driven design of novel halide perovskite alloys

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

36 MATERIALS SCIENCE↗
DOE OSTI2024

Anderson impurity mechanism for a multi-level model in δ-Pu

Abstract Electronic correlations and spin–orbit interactions in plutonium create variations in the bonding behavior of each of its allotropes. In δ -Pu, the 5f electrons lie at the tipping point between itinerant and localized behavior which has made the use of mixed-level models successful in describing its mechanical properties. The mechanism for the emergence of a mixed-level model has not yet been understood. We use a series of density functional theory approximations to understand the interactions that create a mixed-level description of δ -Pu which leads to accurate physical properties. With the intersite interactions present in the hybrid functional, we show that a single 5f electron engages in orbital-selective bonding that can be understood with an Anderson impurity picture. The Anderson model gives us a mechanism to understand how the bonding in δ -Pu evolves as a function of the interactions in the material such that we obtain both the accuracy and physics of the multi-level models from ab initio theory.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Modeling Liquid Water by Climbing up Jacob’s Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials

Within the framework of Kohn–Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob’s ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric–isothermal ensemble. To generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained by using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed by using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water toward the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional. However, some important discrepancies remain. The H-bond network predicted by the SCAN0 functional is still slightly overstructured compared to the experimental results.

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

Minimum conditions for accurate modeling of urea production via co-electrolysis

Co-electrolysis of carbon oxides and nitrogen oxides promise to simultaneously help restore the balance of the C and N cycles while producing valuable chemicals such as urea. However, co-electrolysis processes are still largely inefficient and numerous knowledge voids persist. Here, we provide a solid thermodynamic basis for modelling urea production via co-electrolysis. First, we determine the energetics of aqueous urea produced under electrochemical conditions based on experimental data, which enables an accurate assessment of equilibrium potentials and overpotentials. Next, we use density functional theory (DFT) calculations to model various co-electrolysis reactions producing urea. The calculated reaction free energies deviate significantly from experimental values for well-known GGA, meta-GGA and hybrid functionals. These deviations stem from errors in the DFT-calculated energies of molecular reactants and products. In particular, the error for urea is approximately -0.25 ± 0.10 eV. Finally, we show that all these errors introduce large inconsistencies in the calculated free-energy diagrams of urea production via co-electrolysis, such that gas-phase corrections are strongly advised.

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

Discovery of a Hybrid System for Photocatalytic CO 2 Reduction via Attachment of a Molecular Cobalt-Quaterpyridine Complex to a Crystalline Carbon Nitride

While recent reports have demonstrated the attachment of molecular catalysts to amorphous, graphitic carbon nitrides (g-CN) for light-driven CO 2 reduction, approaches to the utilization of crystalline carbon nitrides have remained undiscovered. Herein, a functional hybrid photocatalyst system has been found using a crystalline carbon nitride semiconductor, poly(triazine imide) lithium chloride (PTI-LiCl), with a surface-attached CoCl 2 (qpy-Ph-COOH) catalyst for CO 2 reduction. The molecular catalyst attaches to PTI-LiCl at concentrations from 0.10 to 4.30 wt % and exhibits ∼96% selectivity for CO production in a CO 2 -saturated, aqueous 0.5 M KHCO 3 solution. Optimal loadings were found to be within 0.42–1.04 wt % with rates between 1,400 and 1,550 μmol CO/g·h at an irradiance of 172 mW/cm 2 (λ = 390 nm) and apparent quantum yields of ∼2%. This optimized loading is postulated to represent a balance between maximal turnover frequency (TOF; 300+ h –1 ) and excess catalyst that can limit excited-electron lifetimes, as probed via transient absorption spectroscopy. An increase in the incident irradiance yields a concomitant increase in the TOFs and CO rates only for the higher catalyst loadings, reaching up to 2,149 μmol CO/g·h with a more efficient use of the catalyst surface capacity. The lower catalyst loadings, by comparison, already function at maximal TOFs. Higher surface loadings are also found to help mitigate deactivation of the molecular catalysts during extended catalytic testing (>24 h) owing to the greater net surface capacity for CO 2 reduction, thus representing an effective strategy to extend lifetime. The hybrid particles can be deposited onto an FTO substrate to yield ∼60% Faradaic efficiency for photoelectrochemical CO production at −1.2 V vs Ag/AgCl bias. In conclusion, these results demonstrate the synergistic combination of a crystalline carbon nitride with a molecular catalyst that achieves among the highest known rates in carbon-nitride systems for the light-driven CO 2 reduction to CO in aqueous solution with >95% selectivity.

CO2 reduction↗
DOE OSTI2021

Hybrid Porous Crystalline Materials from Metal Organic Frameworks and Covalent Organic Frameworks

Abstract Two frontier crystalline porous framework materials, namely, metal‐organic frameworks (MOFs) and covalent organic frameworks (COFs), have been widely explored owing to their outstanding physicochemical properties. While each type of framework has its own intrinsic advantages and shortcomings for specific applications, combining the complementary properties of the two materials allows the engineering of new classes of hybrid porous crystalline materials with properties superior to the individual components. Since the first report of MOF/COF hybrid in 2016, it has rapidly evolved as a novel platform for diverse applications. The state‐of‐art advances in the various synthetic approaches of MOF/COF hybrids are hereby summarized, together with their applications in different areas. Perspectives on the main challenges and future opportunities are also offered in order to inspire a multidisciplinary effort toward the further development of chemically diverse, multi‐functional hybrid porous crystalline materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Delocalization error poisons the density-functional many-body expansion

The many-body expansion is a fragment-based approach to large-scale quantum chemistry that partitions a single monolithic calculation into manageable subsystems. This technique is increasingly being used as a basis for fitting classical force fields to electronic structure data, especially for water and aqueous ions, and for machine learning. Here, we show that the many-body expansion based on semilocal density functional theory affords wild oscillations and runaway error accumulation for ion–water interactions, typified by F − (H 2 O) N with N ≳ 15. We attribute these oscillations to self-interaction error in the density-functional approximation. The effect is minor or negligible in small water clusters, explaining why it has not been noticed previously, but grows to catastrophic proportion in clusters that are only moderately larger. This behavior can be counteracted with hybrid functionals but only if the fraction of exact exchange is ≳50%, whereas modern meta-generalized gradient approximations including ωB97X-V, SCAN, and SCAN0 are insufficient to eliminate divergent behavior. Other mitigation strategies including counterpoise correction, density correction (i.e., exchange–correlation functionals evaluated atop Hartree–Fock densities), and dielectric continuum boundary conditions do little to curtail the problematic oscillations. In contrast, energy-based screening to cull unimportant subsystems can successfully forestall divergent behavior. These results suggest that extreme caution is warranted when the many-body expansion is combined with density functional theory.

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

Many-Body Benchmark of Electronic Charge and Spin Densities for Li 1–x NiO 2

Accurate benchmarks are particularly important for highly correlated oxides as mean-field approximations often fail to describe the subtle balance of charge transfer and magnetism in these materials with an accuracy comparable to experimental needs. Here we present accurate diffusion Monte Carlo (DMC) results of the electronic charge and spin densities for the tunable highly correlated oxide Li 1–x NiO 2 for x = 0, 1/2, and 1. To enable quantitative comparisons, we introduce a robust density-partitioning scheme, extending Voronoi analysis to assign atomic charges from spatially noisy DMC densities. We then benchmark common approximations used in density functional theory (DFT). Comparison against DMC shows that r 2 SCAN delivers the most balanced performance across charge, spin, and radial density descriptors, nearly reproducing DMC results for LiNiO 2 and apical Ni sites in Li 0.5 NiO 2 . Hybrid functionals (PBE0, SCAN0) perform unexpectedly poorly, and PBE + U + V yields inconsistent trends between charge and spin densities. Therefore, the r 2 SCAN functional minimizes errors relative to DMC while capturing the variable valence of the Ni ion and also retaining the computational efficiency of DFT for large-scale simulations of the tunable structural and electronic phases of Li1−xNiO2. Our study highlights the importance of accurate benchmarking of the fundamental quantities involved in DFT to select appropriate DFT approximations in order to advance the predictive modeling of charge-transfer-driven phenomena in correlated electron systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗