Search NASA⌕ Search

SEARCH · Search NASA

Results for “Selection Bias”

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

Estimating irrigation water use from remotely sensed evapotranspiration data: Accuracy and uncertainties at field, water right, and regional scales

Irrigated agriculture is the dominant user of water globally, but most water withdrawals are not monitored or reported. As a result, it is largely unknown when, where, and how much water is used for irrigation. Here, we evaluated the ability of remotely sensed evapotranspiration (ET) data, integrated with other datasets, to calculate irrigation water withdrawals and applications in an intensively irrigated portion of the United States. We compared irrigation calculations based on an ensemble of satellite-driven ET models from OpenET with reported groundwater withdrawals from hundreds of farmer irrigation application records and a statewide flowmeter database at three spatial scales (field, water right group, and management area). At the field scale, we found that ET-based calculations of irrigation agreed best with reported irrigation when the OpenET ensemble mean was aggregated to the growing season timescale (bias = 1.6–4.9%, R 2 = 0.53–0.74), and agreement between calculated and reported irrigation was better for multi-year averages than for individual years. At the water right group scale, linking pumping wells to specific irrigated fields was the primary source of uncertainty. At the management area scale, calculated irrigation exhibited similar temporal patterns as flowmeter data but tended to be positively biased with more interannual variability. Disagreement between calculated and reported irrigation was strongly correlated with annual precipitation, and calculated and reported irrigation agreed more closely after statistically adjusting for annual precipitation. The selection of an ET model was also an important consideration, as variability across ET models was larger than the potential impacts of conservation measures employed in the region. From these results, we suggest key practices for working with ET-based irrigation data that include accurately accounting for changes in soil moisture, deep percolation, and runoff; careful verification of irrigated area and well-field linkages; and conducting application-specific evaluations of uncertainty.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

False Indications of Dose-Response Nonlinearity in Large Epidemiologic Cancer Radiation Cohort Studies; A Simulation Exercise

This study explores the likely prevalence of false indications of dose-response nonlinearity in large epidemiologic cancer radiation cohort studies (A-bomb survivors, INWORKS, Techa River). Reasons: Increasing numbers of tests of nonlinearity are being made in studies. Hypothesized nonlinear dose-response models have been justified to policy makers by analyses that rely in part on isolated findings that could be statistical fluctuations. After removing dose nonlinearity (linearization) by adjusting person-years of observation at each dose category, indications of nonlinearity, necessarily false, were counted in 5,000 randomized replications of six datasets. The average frequency of any false positive for five indicators of nonlinearity tested against a linear null was roughly 25% in Monte Carlo simulations per study, consistent with binomial calculations, increasing to ~50% within 6 studies assessed. Comparable frequencies were found using Akaike's information criterion (AIC) for model selection or multi-model averaging. False above-zero threshold doses were found more than 50% of the time, averaging to 0.05 Gy, consistent with findings in the 6 studies. Such bias, uncorrected, could distort meta-analyses of multiple studies, because meta-analyses can incorporate high P value findings. AIC-based correction for the extra threshold parameter lowered these false occurrences to 8 to 19%. Given the simulation rates, the possibility of false positives might be noted when isolated findings of nonlinearity are discussed in a regulatory context. When reporting a threshold dose with a P value > 0.05, it would be informative to note the expected high false prevalence rate due to bias.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

Sporophyte Stage Genes Exhibit Stronger Selection Than Gametophyte Stage Genes in Haplodiplontic Giant Kelp

Macrocystis pyrifera (giant kelp), a haplodiplontic brown macroalga that alternates between a macroscopic diploid (sporophyte) and a microscopic haploid (gametophyte) phase, provides an ideal system to investigate how ploidy background affects the evolutionary history of a gene. In M. pyrifera , the same genome is subjected to different selective pressures and environments as it alternates between haploid and diploid life stages. We assembled M. pyrifera gene models using available expression data and validated 8,292 genes models using the model alga Ectocarpus siliculosus . Differential expression analysis identified gene models expressed in either or both the haploid and diploid life stages while functional annotation identified processes enriched in each stage. Genes expressed preferentially or exclusively in the gametophyte stage were found to have higher nucleotide diversity (π = 2.3 × 10 –3 and 2.8 × 10 –3 , respectively) than those for sporophytes (π = 1.1 × 10 –3 and 1 × 10 –3 , respectively). While gametophyte-biased genes show faster sequence evolution, the sequence evolution exhibits less signatures of adaptations when compared to sporophyte-biased genes. Our findings contrast the standing masking hypothesis, which predicts higher standing genetic variation at the sporophyte stage, and support the strength of expression theory, which posits that genes expressed more strongly are expected to evolve slower. We argue that the sporophyte stage undergoes more stringent selection compared with the gametophyte stage, which carries a heavy genetic load associated with broadcast spawning. Furthermore, using whole-genome sequencing, we confirm the strong population structure in wild M. pyrifera populations previously established using microsatellite markers, and estimate population genetic parameters, such as pairwise genetic diversity and Tajima’s D , important for conservation and domestication of M. pyrifera .

