SEARCH · Engineering Papers

Results for “parametric uncertainty”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Bayesian Spatial Models for Projecting Corn Yields

Climate change is predicted to impact corn yields. Previous studies analyzing these impacts differ in data and modeling approaches and, consequently, corn yield projections. We analyze the impacts of climate change on corn yields using two statistical models with different approaches for dealing with county-level effects. The first model, which is novel to modeling corn yields, uses a computationally efficient spatial basis function approach. We use a Bayesian framework to incorporate both parametric and climate model structural uncertainty. We find that the statistical models have similar predictive abilities, but the spatial basis function model is faster and hence potentially a useful tool for crop yield projections. We also explore how different gridded temperature datasets affect the statistical model fit and performance. Compared to the dataset with only weather station data, we find that the dataset composed of satellite and weather station data results in a model with a magnified relationship between temperature and corn yields. For all statistical models, we observe a relationship between temperature and corn yields that is broadly similar to previous studies. We use downscaled and bias-corrected CMIP5 climate model projections to obtain detrended corn yield projections for 2020–2049 and 2069–2098. In both periods, we project a decrease in the mean corn yield production, reinforcing the findings of other studies. However, the magnitude of the decrease and the associated uncertainties we obtain differ from previous studies.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2016

Multi-Model Assessment of Global Hydropower and Cooling Water Discharge Potential Under Climate Change

Worldwide, 98% of total electricity is currently produced by thermoelectric power and hydropower. Climate change is expected to directly impact electricity supply, in terms of both water availability for hydropower generation and cooling water usage for thermoelectric power. Improved understanding of how climate change may impact the availability and temperature of water resources is therefore of major importance. Here we use a multi-model ensemble to show the potential impacts of climate change on global hydropower and cooling water discharge potential. For the first time, combined projections of streamflow and water temperature were produced with three global hydrological models (GHMs) to account for uncertainties in the structure and parametrization of these GHMs in both water availability and water temperature. The GHMs were forced with bias-corrected output of five general circulation models (GCMs) for both the lowest and highest representative concentration pathways (RCP2.6 and RCP8.5). The ensemble projections of streamflow and water temperature were then used to quantify impacts on gross hydropower potential and cooling water discharge capacity of rivers worldwide. We show that global gross hydropower potential is expected to increase between +2.4% (GCM-GHM ensemble mean for RCP 2.6) and +6.3% (RCP 8.5) for the 2080s compared to 1971-2000. The strongest increases in hydropower potential are expected for Central Africa, India, central Asia and the northern high-latitudes, with 18-33% of the world population living in these areas by the 2080s. Global mean cooling water discharge capacity is projected to decrease by 4.5-15% (2080s). The largest reductions are found for the United States, Europe, eastern Asia, and southern parts of South America, Africa and Australia, where strong water temperature increases are projected combined with reductions in mean annual streamflow. These regions are expected to affect 11-14% (for RCP2.6 and the shared socioeconomic pathway (SSP)1, SSP2, SSP4) and 41-51% (RCP8.5-SSP3, SSP5) of the world population by the 2080s.

global hydrological models↗

NASA NTRSDate not supplied

Regional and Model-Specific Response Types in A Global Gridded Crop Model Ensemble

