Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ensemble 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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Inference of bipolar neutrino flavor oscillations near a core-collapse supernova based on multiple measurements at Earth

Neutrinos in compact-object environments, such as core-collapse supernovae, can experience various kinds of collective effects in flavor space, engendered by neutrino-neutrino interactions. These include “bipolar” collective oscillations, which are exhibited by neutrino ensembles where different flavors dominate at different energies. Considering the importance of neutrinos in the dynamics and nucleosynthesis in these environments, it is desirable to ascertain whether an Earth-based detection could contain signatures of bipolar oscillations that occurred within a supernova envelope. To that end, we, in this study, continue examining a cost-function formulation of statistical data assimilation (SDA) to infer solutions to a small-scale model of neutrino flavor transformation. SDA is an inference paradigm designed to optimize a model with sparse data. Our model consists of two monoenergetic neutrino beams with different energies emanating from a source and coherently interacting with each other and with a matter background, with radially varying interaction strengths. We attempt to infer flavor transformation histories of these beams using simulated measurements of the flavor content at locations “in vacuum” (that is, far from the source), which could in principle correspond to Earth-based detectors. Within the scope of this small-scale model, we found that: (i) based on such measurements, the SDA procedure is able to infer whether bipolar oscillations had occurred within the protoneutron star envelope, and (ii) if the measurements sample the full amplitude of the neutrino oscillations in vacuum, then the amplitude of the prior bipolar oscillations is well predicted. This result intimates that the inference paradigm can well complement numerical integration codes, via its ability to infer flavor evolution at physically inaccessible locations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Preparation of metrological states in dipolar-interacting spin systems

Abstract Spin systems are an attractive candidate for quantum-enhanced metrology. Here we develop a variational method to generate metrological states in small dipolar-interacting spin ensembles with limited qubit control. For both regular and disordered spatial spin configurations the generated states enable sensing beyond the standard quantum limit (SQL) and, for small spin numbers, approach the Heisenberg limit (HL). Depending on the circuit depth and the level of readout noise, the resulting states resemble Greenberger-Horne-Zeilinger (GHZ) states or Spin Squeezed States (SSS). Sensing beyond the SQL holds in the presence of finite spin polarization and a non-Markovian noise environment. The developed black-box optimization techniques for small spin numbers ( N ≤ 10) are directly applicable to diamond-based nanoscale field sensing, where the sensor size limits N and conventional squeezing approaches fail.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Effective optimization of atomic decoration in giant and superstructurally ordered crystals with machine learning

Crystals with complicated geometry are often observed with mixed chemical occupancy among Wyckoff sites, presenting a unique challenge for accurate atomic modeling. Similar systems possessing exact occupancy on all the sites can exhibit superstructural ordering, dramatically inflating the unit cell size. In this work, a crystal graph convolutional neural network (CGCNN) is used to predict optimal atomic decorations on fixed crystalline geometries. This is achieved with a site permutation search (SPS) optimization algorithm based on Monte Carlo moves combined with simulated annealing and basin-hopping techniques. Our approach relies on the evidence that, for a given chemical composition, a CGCNN estimates the correct energetic ordering of different atomic decorations, as predicted by electronic structure calculations. This provides a suitable energy landscape that can be optimized according to site occupation, allowing the prediction of chemical decoration in crystals exhibiting mixed or disordered occupancy, or superstructural ordering. Verification of the procedure is carried out on several known compounds, including the superstructurally ordered clathrate compound Rb8Ga27Sb16 and vacancy-ordered perovskite Cs2SnI6, neither of which was previously seen during the neural network training. In addition, the critical temperature of an order–disorder phase transition in solid solution CuZn is probed with our SPS routines by sampling site configuration trajectories in the canonical ensemble. This strategy provides an accurate method for determining favorable decoration in complex crystals and analyzing site occupation at unprecedented speed and scale.

Chemistry↗
DOE OSTI2020

