Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Findings”

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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Assessing uncertainty in the timing of energy use during cost-optimal distributed energy technology selection and sizing

This paper empirically derives uncertainty ranges in cost-optimal solar PV and storage sizing by comparing results from the REopt Lite optimization platform from metered data and a set of simulated Department of Energy Commercial Reference Building (CRB) profiles at 65 sites. We find load profile shape alone does not explain a site’s optimal configurations (i.e., PV, Storage, PV and Storage, No System). Still, load profile shape does introduce uncertainty to optimal PV and storage capacities. Across all cases where PV is part of an optimal configuration, we find the average ratio of power capacities derived from metered loads to capacities derived from CRB profiles to be 0.97 (and as high as 1463), where 1 would be a perfect match in system size. For storage, the ratio is 1.6 (and as high as 42). We also assess how, in the absence of complete metered data, a CRB profile can be selected that would be expected to yield the most similar solar PV and storage capacities. From those metrics that can be available from billing data (i.e., peak demand, monthly load totals), we find that uncertainty is most reduced by selecting the CRB’s with an annual peak occurring at the most similar time, or those with the lowest average root mean square error (RMSE) among monthly peak loads. This research can help improve the implementation and interpretation of results derived from simulated load profiles and is an important next step in advancing smart grid solutions.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Structural and thermal properties and the origin of the ultralow thermal conductivity in the defect stannite CuIn 2 Se 4

Phase-pure CuIn 2 Se 4 , a ternary metal chalcogenide that forms in a disordered stannite crystal structure, was synthesized to investigate the structure thermal property relationships as well as reveal the origin of the ultralow thermal conductivity this material possesses over a large temperature range. Modeling of the temperature-dependent heat capacity and thermal conductivity revealed distinctive thermal properties and large lattice anharmonicity. Electron localization function calculations highlight the asymmetric bonding inherent to CuIn 2 Se 4 , which together with lattice anharmonicity directly impacts the thermal properties. Our findings reveal the specific atomic arrangement and bonding governing the thermal properties of this ternary metal chalcogenide. Our findings underscore the specific atomic arrangement and bonding governing the thermal properties in this ternary chalcogenide. This study advances the fundamental understanding of stannites, and our findings can be applied to these and other multinary metal chalcogenides of interest for applications where low thermal conductivity is desirable.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Human in vitro metabolism of an environmental mixture of polycyclic aromatic hydrocarbons (PAH) found at the Portland Harbor Superfund Site

Polycyclic aromatic hydrocarbons (PAHs) are widespread environmental contaminants that pose health risks to humans. Toxicity testing approaches of PAHs have evolved from traditional rodent models to New Approach Methodologies (NAMs), such as high-throughput screening in zebrafish, enabling rapid evaluation of chemical hazards. However, translating toxicity findings from laboratory systems to humans remains difficult due to complexity and species-specific differences. Chemical dosimetry modeling offers a quantitative framework to bridge this gap, but its accuracy depends on robust knowledge of PAH metabolism. The objective of this study was to measure human metabolism rates of Supermix-10, the ten most abundant PAHs found at the Portland Harbor Superfund Site, to support development of human pharmacokinetic models. We incubated individual PAHs from Supermix-10 in pooled human liver microsomes and quantified parent PAH disappearance using high-performance liquid chromatography (HPLC) with UV and florescent detection. To assess the potential of mixture interactions, we also measured metabolism of all 10 compounds in an equimolar mixture and compared rates of parent disappearance to those observed for individual PAHs. All Supermix-10 PAHs demonstrated rapid parent compound disappearance in human hepatic microsomes. PAHs grouped into three metabolism patterns: high metabolism rates and capacity (2-methylnaphthalene, acenaphthylene, fluorene, naphthalene), high affinity metabolism that rapidly achieves low-level saturation (benzo[a]anthracene, chrysene), and moderate metabolism rates and capacity (fluoranthene, pyrene, retene, phenanthrene). Smaller PAHs exhibited faster metabolism, and higher metabolism rates correlated inversely with molecular weight. When incubated in an equimolar mixture, Supermix-10 demonstrated significantly slower metabolism (47–89 %) compared to metabolism of individual PAHs at the same concentration. These findings enhance our understanding of PAH metabolism in humans and demonstrate significant mixture interactions under the conditions tested. Furthermore, our findings offer insights into the metabolic behavior of Supermix-10 and provide critical metabolism rate data to support the development of physiological based pharmacokinetic (PBPK) models. Dosimetry models can translate PAH chemical dosimetry from high-throughput testing platforms, like zebrafish and cellular system assays, to human exposures enhancing the accuracy and reliability of PAH risk assessments.

