Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cost function”

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

Identifying Band Inversions in Topological Materials Using Diffusion Monte Carlo

Topological insulators are characterized by insulating bulk states and robust metallic surface states. Band inversion is a hallmark of topological insulators. At time-reversal invariant points in the Brillouin zone, spin–orbit coupling (SOC) induces a swapping of orbital character at the bulk band edges. Reliably detecting band inversion in solid-state systems with many-body methods would aid in identifying possible candidates for spintronics and quantum computing applications and improve our understanding of the physics behind topologically nontrivial systems. Density functional theory (DFT) methods are a well-established means of investigating these interesting materials due to their favorable balance of computational cost and accuracy but often struggle to accurately model the electron–electron correlations present in the many materials containing heavier elements. In this work, we develop a novel method to detect band inversion within continuum quantum Monte Carlo (QMC) methods that can accurately treat the electron correlation and spin–orbit coupling that are crucial to the physics of topological insulators. Our approach applies a momentum-space-resolved atomic population analysis throughout the first Brillouin zone utilizing the Löwdin method and the one-body reduced density matrix produced with diffusion Monte Carlo (DMC). We integrate this method into QMCPACK, an open source ab initio QMC package, so that these ground-state methods can be used to complement experimental studies and validate prior DFT work on predicting the band structures of correlated topological insulators. Here, we demonstrate this new technique on the topological insulator bismuth telluride, which displays band inversion between its Bi-p and Te-p states at the Γ-point. We show an increase in charge on the bismuth-p orbital and a decrease in charge on the tellurium-p orbital when comparing band structures with and without SOC. Additionally, we use our method to compare the degree of band inversion present in monolayer Bi 2 Te 3 , which has no interlayer van der Waals interactions, to that seen in the bilayer and bulk. The method presented here will enable future many-body studies of band inversion that can shed light on the delicate interplay between correlation and topology in correlated topological materials.

Band structure↗
DOE OSTI2026

Models and Algorithms for Equilibrium Analysis of Mixed-Material Nucleic Acid Systems

Dynamic programming algorithms within the NUPACK software suite enable analysis of equilibrium base-pairing properties for complex and test tube ensembles containing arbitrary numbers of interacting nucleic acid strands. Currently, calculations are limited to single-material systems that are either all-RNA or all-DNA. Here, to enable analysis of mixed-material systems that are critical for modern applications in vitro, in situ, and in vivo, we develop physical models and dynamic programming algorithms that allow the material of the system to be specified at nucleotide resolution. Free energy parameter sets are constructed for both RNA/DNA and RNA/2'OMe-RNA mixed-material systems by combining available empirical mixed-material parameters with single-material parameter sets to enable treatment of the full complex and test tube ensembles. New dynamic programming recursions account for the material of each nucleotide throughout the recursive process. For a complex with N nucleotides, the mixed-material dynamic programming algorithms maintain the O(N 3 ) time complexity of the single-material algorithms, enabling efficient calculation of diverse physical quantities over complex and test tube ensembles (e.g., complex partition function, equilibrium complex concentrations, equilibrium base-pairing probabilities, minimum free energy secondary structure(s), and Boltzmann-sampled secondary structures) at a cost increase of roughly 2.0-3.5×. The results of existing single-material algorithms are exactly reproduced when applying the new mixed-material algorithms to single-material systems. Accuracy is significantly enhanced using mixed-material models and algorithms to predict RNA/DNA and RNA/2'OMe-RNA duplex melting temperatures from the experimental literature as well as RNA/DNA melt profiles from new experiments. In conclusion, mixed-material analyses can be performed online using the NUPACK web app (www.nupack.org) or locally using the NUPACK Python module.

2′OMe-RNA↗
DOE OSTI2026

Data for The Value of Reversible Carbon Storage in a Zero-Emissions World

