Search NASA⌕ Search

SEARCH · Search NASA

Results for “optimized energy density functionals”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

An Analytical Tool to Evaluate Defect Thermodynamics of (La,Ba)Fe1-xMxO3-δ Perovskites for Solid-Oxide Cell Applications

A modeling tool of the defect thermodynamics of (La,Ba)Fe1-xMxO3-δ perovskites which includes energetic information about oxygen vacancy formation, hydration, hydride formation, and charge disproportionation reactions has been developed1. This tool incorporates defect energies and entropies expressed as sixth order polynomial functions to allow refinements of the defect reaction equilibrium constants in the thermodynamic analysis. Calculation of (La,Ba)Fe1-xMxO3-δ Brouwer diagrams as a function of pO2/pH2O and pH2/H2O in a range of temperatures of interest is facilitated by this modeling tool. The results obtained can provide direct guidance how the electronic and ionic defect concentrations of the triple conducting perovskite materials can be used to optimize performance of solid oxide cells for energy applications. The impact of magnetic and electronic structures of the perovskites on the defect reaction energies and entropies as obtained from density function theory modeling and the role played by hydride defect species will also be discussed.

Lee, Yueh-Lin↗
DOE OSTI2023

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Engineering Thermally Resilient and Kinetically Active Reversible Protonic Ceramic Cells via Interfacial Design

Achieving concurrent fast electrode kinetics and long-term thermo-mechanical durability remains a critical challenge for reversible protonic ceramic electrochemical cells (R-PCECs). Herein, we report an interfacial engineering strategy that integrates a perovs.kite-type PrBaRu0.1Co1.9O5+δ (PBRC) nanoparticle layer onto a PrBa0.5Sr0.5Co1.5Fe0.5O5+δ (PBSCF) substrate (PBRC-PBSCF), together with a modified pellet-assisted sintering approach to fabricate dense BaZr0.4Ce0.4Y0.1Yb0.1O3-δ (BZCYYb4411) electrolytes. The in situ reconstructed heterointerface enhances oxygen reduction/evolution reaction (ORR/OER) kinetics, promotes H2O adsorption/dissociation, and improves steam tolerance, as verified by electrochemical measurements and interfacial microstructural analyses. Density functional theory reveals that Ru-induced electronic modulation at the PBRC-PBSCF interface lowers the energy of oxygen vacancy formation and optimizes the position of the O 2p band center, thereby accelerating oxygen redox kinetics and stabilizing the interface. The resulting R-PCECs deliver an excellent peak power density of 1.112 W cm−2 and an electrolysis current density of −1.257 A cm−2 at 1.3 V in 3% H2O wet air at 600°C, with a reasonable faradaic efficiency. Furthermore, the cells demonstrate excellent stability, sustaining 100 h of thermal cycling (400–600°C, 200°C h−1) in both fuel cell and electrolysis modes, with 600 h of stability in electrolysis mode (600°C, −0.5 to −2 A cm−2).

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

Predicting Excitation Energies of Twisted Intramolecular Charge-Transfer States with the Time-Dependent Density Functional Theory: Comparison with Experimental Measurements in the Gas Phase and Solvents Ranging from Hexanes to Acetonitrile

Electronically excited states characterized by intramolecular charge transfer play an essential role in many biological processes and optical devices. The ability to make quantitative ab initio predictions of the relative energetics involved is a challenging yet desirable goal, especially for large molecules in solution. In this work, we present a data set of 61 experimental measurements of absorption and emission processes, both in the gas phase and in solvents representing a broad range of polarities, which involve intramolecular charge transfer mediated by a nonzero, "twisted" dihedral angle between one or more donor and acceptor subunits. Among a variety of density functionals investigated within the framework of linear-response theory, the "optimally tuned" LRC-ωPBE functional, which utilizes a system-specific yet nonempirical procedure to specify the range-separation parameter, emerges as the preferred choice. For the entire set of excitation energies, involving changes in dipole moment ranging from 4 to >20 Debye, the mean signed and absolute errors are 0.02 and 0.18 eV, respectively (compared, e.g., to -0.30 and 0.30 for PBE0, 0.44 and 0.47 for LRC-ωPBEh, 0.83 and 0.83 for ωB97X-V). We analyze the performance of polarizable continuum solvation models available in Q-Chem that partition the solvent response into fast and slow time scales, and clear trends emerge when measurements corresponding to the four small 4-(dimethylamino)benzonitrile (DMABN)-like molecules and a charged species are excluded. We make the case that the large errors found only for small molecules in the gas phase and weak solvents cannot be expected to improve via the optimal tuning procedure, which enforces a condition that is exact only in the well-separated donor-acceptor limit, and present empirical evidence implicating the outsized importance for small donor-acceptor systems of relaxation effects that cannot be accounted for by the linear-response time-dependent density functional theory within the adiabatic approximation. Finally, we demonstrate the utility of the optimally tuned density functional approach by targeting the charge-transfer states of a large biomimetic model system for light-harvesting structures in Photosystem II.

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