Crop models are often employed to project crop yields under changing conditions such as global warming and associated management change for adaptation. Multi-model ensembles are promoted to enhance the robustness of projections, but questions remain on what causes often large differences between projections of individual models. Global Gridded Crop Models (GGCMs) are especially exposed to this question when applied for assessing climate change impacts, adaptation, environmental impacts of agricultural production, because their results are used in downstream analyses, such as in integrated assessment or economic modeling for projecting future land-use change. Even though global gridded crop models are often based on detailed field-scale models or have implemented similar modeling principles in other ecosystem models, global-scale models are subject to substantial uncertainties from both model structure and parametrization as well as from calibration and input data quality. AgMIP’s Global Gridded Crop Model Intercomparison (GGCMI) has thus set out to intercompare GGCMs in order to evaluate model performance, describe model uncertainties, identify inconsistencies within the ensemble and underlying reasons, and to ultimately improve models and modeling capacities. In phase 2 of the GGCMI activities, 12 modeling groups followed a modeling protocol that asked for up to 1404 31-year global simulations at 0.5 arc-degree spatial resolution to assess models’ sensitivities to changes in carbon dioxide (C; 4 different levels) temperature (T; 7 different offset levels), water supply (W; 9 levels), and nitrogen (N; 3 levels), the so-called CTWN experiment (Franke et al. 2020; http://dx.doi.org/10.5194/gmd-13-2315-2020). We here present analyses of model response types using impact response surfaces along the C, T, W, and N dimensions, respectively and collectively. Doing so, we can understand differences in simulated responses per driver rather than aggregated changes in yields. We find that models’ sensitivities to the individual driver dimensions are substantially different and often more different across models than across regions. A cluster analysis finds regional and model-specific patterns. There is some agreement across models with respect to the spatial patterns of response types but strong differences in the distribution of response type clusters across models suggests that models need to undergo further scrutiny. We suggest establishing standards in model process evaluation not only against historical dynamics but also against dedicated experiments across the CTWN dimensions.

crop models↗

DOE OSTI2024

Ensemble Kalman inversion of induced polarization data

SUMMARY This paper explores the applicability of ensemble Kalman inversion (EKI) with level-set parametrization for solving geophysical inverse problems. In particular, we focus on its extension to induced polarization (IP) data with uncertainty quantification. IP data may provide rich information on characteristics of geological materials due to its sensitivity to characteristics of the pore–grain interface. In many IP studies, different geological units are juxtaposed and the goal is to delineate these units and obtain estimates of unit properties with uncertainty bounds. Conventional inversion of IP data does not resolve well sharp interfaces and tends to reduce and smooth resistivity variations, while not readily providing uncertainty estimates. Recently, it has been shown for DC resistivity that EKI is an efficient solver for inverse problems which provides uncertainty quantification, and its combination with level set parametrization can delineate arbitrary interfaces well. In this contribution, we demonstrate the extension of EKI to IP data using a sequential approach, where the mean field obtained from DC resistivity inversion is used as input for a separate phase angle inversion. We illustrate our workflow using a series of synthetic and field examples. Variations with uncertainty bounds in both DC resistivity and phase angles are recovered by EKI, which provides useful information for hydrogeological site characterization. Although phase angles are less well-resolved than DC resistivity, partly due to their smaller range and higher percentage data errors, it complements DC resistivity for site characterization. Overall, EKI with level set parametrization provides a practical approach forward for efficient hydrogeophysical imaging under uncertainty.

Geochemistry & Geophysics↗

NASA NTRS2002

Comparison of Two Methods for Estimating the Sampling-Related Uncertainty of Satellite Rainfall Averages Based on a Large Radar Data Set

The uncertainty of rainfall estimated from averages of discrete samples collected by a satellite is assessed using a multi-year radar data set covering a large portion of the United States. The sampling-related uncertainty of rainfall estimates is evaluated for all combinations of 100 km, 200 km, and 500 km space domains, 1 day, 5 day, and 30 day rainfall accumulations, and regular sampling time intervals of 1 h, 3 h, 6 h, 8 h, and 12 h. These extensive analyses are combined to characterize the sampling uncertainty as a function of space and time domain, sampling frequency, and rainfall characteristics by means of a simple scaling law. Moreover, it is shown that both parametric and non-parametric statistical techniques of estimating the sampling uncertainty produce comparable results. Sampling uncertainty estimates, however, do depend on the choice of technique for obtaining them. They can also vary considerably from case to case, reflecting the great variability of natural rainfall, and should therefore be expressed in probabilistic terms. Rainfall calibration errors are shown to affect comparison of results obtained by studies based on data from different climate regions and/or observation platforms.

Lau, William K. M.↗

NASA NTRS2024

Sizing and Performance Analysis of a MW-Class Electrified Aircraft Propulsion (EAP) System for a Parallel Hybrid Turboprop Concept

Aircraft electrification has emerged as a pivotal research and development focus in the 21st century, gaining international momentum as a strategic approach to curtail fuel consumption and emissions in the civil aviation sector. Through the Electrified Powertrain Flight Demonstration (EPFD) project, NASA is collaborating with industry partners to develop and demonstrate Megawatt (MW) class Electrified Aircraft Propulsion (EAP) systems using regional turboprops, as well as single-aisle commercial transports. With the rapid emergence of novel EAP concepts, current efforts are focused on developing the capabilities to evaluate the impacts of electrification on vehicle-level performance. This requires parametric modeling at the conceptual design stage, using uncertainty propagation techniques to account for the high variability in EAP system architecture. This paper details an integrated approach to parametric sizing of multi-MW EAP systems in a large turboprop freighter concept with a parallel hybrid configuration. The primary goal is to assess vehicle-level performance sensitivities, such as range and fuel burn reduction, to variations in electrical system component-level performance metrics, including specific power, efficiency, and energy density, based on near-term EAP technology levels.

electrified aircraft propulsion↗

DOE OSTI2023

Linear instabilities in the hot-ion regime in a high-field spherical tokamak

Abstract The present operating high-field compact spherical tokamak ST40 is an important step toward an ST-based fusion reactor (Gryaznevich and Asunta 2017 Fusion Eng. Des. 123 177; McNamara 2023 Nucl. Fusion 63 054002). Temperature and density profiles and their uncertainties in recent hot ion ST40 plasmas with central ion temperatures exceeding 8.6 keV, have been inferred using an integrated analysis of several diagnostics including line-of-sight integration and volume average measurements, as well as limited profile information from a charge-exchange-recombination spectrometer. A linear gyrokinetic stability analysis has been carried out to identify the most unstable micro-instabilities in these hot ion plasmas. In one of these plasmas, it is found that linear growth rates of both ion- and electron-scale ( k θ ρ s ⩾ 0.2 where k θ is the poloidal wavenumber and ρ s is the ion gyro-radius with sound speed) modes decrease from the edge toward the core of the plasma (i.e. from ρ = 0.8 to 0.3 where ρ is the square root of the normalized magnetic flux), while the change of linear growth rates at k θ ρ s < 0.2 is non-monotonic. In particular, at ρ = 0.3 no unstable mode was found at k θ ρ s > 5 , and about an order of magnitude reduction in the maximum linear growth rate in the ion-scale wavenumber range of 0.2 ⩽ k θ ρ s < 1 is seen at ρ = 0.3 compared with the ρ = 0.8 location. It is found through parametric scans that while the ion temperature gradient (ITG) mode/trapped electron mode (TEM) and kinetic ballooning mode (KBM)/kinetic shear Alfvén (KSA) mode are more important in the plasma core, ubiquitous mode (UM) is found to be the important ion-scale instability at ρ ⩾ 0.6 . At ρ = 0.8, two instabilities are found at the electron scale, with one being the typical electrostatic electron temperature gradient (ETG) mode at ρ e scale and beyond and the other being an unidentified low-real-frequency instability (at the intermediate scale between electron and ion gyroradii) which can propagate in the ion direction and is driven by ETG. The existence of UM at the ion scale and the low-real-frequency instability and ETG mode at the electron scale at ρ = 0.8 is supported by an extensive linear gyrokinetic stability analysis of another ST40 hot ion plasma in a wider parametric range, taking into account the uncertainties in the experimental profiles.

ubiquitous mode↗

NASA NTRS2003

Modeling and Visualizing Uncertainty in Continuous Variables Predicted using Remotely Sensed Data

The use of remotely sensed images to map continuous biophysical variables, such as those related to terrestrial vegetation amount, sea surface temperature, and many other targets of NASA s Earth Observing System (EOS), includes variable, parametric, positional, spatial support and structural sources of uncertainty. A complete description of uncertainty will lead to a probability distribution at each location, allowing the exploration of the spatial dimension of uncertainty, that is, where the field is not well quantified. To achieve this purpose, convenient visualization tools are required. We have produced such a tool, called PDFVis, that facilitates the display of probability density functions (pdfs) on a per-grid-cell basis. The density estimate from Monte-Carlo generated realizations is interactively displayed as well as parametric and non-parametric summaries of the pdf field (such as mean, median, quartiles, standard deviation, number of modes, and locations of modes). Shaded surface renderings of pdfs along a transect can also be projected onto a plane. This tool will become more useful as richer descriptions of spatial uncertainty become available.

Dungan, Jennifer L.↗

DOE OSTI2022

Bimolecular Peroxy Radical (RO 2 ) Reactions and Their Relevance in Radical Initiated Oxidation of Hydrocarbons

The kinetics of peroxy radical (RO 2 ) reactions have been of long-standing interest in atmospheric and combustion chemistry. Nevertheless, the lack of kinetic studies at higher temperatures for their reactions with other radicals such as OH has precluded the inclusion of this class of reactions in detailed kinetics models developed for combustion applications. In this work, guided by the limited room temperature experimental studies on selected alkyl-peroxy radicals and literature theoretical kinetics on the prototypical CH 3 O 2 + OH system, we have performed parametric studies on the effect of uncertainties in the rate coefficients and branching ratios to potential product channels for RO 2 + OH reactions at higher temperatures. Literature kinetics models were used to simulate autoignition delays, laminar flame speeds, and speciation profiles in flow and stirred reactors for a variety of common combustion-relevant fuels. Inclusion of RO 2 + OH reactions was found to retard autoignition in fuel-lean (φ = 0.5) mixtures of ethane and dimethyl ether in air. The observed effects were noticeably more pronounced in ozone-enriched combustion of ethane and dimethyl ether. The simulations also examined the influence of ozone doping levels, pressures, and equivalence ratios for both ethane and dimethyl ether oxidation. Sensitivity and flux analyses revealed that the RO 2 + OH reaction is a significant sink of RO 2 radicals at the early stage of autoignition, affecting fuel oxidation through RO 2 ↔ QOOH, RO 2 ↔ alkene + HO 2 , or RO 2 + HO 2 ↔ ROOH + O 2 . Additionally, the kinetic stability of the trioxide formed from RO 2 + OH reactions was investigated using master equation analyses. Last, we discuss other bimolecular reactions that are missing in literature kinetics models but are relevant to hydrocarbon oxidation initiated by external radical sources (plasma-enhanced, ozone-enriched combustion, etc.). In conclusion, the present simulations provide a strong motivation for better characterizing the bimolecular kinetics of peroxy radicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Uncertainty propagation in pore water chemical composition calculation using surrogate models

Performance assessment in deep geological nuclear waste repository systems necessitates an extended knowledge of the pore water chemical conditions prevailing in host-rock formations. In the last two decades, important progress has been made in the experimental characterization and thermodynamic modeling of pore water speciation, but the influence of experimental artifacts and uncertainties of thermodynamic input parameters are seldom evaluated. In this respect, we conducted an uncertainty propagation study in a reference geochemical model describing the pore water chemistry of the Callovian-Oxfordian clay formation. Nineteen model input parameters were perturbed, including those associated to experimental characterization (leached anions, exchanged cations, cation exchange selectivity coefficients) and those associated to generic thermodynamic databases (solubilities). A set of 13 quantities of interest were studied by the use of polynomial chaos expansions built non-intrusively with a least-squares forward stepwise regression approach. Training and validation sets of simulations were carried out using the geochemical speciation code PHREEQC. The statistical results explored the marginal distribution of each quantity of interest, their bivariate correlations as well as their global sensitivity indices. The influence of the assumed distributions for input parameters uncertainties was evaluated by considering two parametric domain sizes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Development of a compressive failure model for carbon fiber composites and associated uncertainties.

An approach to increase the value of carbon fiber for wind turbines blades, and other compressive strength driven designs, is to identify pathways to increase its cost-specific compressive strength. A finite element model has been developed to evaluate the predictiveness of current finite element methods and to lay groundwork for future studies that focus on improving the cost-specific compressive strength. Parametric studies are conducted to understand which uncertainties in the model inputs have the greatest impact on compressive strength predictions. Furthermore, a statistical approach is also presented that enables the micromechanical model, which is deterministic, to efficiently account for statistical variability in the fiber misalignment present in composite materials; especially if the results from the hexagonal and square pack models are averaged. The model was found to agree well with experimental results for a Zoltek PX-35 pultrusion. The sensitivity studies suggest that the fiber packing and the interface shear strength have the greatest impact on compressive strength prediction for the fiber reinforced polymer studied here. Based on the performance of the modeling approach presented in this work, it is deemed sufficient for future work which will seek to identify carbon fiber composites with improved cost-specific compressive strength.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Parametric dynamic mode decomposition for reduced order modeling

Dynamic Mode Decomposition (DMD) is a model-order reduction approach, whereby spatial modes of fixed temporal frequencies are extracted from numerical or experimental data sets. The DMD low-rank or reduced operator is typically obtained by singular value decomposition of the temporal data sets. For parameter-dependent models, as found in many multi-query applications such as uncertainty quantification or design optimization, the only parametric DMD technique developed was a stacked approach, with data sets at multiple parameter values were aggregated together, increasing the computational work needed to devise low-rank dynamical reduced-order models. Here in this paper, we present two novel approach to carry out parametric DMD: one based on the interpolation of the reduced-order DMD eigen-pair and the other based on the interpolation of the reduced DMD (Koopman) operator. Numerical results are presented for diffusion-dominated nonlinear dynamical problems, including a multiphysics radiative transfer example. All three parametric DMD approaches are compared.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Analysis of parton distributions in a pion with Bézier parametrizations

We explore the role of parametrizations for nonperturbative QCD functions in global analyses, with a specific application to extending a phenomenological analysis of the parton distribution functions (PDFs) in the charged pion realized in the xFitter fitting framework. The parametrization dependence of PDFs in our pion fits substantially enlarges the uncertainties from the experimental sources estimated in the previous analyses. We systematically explore the parametrization dependence by employing a novel technique to automate generation of polynomial parametrizations for PDFs that makes use of Bézier curves. This technique is implemented in a ++ module that is included in the xFitter program. Our analysis reveals that the sea and gluon distributions in the pion are not well disentangled, even when considering measurements in leading-neutron deep inelastic scattering. For example, the pion PDF solutions with a vanishing gluon and large quark sea are still experimentally allowed, which elevates the importance of ongoing lattice and nonperturbative QCD calculations, together with the planned pion scattering experiments, for conclusive studies of the pion structure. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS2001

Aeroservoelastic Uncertainty Model Identification from Flight Data

Uncertainty modeling is a critical element in the estimation of robust stability margins for stability boundary prediction and robust flight control system development. There has been a serious deficiency to date in aeroservoelastic data analysis with attention to uncertainty modeling. Uncertainty can be estimated from flight data using both parametric and nonparametric identification techniques. The model validation problem addressed in this paper is to identify aeroservoelastic models with associated uncertainty structures from a limited amount of controlled excitation inputs over an extensive flight envelope. The challenge to this problem is to update analytical models from flight data estimates while also deriving non-conservative uncertainty descriptions consistent with the flight data. Multisine control surface command inputs and control system feedbacks are used as signals in a wavelet-based modal parameter estimation procedure for model updates. Transfer function estimates are incorporated in a robust minimax estimation scheme to get input-output parameters and error bounds consistent with the data and model structure. Uncertainty estimates derived from the data in this manner provide an appropriate and relevant representation for model development and robust stability analysis. This model-plus-uncertainty identification procedure is applied to aeroservoelastic flight data from the NASA Dryden Flight Research Center F-18 Systems Research Aircraft.

Brenner, Martin J.↗

DOE OSTI2026

Gaussian processes for inferring parton distributions

The extraction of parton distribution functions (PDFs) from experimental or lattice QCD data is an ill-posed inverse problem, where regularization strongly impacts both systematic uncertainties and the reliability of the results. We study a framework based on Gaussian Process Regression (GPR) to reconstruct PDFs from lattice QCD matrix elements. Within a Bayesian framework, Gaussian processes serve as flexible priors that encode uncertainties, correlations, and constraints without imposing rigid functional forms. We investigate a wide range of kernel choices, mean functions, and hyperparameter treatments. We quantify information gained from the data using the Kullback-Leibler divergence. Synthetic data tests demonstrate the consistency and robustness of the method. Our study establishes GPR as a systematic and non-parametric approach to PDF reconstruction, offering controlled uncertainty estimates and reduced model bias in lattice QCD analyses.

hadronic spectroscopy↗

NASA NTRS2022

Aerosol Models from the AERONET Database: Application to Surface Reflectance Validation Satellite Data

Aerosols play a critical role in radiative transfer within the atmosphere, and they have a significant impact on climate change. In this paper, we propose and implement a framework for developing an aerosol model using their microphysical properties. Such microphysical properties as the size distribution, the complex refractive index, and the percentage of sphericity are derived from the global AERosol RObotic NETwork (AERONET). These measurements, however, are typically retrieved when almucantar measurement procedures are performed (i.e., early mornings and late afternoons with clear sky) and might not have a temporal correspondence to a satellite overpass time, so a valid validation of satellite-derived products cannot be carried out. To address this problem of temporal inconsistency of satellite and ground-based measurements, we developed an approach to retrieve these microphysical properties (and the corresponding aerosol model) using the optical thickness at 440 nm, τ440, and the Ångström coefficient between 440 and 870 nm, α440–870. Such aerosol models were developed for 851 AERONET sites within the last 28 years. Obtained results suggest that empirically microphysical properties can be retrieved with uncertainties of up to 23 %. An exception is the imaginary part of the refractive index ni, for which the derived uncertainties reach up to 38 %. These specific parametric models of aerosol can be used for the studies when retrieval of microphysical properties is required as well as validation of satellite-derived products over land. Specifically, we demonstrate the usefulness of the aerosol models to validate surface reflectance records over land derived from optical remote sensing sensors. We then quantify the propagation of uncertainties in the surface reflectance due to uncertainties with the aerosol model retrieval that is used as a reference from radiative transfer simulations. Results indicate that individual aerosol microphysical properties can impact uncertainties in surface reflectance retrievals between 3.5 × 10(exp -5) to 1 × 10(exp -3) (in reflectance units). The overall impact of microphysical properties combined yields an overall uncertainty in surface reflectance <0.004 (in reflectance units). That corresponds, for example, to 1 to 3 % of the retrieved surface reflectance in the red spectral band (620–670 nm) by the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument. These uncertainty values are well below the specification (0.005 + 0.05ρ; ρ is the retrieved surface reflectance) used for the MODIS atmospheric correction.

Jean-Claude Roger↗

DOE OSTI2022

Learning Stochastic Parametric Differentiable Predictive Control Policies

We present a scalable unsupervised learning-based method for obtaining explicit control policies for model predictive control problems for stochastic linear systems with additive uncertainties subject to nonlinear chance constraints. We call the proposed method stochastic parametric differentiable predictive control (SP-DPC), which extends the recently proposed deterministic DPC policy optimization algorithm. We formulate the SP-DPC as a deterministic approximation to the stochastic parametric constrained optimal control problem via independent sampling of the problem's parameters and uncertainties. This formulation allows us to directly compute the policy gradients via automatic differentiation of the problem's value function, evaluated over sampled parameters and uncertainties. In particular, the computed expectation of the problem's value function is backpropagated through the finite-time closed-loop system rollouts parametrized by a known nominal system dynamics model and neural control policy. We also provide theoretical probabilistic guarantees on closed-loop stability and chance constraints satisfaction for systems controlled by learned neural policies. We demonstrate the computational efficiency and scalability of the proposed policy optimization algorithm in three numerical examples, including systems with a large number of states or subject to nonlinear constraints.

Drgona, Jan↗

DOE OSTI2025

SymbolFit: Automatic Parametric Modeling with Symbolic Regression

We introduce SymbolFit (API: https://github.com/hftsoi/symbolfit), a framework that automates parametric modeling by using symbolic regression to perform a machine-search for functions that fit the data while simultaneously providing uncertainty estimates in a single run. Traditionally, constructing a parametric model to accurately describe binned data has been a manual and iterative process, requiring an adequate functional form to be determined before the fit can be performed. The main challenge arises when the appropriate functional forms cannot be derived from first principles, especially when there is no underlying true closed-form function for the distribution. In this work, we develop a framework that automates and streamlines the process by utilizing symbolic regression, a machine learning technique that explores a vast space of candidate functions without requiring a predefined functional form because the functional form itself is treated as a trainable parameter, making the process far more efficient and effortless than traditional regression methods. We demonstrate the framework in high-energy physics experiments at the CERN Large Hadron Collider (LHC) using five real proton-proton collision datasets from new physics searches, including background modeling in resonance searches for high-mass dijet, trijet, paired-dijet, diphoton, and dimuon events. We show that our framework can flexibly and efficiently generate a wide range of candidate functions that fit a nontrivial distribution well using a simple fit configuration that varies only by random seed, and that the same fit configuration, which defines a vast function space, can also be applied to distributions of different shapes, whereas achieving a comparable result with traditional methods would have required extensive manual effort.

Tsoi, Ho Fung [Univ. of Pennsylvania, Philadelphia↗