Search NASA⌕ Search

SEARCH · Search NASA

Results for “Data fitting”

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 19 records

Automated Resonance Fitting for Nuclear Data Evaluation

Global and national efforts to deliver high-quality nuclear data to users have a wide-ranging impact, affecting applications in national security, reactor operations, basic science, medicine, and more. Cross section evaluation is a major part of this effort, combining theory and experimentation to produce recommended values and uncertainties for reaction probabilities. Resonance region evaluation is a specialized type of nuclear data evaluation that can require significant manual effort and months of time from expert scientists. In this article, non-convex non-linear optimization methods are combined with concepts of inferential statistics to infer a resonance model from experimental data in an automated manner that is not dependent on prior evaluation(s). This methodology aims to enhance the workflow of a resonance evaluator by minimizing time, effort, and the potential for bias from prior assumptions, while enhancing reproducibility and documentation, thereby addressing well-known challenges in the field.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

How robust are estimates of key parameters in standard viral dynamic models?

Mathematical models of viral infection have been developed, fitted to data, and provide insight into disease pathogenesis for multiple agents that cause chronic infection, including HIV, hepatitis C, and B virus. However, for agents that cause acute infections or during the acute stage of agents that cause chronic infections, viral load data are often collected after symptoms develop, usually around or after the peak viral load. Consequently, we frequently lack data in the initial phase of viral growth, i.e., when pre-symptomatic transmission events occur. Missing data may make estimating the time of infection, the infectious period, and parameters in viral dynamic models, such as the cell infection rate, difficult. However, having extra information, such as the average time to peak viral load, may improve the robustness of the estimation. Here, we evaluated the robustness of estimates of key model parameters when viral load data prior to the viral load peak is missing, when we know the values of some parameters and/or the time from infection to peak viral load. Although estimates of the time of infection are sensitive to the quality and amount of available data, particularly pre-peak, other parameters important in understanding disease pathogenesis, such as the loss rate of infected cells, are less sensitive. Viral infectivity and the viral production rate are key parameters affecting the robustness of data fits. Fixing their values to literature values can help estimate the remaining model parameters when pre-peak data is missing or limited. We find a lack of data in the pre-peak growth phase underestimates the time to peak viral load by several days, leading to a shorter predicted growth phase. On the other hand, knowing the time of infection (e.g., from epidemiological data) and fixing it results in good estimates of dynamical parameters even in the absence of early data. While we provide ways to approximate model parameters in the absence of early viral load data, our results also suggest that these data, when available, are needed to estimate model parameters more precisely.

59 BASIC BIOLOGICAL SCIENCES↗

Bottom-up design of actinide materials from molecular clusters: Demonstration of a general-purpose simulation capability leveraging machine-learned atomic potentials

Actinide thin-film coatings such as uranium dioxide (UO 2 ) play an important role in nuclear reactors and other mission-relevant applications, but realization of their potential requires a deep fundamental understanding of the chemical vapor deposition (CVD) processes used for their growth. The slow experimental progress can be attributed, in part, to the standard safety guidelines associated with handling uranium byproducts, which are often corrosive, toxic, and radioactive. Accurate simulation techniques, when used in concert with experiment, can improve laboratory safety, material durability, and deliverable timeframes. However, state-of-the-art computational methods are either insufficiently accurate or intractably expensive. To remedy this situation, in this project we suggested a machine-learning (ML) accelerated workflow for simulating molecular clustering toward deposition. As a benchmark test case, we considered molecular clustering in steam and assessed independent components of our workflow by comparing with measured thermodynamic properties of water. After analyzing each component individually and finding no fundamental barrier to realization of the workflow, we attempted to integrate the ML component, a Sandia-developed tool called FitSNAP. As this was the first application of FitSNAP to atoms and molecules in the gas phase at Sandia, the method required more fitting data than was originally anticipated. Systematic improvements were made by including in the fit data diatomic potentials, molecular single-bond-breaking curves, and symmetry-constrained intermolecular potentials. We concluded that our strategy provides a feasible pathway toward modeling CVD and related processes, but that extensive training data must be generated before it can be of practical use.

36 MATERIALS SCIENCE↗