Predicting synthesizable multi-functional edge reconstructions in two-dimensional transition metal dichalcogenides

Abstract Two-dimensional (2D) transition metal dichalcogenides (TMDCs) have attracted tremendous interest as functional materials due to their exceptionally diverse and tunable properties, especially in their edges. In addition to the conventional armchair and zigzag edges common to hexagonal 2D materials, more complex edge reconstructions can be realized through careful control over the synthesis conditions. However, the whole family of synthesizable, reconstructed edges remains poorly studied. Here, we develop a computational approach integrating ensemble-generation, force-relaxation, and electronic-structure calculations to systematically and efficiently discover additional reconstructed edges and screen their functional properties. Using MoS 2 as a model system, we screened hundreds of edge-reconstruction to discover over 160 reconstructed edges to be more stable than the conventional ones. More excitingly, we discovered nine new synthesizable reconstructred edges with record thermodynamic stability, in addition to successfully reproducing three recently synthesized edges. We also find our predicted reconstructed edges to have multi-functional properties—they show near optimal hydrogen evolution activity over the conventional edges, ideal for catalyzing hydrogen-evolution reaction (HER) and also exhibit half-metallicity with a broad variation in magnetic moments, making them uniquely suitable for nanospintronic applications. Our work reveals the existence of a wide family of synthesizable, reconstructed edges in 2D TMDCs and opens a new materials-by-design paradigm of ‘intrinsic’ edge engineering multifunctionality in 2D materials.

Hu, Guoxiang (ORCID:0000000329428564)↗
DOE OSTI2021

Predicting Thermal Performance of an Enhanced Geothermal System From Tracer Tests in a Data Assimilation Framework

Abstract Predicting the thermal performance of an enhanced geothermal system (EGS) requires a comprehensive characterization of the underlying fracture flow patterns from practically available data such as tracer data. However, due to the inherent complexities of subsurface fractures and the generally insufficient geological/geophysical data, interpreting tracer data for fracture flow characterization and thermal prediction remains a challenging task. The present study aims to tackle the challenge by leveraging a data assimilation method to maximize the utilization of information inherently contained in tracer data, and meanwhile maintain the flexibility to handle various uncertainties. A tracer data interpretation framework was proposed with the following three components integrated: (a) We use principal component analysis (PCA) to reduce the dimensionality of model parameter space. (b) We use ES‐MDA (ensemble smoother with multiple data assimilation) to invert for fracture aperture/flow fields and obtain posterior model ensembles for uncertainty quantification. Various data types are assimilated jointly to improve the predictive ability of the posterior ensemble. (c) The inverted fracture aperture fields are then incorporated into reservoir models to predict thermal performance. We developed a field‐scale EGS model to verify the ability of the framework to characterize highly heterogeneous fracture aperture/flow fields and predicting thermal performance. We also applied the framework to a mesoscale field experiment to demonstrate its potential application in real‐world geothermal reservoirs. The results indicate that the proposed framework can effectively retrieve fracture flow information from tracer data for thermal prediction and uncertainty quantification, and thus provide informative guidance for EGS optimization and risk management.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Mesoscale Modeling to Characterize Eagle Soaring Habitat