Molano, Gary↗
DOE OSTI2021

Etching with electron beam-generated plasmas: Selectivity versus ion energy in silicon-based films

In the ideal case, plasma-enhanced atomic layer etching enables the ability to not only remove one monolayer of material but also leave adjacent layers undamaged. This dual mandate requires fine control over the flux of species to ensure efficacy, while maintaining an often arduously low ion energy. Electron beam-generated plasmas are well-suited for etching at low ion energies as they are generally characterized by highly charged particle densities (10 10 –10 11 cm –3 ) and low electron temperatures (<1.0 eV), which provide the ability to deliver a large flux of ions whose energies are <5 eV. Raising the ion energy with substrate biasing thus enables process control over an energy range that extends down to values commensurate with the bond strength of most material systems. In this work, we discuss silicon nitride etching using pulsed, electron beam-generated plasmas produced in argon-SF6 backgrounds. We pay particular attention to the etch rates and selectivity versus oxidized silicon nitride and polycrystalline silicon as a function of ion energy from a few eV up to 50 eV. Here, we find the blanket etch rate of Si 3 N 4 to be in the range of 1 A/s, with selectivities (versus SiO 2 and poly-Si) greater than 10:1 when ion energies are below 30 eV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Ligand removal of Au 25 nanoclusters by thermal and electrochemical treatments for selective CO 2 electroreduction to CO

Undercoordinated metal nanoclusters have shown great promise for various catalytic applications. However, their activity is often limited by the covalently bonded ligands, which could block the active surface sites. Here, we investigate the ligand removal process for Au 25 nanoclusters using both thermal and electrochemical treatments, as well as its impact on the electroreduction of CO 2 to CO. The Au 25 nanoclusters are synthesized with 2-phenylethanethiol as the capping agent and anchored on sulfur-doped graphene. The thiolate ligands can be readily removed under either thermal annealing at ≥180°C or electrochemical biasing at ≤-0.5 V vs reversible hydrogen electrode, as evidenced by the Cu underpotential deposition surface area measurement, x-ray photoelectron spectroscopy, and extended x-ray absorption fine structure spectroscopy. However, these ligand-removing treatments also trigger the structural evolution of Au 25 nanoclusters concomitantly. The thermally and electrochemically treated Au 25 nanoclusters show enhanced activity and selectivity for the electrochemical CO 2 -to-CO conversion than their pristine counterpart, which is attributed to the exposure of undercoordinated Au sites on the surface after ligand removal. This work provides facile strategies to strip away the staple ligands from metal nanoclusters and highlights its importance in promoting the catalytic performances.

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

Dark Energy Survey Year 3 results: Optimizing the lens sample in a combined galaxy clustering and galaxy-galaxy lensing analysis