Knowledge Graph of RB-Tnseq Data from Fitness Browser (KP-DP1)

Motivation: Predicting microbial gene fitness across environmental conditions remains a central challenge for predictive phenomics and autonomous experimentation. Fitness assays generate large volumes of genotype–phenotype measurements difficult to integrate with experimental metadata and biological function in a form that supports mechanistic reasoning. Knowledge graphs offer a semantic framework for unifying modalities and enabling context-aware inference. Results: We build GIMME (Graph Inference for Microbial Metabolism Exploration), a semantically grounded knowledge graph that unifies gene fitness measurements spanning 10 Pseudomonas species with experimental metadata and biological context. Media are decomposed into chemical components and experiments carry structured links to natural-language descriptions. The resulting graph supports two inference modes: (1) symbolic graph traversal to surface candidate gene–environment and gene–chemical associations, and (2) learned inference using heterogeneous graph neural networks that propagate information across neighborhoods. We formulate link regression over (gene, media, experiment) triplets, combining learned gene embeddings with pretrained LLM sourced text embeddings of node descriptions to predict gene fitness. We then augment a baseline MLP with an auxiliary message-passing encoder (GraphSAGE/GAT) that propagates information over gene–protein–function and media–chemical subgraphs, and fuse the two pathways with a gated residual connection. This approach produces strong agreement with held-out fitness measurements (GraphSAGE Pearson r 0.74) while also highlighting inference challenges in extreme-fitness regimes. We aggregate GAT edge-attention weights by relation type and layer to estimate which biological and environmental relations most influence fitness predictions. Conclusion: This work explores using knowledge graphs as “context graphs” for microbial phenotype prediction. They provide a rich substrate which enables explainable retrieval of supporting evidence, and provides a natural bridge to autonomous workflows that prioritize the next experiment.

59 BASIC BIOLOGICAL SCIENCES↗

Single-channel and single-energy partial-wave analysis with continuity improved through minimal phase constraints

Single-energy partial-wave analysis has often been applied as a way to fit data with minimal model dependence. However, remaining unconstrained, partial waves at neighboring energies will vary discontinuously because the overall amplitude phase cannot be determined through single-channel measurements. This problem can be mitigated through the use of a constraining penalty function based on an associated energy-dependent fit. However, the weight given to this constraint results in a biased fit to the data. In this paper, for the first time, we explore a constraining function which does not influence the fit to data. The constraint comes from the overall phase found in multichannel fits which, in the present study, are the Bonn-Gatchina and Jülich-Bonn multichannel analyses. The data are well reproduced and weighting of the penalty function does not influence the result. The method is applied to K⁢Λ photoproduction data and all observables can be maximally well reproduced. While the employed multichannel analyses display very different multipole amplitudes, we show that the major difference between two sets of multipoles can be related to the different overall phases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

XES_Neo_Public

XES Neo is a fitting software that was based on the already approved EXAFS Neo genetic algorithm fitting software code. Using the principles of genetics, data fitting is done for x-ray emission spectroscopy data. Using the EXAFS Neo open source code, as well as the open source xes_neo code written by other collaborators, a new final repository for XES Neo was created with several necessary changes for general user use.

Humiston, Alaina [Los Alamos National Laboratory]↗

Paired autoencoders for likelihood-free estimation in inverse problems

Abstract We consider the solution of nonlinear inverse problems where the forward problem is a discretization of a partial differential equation. Such problems are notoriously difficult to solve in practice and require minimizing a combination of a data-fit term and a regularization term. The main computational bottleneck of typical algorithms is the direct estimation of the data misfit. Therefore, likelihood-free approaches have become appealing alternatives. Nonetheless, difficulties in generalization and limitations in accuracy have hindered their broader utility and applicability. In this work, we use a paired autoencoder framework as a likelihood-free estimator (LFE) for inverse problems. We show that the use of such an architecture allows us to construct a solution efficiently and to overcome some known open problems when using LFEs. In particular, our framework can assess the quality of the solution and improve on it if needed. We demonstrate the viability of our approach using examples from full waveform inversion and inverse electromagnetic imaging.

Chung, Matthias (ORCID:0000000178224539)↗