2-methylnaphthalene↗
DOE OSTI2020

Microsimulation of Energy and Flow Effects from Optimal Automated Driving in Mixed Traffic

In this paper we study the energy and traffic impact of a proposed Anticipative Cruise Controller in a PTV VISSIM microsimulation environment. We dissect our controller into two parts: 1. the unconnected mode, active when following a human-driven vehicle, and 2. the connected mode, active when following another automated vehicle equipped with connectivity. Probabilistic constraints balance safety considerations with inter-vehicle compactness, and vehicle constraints for acceleration capabilities are expressed through the use of powertrain maps. Emergent highway traffic scenarios are then modeled using time headway distributions from empirical traffic data. To study the impact of automation over a range of demands of free-flow to stop-and-go, we vary vehicle flux from low to high and vary automated vehicle penetration from low to high. When examining all-human driving scenarios, network capacity failed to meet demand in high-volume scenarios, such as rush-hour traffic. We further find that with connected automated vehicles introduced, network capacity was improved to support the high-volume scenarios. Finally, we examine energy efficiencies of the fleet for conventional, electric, and hybrid vehicles. We find that automated vehicles perform at a 10%-20% higher energy efficiency over human drivers when considering conventional powertrains, and find that automated vehicles perform at a 3%-9% higher energy efficiency over human drivers when considering electric and hybrid powertrains. Due to secondary effects of smoothing traffic flow and reducing unnecessary braking, energy benefits also apply to human-driven vehicles that interact with automated ones. Such simulated humans were found to drive up to 10% more energy-efficiently than they did in the baseline all-human scenario.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Funding a Just Transition Away from Coal in the U.S. Considering Avoided Damage from Air Pollution

Abstract Coal is declining in the U.S. as part of the clean energy transition, resulting in remarkable air pollution benefits for the American public and significant costs for the industry. Using the AP3 integrated assessment model, we estimate that fewer emissions of sulfur dioxide, nitrogen oxides, and primary fine particulate matter driven by coal’s decline led to $300 billion in benefits from 2014 to 2019. Conversely, we find that job losses driven by less coal plant and mining activity resulted in $7.84 billion in foregone wages over the same timeframe. While the benefits were greatly distributed (mostly throughout the East), costs were highly concentrated in coal communities. Transferring a small fraction of the benefits to workers could cover these costs while maintaining societal net benefits. Forecasting coal fleet damages from 2020 to 2035, we find that buying out or replacing these plants would result in $589 billion in air quality benefits, which considerably outweigh the costs. The return on investment increases when policy targets the most damaging capacity, and net benefits are maximized when removing just facilities where marginal benefits exceed marginal costs. Evaluating competitive reverse auction policy designs akin to Germany’s Coal Exit Act, we find that adjusting bids based on monetary damages rather than based only on carbon dioxide emissions – the German design – provides a welfare advantage. Our benefit–cost analyses clearly support policies that drive a swift and just transition away from coal, thereby clearing the air while supporting communities needing assistance.

Dennin, Luke R. (ORCID:0000000205405520)↗
DOE OSTI2020

Turbulent exchanges between near-inertial waves and balanced flows