Bayesian optimization of entropy-stabilized metal fluoride conversion cathodes and their synthesis

Fluoride-based conversion cathodes are promising for next-generation Li ion batteries because of their high voltage and energy densities. High entropy earth abundant fluoride cathodes are attractive for sustainable battery because of the potentially higher cycling stability. While Ni and Co containing equimolar alloys stabilized by high entropy were considered previously, Ni and Co free non-equimolar metal alloy compositions might offer new opportunities to increase materials stability and average voltage. This opportunity was pursued using Bayesian optimization of average voltage for the computational prediction of (Fe x1 Cu x2 Zn x3 Mn x4 Mg x5 )F 2 cathode chemistries. Theoretical voltage, room temperature free energy, and electronic structure of potentially high voltage compositions were computed using density functional theory calculations. Here, the experimental realization of non-equimolar quinary Co and Ni free fluorides with stable rutile single phase and improved conversion reaction voltage by up to 28 % compared to the equimolar alloy demonstrated the power of Bayesian optimization methods for the computationally guided experimental synthesis of fluoride cathodes. Finally, non-equimolar cathode chemistries were found with cycling stability superior to known multi-element fluorides.

25 ENERGY STORAGE↗
DOE OSTI2023

High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation

Quantum chemical calculations on atomistic systems have evolved into a standard approach to studying molecular matter. These calculations often involve a significant amount of manual input and expertise, although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput electronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject them to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all components of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow with the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.

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

Shadow energy functionals and potentials in Born–Oppenheimer molecular dynamics

In Born–Oppenheimer molecular dynamics (BOMD) simulations based on the density functional theory (DFT), the potential energy and the interatomic forces are calculated from an electronic ground state density that is determined by an iterative self-consistent field optimization procedure, which, in practice, never is fully converged. The calculated energies and forces are, therefore, only approximate, which may lead to an unphysical energy drift and instabilities. Here, we discuss an alternative shadow BOMD approach that is based on backward error analysis. Instead of calculating approximate solutions for an underlying exact regular Born–Oppenheimer potential, we do the opposite. Instead, we calculate the exact electron density, energies, and forces, but for an underlying approximate shadow Born–Oppenheimer potential energy surface. In this way, the calculated forces are conservative with respect to the approximate shadow potential and generate accurate molecular trajectories with long-term energy stabilities. We show how such shadow Born–Oppenheimer potentials can be constructed at different levels of accuracy as a function of the integration time step, δt, from the constrained minimization of a sequence of systematically improvable, but approximate, shadow energy density functionals. For each energy functional, there is a corresponding ground state Born–Oppenheimer potential. These pairs of shadow energy functionals and potentials are higher-level generalizations of the original “zeroth-level” shadow energy functionals and potentials used in extended Lagrangian BOMD. The proposed shadow energy functionals and potentials are useful only within this extended dynamical framework, where also the electronic degrees of freedom are propagated as dynamical field variables together with the atomic positions and velocities. The theory is quite general and can be applied to MD simulations using approximate DFT, Hartree–Fock, or semi-empirical methods, as well as to coarse-grained flexible charge models.

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

Active Learning Accelerated Discovery of Stable Iridium Oxide Polymorphs for the Oxygen Evolution Reaction