Longitudinal Solid P Longitudinal Solid Polarized Target for CL get for CLAS12 and Study of Spin AS12 and Study of Spin Structure of Nucleons

A suite of experiments measuring target-spin observables in electron-nucleon scattering (dubbed Run Group C) was conducted at Jefferson Lab's Hall B in Newport News, VA with a new polarized nuclear target known as 'APOLLO' (Ammonia POLarized LOngitudinally). This innovative target is engineered to seamlessly integrate with the advanced 12GeV CEBAF (Continuous Electron Beam Accelerator Facility) accelerator and the Hall B CLAS12 (12 GeV CEBAF Large Acceptance Spectrometer) detector array. The 'APOLLO' target harnesses the power of Dynamic Nuclear Polarization (DNP) to achieve longitudinal polarization of solid ammonia, thereby creating a net polarization in both protons (NH3) and deuterons (ND3). These samples are subjected to a 5 Tesla magnetic field produced by the CLAS12 spectrometer central solenoid, cooled to 1 K using helium evaporation, and subsequently exposed to microwave radiation at a frequency of 140 GHz. This project was made possible through a collaborative partnership between the Jefferson Lab Target Group, Old Dominion University, University of Virginia, Christopher Newport University, and the CLAS Collaboration. Integrating the APOLLO target into CLAS12 presented unique challenges because of its specific spatial and dimensional constraints. This dissertation will highlight the innovative solutions developed to address these challenges, including the creation of a new target material transport system, the incorporation of superconducting magnetic correction coils, and the development of a customized Nuclear Magnetic Resonance (NMR) system. In addition to a detailed description of the development and operation of APOLLO, and an overview of RGC, this document also describes the global data fitting of proton asymmetries in the Deep Inelastic scattering (DIS) and the Resonance regions, and it offers insights into the comprehensive analysis of target polarization data.

Pandey, Pushpa↗

Predictive models of the genetic bases underlying budding yeast fitness in multiple environments

Abstract The ability of organisms to adapt and survive depends on the effects of genes and the environment on fitness. However, the multigenic nature of fitness and genotype-by-environment interactions hinder our understanding of the genetic basis of fitness. Here, we established fitness prediction models for 35 environments using machine learning and existing fitness data and different genetic variant types for a Saccharomyces cerevisiae population. Models revealed that the predictive ability of genetic variants varied across environments, with copy number variants explaining the majority of fitness variation in most cases. Model interpretation showed that different variant types identified distinct gene sets associated with predictive variants. These gene sets were significantly enriched in experimentally validated genes affecting fitness in only a subset of environments, indicating that many genes influencing fitness remain unexplored. Notably, non-experimentally validated genes were more important than validated ones for fitness predictions. Gene contributions to predictions were both isolate- and environment-dependent, pointing to gene-by-gene and gene-by-environment interactions. Furthermore, models uncovered experimentally validated and novel candidate genetic interactions for a well-characterized stress, the fungicide benomyl. These findings highlight the feasibility of identifying the genetic basis of fitness by using different genetic variant types and offer novel targets for future functional analysis.

DNA copy number variations↗

Deriving Stable Peak Models to Fit Complex XPS Data From Cu Contaminated Pt Electrocatalysts

X-ray Photoelectron Spectroscopy spectra peak models, designed to partition photoemission signals emanating from different elements or chemical states within an atom, are fitted to data limited to an energy interval over which inelastically scattered photoemission signal can be estimated. While the choice of background approximation and line shapes of components to the peak model requires careful consideration, the energy interval used to define the data to which the peak model is optimized has a significant impact on the final peak model. The relationship between the background intensity and data intensity at the start and end of the energy interval dictates the line shapes used in the peak model. In this work, we devise a method to peak fit a complex overlapping Cu 3p and Pt 4f XPS peak structure to perform the elemental quantification. We first use an Al 2s peak to illustrate how background curves approach data at the limits of the energy interval over which the background is defined, influencing the analysis of XPS spectra. Next, we demonstrate the nature of interactions between specific line shapes (Voigt and pseudo-Voigt profiles) suitable for photoemission peaks and a specific background curve (Shirley) and a peak model is presented that includes components to the peak model that accommodates background intensity during fitting of the peak model to data. The peak model allowed for quantification of the contributions of Pt 4f peaks emanating from the substrate that exhibits strong asymmetry in the presence of the inhomogeneously distributed Cu species, mostly of Lorentzian character.