Wind generated near-inertial waves are ubiquitous in the upper ocean. An improved understanding of near-inertial wave dynamics following their excitation in the ocean and their subsequent interaction with mesoscale geostrophic balanced flows is key to decoding oceanic energy flow pathways. In this regard, multiple oceanic data sets accumulated over the past few decades reveal that the relative strength of near-inertial waves and geostrophic balanced eddy fields is highly variable, both geographically and seasonally. Inspired by these observations, we investigate turbulent interactions and energy exchanges between near-inertial waves and balanced flows using freely evolving numerical simulations of the non-hydrostatic Boussinesq equations. Here, we find accelerated vertical propagation and dissipation of the waves in regimes where balanced and wave fields have comparable strengths. In such regimes we also find that near-inertial waves directly extract energy from balanced flows, with O(10%) being the amount of balanced energy loss. In contrast, we find that near-inertial waves transfer energy to balanced flows in regimes where balance-to-wave energy ratio is small, with the gain in balanced energy being dependent on the relative strength of waves. Furthermore, these regimes are characterized by relatively weaker vertical propagation and dissipation of the near-inertial wave field. One of the key outcomes of this study is the demonstration of the lack of a unique direction for near-inertial wave-balanced flow energy transfers. Depending on the balance-to-wave energy ratio, near-inertial waves can act as an energy sink or energy source for the geostrophic balanced flow.

58 GEOSCIENCES↗
DOE OSTI2024

Self-organization in collisionless, high- β turbulence

The magnetohydrodynamic (MHD) equations, as a collisional fluid model that remains in local thermodynamic equilibrium (LTE), have long been used to describe turbulence in myriad space and astrophysical plasmas. Yet, the vast majority of these plasmas, from the solar wind to the intracluster medium (ICM) of galaxy clusters, are only weakly collisional at best, meaning that significant deviations from LTE are not only possible but common. Recent studies have demonstrated that the kinetic physics inherent to this weakly collisional regime can fundamentally transform the evolution of such plasmas across a wide range of scales. Here, we explore the consequences of pressure anisotropy and Larmor-scale instabilities for collisionless, β $\gg$ 1, turbulence, focusing on the role of a self-organizational effect known as ‘magneto-immutability’. We describe this self-organization analytically through a high-β, reduced ordering of the Chew–Goldberger–Low-MHD (CGL-MHD) equations, finding that it is a robust inertial-range effect that dynamically suppresses magnetic-field-strength fluctuations, anisotropic-pressure stresses and dissipation due to heat fluxes. As a result, the turbulent cascade of Alfvénic fluctuations continues below the putative viscous scale to form a robust, nearly conservative, MHD-like inertial range. These findings are confirmed numerically via Landau-fluid CGL-MHD turbulence simulations that employ a collisional closure to mimic the effects of microinstabilities. We find that microinstabilities occupy a small (~5%) volume-filling fraction of the plasma, even when the pressure anisotropy is driven strongly towards its instability thresholds. We discuss these results in the context of recent predictions for ion-vs-electron heating in low-luminosity accretion flows and observations implying suppressed viscosity in ICM turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

On reduced modelling of the modulational dynamics in magnetohydrodynamics

This paper explores structure formation in two-dimensional magnetohydrodynamic (MHD) turbulence as a modulational instability (MI) of turbulent fluctuations. We focus on the early stages of structure formation and consider simple backgrounds that allow for a tractable model of the MI while retaining the full chain of modulational harmonics. This approach allows us to systematically examine the validity of popular closures such as the quasilinear approximation and other low-order truncations. We find that, although such simple closures can provide quantitatively accurate approximations of the MI growth rates in some regimes, they can fail to capture the modulational dynamics in adjacent regimes even qualitatively, falsely predicting MI when the system is actually stable. We find that this discrepancy is due to the excitation of propagating spectral waves (PSWs) which can ballistically transport energy along the modulational spectrum, unimpeded until dissipative scales, thereby breaking the feedback loops that would otherwise sustain MIs. The PSWs can be self-maintained as global modes with real frequencies and drain energy from the primary structure at a constant rate until the primary structure is depleted. To describe these waves within a reduced model, we propose an approximate spectral closure that captures them and MIs on the same footing. We also find that introducing corrections to ideal MHD, conservative or dissipative, can suppress PSWs and reinstate the accuracy of the quasilinear approximation. In this sense, ideal MHD is a ‘singular’ system that is particularly sensitive to the accuracy of the closure within mean-field models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Design Principles for Efficient and Stable Water Splitting Photoelectrocatalysts