Atmospheric carbon dioxide removal (CDR) is required to stabilize global temperature. CDR can be achieved via ecosystem-based approaches that are cost-effective but reversible (e.g., soil and forest management) or by more durable but expensive approaches (e.g., direct air capture coupled with geologic storage). Here, we examine trade-offs between these approaches, focusing on timing, climate impacts, and cost. We simulated reversible carbon accrual for a range of CDR contract structures using a general minimalist model of ecosystem carbon cycling, and parameterized it to simulate US agricultural soil management─specifically cover cropping─as a case study. We then quantified the resulting impact on atmospheric carbon and global temperature using a climate model emulator. We find that maintaining a patchwork of reversible CDR projects by replacing lapsed projects with new projects can reduce warming by 22–195 μ°C in 2100 and that the magnitude of this cooling effect depends on how effectively the patchwork is maintained. Long-term maintenance of reversible CDR projects requires institutional stability that cannot be guaranteed over multiple decades. Consequently, effective CDR ultimately requires replacing reversible projects with durable projects. To address this problem, we modeled the cost of replacing reversible agricultural soil CDR with geologic CDR. We found that using reversible CDR as a bridge to durable CDR is potentially more cost-effective as a global cooling strategy (0.20–0.81 billion USD per μ°C avoided) than perpetual maintenance of reversible CDR (0.32–1.31 billion USD per μ°C avoided) or an immediate transition to durable CDR (1.37–2.19 billion USD per μ°C avoided). However, we emphasize that institutional commitments to maintain reversible CDR projects cannot be guaranteed. Reliance on reversible CDR as a bridge to durable CDR therefore carries an unknown amount of risk and will only function if efforts to maintain reversible CDR are robust.

Carbon↗
DOE OSTI2026

Optimization of Desalination Systems with Detailed Water Chemistry through Integration of Reaktoro in WaterTAP

Chemistry predictions are critical for an accurate estimation of performance and costs in desalination process models, which allows for the estimation of the value of new technologies and the viability of treating new water sources. Herein, we present how an implicit function formulation can be used to integrate the chemical modeling package, Reaktoro, into the techno-economic assessment and modeling platform, WaterTAP. This approach resolves the critical issues of integrating large-scale thermodynamic models and databases into equation-oriented process models while allowing more flexibility relative to previously presented surrogate-based methods. We describe how this integration into Pyomo and WaterTAP models is implemented and used through the open-source package Reaktoro-PSE . We first validate this integration approach by performing optimization on a previously presented desalination treatment train with softening and acid addition as the pretreatment steps. Then, to demonstrate the value of this approach, we extend the cost-optimization problem to include the simultaneous addition of lime and soda ash for softening, and HCl and H 2 SO 4 in the acidification steps. Finally, we were able to confirm the previously established results that were obtained by using surrogate models and demonstrate that the implicit function approach enables exploration of different feedwater compositions and a larger number of chemicals and their combinations.

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

NSFnets (Navier-Stokes flow nets): Physics-informed neural networks for the incompressible Navier-Stokes equations