The discovery of high-performing and stable materials for sustainable energy applications is a pressing goal in catalysis and materials science. Understanding the relationship between a material’s structure and functionality is an important step in the process, such that viable polymorphs for a given chemical composition need to be identified. Machine-learning-based surrogate models have the potential to accelerate the search for polymorphs that target specific applications. Herein, we report a readily generalizable active-learning (AL) accelerated algorithm for identification of electrochemically stable iridium oxide polymorphs of IrO 2 and IrO 3 . The search is coupled to a subsequent analysis of the electrochemical stability of the discovered structures for the acidic oxygen evolution reaction (OER). Structural candidates are generated by identifying all 956 structurally unique AB2 and AB3 prototypes in existing materials databases (more than 38000). Next, using an active learning approach, we find 196 IrO 2 polymorphs within the thermodynamic amorphous synthesizability limit and reaffirm the global stability of the rutile structure. We find 75 synthesizable IrO 3 polymorphs and report a previously unknown FeF 3 -type structure as the most stable, termed α-IrO 3 . To test the algorithms performance, we compare to a random search of the candidate space and report at least a 2-fold increase in the rate of discovery. Additionally, the AL approach can acquire the most stable polymorphs of IrO 2 and IrO 3 with fewer than 30 density functional theory optimizations. Analysis of the structural properties of the discovered polymorphs reveals that octahedral local coordination environments are preferred for nearly all low-energy structures. Subsequent Pourbaix Ir–H 2 O analysis shows that α-IrO3 is the globally stable solid phase under acidic OER conditions and supersedes the stability of rutile IrO 2 . Calculation of theoretical OER surface activities reveal ideal weaker binding of the OER intermediates on α-IrO 3 than on any other considered iridium oxide. We emphasize that the proposed AL algorithm can be easily generalized to search for any binary metal oxide structure with a defined stoichiometry.

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

High-shade dryland agrivoltaic conditions enhanced carbon uptake and water-use efficiency in zucchini ( Cucurbita pepo )

Introduction: The increasing global demand for food and energy is intensifying land-use competition. Agrivoltaic systems are a multifunctional land-use approach that vertically integrates the production of agricultural crops and solar power on the same land area. Most food crops are adapted to full-sun conditions, and the physiological responses of these crops to the novel microclimate under solar panels remain poorly understood. We hypothesized that the microclimate beneath the high-density photovoltaic system would influence carbon uptake, water use, and yield outcomes of zucchini summer squash.Methods: We conducted a field experiment in a hot, semi-arid climate on zucchini (Cucurbita pepo). Plants were grown under an agrivoltaic system with a 75% ground cover ratio (GCR) and in a full-sun control plot, each with two irrigation regimes (100 and 50%). We measured leaf-level photosynthesis, microclimate variables, and fruit yield at plant maturity and throughout the growing season.Results: The agrivoltaic array reduced photosynthetically active radiation (PAR) by ~79%, resulting in a cooler (−1.1 °C), more humid environment with higher soil moisture. These microclimatic conditions enhanced midday photosynthesis and daily cumulative carbon uptake. However, fruit yield was consistently lower under the panels, indicating a shift in carbon allocation toward vegetative growth. Photosynthesis was primarily driven by PAR across treatments, while soil moisture significantly influenced photosynthesis only in the control plots, suggesting water limitation was alleviated under the panels.Discussion: These findings highlight a trade-off between improved physiological performance and reduced yield under high-density agrivoltaics. While the system buffered heat and drought stress and improved overall plant function, excessive shade reduced reproductive output. Optimizing panel density or selecting crops cultivated for non-fruit yields will be essential for balancing food production and energy generation in dryland agrivoltaic settings.

14 SOLAR ENERGY↗
DOE OSTI2020

A fundamental look at electrocatalytic sulfur reduction reaction

The fundamental kinetics of the electrocatalytic sulfur reduction reaction (SRR), a complex 16-electron conversion process in lithium–sulfur batteries, is so far insufficiently explored. In this study, by directly profiling the activation energies in the multistep SRR, we reveal that the initial reduction of sulfur to the soluble polysulfides is relatively easy owing to the low activation energy, whereas the subsequent conversion of the polysulfides into the insoluble Li 2 S 2 /Li 2 S has a much higher activation energy, contributing to the accumulation of polysulfides and exacerbating the polysulfide shuttling effect. We use heteroatom-doped graphene as a model system to explore electrocatalytic SRR. We show that nitrogen and sulfur dual-doped graphene considerably reduces the activation energy to improve SRR kinetics. Density functional calculations confirm that the doping tunes the p-band centre of the active carbons for an optimal adsorption strength of intermediates and electroactivity. This study establishes electrocatalysis as a promising pathway to tackle the fundamental challenges facing lithium–sulfur batteries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1995

The Structure and Stability of Bn(+) Clusters

The geometries of B+n clusters for n less than 14 have been optimized using density functional theory with the B3LYP functional. The most stable structure for each cluster is planar or quasi-planar. The B3LYP fragmentation energies are calibrated using coupled cluster theory. Overall, our corrected fragmentation energies are in reasonable agreement with experiment. Our results are compared with previous theoretical results.