We investigate potential gains in cosmological constraints from the combination of galaxy clustering and galaxy-galaxy lensing by optimizing the lens galaxy sample selection using information from Dark Energy Survey (DES) Year 3 data and assuming the DES Year 1 metacalibration sample for the sources. We explore easily reproducible selections based on magnitude cuts in i-band as a function of (photometric) redshift, zphot, and benchmark the potential gains against those using the well-established redMaGiC [E. Rozo et al., Mon. Not. R. Astron. Soc. 461, 1431 (2016)MNRAA40035-871110.1093/mnras/stw1281] sample. We focus on the balance between density and photometric redshift accuracy, while marginalizing over a realistic set of cosmological and systematic parameters. Our optimal selection, the MagLim sample, satisfies i<4zphot+18 and has ∼30% wider redshift distributions but ∼3.5 times more galaxies than redMaGiC. Assuming a wCDM model (i.e. with a free parameter for the dark energy equation of state) and equivalent scale cuts to mitigate nonlinear effects, this leads to 40% increase in the figure of merit for the pair combinations of Ωm, w, and σ8, and gains of 16% in σ8, 10% in Ωm, and 12% in w. Similarly, in ΛCDM, we find an improvement of 19% and 27% on σ8 and Ωm, respectively. We also explore flux-limited samples with a flat magnitude cut finding that the optimal selection, i<22.2, has ∼7 times more galaxies and ∼20% wider redshift distributions compared to MagLim, but slightly worse constraints. We show that our results are robust with respect to the assumed galaxy bias and photometric redshift uncertainties with only moderate further gains from increased number of tomographic bins or the inclusion of bin cross-correlations, except in the case of the flux-limited sample, for which these gains are more significant.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Deep Learning Based Frequency Stability Assessment in Power Grid with High Renewables

Frequency stability assessment is one critical aspect of power system security assessment. Traditional N-1 screening method is based on the simulations of a few typical daily and seasonal operation scenarios. However, the increasing integration of inverter-based renewables and the retirement of conventional synchronous generators result in decreasing system inertia and growing complexity of system operating conditions. Selecting a few typical operation scenarios cannot cover all operating conditions, and the time-domain simulation of all operation conditions requires tremendous time. This paper proposes a more efficient frequency stability assessment method based on deep learning. The affinity propagation clustering algorithm is used to divide the dataset into different clusters, so the selected dataset for training can cover the diversified operating conditions as much as possible. Also, feature normalization is applied to both the training dataset and testing dataset in order to remove any unnecessary bias. Especially, trained model based on full dataset normalization has bounded error in the prediction. The case study on the reduced 240-bus WECC system demonstrates that the proposed method can predict accurate frequency nadir with limited training dataset. The deep learning model using the revised feature normalization can predict more accurate frequency nadir than that using the traditional feature normalization and has very small maximum prediction error.

affinity propagation↗
DOE OSTI2023

Inferring the Energy and Distance Distributions of Fast Radio Bursts Using the First CHIME/FRB Catalog

Abstract Fast radio bursts (FRBs) are brief, energetic, typically extragalactic flashes of radio emission whose progenitors are largely unknown. Although studying the FRB population is essential for understanding how these astrophysical phenomena occur, such studies have been difficult to conduct without large numbers of FRBs and characterizable observational biases. Using the recently released catalog of 536 FRBs published by the Canadian Hydrogen Intensity Mapping Experiment/Fast Radio Burst (CHIME/FRB) collaboration, we present a study of the FRB population that also calibrates for selection effects. Assuming a Schechter function, we infer a characteristic energy cut-off of E char = 2.38 − 1.64 + 5.35 × 10 41 erg and a differential power-law index of γ = − 1.3 − 0.4 + 0.7 . Simultaneously, we infer a volumetric rate of [ 7.3 − 3.8 + 8.8 (stat.) − 1.8 + 2.0 ( sys . ) ] × 10 4 Gpc −3 yr −1 above a pivot energy of 10 39 erg and below a scattering timescale of 10 ms at 600 MHz, and find we cannot significantly constrain the cosmic evolution of the FRB population with star-formation rate. Modeling the host’s dispersion measure (DM) contribution as a log-normal distribution and assuming a total Galactic contribution of 80 pc cm −3 , we find a median value of DM host = 84 − 49 + 69 pc cm −3 , comparable with values typically used in the literature. Proposed models for FRB progenitors should be consistent with the energetics and abundances of the full FRB population predicted by our results. Finally, we infer the redshift distribution of FRBs detected with CHIME, which will be tested with the localizations and redshifts enabled by the upcoming CHIME/FRB Outriggers project.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Rapid Evaluation Framework for the CMIP7 Assessment Fast Track