Photoelectrochemical water splitting is a promising avenue for sustainable production of hydrogen used in the chemical industry and hydrogen fuel cells. The basic components of most photoelectrochemical water splitting systems are semiconductor light absorbers coupled to electrocatalysts, which perform the desired chemical reactions. A critical challenge for the design of these systems is the lack of stability for the majority of desired semiconductors under operating water splitting conditions. One strategy to address this issue is to protect the semiconductor by covering it with a stabilizing insulator layer, creating a metal–insulator–semiconductor (MIS) architecture, which has demonstrated improved stability. In addition to enhanced stability, the insulator layer may significantly affect the electron and hole transfer, which governs the recombination rates. Furthermore, the insertion of an insulator layer leads to the introduction of additional insulator/electrocatalyst and insulator/semiconductor interfaces. These interfaces can impact the system’s performance significantly, and they need to be carefully engineered to optimize the efficiencies of MIS systems. In this Account, we describe our recent progress in shedding light on the critical role of the insulator and the interfaces on the performance of MIS systems. We discuss our findings by focusing on the concrete example of planar n-type Si protected by a HfO 2 insulator layer and coupled to a Ni or Ir electrocatalyst that performs the oxygen evolution reaction, one of the water splitting half-reactions. To improve our fundamental understanding of the insulator layer, we precisely control the HfO 2 insulator thickness using atomic layer deposition (ALD), and we perform a series of rigorous electrochemical experiments coupled with theory and modeling. We demonstrate that by tuning the insulator thickness, we can control the flux and recombination of photogenerated electrons and holes to optimize the generated photovoltage. Despite optimizing the thickness, we find that the maximum generated photovoltage in MIS systems is often significantly lower than the upper performance limit, i.e., there are additional losses in the system that could not be addressed by optimizing the insulator thickness. We identify the sources of these losses and describe strategies to minimize them by a combination of improving the semiconductor light absorption, removing nonidealities associated with interfacial defects, and finding alternative insulators with improved charge carrier selectivity. Finally, we quantify the improvements that can be obtained by implementing these specific strategies. Our collective work outlines strategies to analyze MIS systems, identify the sources of efficiency losses, and optimize the design to approach the fundamental performance limits. Furthermore, these general approaches are broadly applicable to photoelectrochemical materials that utilize sunlight to produce value-added chemicals.

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

Site Disorder Drives Cyanide Dynamics and Fast Ion Transport in Li 6 PS 5 CN

Halide argyrodite solid-state electrolytes of the general formula Li 6 PS 5 X exhibit complex static and dynamic disorder that plays a crucial role in ion transport processes. Here, we unravel the rich interplay between site disorder and dynamics in the plastic crystal argyrodite Li 6 PS 5 CN and the impact on ion diffusion processes through a suite of experimental and computational methodologies, including temperature-dependent synchrotron powder X-ray diffraction, AC electrochemical impedance spectroscopy, 7 Li solid-state NMR, and machine learning-assisted molecular dynamics simulations. Sulfide and (pseudo)halide site disorder between the two anion sublattices unilaterally improves long-range lithium diffusion irrespective of the (pseudo)halide identity, which demonstrates the importance of site disorder in dictating bulk ionic conductivity in the argyrodite family. Furthermore, we find that anion site disorder modulates the presence and time scales of cyanide rotational dynamics. Ordered configurations of anions enable fast, quasi-free rotations of cyanides that occur on time scales of 10 11 Hz at T = 300 K. In contrast, we find that cyanide dynamics are slow or frozen in Li 6 PS 5 CN when site disorder between the cyanide and sulfide sublattices is present at T = 300 K. We rationalize the observed differences in cyanide dynamics in the context of elastic dipole interactions between neighboring cyanide anions and local strain induced by the configurations of site disorder that may impact the energetic landscape for cyanide rotational dynamics. Through this study, we find that anion disorder plays a decisive role in dictating the extent and time scales of both lithium ion and cyanide dynamics in Li 6 PS 5 CN.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Advancing Porous Carbons: Understanding the Importance of Surface Chemistry for the Energy–Environment Nexus