In the last 50 years there has been a tremendous progress in solving numerically the Navier-Stokes equations using finite differences, finite elements, spectral, and even meshless methods. Yet, in many real cases, we still cannot incorporate seamlessly (multi-fidelity) data into existing algorithms, and for industrial-complexity applications the mesh generation is time consuming and still an art. Moreover, solving ill-posed problems (e.g., lacking boundary conditions) or inverse problems is often prohibitively expensive and requires different formulations and new computer codes. Here, we employ physics-informed neural networks (PINNs), encoding the governing equations directly into the deep neural network via automatic differentiation, to overcome some of the aforementioned limitations for simulating incompressible laminar and turbulent flows. We develop the Navier-Stokes flow nets (NSFnets) by considering two different mathematical formulations of the Navier-Stokes equations: the velocity-pressure (VP) formulation and the vorticity-velocity (VV) formulation. Since this is a new approach, we first select some standard benchmark problems to assess the accuracy, convergence rate, computational cost and flexibility of NSFnets; analytical solutions and direct numerical simulation (DNS) databases provide proper initial and boundary conditions for the NSFnet simulations. The spatial and temporal coordinates are the inputs of the NSFnets, while the instantaneous velocity and pressure fields are the outputs for the VP-NSFnet, and the instantaneous velocity and vorticity fields are the outputs for the VV-NSFnet. This is unsupervised learning and, hence, no labeled data are required beyond boundary and initial conditions and the fluid properties. The residuals of the VP or VV governing equations, together with the initial and boundary conditions, are embedded into the loss function of the NSFnets. No data is provided for the pressure to the VP-NSFnet, which is a hidden state and is obtained via the incompressibility constraint without extra computational cost. Unlike the traditional numerical methods, NSFnets inherit the properties of neural networks (NNs), hence the total error is composed of the approximation, the optimization, and the generalization errors. Here, we empirically attempt to quantify these errors by varying the sampling (“residual”) points, the iterative solvers, and the size of the NN architecture. For the laminar flow solutions, we show that both the VP and the VV formulations are comparable in accuracy but their best performance corresponds to different NN architectures. The initial convergence rate is fast but the error eventually saturates to a plateau due to the dominance of the optimization error. For the turbulent channel flow, we show that NSFnets can sustain turbulence at , but due to expensive training we only consider part of the channel domain and enforce velocity boundary conditions on the subdomain boundaries provided by the DNS data base. We also perform a systematic study on the weights used in the loss function for balancing the data and physics components, and investigate a new way of computing the weights dynamically to accelerate training and enhance accuracy. In the last part, we demonstrate how NSFnets should be used in practice, namely for ill-posed problems with incomplete or noisy boundary conditions as well as for inverse problems. We obtain reasonably accurate solutions for such cases as well without the need to change the NSFnets and at the same computational cost as in the forward well-posed problems. As a result, we also present a simple example of transfer learning that will aid in accelerating the training of NSFnets for different parameter settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Exploring the Functional Performance of a Commercial High-Temperature Photopolymer Resin for Vat Photopolymerized Injection Molding Tools

Additive manufacturing (AM) presents a promising approach for rapid tooling, offering significant reductions in production costs and lead times compared to that of conventional moldmaking processes. This work examines the fabrication of molds for injection molding through vat photopolymerization (VPP) AM, using a commercially available photopolymer resin intended for high temperature applications (High Temp V2 resin). The photorheological behavior of the resin is examined to evaluate its curing kinetics and mechanical development during polymerization, ensuring suitability for mold fabrication. The thermal and mechanical properties of the printed photopolymer resin are characterized, alongside assessments of dimensional accuracy and surface finish, which are critical factors influencing the performance and quality of injection molded parts. Additionally, specimens produced using VPP-fabricated molds are evaluated and benchmarked against those molded with traditional metal molds.