XPS↗

High-throughput genetics enables identification of nutrient utilization and accessory energy metabolism genes in a model methanogen

Archaea are widespread in the environment and play fundamental roles in diverse ecosystems; however, characterization of their unique biology requires advanced tools. This is particularly challenging when characterizing gene function. Here, we generate randomly barcoded transposon libraries in the model methanogenic archaeon Methanococcus maripaludis and use high-throughput growth methods to conduct fitness assays (RB-TnSeq) across over 100 unique growth conditions. Using our approach, we identified new genes involved in nutrient utilization and response to oxidative stress. We identified novel genes for the usage of diverse nitrogen sources in M. maripaludis including a putative regulator of alanine deamination and molybdate transporters important for nitrogen fixation. Furthermore, leveraging the fitness data, we inferred that M. maripaludis can utilize additional nitrogen sources including $\tiny{L}$-glutamine, $\tiny{D}$-glucuronamide, and adenosine. Under autotrophic growth conditions, we identified a gene encoding a domain of unknown function (DUF166) that is important for fitness and hypothesize that it has an accessory role in carbon dioxide assimilation. Finally, comparing fitness costs of oxygen versus sulfite stress, we identified a previously uncharacterized class of dissimilatory sulfite reductase-like proteins (Dsr-LP; group IIId) that is important during growth in the presence of sulfite. When overexpressed, Dsr-LP conferred sulfite resistance and enabled use of sulfite as the sole sulfur source. The high-throughput approach employed here allowed for generation of a large-scale data set that can be used as a resource to further understand gene function and metabolism in the archaeal domain.

59 BASIC BIOLOGICAL SCIENCES↗

Models for Single–Site Heterogeneous Catalysts on Carbon: MoO 2 Epoxidation Catalyst Anchored to a Fullerene

Single-site molybdenum dioxo catalysts, fullerenol/MoO 2 , are prepared via grafting precursor (DME)MoO 2 Cl 2 onto a highly polyhydroxylated fullerene (ful) and an isomerically-pure and well-defined fullerene (ful*). These catalyst structures are characterized by ICP-OES, XPS, XANES, EXAFS, DRIFT, Raman, and NMR spectroscopy, and DFT. Mo 3d 5/2 XPS and Mo K-edge XANES assign the oxidation state as Mo(VI). Mo EXAFS data fitting reveals two Mo=O double and two Mo–O single bonds at distances of 1.7 and 1.9 Å, respectively, while an Mo=O stretchingl mode is observed at ~950 cm –1 by DRIFT and Raman spectroscopy. These data align well with DFT computational results, supporting the proposed catalyst structure as Fullerene(-μ-O-) 2 M(=O) 2 . Additionally, DFT provides insight into the energetically favorable grafting sites for an isomerically pure fullerenol. The scope of fullerenol/MoO 2 mediated alkene epoxidation includes abiotic alkenes, natural occurring terpenes, and conjugated olefins. For cyclooctene the rate law is first-order in [Mo], near first order in [olefin] and zero-order in [t-butyl hydroperoxide]. A plausible reaction mechanism involves peroxide addition first and then cyclooctene addition directly across the peroxo bond forming the epoxide product, consistent with DFT computation. Overall, fullerenol/MoO 2 shows promise as a sustainable and structurally well-defined system with versatile catalytic activity and good epoxidation recyclability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

RB-TnSeq elucidates dicarboxylic-acid-specific catabolism in β-proteobacteria for improved plastic monomer upcycling