This review intends, in a critical way, the comprehensive view of the importance of porous carbons surface chemistry for their applications in an energy− environment nexus. Surface chemistry is presented as a combination of functional heteroatom-containing groups, dopants, and structural defects. First, we briefly address carbon surface chemical environment and the methods of its modification and characterization, indicating their practical limitations. Then, the effects of surface chemistry on separation, catalysis, energy storage, sensing and microwave absorption are introduced. Besides a critical analysis of published findings on these topics, we also include our views on the advancement in the processes which rely on porous carbons surface chemistry, and identify strategic areas and directions that should deserve further attention. We focus on new findings and important original contributions to the field. Since the community of carbon researchers grows following the strategic application of these materials, the role of functional groups, dopants and structural defects in various cutting-edge applications is emphasized, showing the progress in the field and the evolution of findings. A clear determination of the effects of carbon surface is often a challenge since carbons porosity and the locations of specific bonds/sites/ defects in the carbon texture provide nanoconfinement effects.

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

Temporal Variability of Emissions Revealed by Continuous, Long-Term Monitoring of an Underground Natural Gas Storage Facility

Temporal variability contributes to uncertainty in inventories of methane emissions from the natural gas supply chain. Extrapolation of instantaneous, “snapshot-in-time” measurements, for example, can miss temporal intermittency and confound bottom-up/top-down comparisons. Importantly, no continuous long-term datasets record emission variability from underground natural gas storage facilities despite substantial contributions to sector-wide emissions. We present 11 months of continuous observations on a section of a storage site using dual-frequency comb spectroscopy (DCS observing system) and aircraft measurements. Here, we find high emission variability and a skewed distribution in which the 10% highest 3 h emission periods observed by the continuous DCS observing system comprise 41% of the total observed 3-hourly emissions. Monthly emission rates differ by >12×, and 3-hourly rates vary by 17× in 24 h. We find links to the operating phase of the facility—emission rates, including as a percentage of the total gas flow rate, are significantly higher during periods of injection compared to those of withdrawal. We find that if a high frequency of aircraft flights can occur, then the ground- and aircraft-based approaches show excellent agreement in emission distributions. A better understanding of emission variability at underground natural gas storage sites will improve inventories and models of methane emissions and clarify pathways toward mitigation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Simulating Droplet-Resolved Haze and Cloud Chemistry Forming Secondary Organic Aerosols in Turbulent Conditions within Laboratory and Cloud Parcels

Most of our existing knowledge of cloud chemistry in regards to forming secondary organic aerosols (SOA) is based on measurements in bulk aqueous solutions. However, SOA reaction kinetics derived from bulk solution measurements might differ from the kinetics in actual cloud droplets, since turbulent mixing and ionic strengths, and ratio of surface area to volume in individual cloud droplets might vary substantially from bulk solutions in the real atmosphere. Three-dimensional models at various scales have been used to simulate aqueous chemistry. However, most of these models do not resolve turbulence down to the smallest length scales of 1 mm and do not simulate cloud chemistry in individual cloud droplets due to large computational costs. Here we incorporate the formation of isoprene epoxydiol SOA (IEPOX-SOA) in individual droplets within a one-dimensional explicit mixing parcel model (EMPM-Chem). We apply EMPM-Chem to simulate turbulence and droplet-resolved IEPOX-SOA formation using a configuration based on the Michigan Tech Pi chamber. We find that the dissolution of IEPOX gases is weighted more towards larger cloud droplets due to their large liquid water content (compared to smaller droplets), while the conversion of dissolved IEPOX to IEPOX-SOA is much greater within smaller deliquesced haze particles due to their higher acidity and ionic strengths compared to cloud droplets. We also find that as droplet residence times increase in the chamber, e.g., due to increasing ammonium bisulfate seed aerosol injection rates and/or increasing heights of the chamber, formation of IEPOX-SOA increases substantially. Thus, our EMPM-Chem model could be used to design future cloud chambers to maximize SOA production from cloud chemistry. We also apply the EMPM-Chem model to simulate how IEPOX-SOA formation evolves in individual cloud droplets within rising cloudy parcels in the atmosphere. We find that as subsaturated air is entrained into and turbulently mixed with the cloud parcel, evaporation causes a reduction in droplet sizes, which leads to corresponding increases in per droplet ionic strength and acidity. Increased droplet acidity in turn greatly accelerates the kinetics of IEPOX-SOA formation. In conclusion, our results provide key insights into single-cloud-droplet chemistry, suggesting that entrainment mixing may be an important process that increases SOA formation in the real atmosphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Exploring the Limits of Dative Boratrane Bonding: Iron as a Strong Lewis Base in Low-Valent Non-Heme Iron-Nitrosyl Complexes

