Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lead optimization”

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,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

HydraGNN_Predictive_GFM_2026 - Ensemble of predictive graph foundation models for atomistic materials modeling

This release contains data and parameters of HydraGNN-based graph foundation models trained as a result of the work published in the pre-print "Exascale Multi-Task Graph Foundation Models for Imbalanced, Multi-Fidelity Atomistic Data" by M. Lupo Pasini et al. (https://arxiv.org/abs/2604.15380). We jointly train on 16 open first-principles datasets (544+ million structures covering 85+ elements) using a multi-task architecture with per-dataset heads and a scalable ADIOS2/DDStore data pipeline. On Frontier, we execute six large-scale DeepHyper hyperparameter optimization campaigns in FP64 and promote the top-performing message-passing models to sustained 2,048-node training, yielding a PaiNN-based lead model. The version of HydraGNN used to generate the outputs provided in this release is HydraGNN v5.0 (https://github.com/ORNL/HydraGNN/releases/tag/v5.0) The list of datasets used for the training of the graph foundation model is the following: 1) Alexandria [1] 2) ANI1x [2] 3) MPTrj [3] 4) Open Catalyst 2020 (OC20) [4] 5) Open Catalyst 2022 (OC22) [5] 6) Open Catalyst 2025 (OC25) [6] 7) Open Direct ir Capture 2023 (ODAC23) [7] 8) Open Materials 2024 (OMat24) [8] 9) Open Molecules 2025 (OMol25) [9] 10) OMol25-neutral (subset of OMol25 that contains only molecules with zero total charge) 11) OMol25-non-neutral (subset of OMol25 that contains only molecules with non-zero total charge) 12) Open Polymers 2026 (OPoly2026) [10] 13) Nabla2DFT [11] 14) QCML [12] 15) QM7X [reference 13] 16) transition1x [14] Dataset references: [1] J. Schmidt et al., “A dataset of 175k stable and metastable materials calculated with the PBEsol and SCAN functionals,” Scientific Data, vol. 9, p. 64, 2022. [2] J. S. Smith et al., “The ANI-1ccx and ANI-1x data sets, coupled-cluster and density functional theory properties for molecules,” Scientific Data, vol. 7, p. 134, 2020. [Online]. Available: https: //www.nature.com/articles/s41597-020-0473-z [3] A. Jain et al., “Commentary: The Materials Project: A materials genome approach to accelerating materials innovation,” APL Materials, vol. 1, no. 1, p. 011002, 07 2013. [Online]. Available: https://doi.org/10.1063/1.4812323 [4] L. Chanussot et al., “Open catalyst 2020 (oc20) dataset and community challenges,” ACS Catalysis, vol. 11, no. 10, pp. 6059–6072, 2021. [Online]. Available: https://doi.org/10.1021/acscatal.0c04525 [5] K. Tran et al., “Open catalyst 2022 (oc22) dataset and challenges for oxidation electrocatalysts,” ACS Catalysis, vol. 13, no. 5, pp. 3066–3084, 2023. [Online]. Available: https://doi.org/10.1021/acscatal.2c05426 [6] S. J. Sahoo et al., “The open catalyst 2025 (oc25) dataset and models for solid-liquid interfaces,” arXiv preprint arXiv:2509.17862, 2025. [Online]. Available: https://arxiv.org/abs/2509.17862 [7] A. Sriram et al., “The open DAC 2023 dataset and challenges for sorbent discovery in direct air capture,” ACS Central Science, vol. 10, no. 5, pp. 923–941, 2024. [8] L. Barroso-Luque et al., “Open materials 2024 (omat24) inorganic materials dataset and models,” 2024. [Online]. Available: https://arxiv.org/abs/2410.12771 [9] D. S. Levine et al., “The open molecules 2025 (OMol25) dataset, evaluations, and models,” 2025. [Online]. Available: https://arxiv.org/abs/2505.08762 [10] D. S. Levine et al., The open polymers 2026 (OPoly26) dataset and evaluations,” arXiv preprint arXiv:2512.23117, 2025. [Online]. Available: https://arxiv.org/abs/2512.23117 [11] K. Khrabrov et al., “Nabla2dft: A universal quantum chemistry dataset of drug-like molecules and a benchmark for neural network potentials,” in NeurIPS 2024 Datasets and Benchmarks Track, 2024. [Online]. Available: https://openreview.net/forum?id=ElUrNM9U8c [12] S. Ganscha et al., “The QCML dataset, quantum chemistry reference data from 33.5M DFT and 14.7B semi-empirical calculations,” Scientific Data, vol. 12, p. 406, 2025. [13] J. Hoja et al., “QM7-X, a comprehensive dataset of quantum-mechanical properties spanning the chemical space of small organic molecules,” Scientific Data, vol. 8, p. 43, 2021. [Online]. Available: https://www.nature.com/articles/s41597-021-00812-2 [14] M. Schreiner et al., “Transition1x - a dataset for building generalizable reactive machine learning potentials,” Scientific Data, vol. 9, p. 779, 2022. The folder "datasets_ADIOS2_format" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "datasets_ADIOS2_format" directory contains 2 sub-directories, one for the version "v1" of the datasets and one for the version "v2" of the datasets. The version "v1" of the datasets provides values of the total energy as they are extracted from the original data as it was released by the respective institutions. The version "v2" of the datasets provides values of the energy that have been realigned. The realignment was performed by training a linear regression model that predicts the total energy as a function of the chemical composition of the atomistic structure, and then subtract such prediction from the original value of the total energy. Both folders "v1" and "v2" contain 16 sub-directories, each corresponding to an ADIOS2-formatted dataset The folder "DeepHyper-results" contains the configurational files and model's parameters for all the 186 HPO trials that were successfully completed by the scalable hyperparameter optimization (HPO) runs on Frontier. The content of the folder "DeepHyper-results" I structured as follows: 1) task-list.txt: list of mpnn name, jobid, and deephyper task id 2) gfm_${MPNN}_${JOBID}_0.${TASKID}: run directory with checkpoint files 3) gfm_${MPNN}: deephyper summary directory (*.csv) for each specific MPNN type 4) deephyper-experiment-${JOBID}: output and error logs for each job The file "deephyper-sorted.csv" contains the details of each HydraGNN model built and tested by HPO, obtained by merging the (*.csv) filed from each HPO run executed. Out of all the HPO trials, we selected 10 to continue the training of the respective HydraGNN models. Due to limited computational budget available in the LRN070 allocation we could not complete the training till convergence for all these 10 selected models. The folder "models" contains multiple sub-folders, one per each HydraGNN model trained. Each model sub-folder contains the parameters of each HydraGNN model, with multiple checkpoint-restarts. The list of sub-folders are as follows: 1) multidataset_hpo-BEST1-fp64 2) multidataset_hpo-BEST2-fp64 3) multidataset_hpo-BEST3-fp64 4) multidataset_hpo-BEST4-fp64 5) multidataset_hpo-BEST5-fp64 6) multidataset_hpo-BEST6-fp64 7) multidataset_hpo-BEST7-fp64 8) multidataset_hpo-BEST8-fp64 9) multidataset_hpo-BEST9-fp64 10) multidataset_hpo-BEST10-fp64 Within each one of these folders, additional auxiliary log files are provided with descriptions about how the training proceeded. The lead PaiNN-model is contained inside "multidataset_hpo-BEST6-fp64". The file "mlp_branch_weights" contains the parameters of the multi-layer perceptron (MLP) used to reconcile the predictions of the 16 output decoding heads of the HydragNN architectures. The MLP takes in input the chemical composition of the atomistic structure and predicts averaging weights to linearly mix the predictions of each output decoding head toward consolidating them into a single one. The folder "1.1billion-structure-inference" contains 1.1 billion atomistic structures randomly generated. Each structures is associated with energy and forces predicted with the lead-PaiNN model combined with the MLP model for reconciliation of the multi-branch predictions generated by the 16 output decoding heads. The folder "1.1billion-structure-inference" contains 9,300 (*.tar.gz) subdirectories, one per Frontier compute node used to execute the inference at exascale. Once uncompressed, each (*.tar.gz) subdirectory contains an ADIOS2 (*.bp) file container, where each atomistic structure is stored as a PyTorch-Geometric Data object. The file "export_dataset_environment_variables.sh" contains the environment variables that need to be set before running the HydraGNN code to reproduce the results provided in this dataset release. The code that can be used to load the ADIOS2 files, load HydraGNN models, and run inference is available at: https://github.com/ORNL/HydraGNN/releases/tag/v5.0

36 MATERIALS SCIENCE↗
DOE OSTI2023

Optimization of operational strategies for rich gas enhanced oil recovery based on a pilot test in the Bakken tight oil reservoir

Horizontal well drilling and multistage hydraulic fracturing have been demonstrated as effective approaches for stimulating oil production in the Bakken tight oil reservoir. However, after multiple years of production, primary oil recovery in the Bakken is generally less than 10% of the estimated original oil in place. Gas huff ‘n’ puff (HnP) has been tested in the Bakken Formation as an enhanced oil recovery (EOR) method; however, most field pilot test results showed no significant incremental oil production. One of the factors affecting HnP EOR performance is premature gas breakthrough, which is one of the most critical issues observed in the field because of the presence of interwell fractures. Consequently, injected gas rapidly reaches adjacent production wells without contacting reservoir rock and increasing oil recovery. Proper conformance control is therefore needed to avoid early gas breakthrough and improve EOR performance. In this study, a rich gas EOR pilot in the Bakken was carefully analyzed to collect the essential reservoir and operational data. A simulation model with 16 wells was then developed to reproduce the production history and predict the EOR performance with and without conformance control. EOR operational strategies, including single- and multiple-well HnP, with different gas injection constraints were investigated. The simulation results of single-well HnP without conformance control showed that a rich gas injection rate of at least 10 MMscfd was needed to yield meaningful incremental oil production. The strategy of conformance control via water injection could significantly improve oil production in the HnP well, but injecting an excessive amount of water also leads to water breakthrough and loss of oil production in the offset wells. By analyzing the production performance of the wells individually, the arrangement of wells was optimized for multiple-well HnP EOR. The multiwell results showed that rich gas EOR could improve oil production up to 7.4% by employing conformance control strategies. Furthermore, replacing rich gas with propane as the injection gas could result in 14% of incremental oil production.

02 PETROLEUM↗
DOE OSTI2022

Mitigating the Impacts of Uncertain Geomagnetic Disturbances on Electric Grids: A Distributionally Robust Optimization Approach

Severe geomagnetic disturbances (GMDs) increase the magnitude of the electric field on the Earth’s surface (E-field) and drive geomagnetically-induced currents (GICs) along the transmission lines in electric grids. These additional currents can pose severe risks, such as current distortions, transformer saturation and increased reactive power losses, each of which can lead to system unreliability. Today several mitigation actions (e.g., changing grid topology) exist that can reduce the harmful GIC effects on the grids. Making such decisions can be challenging, however, because the magnitude and direction of the E-field are uncertain and non-stationary. In this paper, we model uncertain E-fields using the distributionally robust optimization (DRO) approach that determines optimal transmission grid operations such that the worst-case expectation of the system cost is minimized. We also capture the effect of GICs on the nonlinear AC power flow equations. For solution approaches, we develop an accelerated column-and-constraint generation (CCG) algorithm by exploiting a special structure of the support set of uncertain parameters representing the E-field. Extensive numerical experiments based on “epri-21” and “uiuc-150” systems, designed for GMD studies, demonstrate (i) the computational performance of the accelerated CCG algorithm, (ii) the superior performance of distributionally robust grid operations that satisfy nonlinear, nonconvex AC power flow equations and GIC constraints, in comparison with standard stochastic programming-based methods during the out-of-sample testing.

42 ENGINEERING↗
DOE OSTI2024

Optimizing the sextupole configuration for simultaneous correction of third order resonances at the recycler ring

For the Recycler Ring at Fermilab, space charge tune shifts of almost 0.1 will have to be dealt with under the Proton Improvement Plan (PIP-II) framework. This will lead to the excitation of third order resonances. The minimization of Resonance Driving Terms (RDTs) allows to mitigate the harmful effect of these betatron resonances. Past work has shown that previously-installed sextupoles can compensate the RDTs of individual third order resonance lines, thus reducing particle losses in these operational regimes. Nevertheless, trying to compensate multiple resonances of the same order simultaneously with these existing sextupoles is limited due to current constraints in the magnets. The following study showcases the procedure to install additional sextupoles in order to aid the compensation of multiple resonances. This includes the optimization of the new sextupoles' locations in order to cancel out multiple RDTs while minimizing the currents needed. This is followed by a verification of their effectiveness by means of the RDT response matrix.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Variational Monte Carlo Calculations of A ≤ 4 Nuclei with an Artificial Neural-Network Correlator Ansatz

Here, the complexity of many-body quantum wave functions is a central aspect of several fields of physics and chemistry where nonperturbative interactions are prominent. Artificial neural networks (ANNs) have proven to be a flexible tool to approximate quantum many-body states in condensed matter and chemistry problems. In this work we introduce a neural-network quantum state ansatz to model the ground-state wave function of light nuclei, and approximately solve the nuclear many-body Schrodinger equation. Using efficient stochastic sampling and optimization schemes, our approach extends pioneering applications of ANNs in the field, which present exponentially scaling algorithmic complexity. We compute the binding energies and point-nucleon densities of A ≤ 4 nuclei as emerging from a leading-order pionless effective field theory Hamiltonian. We successfully benchmark the ANN wave function against more conventional parametrizations based on two- and three-body Jastrow functions, and virtually exact Green's function Monte Carlo results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1999

Ames Optimized TCA Configuration

Configuration design at Ames was carried out with the SYN87-SB (single block) Euler code using a 193 x 49 x 65 C-H grid. The Euler solver is coupled to the constrained (NPSOL) and the unconstrained (QNMDIF) optimization packages. Since the single block grid is able to model only wing-body configurations, the nacelle/diverter effects were included in the optimization process by SYN87's option to superimpose the nacelle/diverter interference pressures on the wing. These interference pressures were calculated using the AIRPLANE code. AIRPLANE is an Euler solver that uses a unstructured tetrahedral mesh and is capable of computations about arbitrary complete configurations. In addition, the buoyancy effects of the nacelle/diverters were also included in the design process by imposing the pressure field obtained during the design process onto the triangulated surfaces of the nacelle/diverter mesh generated by AIRPLANE. The interference pressures and nacelle buoyancy effects are added to the final forces after each flow field calculation. Full details of the (recently enhanced) ghost nacelle capability are given in a related talk. The pseudo nacelle corrections were greatly improved during this design cycle. During the Ref H and Cycle 1 design activities, the nacelles were only translated and pitched. In the cycle 2 design effort the nacelles can translate vertically, and pitch to accommodate the changes in the lower surface geometry. The diverter heights (between their leading and trailing edges) were modified during design as the shape of the lower wing changed, with the drag of the diverter changing accordingly. Both adjoint and finite difference gradients were used during optimization. The adjoint-based gradients were found to give good direction in the design space for configurations near the starting point, but as the design approached a minimum, the finite difference gradients were found to be more accurate. Use of finite difference gradients was limited by the CPU time limit available on the Cray machines. A typical optimization run using finite difference gradients can use only 30 to 40 design variables and one optimization iteration within the 8 hour queue limit for the chosen grid size and convergence level. The efficiency afforded by the adjoint method allowed for 50-120 design variables and 5-10 optimization iterations in the 8 hour queue. Geometric perturbations to the wing and fuselage were made using the Hicks/Henne (HH) shape functions. The HH functions were distributed uniformly along the chords of the wing defining sections and lofted linearly. During single-surface design, constraints on thickness and volume at selected wing stations were imposed. Both fuselage camber and cross-sectional area distributions were permitted to change during design. The major disadvantage to the use of these functions is the inherent surface waviness produced by repeated use of such functions. Many smoothing operations were required following optimization runs to produce a configuration with reasonable smoothness. Wagner functions were also used on the wing sections but were never used on the fuselage. The Wagner functions are a family of increasingly oscillatory functions that have also been used extensively in airfoil design. The leading and trailing edge regions of the wing were designed by use of polynomial and monomial functions respectively. Twist was attempted but was abandoned because of little performance improvement available from changing the baseline twist.

Cliff, Susan E.↗
DOE OSTI2022

Bounds on Optimal Revisit Times in Persistent Monitoring Missions With a Distinct and Remote Service Station</strong

We report persistent monitoring missions require an up-to-date knowledge of the changing state of the underlying environment. Unmannned aerial vehicles (UAVs) can be gainfully employed to continually visit a set of targets representing tasks (and locations) in the environment and collect data therein for long time periods. The enduring nature of these missions requires the UAV to be regularly recharged at a service station. In this article, we consider the case in which the service station is not colocated with any of the targets. An efficient monitoring requires the revisit time, defined as the maximum of the time elapsed between successive revisits to targets, to be minimized. Here, we consider the problem of determining UAV routes that lead to the minimum revisit time. The problem is NP-hard, and its computational difficulty increases with the fuel capacity of the UAV. We develop an algorithm to construct near-optimal solutions to the problem quickly when the fuel capacity exceeds a threshold. We also develop lower bounds to the optimal revisit time and use these bounds to demonstrate (through numerical simulations) that the constructed solutions are, on an average, at most 0.01% away from the optimum.

42 ENGINEERING↗
DOE OSTI2025

Quantum Monte Carlo Approaches to Na Intercalation on Bilayer Graphene

We have performed Quantum Monte Carlo (QMC) simulations on Na-intercalated bilayer graphene to study the evolution of electronic and optical properties upon Na intercalation into hard carbon layers. The objective was to model the optimal configuration of Na intercalation into a hard carbon matrix containing graphene regions. Our study showed that Na intercalation can be energetically stabilized at large interlayer distances (over 6 Å) in both AA- and AB-stacked bilayer graphene. In the QMC results, we found a significant band gap opening at the equilibrium interlayer distance of Na-intercalated bilayer graphene, while corresponding density functional theory (DFT) results showed no gap. This difference between DFT and QMC results indicates that the gap opening induced by Na intercalation into a hard carbon is underestimated within the DFT framework. In addition, a zigzag configuration of Na atoms was found to be energetically stable at interlayer distances up to 10 Å, leading us to predict the existence of a local minimum of Na intercalation at large interlayer distance. These computation and modeling results can provide guidance on how to synthesize and optimize hard carbon with bilayer graphene regions that permit a zigzag intercalation configuration that will maximize and stabilize sodium hosting.

Binding energy↗
DOE OSTI2023

Face‐on Orientation Matches Vertical Organic Electrochemical Transistors for High Transconductance and Superior Non‐Volatility

Abstract The recently developed vertical structure of organic electrochemical transistors (OECTs) can integrate volatile and non‐volatile functions into one reconfigurable device, making it highly promising. However, comparing with the conventional planar OECT (c‐OECT), the understanding of vertical OECT (v‐OECT) working principles and device engineering strategies is still lacking, impeding rational optimization. Since a major difference between c‐ and v‐OECTs is their charge transport directionality, which is highly influenced by crystallite orientations, the orientation–device structure match thus becomes an important yet outstanding topic for OECTs. Herein, using ann‐type small molecule IDIC‐MEG, investigate how much impact such match can have on OECT performance. The IDIC‐MEG c‐OECT fails to work due to the seriously hindered in‐plane electron transport by face‐on orientation. Surprisingly, simply changing the device structure from planar to vertical allows the resultant v‐OECT to exhibit the highest reported transconductance (46.3 mS) among all small‐molecule OECTs, thanks to the match between face‐on orientation and the vertical structure. Such match also leads to excellent non‐volatility, including highly predictable programmability and good operational stability. This work, for the first time, explicitly demonstrates the significance of orientation–device structure match for OECT optimization, establishing new guidelines for achieving high‐performance volatile and non‐volatile OECTs.

Chemistry↗
DOE OSTI2025

Organic Electrochemical Transistor Channel Materials: Copolymerization Versus Physical Mixing of Glycolated and Alkoxylated Polymers

Organic electrochemical transistors (OECTs) feature a polymer channel capable of conducting both ions and electronic charges. The choice of the channel material is critical for OECT performance. Many efforts have focused on improving performance via the chemical tunability of conjugated polymers – through backbone, side chain, and molar mass engineering – leading to useful design principles for accumulation-mode OECT materials. However, tuning the chemical structure of conjugated polymers often requires time-consuming optimization of the synthesis route. Meanwhile, variations in molar mass, dispersity, structural defects, and metal content present challenges when attempting to analyze the detailed effects of structural modifications, as multiple performance-determining factors are often (unintentionally) changed at the same time. Therefore, this study explores blended channel materials obtained by physically mixing glycolated and alkoxylated polymers in different ratios, and compares their OECT performance with the corresponding statistical copolymers. It is shown that mixing two well-performing materials creates blends that enable rational tuning of the transistor properties without compromising on performance. Thus, channels based on blends of alkoxylated and glycolated polymers hold promise for OECT technology with tailored response, as only two materials are needed to achieve any desired side chain ratio, simplifying the optimization of OECT characteristics.

copolymerization↗
DOE OSTI2024

Physics-Informed Active Learning With Simultaneous Weak-Form Latent Space Dynamics Identification

The parametric greedy latent space dynamics identification (gLaSDI) framework has demonstrated promising potential for accurate and efficient modeling of high-dimensional nonlinear physical systems. However, it remains challenging to handle noisy data. Here, to enhance robustness against noise, we incorporate the weak-form estimation of nonlinear dynamics (WENDy) into gLaSDI. In the proposed weak-form gLaSDI (WgLaSDI) framework, an autoencoder and WENDy are trained simultaneously to discover intrinsic nonlinear latent-space dynamics of high-dimensional data. Compared with the standard sparse identification of nonlinear dynamics (SINDy) employed in gLaSDI, WENDy enables variance reduction and robust latent space discovery, therefore leading to more accurate and efficient reduced-order modeling. Furthermore, the greedy physics-informed active learning in WgLaSDI enables adaptive sampling of optimal training data on the fly for enhanced modeling accuracy. The effectiveness of the proposed framework is demonstrated by modeling various nonlinear dynamical problems, including viscous and inviscid Burgers' equations, time-dependent radial advection, and the Vlasov equation for plasma physics. With data that contains 5%–10% Gaussian white noise, WgLaSDI outperforms gLaSDI by orders of magnitude, achieving 1%–7% relative errors. Compared with the high-fidelity models, WgLaSDI achieves 121 to 1779x speed-up.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Carbon fiber synthesis from pitch: Insights from ReaxFF based molecular dynamics simulations

Despite the potential to lowering fabrication costs of pitch-based carbon fiber (CF) compared to polyacrylonitrile-based CF, large scale utilization is currently limited by the complexity in manufacturing. This is mainly due to the lack of consistency between the pitch structure/chemistry and CF properties under different processing conditions. Here, we employed full atomistic simulations based on reactive force fields to investigate the conversion of pitch precursors to CFs. Additionally, it was found that at the carbonization stage, the consumption rates of precursors, gas formations, and solid yields (carbonized CFs) are all dependent on the composition and arrangement of molecular moieties in the precursor. Specifically, precursors with sp 3 carbons in the cata-condensed core area have limited decomposition pathways (i.e., limited reaction sites), which leads to the formation of carbonized CFs with similar compositions. Furthermore, carbonized CFs possess distinct structural characteristics, and finally lead to the formation of CFs with different moduli. By analyzing CF crystallinity and elemental composition, general rationales are proposed for the selection of optimal precursors and processing conditions, and thus shed lights on parameter screening during industrial CF synthesis. For instance, inclusion of sp 3 carbons in cata- condensed rings can be preferred over the presence of sp 3 carbons on side chains.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

A high-performance implementation of atomistic spin dynamics simulations on x86 CPUs

Atomistic spin dynamics simulations provide valuable information about the energy spectrum of magnetic materials in different phases, allowing one to identify instabilities and the nature of their excitations. However, the time cost of evaluating the dynamical correlation function S(q,t) increases quadratically as the number of spins N, leading to significant computational effort, making the simulation of large spin systems very challenging. In this work, here, we propose to use a highly optimized general matrix multiply (GEMM) subroutine to calculate the dynamical spin-spin correlation function that can achieve near-optimal hardware utilization. Furthermore, we fuse the element-wise operations in the calculation of S(q,t) into the in-house GEMM kernel, which results in further performance improvements of 44\% - 71\% on several relatively large lattice sizes when compared to the implementation that uses the GEMM subroutine in OpenBLAS, which is the state-of-the-art open source library for Basic Linear Algebra Subroutine (BLAS).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Strangers in a foreign land: ‘Yeastizing’ plant enzymes

Abstract Expressing plant metabolic pathways in microbial platforms is an efficient, cost‐effective solution for producing many desired plant compounds. As eukaryotic organisms, yeasts are often the preferred platform. However, expression of plant enzymes in a yeast frequently leads to failure because the enzymes are poorly adapted to the foreign yeast cellular environment. Here, we first summarize the current engineering approaches for optimizing performance of plant enzymes in yeast. A critical limitation of these approaches is that they are labour‐intensive and must be customized for each individual enzyme, which significantly hinders the establishment of plant pathways in cellular factories. In response to this challenge, we propose the development of a cost‐effective computational pipeline to redesign plant enzymes for better adaptation to the yeast cellular milieu. This proposition is underpinned by compelling evidence that plant and yeast enzymes exhibit distinct sequence features that are generalizable across enzyme families. Consequently, we introduce a data‐driven machine learning framework designed to extract ‘yeastizing’ rules from natural protein sequence variations, which can be broadly applied to all enzymes. Additionally, we discuss the potential to integrate the machine learning model into a full design‐build‐test cycle.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Efficient Training of Deep Neural Operator Networks via Randomized Sampling

Neural operators (NOs) employ deep neural networks to learn the mappings between infinitedimensional function spaces. Deep operator network (DeepONet), a popular NO architecture, has demonstrated success in the real-time prediction of complex dynamics across various scientific and engineering applications. In this work, we introduce a random sampling technique to be adopted during the training of DeepONet, aimed at improving the generalization ability of the model, while significantly reducing the computational time. The proposed approach targets the trunk network of the DeepONet model that outputs the basis functions corresponding to the spatiotemporal locations of the bounded domain on which the physical system is defined. While constructing the loss function, DeepONet training traditionally considers a uniform grid of spatiotemporal points at which all the output functions are evaluated for each iteration. This approach leads to a larger batch size, resulting in poor generalization and increased memory demands, due to the limitations of the stochastic gradient descent (SGD) optimizer. The proposed random sampling over the inputs of the trunk net mitigates these challenges, improving generalization and reducing the memory requirements during training, resulting in significant computational gains. We validate our hypothesis through three benchmark examples, demonstrating substantial reductions in training time while achieving comparable or lower overall test errors relative to the traditional training approach. Our results indicate that incorporating randomization in the trunk network inputs during training enhances the efficiency and robustness of DeepONet, offering a promising avenue for improving the framework’s performance in modeling complex physical systems.

Karumuri, Sharmila [Department of Civil & Systems ↗
DOE OSTI2022

The Critical Stack Pressure to Alter Void Generation at Li/Solid-Electrolyte Interfaces during Stripping

The lithium stripping process generates vacancies, which may accumulate as voids and lead to uneven current distribution and dendrite growth in the following plating cycles. A stack pressure is typically required during stripping, but how to optimize the stack pressure is not clear. In this work, extremely lithiophilic Li/Li 2 O and lithiophobic Li/LiF interfaces were used to reveal the combining effect of interface interaction and stack pressure induced lithium creep on the stripping critical current density (CCD). A multiscale simulation scheme with Density Functional Theory (DFT), kinetic Monte Carlo (KMC) simulations, and an analytical model was developed. The analytical model predicted lithiophobic interfaces require a higher stack pressure than lithiophilic interfaces to reach the same CCD. The KMC simulations also showed higher stack pressure is needed at lithiophobic interfaces to accelerate Li vacancy diffusion into the bulk and maintain a flat surface. This stack pressure needs to be high enough to alter the Li forward-and-backward hopping barriers at the interface. This multiscale simulation scheme illustrates the importance to include the chemical-mechanical effects during Li stripping morphology evolution. It can be used to design ideal interlayer coating materials to maintain a flat Li surface during cycling.

25 ENERGY STORAGE↗
DOE OSTI2025

Optimal production of Phanerochaete chrysosporium manganese peroxidases and Trametes sp. C30 laccase hybrid Lac131 in Aspergillus niger for lignin bioconversion

Background Incorporating the production of related ligninolytic enzymes into industrial filamentous fungus Aspergillus niger will enhance the bioconversion of lignocelluloses to various chemical products. Results In this study, transgenic expression of Phanerochaete chrysosporium manganese peroxidases (mnps) and Trametes sp. C30 laccase hybrid Lac131 (lac131) were examined and optimized in A. niger 11414 prtT∆ strain. Five mnps (mnp1, mnp2, mnp3, mnp4, and mnp5) and lac131 genes were expressed separately or in combination. The transgenic strain containing the entire mnp2 genomic coding sequence (gmnp2) exhibited the highest mnP activity among the five mnp over-expression strains in the modified minimal medium (mMM) with addition of 5 g/L bovine hemoglobin (bHg). We examined the effects of hemin and bHg on mnP production in the gmnp2 strain cultures and found that at least 1 g/L bHg was required, while hemin was not. Culture conditions for mnP production were further optimized for the gmnp2 strain and the highest mnP activities were detected in the cultures grown at 25 °C and 200 rpm with an initial pH of 4.5. Effects of soy protein, skim milk, and bovine serum albumin on mnP production were investigated; 5 g/L of soy proteins or skim milk had comparable effects to 2.5 g/L bHg, while cultures with bovine serum albumin had diminished mnP activity. Disruption of both prtT and vsm1 substantially augmented the mnP production and its activity reached 575 U/L. Trametes sp. C30 laccase hybrid lac131 was strongly expressed in either A. niger gmnp2 (1975 U/L) or 11414prtT∆ (3895 U/L) strain. Both mnP and laccase in the culture supernatants effectively decolorized selected phenolic compounds (dyes) and cleaved tagged model lignin dimers. Conclusion The mnP was successfully produced in A. niger by optimizing the culture conditions and host strain. Co-expression of all four mnp genes in the same expression host by multiplex CRISPR will lead to the mnP production reaching levels comparable to P. chrysosporium, while only requiring 36 h at 25 °C. The Lac131 activity in transgenic A. niger strain is 4- to 7-times higher than that in previous studies. Co-production of mnP and laccase in A. niger will enhance the lignin bioconversion efficiency.

Aspergillus niger↗
DOE OSTI2023

Life cycle cost assessment of material recovery from perovskite solar cells

Perovskite photovoltaic (PV) cells have created a significant interest over the last few years due to their low-cost and high-power conversion efficiencies. While perovskite PVs are still under development, we analyze cost of alternative sustainable end-of-life management for perovskite PV cells; this allows to reuse economically valuable materials and prevent environmental contamination. Here, we studied the life cycle cost assessment of recovering metals such as lead, aluminum, gold, nickel, and silver and other valuable materials such as glass from waste perovskite PVs. We also developed a recycling scheme to optimize the material recovery and determine the economically feasible ways to separate the cells layer-by-layer. Further, we assessed the cost and the net cost of the perovskite PV recycling processes. We found that the most feasible recycling technique’ cost is $\$$10.70 and its net cost is $\$$– 2.95 per 1 m 2 module. These values could be further improved by optimizing and reusing chemicals involved in process. These results indicate that the cost of perovskite PV recycling is economically feasible and there is potentional to gain benefit from sustainable end-of-life management of perovskite PVs.

14 SOLAR ENERGY↗