Dicarboxylic acids are key components of many polymers and plastics, making them a target for both engineered microbial degradation and sustainable bioproduction. In this study, we generated a comprehensive data set of functional evidence for the genetic basis of dicarboxylic and fatty acid metabolism using randomly barcoded transposon sequencing (RB-TnSeq). We identified four β-proteobacteria that displayed robust growth with dicarboxylic acid sole carbon source and cultured their mutant libraries with dicarboxylic and fatty acids with carbon chain lengths from C3 to C12. The resulting fitness data suggested that dicarboxylic and fatty acid metabolisms are largely distinct, and different sets of β-oxidation genes are required for catabolizing dicarboxylic versus fatty acids of the same carbon chain lengths. In addition, we identified transcriptional regulators and transporters with strong fitness phenotypes related to dicarboxylic acid utilization. In Ralstonia sp. UNC404CL21Col (R. CL21), we deleted two transcriptional repressors to improve its utilization of short-chain dicarboxylic acids. We exploited the diacid-utilizing catabolism of R. CL21 to upcycle a mock mixture of the dicarboxylic acids produced when polyethylene is oxidized. After introducing a heterologous indigoidine production pathway, this engineered Ralstonia produced 0.56 ± 0.02 g/L indigoidine from a mixture of dicarboxylic acids as a carbon source, demonstrating the potential of R. CL21 to upcycle plastic wastes to products derived from tricarboxylic acid (TCA) cycle intermediates. IMPORTANCE: Upcycling the carbon in plastic wastes to value-added products is a promising approach to address the plastic waste and climate crises, and dicarboxylic acid metabolism is an important facet of several approaches. Improving our understanding of the genetic basis of this metabolism has the potential to uncover new enzymes and genetic parts for engineered pathways involving dicarboxylic acids. Our data set is the most comprehensive interrogation of dicarboxylic acid catabolism to date, and this work will be of utility to researchers interested in both plastics bioproduction and upcycling applications.

Pearson, Allison N↗

Information theory optimization of signals from small-angle scattering measurements

Small-angle X-ray scattering (SAXS) of particles in solution informs on the conformational states and assemblies of biological macromolecules (bioSAXS) outside of cryo- and solid-state conditions. In bioSAXS, the SAXS measurement under dilute conditions is resolution limited, and through an inverse Fourier transform, the measured SAXS intensities directly relate to the physical space occupied by the particles via the P (r)-distribution. Yet, this inverse transform of SAXS data has been historically cast as an ill-posed, ill-conditioned problem requiring an indirect approach. Here, we show that through the applications of matrix and information theories, the inverse transform of SAXS intensity data is a well-conditioned problem. The so-called ill-conditioning of the inverse problem is directly related to the Shannon number. By exploiting the oversampling enabled by modern detectors, a direct inverse Fourier transform of the SAXS data is possible, provided the recovered information does not exceed the Shannon number. The Shannon limit corresponds to the maximum number of significant singular values that can be recovered in a SAXS experiment, suggesting this relationship is a fundamental property of band-limited inverse integral transform problems. This correspondence reduces the complexity of the inverse problem to the Shannon limit and maximum dimension. We propose a hybrid scoring function using an information theory framework that assesses both the quality of the model-data fit as well as the quality of the recovered P (r)-distribution. The hybrid score utilizes the Akaike information criteria and Durbin-Watson statistic that considers parameter-model complexity, i.e., degrees of freedom, and the randomness of the model-data residuals. The described tests and findings extend the boundaries for bioSAXS by completing the information theory formalism initiated by Peter B. Moore to enable a quantitative measure of resolution in SAXS, robustly determine maximum dimension, and more precisely define the best parameter model appropriately representing the observed scattering data.