We previously reported the synthesis and preliminary characterization of a unique series of low-spin (ls) {FeNO} 8–10 complexes supported by an ambiphilic trisphosphineborane ligand, [Fe(TPB)(NO)] +/0/– . In this study, we use advanced spectroscopic techniques and density functional theory (DFT) calculations to extract detailed information as to how the bonding changes across the redox series. We find that, in spite of the highly reduced nature of these complexes, they feature an NO + ligand throughout with strong Fe–NO π-backbonding and essentially closed-shell electronic structures of their FeNO units. This is enabled by an Fe–B interaction that is present throughout the series. In particular, the most reduced [Fe(TPB)(NO)] – complex, an example of a ls-{FeNO} 10 species, features a true reverse dative Fe → B bond where the Fe center acts as a strong Lewis-base. Hence, this complex is in fact electronically similar to the ls-{FeNO} 8 system, with two additional electrons “stored” on site in an Fe–B single bond. The outlier in this series is the ls-{FeNO} 9 complex, due to spin polarization (quantified by pulse EPR spectroscopy), which weakens the Fe–NO bond. These data are further contextualized by comparison with a related N 2 complex, [Fe(TPB)(N 2 )] – , which is a key intermediate in Fe(TPB)-catalyzed N 2 fixation. Our present study finds that the Fe → B interaction is key for storing the electrons needed to achieve a highly reduced state in these systems, and highlights the pitfalls associated with using geometric parameters to try to evaluate reverse dative interactions, a finding with broader implications to the study of transition metal complexes with boratrane and related ligands.

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

Synthesis of Solid-Solution Mn x Zn 1– x O Nanoparticles and their Electrochemical Oxidation of Furfural

Electrochemical valorization of biomass-derived substrates has become a prominent area of research due to its potential to produce value-added products from renewable feedstocks in a more sustainable way. First-row transition metal electrodes are compelling candidates for these conversions due to their stability, abundance, and cost-effectiveness. Herein, we report on the colloidal synthesis of Mn x Zn 1– x O (x = 0.3-0.7) nanoparticles and their electrocatalytic activity towards furfural oxidation. We find that the hindrance of a MnO impurity can be achieved by leveraging the oxidation state of the Mn precursor. The Mn x Zn 1– x O composition closely follows the ratio of precursors, with all the nanoparticles having a wurtzite structure as determined by ICP-MS and PXRD, respectively. XANES and XPS revealed the presence of Mn in different oxidation states with the ratio of these varying based on the composition. When comparing the electrocatalytic activity of the monometallic and bimetallic oxides for furfural oxidation, a decrease in current density was observed with increasing Zn content. We find the Mn x Zn 1– x O nanoparticles favor the formation of the 6 e - oxidation product 5-hydroxy-2(5H)-furanone, while both monometallic oxides primarily yield CO 2 and other deeply oxidized products as the majority pathway. Furthermore, these findings can contribute towards the design and synthesis of more active and selective electrocatalyst.

Aldehydes↗
DOE OSTI2020

Spin Splitting Energy of Transition Metals: A New, More Affordable Wave Function Benchmark Method and Its Use to Test Density Functional Theory