Ricca, Alessandra↗
DOE OSTI2022

Combined Experimental and Computational Study of Molybdenum and Niobium for Nuclear Sensor Application

Due to their novel electromagnetic and thermal properties, molybdenum (Mo) and niobium (Nb) become optimal temperature sensor materials for nuclear energy applications. We leveraged voltage recorded during a heat ramp to tune a computational method to predict the Seebeck electromotive force (EMF) of Mo and Nb. Using a combined Density Functional Theory (DFT) and Boltzmann Transport Equations (BTE) method the voltage was predicted but did not include the effects of temperature on atomic structure. Combining Ab Initio Molecular Dynamics (AIMD) and BTE included temperature effects on structure optimization and yielded voltages in a good agreement with experiment. Lanthanum (La) and Phosphorus (P) additives in Mo and Nb, respectively, could increase the EMF compared to those of the pure metals. The presence of oxygen (O) in Mo increases the EMF while O in Nb slightly reduces the EMF. Furthermore, our studies suggested that heat treatment-induced structural changes that lead to a reduction in voltage occur not only at the mesoscale as previously understood but also at the atomic scale.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Structured for success: conjugated polymer binders with tailored composition and architecture for lithium-ion batteries

Conjugated polymer binders are replacing conventional binders in lithium-ion batteries. Herein, we examine how molecular engineering and hierarchical nanostructuring govern binder functionality and electrochemical performance. Lithium-ion batteries (LIBs) are the leading energy storage technology, yet enhancing their energy density and cycle life remains critical. Significant progress has been made in high-capacity anodes and high-voltage cathodes, but their performance is hindered by electrode degradation, where it is related to the behaviors of binders at the surface and interface. Conventional non-conductive binders like poly(vinylidene difluoride) (PVDF), combined with conductive additives, often fail to maintain electrical pathways under repeated volume changes. Alternatively, conjugated polymer binders have emerged as a superior alternative, simultaneously offering intrinsic conductivity, mechanical flexibility, and strong adhesion through π-conjugated backbones and functional groups. Their tunable molecular structure enables efficient electron/ion transport while mitigating electrode cracking. Additionally, the development of hierarchically ordered nanostructures in conjugated polymer binder can further enhance their electrochemical performance. This review examines the design principles of conjugated polymer binders, focusing on molecular engineering and nanostructural control to optimize their performance in high-loading electrodes, such as silicon-based anodes. By addressing key challenges in binder functionality, these advanced materials pave the way for next-generation high-energy-density LIBs.

Jin, Xiuyu↗
DOE OSTI2023

Density functional modeling of the binding energies between aluminosilicate oligomers and different metal cations

Interactions between negatively charged aluminosilicate species and positively charged metal cations are critical to many important engineering processes and applications, including sustainable cements and aluminosilicate glasses. In an effort to probe these interactions, here we have calculated the pair-wise interaction energies (i.e., binding energies) between aluminosilicate dimer/trimer and 17 different metal cations M n+ (M n+ = Li + , Na + , K + , Cu + , Cu 2+ , Co 2+ , Zn 2+ , Ni 2+ , Mg 2+ , Ca 2+ , Ti 2+ , Fe 2+ , Fe 3+ , Co 3+ , Cr 3+ , Ti 4+ and Cr 6+ ) using a density functional theory (DFT) approach. Analysis of the DFT-optimized structural representations for the clusters (dimer/trimer + M n+ ) shows that their structural attributes (e.g., interatomic distances) are generally consistent with literature observations on aluminosilicate glasses. The DFT-derived binding energies are seen to vary considerably depending on the type of cations (i.e., charge and ionic radii) and aluminosilicate species (i.e., dimer or trimer). A survey of the literature reveals that the difference in the calculated binding energies between different M n+ can be used to explain many literature observations associated with the impact of metal cations on materials properties (e.g., glass corrosion, mineral dissolution, and ionic transport). Analysis of all the DFT-derived binding energies reveals that the correlation between these energy values and the ionic potential and field strength of the metal cations are well captured by 2nd order polynomial functions ( R 2 values of 0.99–1.00 are achieved for regressions). Given that the ionic potential and field strength of a given metal cation can be readily estimated using well-tabulated ionic radii available in the literature, these simple polynomial functions would enable rapid estimation of the binding energies of a much wider range of cations with the aluminosilicate dimer/trimer, providing guidance on the design and optimization of sustainable cements and aluminosilicate glasses and their associated applications. Finally, the limitations associated with using these simple model systems to model complex interactions are also discussed.