As Earth system models (ESMs) grow in complexity and in volume of output data, there is an increasing need for rapid, comprehensive evaluation of their scientific performance. The upcoming Assessment Fast Track for the Seventh Phase of the Coupled Model Intercomparison Project (CMIP7) will require expeditious response for model analyses designed to inform and drive integrated Earth system assessments. To meet this challenge, the Rapid Evaluation Framework (REF), a community-driven platform for benchmarking and performance assessment of ESMs, was designed and developed. The initial implementation of the REF, constructed to meet the near-term needs of the CMIP7 Assessment Fast Track, builds upon four disparate community evaluation and benchmarking tools that are coupled together using the Coordinated Model Evaluation Capabilities (CMEC) framework. The REF runs within a containerized workflow for portability and reproducibility and is aimed at generating and organizing diagnostics covering a variety of model variables. The REF leverages well documented observational datasets to provide assessments of model fidelity across a collection of diagnostics. All diagnostics were identified and selected with community involvement and consultation. Operational integration with the Earth System Grid Federation (ESGF) will permit automated execution of the REF for selected diagnostics as soon as model output data are published on ESGF by the originating modeling centers. The REF is designed to be portable across a range of current computational platforms to facilitate use by modeling centers for assessing the evolution of model versions or gauging the relative performance of CMIP simulations before being published on ESGF. When integrated into production simulation workflows, results from the REF provide immediate quantitative feedback that allows model developers and scientists to quickly identify model biases and performance issues. After the REF is released to the community, its subsequent development and support will be prioritized by an international consortium of scientists and engineers, enabling a broader impact across Earth science disciplines. For instance, the REF will facilitate improvements to models and will enhance confidence in model projections through process-based selection of models based on their performance with respect to observations. Production of reproducible diagnostics and community-based assessments are key features of the REF. Furthermore, providing interoperability with existing evaluation packages assures that contributions from previous community efforts will be available for use in future model intercomparison projects.

Hoffman, Forrest [ORNL] (ORCID:0000000158024134)↗
DOE OSTI2021

Planet Occurrence Rate Correlated to Stellar Dynamical History: Evidence from Kepler and Gaia

The dynamical history of stars influences the formation and evolution of planets significantly. To explore the influence of dynamical history on the planet formation and evolution using observations, we assume stars that experienced significantly different dynamical histories tend to have different relative velocities. Utilizing the accurate Gaia–Kepler Stellar Properties Catalog, we select single main-sequence stars and divide these stars into three groups according to their relative velocities, i.e., high-V, medium-V, and low-V stars. After considering the known biases from Kepler data and adopting prior and posterior correction to minimize the influence of stellar properties on planet occurrence rate, we find that high-V stars have a lower occurrence rate of super-Earths and sub-Neptunes (1–4 R {sub ⊕}, P < 100 days) and a higher occurrence rate of sub-Earth (0.5–1 R {sub ⊕}, P < 30 days) than low-V stars. Additionally, high-V stars have a lower occurrence rate of hot Jupiter-sized planets (4–20 R {sub ⊕}, P < 10 days) and a slightly higher occurrence rate of warm or cold Jupiter-sized planets (4–20 R {sub ⊕}, 10 < P < 400 days). After investigating multiplicity and eccentricity, we find that high-V planet hosts prefer a higher fraction of multiplanet systems and lower average eccentricity, which are consistent with the eccentricity–multiplicity dichotomy of Kepler planetary systems. All of these statistical results favor the scenario that high-V stars with large relative velocity may experience fewer gravitational events, while low-V stars may be influenced by stellar clustering significantly.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Uncovering hidden bias in neutron diffraction residual strain measurements

When calculating residual strain via neutron or X-ray diffraction, uncertainties propagated from the peak fit are often inadequate to describe the true scatter of measurements about a singular strain state, such as one that should describe a macroscopic continuum. Because diffraction is inherently a selective process, orientation-dependent scatter arises from the sub-sampling of strong microstructure and strain gradients. This paper investigates the appropriateness of propagated uncertainties with reference to their original intention, i.e. noise about a mean value. Thirty-six unique orientations of strain measurements are taken at multiple locations within an additive friction-stir deposition component with fine-scale gradients (∼200 µm) of plastic strain, texture and residual elastic strain. Multiple strain and stress calculation pathways are compared: direct substitution of three measurements into Hooke's law, direct inversion of any six unique orientations into the strain state tensor and thirty-six measurement least-squares estimation. For the last two cases, the appropriateness of the uncertainty interval is statistically evaluated on the basis of a physical constraint: common agreement under the strain transformation law. For this sample, the direct inversion of six measurements retains a conservative estimate of the uncertainty. However, propagated uncertainties in the least-squares solution greatly underestimate the true experimental scatter. A simple pathway to estimate appropriate uncertainty intervals is suggested. These results demonstrate that the interpretation of uncertainty in residual strain is strongly dependent on intrinsic sample-dependent effects, and that oversampling orientations and statistical analysis can give more accurate results with realistic uncertainties.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Systematic Benchmarking of Climate Models: Methodologies, Applications, and New Directions