Jones, Haley W. [Savannah River National Laborator↗
DOE OSTI2022

Transformational Solid Oxide Fuel Cell (SOFC) Technology

This project was conducted under the Co-operative Agreement No. DE-FE0027584 with the US Department Energy to developed advanced Solid Oxide Fuel Cell (SOFC) Technologies. The overall objective of this project was to advance SOFC technology at the cell and stack level to enhance cell robustness and durability, increase performance, and reduce balance-of-plant (BOP) requirements. By reducing system complexity combined with the increases in power density and efficiency, the ultimate goal of the project was to increased reliability and to reduce capital and operating costs of installed systems. The project was focused on pathways that will reduce the cost of the SOFC cell and stack, including the following areas: Robust, redox tolerant cell technology Lower cost cell manufacturing through advances in cell design, which will reduce the amount of material, energy and time used in the fabrication of SOFCs High performance, low temperature electrolyte based on improvement of established materials Innovative SOFC stack architecture which truly integrates Balance of Plant functionality into the stack level design Thermal management of the fuel cell stack for increased durability and expanded window of operation Novel stack design amenable for use in sub-MW to multi-MW-scale power plants and having low replacement cost The incorporation of balance-of-plant (BOP) equipment into the stack platform increased the economic viability of smaller scale systems. The project objectives were met by a multi-prong approach, including new cell design complemented with modifications to existing cell technology, as well as a new stack design incorporating components typically included in the BOP, such as heat exchangers, oxidizer, fuel reformers, and recycle systems. The project culminated with demonstration of a stack test validating the viability of the cell and stack improvements. A cost model was also developed to estimate costs for the advanced stack technology at high volume manufacturing levels. The net outcome of the project is SOFC cell and stack technology with costs significantly below current DOE targets without compromising and, in some cases, improving on the performance and degradation rate demonstrated with the current state-of-the-art stack design. The results of this project advanced the reliability, robustness, and endurance of low-cost SOFC technology that ultimately are ready to be deployed in coal power systems with greater than 60 percent efficiency (based on higher heating value of fuel) and the capability for ≥97% CO 2 capture at a cost-of-electricity that is approximately 40 percent below presently available Integrated Gasification Combined Cycle systems.

03 NATURAL GAS↗
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 OSTI2025

Open Architecture for Cost Savings in Advanced Nuclear Reactors

Recently, nuclear power plant build projects in the West have run over budget due to high capital costs and schedule overruns. Compared to other sources of energy, nuclear power plants have higher capital costs. Reactors are often different at every site, resulting in a lack of standardization. Nuclear is expected to compete with other low carbon sources of energy which have lower capital costs making it essential for nuclear to develop ways of reducing costs. Strategies such as standardization, learning rates, modularization, and schedule reduction in advanced reactors can reduce nuclear costs by about 40%. Standardization as a way of cutting capital costs has been explored even in large nuclear power plants. Standardization of certain plant components can result in lower component and installation costs and higher learning from experience. Standardization can be achieved by adopting a criterion of key performance indicators and general design principles for a specific system or component such as the balance of plant. Modularization allows the construction of certain components of SMRs in a factory, which saves time, increases productivity, and encourages higher learning rates. Production learning decreases the time and the cost related to an activity. The potential for modularized components of advanced reactors to be manufactured in factories makes it conducive to achieving higher learning rates. Developing large-capacity nuclear programs through sequential builds cultivates a higher learning rate, which in effect may reduce schedule overruns. Open architecture has been identified as a way to drive standardization among advanced reactor designs and result in cost savings. Open architecture (OA) is defined as a design enabling a diverse supply chain by defining and publishing requirements of systems or equipment in functional and/or interface terms, utilizing technical standards in widespread use. Currently, the nuclear industry’s approach is to use closed architecture, making most designs proprietary. However, collaboration between various advanced reactor vendors and suppliers utilizing the concept of open architecture can result in modular and standardized architecture of subsystems or subcomponents of a nuclear power plant. Completely standardizing nuclear power plants may be impossible, however, certain common subsystems amongst the various reactor designs could be standardized and/or access a wider supply chain and leverage existing learning from other sectors. Open architecture will save time and allocate resources to the parts of the plants that have the most unique features. A key advantage of open architecture is its ability to improve production learning across advanced reactors (AR) types in the industry, by providing and utilizing the same kind of component. Sodium fast reactor (SFR), High Temperature Gas Reactor (HTGR) and Molten Salt Reactor (MSR) are the advanced reactors considered for this project. This paper aims to determine the cost savings in advanced reactor programs due to open architecture learning rate. This work is an extension of work done on light water reactor small modular reactors; the cost methodology was utilized to investigate the impact of open architecture on advanced reactors with a particular focus on sodium fast reactors. The cost data on sodium fast reactors used in the model presented the most adequate information required for the analysis.

Advanced Nuclear Reactors↗
DOE OSTI2022

Computing x-ray absorption spectra from linear-response particles atop optimized holes

State specific orbital optimized density functional theory (OO-DFT) methods, such as restricted open-shell Kohn–Sham (ROKS), can attain semiquantitative accuracy for predicting x-ray absorption spectra of closed-shell molecules. OO-DFT methods, however, require that each state be individually optimized. In this Communication, we present an approach to generate an approximate core-excited state density for use with the ROKS energy ansatz, which is capable of giving reasonable accuracy without requiring state-specific optimization. Herein, this is achieved by fully optimizing the core-hole through the core-ionized state, followed by the use of electron-addition configuration interaction singles to obtain the particle level. This hybrid approach can be viewed as a DFT generalization of the static-exchange (STEX) method and can attain ~0.6 eV rms error for the K-edges of C–F through the use of local functionals, such as PBE and OLYP. This ROKS(STEX) approach can also be used to identify important transitions for full OO ROKS treatment and can thus help reduce the computational cost of obtaining OO-DFT quality spectra. ROKS(STEX), therefore, appears to be a useful technique for the efficient prediction of x-ray absorption spectra.

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

Analytical harmonic vibrational frequencies with VV10-containing density functionals: Theory, efficient implementation, and benchmark assessments

VV10 is a powerful nonlocal density functional for long-range correlation that is used to include dispersion effects in many modern density functionals, such as the meta-generalized gradient approximation (mGGA), B97M-V, the hybrid GGA, ωB97X-V, and the hybrid mGGA, ωB97M-V. While energies and analytical gradients for VV10 are already widely available, this study reports the first derivation and efficient implementation of the analytical second derivatives of the VV10 energy. The additional compute cost of the VV10 contributions to analytical frequencies is shown to be small in all but the smallest basis sets for recommended grid sizes. Here, this study also reports the assessment of VV10-containing functionals for predicting harmonic frequencies using the analytical second derivative code. The contribution of VV10 to simulating harmonic frequencies is shown to be small for small molecules but important for systems where weak interactions are important, such as water clusters. In the latter cases, B97M-V, ωB97M-V, and ωB97X-V perform very well. The convergence of frequencies with respect to the grid size and atomic orbital basis set size is studied, and recommendations are reported. Finally, scaling factors to allow comparison of scaled harmonic frequencies with experimental fundamental frequencies and to predict zero-point vibrational energy are presented for some recently developed functionals (including r2SCAN, B97M-V, ωB97X-V, M06-SX, and ωB97M-V).

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

Pushing Cu uphill of the volcano curve: Impact of a WC support on the catalytic activity of copper toward the hydrogen evolution reaction

Here, the adsorption of atomic H and H 2 on copper mono- and submonolayers supported on hexagonal WC(0001) surfaces has been investigated using density functional theory with the Perdew–Burke–Ernzerhof exchange correlation functional and D2 van der Waals corrections. Results evidence the impact of the termination of the carbide substrate on fundamental properties of Cu adatoms, and, hence, on the stability of molecular and atomic hydrogen, defining copper's catalytic activity for hydrogen evolution reaction. Using H adsorption energy as a descriptor, catalytic activity of Cu adlayers for hydrogen evolution reaction was estimated using traditional volcano curves and a curve, obtained at low hydrogen coverage. Obtained results evidence that copper adlayers supported on the WC may present a viable low-cost alternative to noble metal-based catalysts, with improved catalytic activity compared to that of copper. This, potentially, can be a useful basis for designing and developing novel functional materials with predetermined catalytic properties.

08 HYDROGEN↗
DOE OSTI2022

Muconic Acid Production from P. putida Using High Protein Algae Hydrolysate

The composition of algal biomass is highly dynamic, with protein, lipid, and carbohydrate contents varying in response to nutrient and environmental conditions during cultivation. Because shifts in biomass composition are often associated with reduced biomass productivity, production costs can often increase if targeting higher biomass compositional quality (enriched in carbohydrates or lipids at reduced protein content) as input for the algal biorefinery. The optimal algal biorefinery configuration is thus a function of many factors. One of the key strengths of the Combined Algal Processing (CAP) process is the versatility of feedstocks and products produced. The concept has been demonstrated with ethanol and a variety of carboxylic acids (succinic, butyric, muconic) as coproducts along with lipid upgrading to biofuel. Modification of the approaches, processes and downstream upgrading to fuels has allowed the CAP process to reduce costs and improve efficiency. Muconic acid is a high-value, potential fermentation coproduct of interest because it can be easily converted to adipic acid, a high-volume monomer for the production of nylon and other valuable consumer plastics. As such, the production of muconic acid through CAP was explored to expand the suite of products from algal biomass and to begin exploring the valorization of high protein content biomass from rapidly grown algae biomass. We have established initial performance parameters and shown that the range of substrates consumed by the muconic acid-producing microbe, Pseudomonas putida, includes at least glucose, mannose, glycerol, and lactic acid. We achieved complete utilization of these four major hydrolysate substrates achieving productivities of 0.037 (g/L/h) from Scenedesmus obliquus and 0.029 (g/L/h) from Monoraphidium minutum hydrolysates. Final titer and process yield (mol of muconic acid per mol of substrate (molP/molS)) were 0.99 g/L and 0.42 molP/molS from S. obliquus hydrolysate and 0.75 g/L and 0.23 molP/molS from M. minutum hydrolysate, respectively.

BIOMASS FUELS↗
DOE OSTI2023

Transportation and electricity systems integration via electric vehicle charging-as-a-service: A review of techno-economic and societal benefits

Here, in support of global decarbonization efforts, the adoption of electric vehicles is proceeding rapidly across transportation sectors, population groups, and regions. However, access to electric vehicle charging infrastructure remains sparse and inequitable, and technical, economic, and procedural challenges have encumbered its expansion. Charging-as-a-service, which can mitigate cost and effort risks of charging equipment ownership, adds a means of meeting many of these challenges. As charging-as-a-service has not been analyzed or addressed in the scholarly literature, we synthesize media, marketing, and analogous scholarly publications to categorize and describe the functions asserted by existing charging-as-a-service vendors. We demonstrate and quantify the benefits of many of these functions via four example business cases, sited at a detached residence, a multi-unit dwelling, a consumer-facing business, and a commercial fleet depot. Outcomes show that savings and revenues realized via charging-as-a-service can reduce charging costs, beneficially shape grid loads, and make charging installation economical where it otherwise would not be. Discussion highlights the need to deploy charging-as-a-service in cooperation with initiatives promoting societal benefits, including charging access equity.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

3D printing of metal-based materials for renewable energy applications

Large-scale renewable energy must overcome conversion and storage challenges before it can replace fossil fuels due to its intermittent nature. However, current sustainable energy devices still suffer from high cost, low efficiency, and poor service life problems. Recently, porous metal-based materials have been widely used as desirable cross-functional platforms for electrochemical and photochemical energy systems for their unique electrical conductivity, catalytic activity, and chemical stability. To tailor the porosity length scale, ordering, and compositions, 3D printing has been applied as a disruptive manufacturing revolution to create complex architected components by directly joining sequential layers into designed structures. In conclusion, this article intends to summarize cutting-edge advances of metal-based materials for renewable energy devices (e.g., fuel cells, solar cells, supercapacitors, and batteries) over the past decade.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Application of machine learning interatomic potentials in heterogeneous catalysis

Heterogeneous catalysts are crucial in modern societies as they promote sustainability by enabling lower-energy pathways for various chemical reactions. While Density Functional Theory (DFT) computations can provide critical insights into how heterogeneous catalysts operate at the atomic level, they are limited by computational costs and unfavorable scaling with system size. Recently, machine learning interatomic potentials (MLIPs) have emerged as a promising alternative to DFT, offering near-DFT accuracy at significantly reduced cost. Here, in this perspective, we discuss the application of MLIPs in heterogeneous catalyst modeling as a surrogate for DFT. We detail how MLIPs have been applied in thermal catalysis to probe active sites, enable studying complex metallic and nanoporous catalysts, and investigate the reconstruction of catalytic surfaces. We review the use of MLIPs in electrocatalysis and photocatalysis, emphasizing their capabilities in studying transition metal oxide surfaces and solid–liquid interfaces. We also discuss the current limitations of MLIPs, particularly their challenges with transferability and description of non-local interactions. Finally, we conclude by identifying promising and underexplored domains in which MLIPs can further advance our understanding of heterogeneous catalysts.

Catalytic surfaces↗
DOE OSTI2022

Using soil library hyperspectral reflectance and machine learning to predict soil organic carbon: Assessing potential of airborne and spaceborne optical soil sensing

Soil organic carbon (SOC) is a key variable to determine soil functioning, ecosystem services, and global carbon cycles. Spectroscopy, particularly optical hyperspectral reflectance coupled with machine learning, can provide rapid, efficient, and cost-effective quantification of SOC. However, how to exploit soil hyperspectral reflectance to predict SOC concentration, and the potential performance of airborne and satellite data for predicting surface SOC at large scales remain relatively underknown. Here, this study utilized a continental-scale soil laboratory spectral library (37,540 full-pedon 350–2500 nm reflectance spectra with SOC concentration of 0–780 g·kg –1 across the US) to thoroughly evaluate seven machine learning algorithms including Partial-Least Squares Regression (PLSR), Random Forest (RF), K-Nearest Neighbors (KNN), Ridge, Artificial Neural Networks (ANN), Convolutional Neural Networks (CNN), and Long Short-Term Memory (LSTM) along with four preprocessed spectra, i.e. original, vector normalization, continuum removal, and first-order derivative, to quantify SOC concentration. Furthermore, by using the coupled soil-vegetation-atmosphere radiative transfer model, we simulated twelve airborne and spaceborne hyper/multi-spectral remote sensing data from surface bare soil laboratory spectra to evaluate their potential for estimating SOC concentration of surface bare soils. Results show that LSTM achieved best predictive performance of quantifying SOC concentration for the whole data sets (R 2 = 0.96, RMSE = 30.81 g·kg –1 ), mineral soils (SOC ≤ 120 g·kg –1 , R 2 = 0.71, RMSE = 10.60 g·kg –1 ), and organic soils (SOC > 120 g·kg –1 , R 2 = 0.78, RMSE = 62.31 g·kg –1 ). Spectral data preprocessing, particularly the first-order derivative, improved the performance of PLSR, RF, Ridge, KNN, and ANN, but not LSTM or CNN. We found that the SOC models of mineral and organic soils should be distinguished given their distinct spectral signatures. Finally, we identified that the shortwave infrared is vital for airborne and spaceborne hyperspectral sensors to monitor surface SOC. This study highlights the high accuracy of LSTM with hyperspectral/multispectral data to mitigate a certain level of noise (soil moisture <0.4 m 3 ·m –3 , green leaf area < 0.3 m 2 ·m –2 , plant residue <0.4 m 2 ·m –2 ) for quantifying surface SOC concentration. Forthcoming satellite hyperspectral missions like Surface Biology and Geology (SBG) have a high potential for future global soil carbon monitoring, while high-resolution satellite multispectral fusion data can be an alternative.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Ethylene oxidation on unpromoted silver catalysts: Reaction pathway and selectivity analysis using DFT calculations

Chemical conversions in catalytic partial oxidation processes of light hydrocarbons are responsible for the production of numerous industrial chemicals, plastics, and intermediates. These processes are relatively expensive to perform, and are typically operated at high thermodynamic inefficiency, so the development of novel, highly efficient catalysts would prove to be very cost effective. Herein, our study focused on surface catalytic mechanisms of the ethylene oxide (EO) formation process. Periodic plane-wave Density Functional Theory (DFT) methods were used to analyze related reaction mechanisms on the Ag(111) surface facet with low coverage. Energetic changes of related species and pathways were calculated. Key surface species are identified to suggest the factors for the observed selectivity during EO formation. Our results are consistent with previous kinetic modeling efforts in the literature which did not employ DFT analysis. Lastly, our study demonstrates how fundamental theoretical investigations and multi-scale modeling techniques are currently impacting the advancement of rational catalyst design and microkinetic modeling techniques in the light hydrocarbon processing industry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