Gong, Kai↗
DOE OSTI2021

Nitromethane Decomposition via Automated Reaction Discovery and an Ab Initio Corrected Kinetic Model

In the explore the systematic construction of kinetic models from in silico reaction data for the decomposition of nitromethane. Our models are constructed in a computationally affordable manner by using reactions discovered through accelerated molecular dynamics simulations using the ReaxFF reactive force field. The reaction paths are then optimized to determine reaction rate parameters. We introduce a reaction barrier correction scheme that combines accurate thermochemical data from density functional theory with ReaxFF minimal energy paths. We validate our models across different thermodynamic regimes, showing predictions of gas phase CO and NO concentrations and high-pressure induction times that are similar to experimental data. The kinetic models are analyzed to find fundamental decomposition reactions in different thermodynamic regimes.

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

Binding of Sulfates and Water to Monovalent Cations

The binding of the sulfate ligand group to monovalent cations in the presence of water is important for many systems. To understand the structure and energetics of sulfate complexes, we use density functional theory to study ethyl sulfate binding to the monovalent cations Li + , Na + , and K + , and to water. The free energies of binding and optimal structures are calculated for a range of the number of ethyl sulfates and waters. Without water, the most optimal structure for all the cations is bidentate binding by two ethyl sulfates, yielding a 4-fold coordination. With water, the lowest free energy structures also have two ethyl sulfates, but the coordination varies with cations. For complexes with water, the four oxygen atoms in the sulfate group enable multiple binding geometries for the cations and for hydrogen bonding with water. Many of these geometries differ in free energy by only a small amount (1–2 kcal/mol), meaning there will be multiple binding configurations in bulk solution. In comparison to the optimal structures for binding to the carboxylate group, there is more variation for binding to the sulfate group as a function of cation type and the number of waters. Further, the polarization of the atoms is significant and varies among the sulfate oxygen atoms. The water oxygen charge is often larger than that of sulfate oxygen, which plays a role in the preference for monodentate ligand binding to cations in the presence of water.

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

Simulating excited-state complex ensembles: Fluorescence and solvatochromism in amine-arene exciplexes

Exciplexes are excited-state complexes formed as a result of partial charge transfer from the donor to the acceptor species when one moiety of the donor–acceptor pair is electronically excited. The arene–amine exciplex formed between oligo-(p-phenylene) (OPP) and triethylamine (TEA) is of interest in the catalytic photoreduction of CO2 because it can compete with complete electron transfer to the OPP catalyst. Therefore, formation of the exciplex can hinder the generation of a radical anion OPP·− necessary for subsequent CO2 reduction. We report an implementation of a workflow automating quantum-chemistry calculations that generate and characterize an ensemble of structures to represent this exciplex state. We use FireWorks, Pymatgen, and Custodian Python packages for high-throughput ensemble generation. The workflow includes time-dependent density functional theory optimization, verification of excited-state minima, and exciplex characterization with natural transition orbitals, exciton analysis, excited-state Mulliken charges, and energy decomposition analysis. Fluorescence spectra computed for these ensembles using Boltzmann-weighted contributions of each structure agree better with experiment than our previous calculations based on a single representative exciplex structure [Kron et al., J. Phys. Chem. A 126, 2319–2329 (2022)]. The ensemble description of the exciplex state also reproduces an experimentally observed red shift of the emission spectrum of [OPP-4–TEA]* relative to [OPP-3–TEA]*. The workflow developed here streamlines otherwise labor-intensive calculations that would require significant user involvement and intervention.

Chemistry↗
DOE OSTI2021

Analytical gradients for molecular-orbital-based machine learning

We report molecular-orbital-based machine learning (MOB-ML) enables the prediction of accurate correlation energies at the cost of obtaining molecular orbitals. Here, we present the derivation, implementation, and numerical demonstration of MOB-ML analytical nuclear gradients, which are formulated in a general Lagrangian framework to enforce orthogonality, localization, and Brillouin constraints on the molecular orbitals. The MOB-ML gradient framework is general with respect to the regression technique (e.g., Gaussian process regression or neural networks) and the MOB feature design. We show that MOB-ML gradients are highly accurate compared to other ML methods on the ISO17 dataset while only being trained on energies for hundreds of molecules compared to energies and gradients for hundreds of thousands of molecules for the other ML methods. The MOB-ML gradients are also shown to yield accurate optimized structures at a computational cost for the gradient evaluation that is comparable to a density-corrected density functional theory calculation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