As climate models become increasingly complex, there is a growing need to comprehensively and systematically assess model performance with respect to observations. Given the increasing number and diversity of climate model simulations in use, the community has moved beyond simple model intercomparison and toward developing methods capable of benchmarking a large number of simulations against a suite of climate metrics. Here, we present a detailed review of evaluation and benchmarking methods and approaches developed in the last decade, focusing primarily on scientific implications for Coupled Model Intercomparison Project (CMIP) simulations and CMIP6 results that contributed to the Intergovernmental Panel on Climate Change (IPCC) Sixth Assessment Report (AR6). Based on this review, we explain the resulting contemporary philosophy of model benchmarking, and provide clear distinctions and definitions of the terms model verification, process validation, evaluation, and benchmarking. While significant progress has been made in model development based on systematic evaluation and benchmarking efforts, some climate system biases still remain. The development of open‐source community software packages has played a fundamental role in identifying areas of significant model improvement and bias reduction. We review the key features of several software packages that have been commonly used over the past decade to evaluate and benchmark global and regional climate models. Additionally, we discuss best practices for the selection of evaluation and benchmarking metrics and for interpreting the obtained results, the importance of selecting suitable sources of reference data and accurate uncertainty quantification.

Environmental sciences↗
DOE OSTI2021

Direct determination of cellulosic glucan content in starch-containing samples

A simple and highly selective analytical procedure is presented for the determination of cellulosic glucan content in samples that contain both cellulose and starch. This method eliminates the unacceptably large compounding errors of current two-measurement methods. If both starch and cellulose are present before analytical hydrolysis, both will be hydrolyzed to glucose causing bias and inaccuracy in the method. To prevent this interference, the removal of starch prior to cellulosic quantification is crucial. The method presented here is a concise in-series procedure with minimal measurements, eliminating large compounding errors. Sample preparation consists of a starch extraction employing enzymatic hydrolysis followed by a simple filtration and wash. The samples are then subjected to a two-stage acid hydrolysis. The concentration of glucose is determined by ion exchange high-performance liquid chromatography with a Pb 2+ column and a refractive index detector. The cellulosic glucan content is calculated based on the initial dry weight of the starting material. Data for the native biomass materials studied show excellent reproducibility, with coefficients of variance of 3.0% or less associated with the method. This selectivity for cellulosic glucan by the procedure was validated with several analytical techniques such as liquid chromatography coupled with mass spectrometry (LC–MS), Raman spectroscopy, and nuclear magnetic resonance.

09 BIOMASS FUELS↗
DOE OSTI2023

Discovery of a Hybrid System for Photocatalytic CO 2 Reduction via Attachment of a Molecular Cobalt-Quaterpyridine Complex to a Crystalline Carbon Nitride

While recent reports have demonstrated the attachment of molecular catalysts to amorphous, graphitic carbon nitrides (g-CN) for light-driven CO 2 reduction, approaches to the utilization of crystalline carbon nitrides have remained undiscovered. Herein, a functional hybrid photocatalyst system has been found using a crystalline carbon nitride semiconductor, poly(triazine imide) lithium chloride (PTI-LiCl), with a surface-attached CoCl 2 (qpy-Ph-COOH) catalyst for CO 2 reduction. The molecular catalyst attaches to PTI-LiCl at concentrations from 0.10 to 4.30 wt % and exhibits ∼96% selectivity for CO production in a CO 2 -saturated, aqueous 0.5 M KHCO 3 solution. Optimal loadings were found to be within 0.42–1.04 wt % with rates between 1,400 and 1,550 μmol CO/g·h at an irradiance of 172 mW/cm 2 (λ = 390 nm) and apparent quantum yields of ∼2%. This optimized loading is postulated to represent a balance between maximal turnover frequency (TOF; 300+ h –1 ) and excess catalyst that can limit excited-electron lifetimes, as probed via transient absorption spectroscopy. An increase in the incident irradiance yields a concomitant increase in the TOFs and CO rates only for the higher catalyst loadings, reaching up to 2,149 μmol CO/g·h with a more efficient use of the catalyst surface capacity. The lower catalyst loadings, by comparison, already function at maximal TOFs. Higher surface loadings are also found to help mitigate deactivation of the molecular catalysts during extended catalytic testing (>24 h) owing to the greater net surface capacity for CO 2 reduction, thus representing an effective strategy to extend lifetime. The hybrid particles can be deposited onto an FTO substrate to yield ∼60% Faradaic efficiency for photoelectrochemical CO production at −1.2 V vs Ag/AgCl bias. In conclusion, these results demonstrate the synergistic combination of a crystalline carbon nitride with a molecular catalyst that achieves among the highest known rates in carbon-nitride systems for the light-driven CO 2 reduction to CO in aqueous solution with >95% selectivity.