Rambo, Robert P. [Science and Technology Facilitie↗

Initial measurement of reactor antineutrino oscillation at SNO+

The SNO+ collaboration reports its first spectral analysis of long-baseline reactor antineutrino oscillation using 114 tonne-years of data. Fitting the neutrino oscillation probability to the observed energy spectrum yields constraints on the neutrino mass-squared difference Δ$m^2_{21}$. In the ranges allowed by previous measurements, the best-fit Δ$m^2_{21}$ is ($8.85^{+1.10}_{-1.33}$) × 10 -5 eV 2 . This measurement is continuing in the next phases of SNO+ and is expected to surpass the present global precision on Δ$m^2_{21}$ with about three years of data. In the following, we first describe the configuration of the SNO+ detector when it was partially filled with scintillator, and a characterization of the detector response using intrinsic radioactivity. Next, we detail the event selection and expectations for reactor IBDs and (α, n) reactions. Then, we present the results of an energy spectrum analysis using 114 tonne-years of data. We conclude with prospects of future results from the SNO+ detector, which has been operating fully-filled with 780 tonnes of scintillator

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Elucidation of Marcus Relationships for Hydride Transfer Reactions Involving Transition Metal Hydrides

The rate of hydride transfer from three Ir hydride complexes of the type Cp*Ir( R bpy)H + (Cp* = C 5 Me 5 ; R bpy = 4,4′-R-2,2′-bipyridine, R = OMe, H, CO 2 Me) to six N-methylacridinium ( R Acr + ) acceptors with electronically different substituents in the 2- or 2,7-positions were measured. Using the thermodynamic hydricity of the donors and the hydride affinity of the acceptors the thermodynamic driving forces for hydride transfer were determined. Brønsted plots, which correlate kinetic and thermodynamic hydricity, demonstrate distinct linear free energy relationships for each complex, with different Brønsted α values. Thus, at the same driving force hydride transfer from Cp*Ir( OMe bpy)H + is faster than for Cp*Ir(bpy)H + or Cp*Ir( CO2Me bpy)H + . Experimental and computational analyses are consistent with a concerted hydride transfer mechanism for all Ir complexes. As the thermodynamic driving force increases an earlier transition state is observed and all transition states also include π-stacking interactions between the donor and acceptor, which likely contribute to the different α values. The experimental data fits well to the Marcus model, enabling the determination of reorganization energies (λ) that range from 58 to 69 kcal mol -1 . These are lower than λ values for hydride transfer reactions involving organic donors and acceptors. This work provides a rare example of the correlation of kinetic and thermodynamic hydricity using only experimental data and shows that hydride transfer reactions involving metal hydrides can follow Marcus theory. Furthermore, the findings offer insight into controlling metal-catalyzed hydride transfer reactions, which is valuable for designing improved systems for a range of transformations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Polymer Size–Catalytic Activity Relationships in Solution by Fluorescence Correlation Spectroscopy

Measuring the catalytic activity of specific sizes of polymers with active catalysts in solution is typically challenging, due to limited instrument detection sensitivity and/or dynamic range. Here, a fluorescence correlation spectroscopy (FCS) method is developed to determine the catalytic activity of living polymers of a specific apparent size in solution. Deviation from a single-component FCS data fitting, as assessed by χ2, is also introduced and developed as a “speciation index”—a method to evaluate and track changes in the relative amount of distinct polymer sizes with reaction progress. These methods are enabled by incorporating a selectively reactive fluorescent monomer into growing polydicyclopentadiene or polynorbornene during ring-opening metathesis polymerization (ROMP). Compared to polynorbornene, data showed that catalysts in aggregates of polyDCPD retained higher activity for longer—outcomes not directly inferable from simple diffusional-access predictions. Here, the ability to assign catalytic activity to polymers of specific sizes, and then to determine how this activity evolves with reaction progress, support long-term goals in the development and measurement of nano-objects that possess size-dependent catalytic activity.

Active catalyst↗

Fit-free optical determination of electronic thermalization time in nematic iron-based superconductors

We present a nematic response function model (NRFM) for fit-free direct extraction of the characteristic time of ultrafast electronic thermalization in iron-based superconductors, materials with electronic nematicity. By combining the NRFM for polarization-dependent pump–probe measurements of electronic nematic response with the two-temperature model (TTM) for sub-picosecond quasiparticle relaxation, we quantify the electronic thermalization timescales and their anisotropy. The nematic response function is modeled as the difference in normalized reflectivity signals, revealing a pronounced sub-picosecond extremum in signal evolution that directly yields the characteristic electronic thermalization time. This method demonstrates that the NRFM is consistent with TTM fits of transient optical response, yielding electronic thermalization time constants on the order of 110–230 fs for the FeSe 1−x Te x and Ba(Fe 0.92 Co 0.08 ) 2 As 2 thin films. Here, the proposed approach can be applied to any material that exhibits electronic nematicity, providing a powerful tool for direct mapping of the relaxation time in nematic materials, avoiding complex experimental data-fitting procedures.

Bartenev, Alexander [University of Puerto Rico, Ma↗