Accurately predicting the spin splitting energy of chemical species is important for understanding their reactivity and magnetic properties, but it is very challenging, especially for molecules containing transition metals. One impediment to progress is the scarcity of accurate benchmark data. Here we report a set of calculations designed to yield reliable benchmarks for simple transition-metal complexes that can be used to test density functional methods that are affordable for large systems of more practical interest. Various wave function methods are tested against experiment for Fe 2+ , Fe 3+ , and Co 3+ , including CASSCF, CASPT2, CASPT3, MRCISD, MRCISD+Q, ACPF, AQCC, CCSD(T), and CASPT2/CCSD(T) and also a new method called CASPT2.5, which is performed by taking the average of the CASPT2 and CASPT3 energies. Furthermore, we find that MRCISD+Q, ACPF, and AQCC require smaller active spaces for good accuracy than are required by CASPT2 and CASPT3, and this aspect may be important for calculations on larger molecules; here we find that CASPT2.5 extrapolated to a complete basis set is the most suitable method—in terms of computational cost and in terms of accuracy on monatomic systems—and therefore we chose this method for molecular benchmarks. Then Kohn–Sham density functional calculations with 60 exchange-correlation functionals are tested for FeF 2 , FeCl 2 , and CoF 2 . We find that MN15-L, M06-SX, and revM06 have very good agreement with CASPT2.5 benchmarks in terms of both the spin splitting energy and the optimized geometry for each spin state. In addition, we recommend def2-TZVP as the most suitable basis set to perform density functional calculations for molecular spin splitting energies; extra polarization functions in the basis set do not help to increase the accuracy of the spin splitting energy in KS calculations.

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

Large Scale Study of Ligand–Protein Relative Binding Free Energy Calculations: Actionable Predictions from Statistically Robust Protocols

The accurate and reliable prediction of protein–ligand binding affinities can play a central role in the drug discovery process as well as in personalized medicine. Of considerable importance during lead optimization are the alchemical free energy methods that furnish an estimation of relative binding free energies (RBFE) of similar molecules. Recent advances in these methods have increased their speed, accuracy, and precision. This is evident from the increasing number of retrospective as well as prospective studies employing them. However, such methods still have limited applicability in real-world scenarios due to a number of important yet unresolved issues. Here, we report the findings from a large data set comprising over 500 ligand transformations spanning over 300 ligands binding to a diverse set of 14 different protein targets which furnish statistically robust results on the accuracy, precision, and reproducibility of RBFE calculations. We use ensemble-based methods which are the only way to provide reliable uncertainty quantification given that the underlying molecular dynamics is chaotic. These are implemented using TIES (Thermodynamic Integration with Enhanced Sampling). Results achieve chemical accuracy in all cases. Ensemble simulations also furnish information on the statistical distributions of the free energy calculations which exhibit non-normal behavior. We find that the “enhanced sampling” method known as replica exchange with solute tempering degrades RBFE predictions. We also report definitively on numerous associated alchemical factors including the choice of ligand charge method, flexibility in ligand structure, and the size of the alchemical region including the number of atoms involved in transforming one ligand into another. Our findings provide a key set of recommendations that should be adopted for the reliable application of RBFE methods.

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

Toward a Machine Learning Approach to Interpreting X-ray Spectra of Trace Impurities by Converting XANES to EXAFS

The fact that the photoabsorption spectrum of a material contains information about the atomic structure, commonly understood in terms of multiple scattering theory, is the basis of the popular extended X-ray absorption spectroscopy (EXAFS) technique. How much of the same structural information is present in other complementary spectroscopic signals is not obvious. Here we use a machine learning approach to demonstrate that within theoretical models that accurately predict the EXAFS signal, the extended near-edge region does indeed contain the EXAFS-accessible structural information. We do this by exhibiting deep operator neural networks (DeepONets) that have learned the relationship between the extended and near edge portions of the X-ray absorption spectrum to predict the former from the latter. We find that we can accurately predict the EXAFS spectrum between 6 and 14 Å –1 from the first 6 Å –1 (≈100 eV) of the absorption spectrum of Cu 2 + substitutional defects in the Fe 3+ mineral hematite (α-Fe 2 O 3 ). This surprising finding implies that theoretical analyses of X-ray absorption spectra could be implemented that extract the same conclusions as high-quality EXAFS studies from spectra collected over a much smaller range of photon energies. This relaxes a host of experimental limitations related to the X-ray source and measurement sample, including collection time, minimum dopant concentration, source brilliance, and energy range. We describe the theoretical data sets and DeepONet construction and show that the resulting DeepONets produce EXAFS that recovers linear combination fits to experimental data with accuracy approaching the original ab initio calculations. We discuss the implications of our findings for minor constituent characterization and for understanding the information content of spectroscopic data more broadly, including how this approach might be applied to measured experimental spectra. In conclusion, to encourage similar efforts, the simulated X-ray spectra, machine learning, and fitting code are publicly available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