CO2 reduction↗
DOE OSTI2025

An Antibody Cocktail-Based Immunoaffinity-LC-MS Method Enabled Ultra-Sensitive and Robust Quantification of Circulating Proinsulin Proteoforms and C-Peptide

Accurately measuring circulating proinsulin proteoforms is crucial for clinical investigation of diabetes, but was previously not feasible owing to limited assay specificity/sensitivity. Here, in this study, we devised a highly sensitive LC-MS-based strategy to quantify intact proinsulin, des-31,32 and des-64,65 proinsulin, and C-peptide in circulation. The method involves: (i) quantitative, robust affinity capture using an optimized antibody cocktail, eliminating the severe quantitative bias across multiple proteoforms typically introduced when using a single antibody; (ii) Lys-C digestion producing unique signature peptides for each proteoform, and (iii) trapping-nano-LC coupled with FAIMS/dCV-MS for an ultrasensitive analysis. The selective trapping/delivery ensured sensitive/selective analysis of the targets while achieving excellent analytical robustness that is critical for clinical assays, and the FAIMS/dCV substantially reduces baseline noise/interferences, further enhancing S/N. The assay achieved exceptional sensitivity, with serum LOQs of 1.7, 2.3, and 3.6 pg/mL respectively for intact-proinsulin, des-31,32 and des-64,65, representing the first assay capable of sensitively quantifying these major circulating proinsulin proteoforms. We applied this assay to 78 subjects, including autoantibody positive (n = 20) and new-onset type 1 diabetes (T1D, n = 19) with respective age/sex/BMI-matched controls, enabling the first accurate profiling of proinsulin proteoforms in clinical groups. The assay results demonstrated a clear separation of control and new-onset T1D groups that a parallel total-proinsulin ELISA assay fails to capture. Furthermore, distinct expression patterns in relative abundance ratios among proteoforms were observed across clinical groups. This assay may provide valuable insights into the β-cell functions and the onset/progression of diabetes and other associated conditions. Moreover, the strategy is broadly applicable to targeted measurement of other biomarker proteoforms.

Shen, Qingqing [University at Buffalo, NY (United ↗
DOE OSTI2026

Hierarchical Multi-agent Large Language Model Reasoning for Autonomous Heterogeneous Catalyst Discovery

Artificial intelligence is reshaping scientific exploration, but most methods automate procedural tasks without engaging in scientific reasoning, limiting autonomy in discovery. We demonstrate that hierarchical agentic large language model reasoning can efficiently drive simulation and scientific exploration. Across two chemical applications, CO adsorption on Cu surface transition metal adatoms and on M–N–C catalysts, reasoning-guided exploration reduces required atomistic simulations by up to 90% relative to heuristic or random selection. Comparisons across single-agent, multi-agent, and stochastic baselines show that hierarchical strategies yield more coherent and information-efficient search trajectories. Reasoning traces reveal chemically grounded decisions that cannot be explained by semantic bias or stochastic sampling. We realize these agentic reasoning strategies in Materials Agents for Simulation and Theory in Electronic-structure Reasoning (MASTER), a multimodal system that translates natural language into density functional theory workflows. Altogether, multi-agent collaboration accelerates heterogeneous catalyst discovery and marks a step toward more autonomous, reasoning-guided scientific exploration.

30 DIRECT ENERGY CONVERSION↗