Uncovering drivers of risk is crucial to understanding interactions between wildlife and wind turbines, and identifying options for impact minimization. These drivers tie to co-variates linked to behavior and movement patterns that allow us to estimate locations and periods of risk. For volant species, atmospheric flow can have significant influence on flight patterns. For obligate soaring birds, like golden eagles, updraft velocities can inform where eagles are likely to travel, at what altitude, and where conditions are not likely sufficient to sustain soaring flight. This has been an active area of study in recent years, using relatively coarse atmospheric data generally at the 20km x 20km scale or larger. Leveraging a 20-year dataset from the Weather Research and Forecasting Model (WRF) (https://www.mmm.ucar.edu/weather-research-and-forecasting-model), we are quantifying vertical velocities across the continental United States at a 2km x 2km resolution. Wind resource data sets originally were static maps showing the mean annual wind speed over an area. However, for these data sets to be optimally used for various applications they must be high-resolution time series, seamlessly span large geographic contexts, and account for uncertainty in wind speed. The National Renewable Energy Laboratory is producing public available datasets that meet these criteria and will be bias corrected to yield the most accurate wind resource data. This effort is an augment to the current WIND Toolkit which houses a high resolution data set. In the new iteration of the WIND Toolkit, a 20-year dataset will be used to improve the accuracy and estimate uncertainty using ensemble and machine-learning techniques. The resulting product will be the most accurate dataset of its size and at a 2km x 2km spatial and 5-minute temporal resolution. Through this work, a mesoscale vertical velocity layer will be produced by calculating the likelihood of orographic updraft and thermal updraft conditions across the continental United States. Specifically, we will use WRF model output combined with digital elevation maps to predict updrafts and then determine if vertical velocities are sufficient to support Golden Eagle soaring and gliding. Ultimately this data layer will be made available as a GIS layer in the Wind Prospector (maps.nrel.gov/wind-prospector/) tool or a similar framework. Data that will be incorporated include wind speed, direction temperature, relative humidity, barometric pressure, air density, precipitation rate, solar radiation, atmospheric stability, skin temperature, and upward heat flux. These products will advance research on interactions between volant species and wind energy by providing open access to highly resolved data with uncertainty quantification not previously available at this scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Variable domain mutational analysis to probe the molecular mechanisms of high viscosity of an IgG 1 antibody

Subcutaneous injection is the preferred route of administration for many antibody therapeutics for reasons that include its speed and convenience. However, the small volume limit (typically ≤2 mL) for subcutaneous delivery often necessitates antibody formulations at high concentrations (commonly ≥100 mg/mL), which may lead to physicochemical problems. For example, antibodies with large hydrophobic or charged patches can be prone to self-interaction giving rise to high viscosity. Here, we combined X-ray crystallography with computational modeling to predict regions of an anti-glucagon receptor (GCGR) IgG 1 antibody prone to self-interaction. An extensive mutational analysis was undertaken of the complementarity-determining region residues residing in hydrophobic surface patches predicted by spatial aggregation propensity, in conjunction with residue-level solvent accessibility, averaged over conformational ensembles from molecular dynamics simulations. Dynamic light scattering (DLS) was used as a medium throughput screen for self-interaction of ~200 anti-GCGR IgG 1 variants. A negative correlation was found between the viscosity determined at high concentration (180 mg/mL) and the DLS interaction parameter measured at low concentration (2–10 mg/mL). Additionally, anti-GCGR variants were readily identified with reduced viscosity and antigen-binding affinity within a few fold of the parent antibody, with no identified impact on overall developability. The methods described here may be useful in the optimization of other antibodies to facilitate their therapeutic administration at high concentration.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Development of uncertainty-aware equation-of-state models: Application to copper

Sophisticated hydrodynamic codes are commonly used to understand and predict events relevant to natural and applied sciences. The degree to which these simulations reflect reality, however, is dependent on how well we understand the materials and underlying physics involved. These research communities need material models that communicate the uncertainty in the physical properties, which at their basest form comes from the uncertainty in the underlying experimental measurements. Additionally, we have constructed a new framework for using both experimental measurements and the associated experimental uncertainties to build equation-of-state models that reflect not only current best measurements but also the accuracy of those measurements. This method had been used to construct an ensemble of equation-of-state models for copper that communicates the experimental uncertainties in the data through the equation-of-state model, which is available for application in any simulation metric of interest.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Single-Molecule Fluorescence Microscopy Reveals Energy Transfer Active versus Inactive Nanocrystal/Dye Conjugate Pairs

Defect-mediated energy transfer (EnT) is a radiative process that occurs between donor defect states in the forbidden bandgap of semiconductor nanocrystals (NCs) and dye molecules bound to their surfaces. The EnT efficiency depends on the number of dye molecules attached to each NC, the donor–acceptor distance, and the dipole orientation factor between the donor and acceptor, all of which vary across all individual NCs in a sample. While ensemble-level fluorescence spectroscopy measurements have provided average values for donor–acceptor distances, dye-to-NC ratios, and EnT rate constants, questions remain about the impact of donor/acceptor heterogeneity on observed EnT efficiencies. Notably, ensemble-level measurements cannot distinguish between bare NCs and EnT-active versus inactive NC/dye pairs in the same sample batch, limiting the ability to design systems with 100% EnT efficiency. To address this, we studied defect-mediated EnT between AlexaFluor 555 dye acceptors chemically bound to ZnO NC donors at the level of single molecules and single NCs. Interestingly, 20% of bound NC/dye pairs are EnT-inactive, likely contributing to residual defect photoluminescence (PL) observed in ensemble-level measurements and reducing overall EnT efficiency. Single particle-level ZnO defect PL and acceptor fluorescence trajectories exhibited distinct microfluctuations, which are absent in bare ZnO NCs. We hypothesized that our observations can be explained with a competitive dye fluorescence quenching pathway, possibly due to charge transfer between the excited state dye and the ZnO NC. Numerical simulations of single-molecule PL traces for this scenario produced microfluctuations consistent with the experimental results. These findings highlight the impact of sample heterogeneity on EnT processes and provide insights for designing light-harvesting systems with optimized EnT efficiency.

25 ENERGY STORAGE↗
DOE OSTI2023

Agilent AgileBioFoundry CRADA (Final Report)

The mission of this CRADA with Agilent was to couple powerful MS platforms (QQQ, IM-QTOF-MS) with Agilent’s novel Ultra-High-Performance Liquid Chromatography (UHPLC) fast metabolomic workflows and perform ABF Machine Learning (ML) to generated datasets. Agilent transferred UHPLC methods to PNNL and LBNL and methods were implemented and demonstrated in both labs, achieving total acquisition times of < 10 min. Metabolites analyzed using Agilent’s shared methods included metabolites from central carbon metabolism, common across hosts, and metabolites unique to engineered strains. Standards were acquired in an UHPLC-Drift Tube Ion Mobility Mass Spectrometer (DTIMS) system for the first time within the context of ABF and methods were optimized based on Agilent’s protocols. Samples from ABF hosts Pseudomonas putida, Aspergillus pseudoterreus, Aspergillus niger and Rhodosporidium toruloides were analyzed using the UHPLC-DTIMS platform for a total of 276 runs. A data analysis workflow compatible with the Experimental Data Depot (EDD) and completely shareable was developed for the acquired UHPLC-DTIMS data. Samples were analyzed using a Data Independent Acquisition Approach (DIA), which for most of the standards provided more transitions therefore increasing detection confidence. Using the data acquired by PNNL, LBNL, and Agilent’s specifications from previous ML projects, SNL applied an ensemble ML strategy to pick the best performing model for automated LC-method selection. Finally, with the contribution of the participant labs and Agilent, SNL developed an Automated Method Selection (AMS) software tool to predict the best liquid chromatography method for analysis of any new molecules of interest. Samples with novel pathways and new metabolite targets of interest are generated at a high pace in the ABF. Therefore, our accomplishment in this CRADA improved the efficiency and accuracy of strain testing by developing and implementing fast analytical methods, robust processing tools, and software for predicting the best methods for UHPLC analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Development and implementation of high-throughput proteomic and metabolomics assays by using advanced chromatographic and mass spectrometric systems (CRADA Final Report)

The mission of this CRADA with Agilent was to couple powerful MS platforms (QQQ, IM-QTOFMS) with Agilent’s novel Ultra-High-Performance Liquid Chromatography (UHPLC) fast metabolomic workflows and perform ABF Machine Learning (ML) to generated datasets. Agilent transferred UHPLC methods to PNNL and LBNL and methods were implemented and demonstrated in both labs, achieving total acquisition times of < 10 min. Metabolites analyzed using Agilent’s shared methods included metabolites from central carbon metabolism, common across hosts, and metabolites unique to engineered strains. Standards were acquired in an UHPLC-Drift Tube Ion Mobility Mass Spectrometer (DTIMS) system for the first time within the context of ABF and methods were optimized based on Agilent’s protocols. Samples from ABF hosts Pseudomonas putida, Aspergillus pseudoterreus, Aspergillus niger and Rhodosporidium toruloides were analyzed using the UHPLC-DTIMS platform for a total of 276 runs. A data analysis workflow compatible with the Experimental Data Depot (EDD) and completely shareable was developed for the acquired UHPLC-DTIMS data. Samples were analyzed using a Data Independent Acquisition Approach (DIA), which for most of the standards provided more transitions therefore increasing detection confidence. Using the data acquired by PNNL, LBNL, and Agilent’s specifications from previous ML projects, SNL applied an ensemble ML strategy to pick the best performing model for automated LC-method selection. Finally, with the contribution of the participant labs and Agilent, SNL developed an Automated Method Selection (AMS) software tool to predict the best liquid chromatography method for analysis of any new molecules of interest. Samples with novel pathways and new metabolite targets of interest are generated at a high pace in the ABF. Overall, the project advanced rapid metabolomics by combining liquid chromatography, ion mobility spectrometry, and data-independent mass spectrometry with machine learning. This multidimensional approach uses retention time, collision cross-section, precursor mass, and fragment-ion information to distinguish chemically similar metabolites that can be difficult to resolve using conventional liquid- or gas-chromatography methods. The resulting workflow also provided automated metabolite-identification error estimates, addressing a recognized need for statistical confidence measures in metabolomics.

Petzold, Christopher [Lawrence Berkeley National L↗
DOE OSTI2024

Rational Optimization of Microbial Processing for High Yield CO 2 -to-Isopropanol Conversion: Cooperative Research and Development Final Report, CRADA Number CRD-20-17114

This project focuses on the production of the fuel blendstock isopropanol using a CO 2 -fixing Clostridium by metabolic engineering and process optimizations. The project will initiate from a baseline isopropanol producer and pursue isopropanol production at high carbon-conversion efficiency. We will lead engineering work by in-depth pathway analyses including thermodynamics optimization, enzyme expense analysis, metabolic robustness analysis, and -omics analysis. The isopropanol production will be optimized via genome editing followed by fermentation optimizations. This project will deliver a novel microbial process that efficiently converts waste CO 2 to isopropanol at ~g/L titer level within 18-months. This project will layout a solid knowledge basis and technology platform for renewable CO 2 valorization to bio-blendstock that help achieve Co-Optima and Shell's goals.

09 BIOMASS FUELS↗
DOE OSTI2022

Size-Dependent Photon Avalanching in Tm 3+ Doped LiYF 4 Nano, Micro, and Bulk Crystals

Photon avalanche (PA) is a highly nonlinear mode of upconversion that is characterized by 100–1000-fold increase in luminescence intensity upon minute increments of pumping power. The practical realization of numerous possible nano-bio-technology applications utilizing the PA phenomenon will require information on its susceptibility to the material volume and surface. In this report these parameters are investigated via experimental and theoretical PA. The two-color, highly nonlinear PA emission at 475 and 800 nm is clearly observed in bulk single crystal, individual microcrystals, and ensembles of colloidal core and core–shell nanoparticles of LiYF 4 host doped with either 3 or 8% of thulium ions. The properties of PA emission, such as PA nonlinearity, PA gain, PA intensity, and luminescence kinetics in these materials show dependence on crystal volume and surface quenching. Theoretical simulations provide understanding of key physical processes that influence PA performance. Moreover, photon avalanche single beam super-resolution imaging is realized for the first time in 3% Tm 3+ doped LiYF 4 core–shell nanoparticles. The obtained insights and predictions form a solid background for further development and applications of new optimized PA materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Simultaneous quantification of uranium( VI ), samarium, nitric acid, and temperature with combined ensemble learning, laser fluorescence, and Raman scattering for real-time monitoring

In this work, laser-induced fluorescence spectroscopy (LIFS), Raman spectroscopy, and a stacked regression ensemble was developed for near real-time quantification of uranium(VI) (1–100 μg mL –1 ), samarium (0–200 μg mL –1 ) and nitric acid (0.1–4 M) with varying temperature (20 °C–45 °C). LIFS applications range from fundamental lab-scale studies to real-time process monitoring at industrial levels, such as nuclear reprocessing applications, provided the phenomena affecting the fluorescence spectrum are accounted for (e.g., absorption, quenching, complexation). Multiple chemometric models were examined and compared to a more traditional multivariate regression approach called partial least squares (PLS). Results obtained on synthetic samples selected using D-optimal experimental design indicated that a stacked regression method, which included ridge regression, random forest, PLS, and an eXtreme gradient boost algorithm, successfully measured uranium(VI) concentrations directly in nitric acid without measuring luminescence lifetimes or standard addition. The top model resulted in percent root-mean-square error of prediction values of 5.2, 1.9, 3.0, and 2.3% for U(VI), Sm 3+ , HNO 3 , and temperature, respectively. The approach may be useful for quantifying fluorescent fission products (e.g., Sm 3+ ) to provide information on burnup of irradiated nuclear fuel. This novel framework reinforces the applicability of LIFS for real-time applications in nuclear fuel cycle applications.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

Uncertainty-Aware Machine Learning for Small-Angle X-ray Scattering Analysis in Autonomous Experimentation

Small-angle X-ray scattering (SAXS) is a powerful high-throughput characterization tool for probing nanoscale structure in native sample environments, providing real-time morphological information such as nanoparticle size and shape during synthesis. However, automated SAXS data analysis for extracting meaningful structural parameters is non-trivial and remains a bottleneck in closed-loop experimentation towards autonomous materials discovery, which demands fast, reliable, and uncertainty-aware data analysis. Here, we develop a machine-learning approach for automated SAXS analysis tailored to closed-loop nanoparticle synthesis. A Random Forest (RF) regression model is trained on 100,000 synthetic SAXS curves generated from polydisperse spherical nanoparticles with realistic background contributions. Using normalized one-dimensional SAXS intensity profiles as input, the RF model directly predicts nanoparticle radius, size polydispersity, and background parameters, while the ensemble standard deviation across trees provides built-in uncertainty quantification (UQ). On synthetic data, we show that combining fit-quality metrics (R 2 , MAE) with thresholds on prediction uncertainty reliably identifies accurate parameter estimates without access to ground truth. We then apply the trained model to 365 experimental SAXS profiles of citrate-reduced gold nanoparticles synthesized using an automated droplet-flow microreactor with in situ SAXS at a synchrotron beamline, classifying the results into high- and low-confidence subsets based on UQ metrics. Finally, we integrate RF-based SAXS analysis into a simulated closed-loop optimization campaign using Gaussian process Bayesian optimization to minimize nanoparticle polydispersity, benchmarking against conventional automated Levenberg–Marquardt fitting. The RF-guided campaign exhibits substantially faster convergence and lower relative opportunity cost (∼0.07 vs ∼0.3), demonstrating that uncertainty-aware machine-learning SAXS analysis significantly enhances the efficiency and robustness of autonomous nanomaterials synthesis workflows.

Bayesian optimization↗
DOE OSTI2026

Active operator learning with predictive uncertainty quantification for partial differential equations

With the increased prevalence of neural operators being used to provide rapid solutions to partial differential equations (PDEs), understanding the accuracy of model predictions and the associated error levels is necessary for deploying reliable surrogate models in scientific applications. Existing uncertainty quantification (UQ) frameworks employ ensembles or Bayesian methods, which can incur substantial computational costs during both training and inference. Here, we propose a lightweight predictive UQ method tailored for Deep operator networks (DeepONets) that also generalizes to other operator networks. Numerical experiments on linear and nonlinear PDEs demonstrate that the framework’s uncertainty estimates are unbiased and provide accurate out-of-distribution uncertainty predictions with a sufficiently large training dataset. Our framework provides fast inference and uncertainty estimates that can efficiently drive outer-loop analyses that would be prohibitively expensive with conventional solvers. We demonstrate how predictive uncertainties can be used in the context of Bayesian optimization and active learning problems to yield improvements in accuracy and data-efficiency for outer-loop optimization procedures. In the active learning setup, we extend the framework to Fourier Neural Operators (FNO) and describe a generalized method for other operator networks. To enable real-time deployment, we introduce an inference strategy based on precomputed trunk outputs and a sparse placement matrix, reducing evaluation time by more than a factor of five. Our method provides a practical route to uncertainty-aware operator learning in time-sensitive settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Mapping the global distribution of C 4 vegetation using observations and optimality theory

Plants with the C 4 photosynthesis pathway typically respond to climate change differently from more common C 3 -type plants, due to their distinct anatomical and biochemical characteristics. These different responses are expected to drive changes in global C 4 and C 3 vegetation distributions. However, current C 4 vegetation distribution models may not predict this response as they do not capture multiple interacting factors and often lack observational constraints. Here, we used global observations of plant photosynthetic pathways, satellite remote sensing, and photosynthetic optimality theory to produce an observation-constrained global map of C 4 vegetation. We find that global C 4 vegetation coverage decreased from 17.7% to 17.1% of the land surface during 2001 to 2019. This was the net result of a reduction in C 4 natural grass cover due to elevated CO 2 favoring C 3 -type photosynthesis, and an increase in C 4 crop cover, mainly from corn (maize) expansion. Using an emergent constraint approach, we estimated that C 4 vegetation contributed 19.5% of global photosynthetic carbon assimilation, a value within the range of previous estimates (18–23%) but higher than the ensemble mean of dynamic global vegetation models (14 ± 13%; mean ± one standard deviation). Our study sheds insight on the critical and underappreciated role of C 4 plants in the contemporary global carbon cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Basal Plane Doping to Activate Colloidal MoS 2 Nanosheets for Catalytic Hydrodeoxygenation of para- Cresol

The valorization of biomass into biofuels is a critical process for producing renewable fuels. Hydrodeoxygenation (HDO), particularly over doped molybdenum disulfide (MoS 2 ), a transition metal dichalcogenide (TMD) material, is a common representative catalytic reaction system for converting biomass-derived materials into useful hydrocarbons. However, the location and role of dopants, such as Co, in HDO is not fully understood. The effects of dopant location and oxidation state are often precluded by inhomogeneity in the ensemble properties of nanosheet size and dopant dispersion, as well as difficulty in observing the behavior of atomic site behavior directly. Using a colloidal approach to synthesize cobalt-doped MoS 2 nanosheets with controlled dopant concentration, combined with X-ray absorption spectroscopy (XAS) and density functional theory (DFT) calculations, we determine that basal plane doped Co (25% Co:Mo mole ratio) shows peak catalytic activity in HDO of para-cresol, a model biomass-derived compound, and that basal Co sites are demonstrably more active than edge sites. By observing these doping effects in MoS 2 catalysts for HDO, we can further optimize not only the production of carbon-neutral fuels but also direct the tailoring of doped TMD catalysts toward their intended applications.

36 MATERIALS SCIENCE↗