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 55 records · Page 3

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↗

Extracting Neutrino Oscillation Parameters using a Simultaneous Fit of the $\nu_{e}$ Appearance and $\nu_{\mu}$ Disappearance Data in the NOvA Experiment

NOvA is a two detector, long baseline neutrino oscillation experiment designed to study $\nu_{e}$ ($\bar{\nu}_{e}$) appearance and $\nu_{\mu}$ ($\bar{\nu}_{\mu}$) disappearance in a $\nu_{\mu}$ ($\bar{\nu}_{\mu}$) beam produced at Fermilab. The near detector (ND) is located 100 meters underground at a distance of 1 km from the target whereas the far detector (FD) is located on the surface, 810 km away from the beam source in Ash River, MN. The ND is used to measure the beam before oscillations and the FD measures the oscillated spectrum. The ND and the FD are functionally identical detectors and the ND spectra are extrapolated to the FD to predict the signal and background spectra expected in the FD. The extrapolation and data fitting techniques developed for these analyses within NOvA are presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An improved method to increase the predictive accuracy of the ECR technique

The technique of electrical conductivity relaxation (ECR) is a widely utilized method for obtaining oxygen chemical surface exchange (k) and bulk diffusion (D) coefficients for mixed ionic electronic conductors. Although the technique is conceptually simple, its practical application does not frequently produce reliably accurate values of k and D, and many examples of poor agreement can be found in the published literature. In this present paper, methods for enhancing the accuracy of the ECR technique relative to data processing are described in detail. Physically accurate determination of the actual starting time (t{sub 0}) and improved methods for data fitting of ECR measurements are discussed. To test the methods described, generated relaxation data sets are analyzed by using non-linear-least-square (NLLS) methods. The results showed that fixing t{sub 0} based on an ideally stirred reactor model improves the fitted results compared to treating t{sub 0} as a fitted parameter. Additionally, evaluation of the global minimum in the fitted error is made by analyzing multiple samples of varying thickness, which improves parameter evaluation compared to analysis of a single relaxation data set. The improved accuracy of the computed values of k and D is described in comparison to conventional methods.

Bulk diffusion coefficient↗

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↗

Probing surface recombination velocities in semiconductors using two-photon microscopy

The determination of minority-carrier lifetimes and surface recombination velocities is essential for the development of semiconductor technologies such as solar cells. The recent development of two-photon time-resolved microscopy allows for better measurements of bulk and subsurface interfaces properties. Here, we analyze the diffusion problem related to this optical technique. Our three-dimensional treatment enables us to separate lifetime (recombination) from transport effects (diffusion) in the photoluminescence intensity. It also allows us to consider surface recombination occurring at a variety of geometries: a single plane (representing an isolated exposed or buried interface), a two parallel planes (representing two inequivalent interfaces), and a spherical surface (representing the enclosing surface of a grain boundary). We provide fully analytical results and scalings directly amenable to data fitting and apply those to experimental data collected on heteroepitaxial CdTe/ZnTe/Si.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

evSeq: Cost-Effective Amplicon Sequencing of Every Variant in a Protein Library

Widespread availability of protein sequence-fitness data would revolutionize both our biochemical understanding of proteins and our ability to engineer them. Unfortunately, even though thousands of protein variants are generated and evaluated for fitness during a typical protein engineering campaign, most are never sequenced, leaving a wealth of potential sequence-fitness information untapped. Primarily, this is because sequencing is unnecessary for many protein engineering strategies; the added cost and effort of sequencing is thus unjustified. It also results from the fact that, even though many lower cost sequencing strategies have been developed, they often require at least some sequencing or computational resources, both of which can be barriers to access. In this work, we present every variant sequencing (evSeq), a method and collection of tools/standardized components for sequencing a variable region within every variant gene produced during a protein engineering campaign at a cost of cents per variant. evSeq was designed to democratize low-cost sequencing for protein engineers and, indeed, anyone interested in engineering biological systems. Execution of its wet-lab component is simple, requires no sequencing experience to perform, relies only on resources and services typically available to biology labs, and slots neatly into existing protein engineering workflows. Analysis of evSeq data is likewise made simple by its accompanying software (found at github.com/fhalab/evSeq, documentation at fhalab.github.io/evSeq), which can be run on a personal laptop and was designed to be accessible to users with no computational experience. Here, low-cost and easy to use, evSeq makes collection of extensive protein variant sequence-fitness data practical.

59 BASIC BIOLOGICAL SCIENCES↗

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]↗

RISE: robust iterative surface extension for sub-nanometer X-ray mirror fabrication

Precision optics have been widely required in many advanced technological applications. X-ray mirrors, as an example, serve as the key optical components at synchrotron radiation and free electron laser facilities. They are rectangular silicon or glass substrates where a rectangular Clear Aperture (CA) needs to be polished to sub-nanometer Root Mean Squared (RMS) to keep the imaging capability of the incoming X-ray wavefront at the diffraction limit. The convolutional polishing model requires a CA to be extended with extra data, from which the dwell time is calculated via deconvolution. However, since deconvolution is very sensitive to boundary errors and noise, the existing surface extension methods can hardly fulfill the sub-nanometer requirement. On one hand, the figure errors in a CA were improperly modeled during the extension, leading to continuity issues along the boundary. On the other hand, uncorrectable high-frequency errors and noise were also extended. In this study, we propose a novel Robust Iterative Surface Extension (RISE) method that resolves these problems with a data fitting strategy. RISE models the figure errors in a CA with orthogonal polynomials and ensures that only correctable errors are fit and extended. Combined with boundary conditions, an iterative refinement of dwell time is then proposed to compensate the errors brought by the extension and deconvolution, which drastically reduces the estimated figure error residuals in a CA while the increase of total dwell time is negligible. To our best knowledge, RISE is the first data fitting-based surface extension method and is the first to optimize dwell time based on iterative extension. An experimental verification of RISE is given by fabricating two elliptic cylinders (10 mm × 80 mm CAs) starting from a sphere with a radius of curvature around 173 m using ion beam figuring. The figure errors in the two CAs greatly improved from 204.96 nm RMS and 190.28 nm RMS to 0.62 nm RMS and 0.71 nm RMS, respectively, which proves that RISE is an effective method for sub-nanometer level X-ray mirror fabrication.

36 MATERIALS SCIENCE↗

Applying uncertainty quantification to multiphase flow computational fluid dynamics

Multiphase computational fluid dynamics plays a major role in design and optimization of fossil fuel based reactors. There is a growing interest in accounting for the influence of uncertainties associated with physical systems to increase the reliability of computational simulation based engineering analysis. The U.S. Department of Energy's National Energy Technology Laboratory (NETL) has recently undertaken an initiative to characterize uncertainties associated with computer simulation of reacting multiphase flows encountered in energy producing systems such as a coal gasifier. The current work presents the preliminary results in applying non-intrusive parametric uncertainty quantification and propagation techniques with NETL's open-source multiphase computational fluid dynamics software MFIX. For this purpose an open-source uncertainty quantification toolkit, PSUADE developed at the Lawrence Livermore National Laboratory (LLNL) has been interfaced with MFIX software. In this study, the sources of uncertainty associated with numerical approximation and model form have been neglected, and only the model input parametric uncertainty with forward propagation has been investigated by constructing a surrogate model based on data-fitted response surface for a multiphase flow demonstration problem. Monte Carlo simulation was employed for forward propagation of the aleatory type input uncertainties. Several insights gained based on the outcome of these simulations are presented such as how inadequate characterization of uncertainties can affect the reliability of the prediction results. Also a global sensitivity study using Sobol' indices was performed to better understand the contribution of input parameters to the variability observed in response variable.

42 ENGINEERING↗

NMR_DIFFSIM v. 1.0

The software generates and analyzes nuclear magnetic resonance spin diffusion simulations for the purpose of measuring domain sizes or spin diffusion constants in a heterogenous sample. Its primary capabilities are: 1. It solves the spin magnetization diffusion equation to simulate the results of a spin diffusion experiment using any structural input geometry, including structure models generated by molecular dynamics or other computational or experimental methods. 2. It compares the results of simulations on an input structure to experimental data, providing a data fitting algorithm to extract physical information from the experimental data via the simulation. 3. It can provide statistical measures to compare correlations between structural models, to determine domain clustering and connectedness, and to provide goodness-of-fit measures including chord length distributions, chi-squared values, Pearson coefficients and pair correlation functions. 4. It provides output reports, data files, and movies to visualize simulation results.

Sorte, Eric↗

Rate-equation modelling and ensemble approach to extraction of parameters for viral infection-induced cell apoptosis and necrosis

We develop a theoretical approach that uses physiochemical kinetics modelling to describe cell population dynamics upon progression of viral infection in cell culture, which results in cell apoptosis (programmed cell death) and necrosis (direct cell death). Several model parameters necessary for computer simulation were determined by reviewing and analyzing available published experimental data. By comparing experimental data to computer modelling results, we identify the parameters that are the most sensitive to the measured system properties and allow for the best data fitting. Our model allows extraction of parameters from experimental data and also has predictive power. Using the model we describe interesting time-dependent quantities that were not directly measured in the experiment and identify correlations among the fitted parameter values. Numerical simulation of viral infection progression is done by a rate-equation approach resulting in a system of “stiff” equations, which are solved by using a novel variant of the stochastic ensemble modelling approach. The latter was originally developed for coupled chemical reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Chitosan–Gelatin @ Tin (IV) Tungstatophosphate Nanocomposite Ion Exchanger: Synthesis, Characterization and Applications in Environmental Remediation

The present study reported the synthesis of novel organic–inorganic hybrid nanocomposite by incorporating tin (IV) based ion exchanger into the hybrid polymer network of chitosan and gelatin prepared under vacuum for the efficient removal of heavy metal ions and toxic dyes from an aqueous fluid. The physicochemical studies such as ion exchange capacity (IEC), chemical stability, thermal stability, pH titration and distribution behaviour studies were also carried out to determine the cation exchange behaviour of the material. The surface morphology and structural properties were studied by the techniques such as FTIR, FESEM, EDS, TEM and XRD. Distribution studies confirmed the synthesized CG/STPNC had the highest selectivity for Pb{sup 2+} ions (85.3 mL/g). Maximum adsorption of methylene blue (82%) was achieved within 240 min at 500 mg of adsorbent dose, 10 mg/L of the initial concentration of dye, pH of 7 and 30 °C of temperature. Adsorption kinetic data fitted well with pseudo-second order rate model with R{sup 2} = 0.995. The correlation value 0.95 and favourable R{sub L} = 0.21 of adsorption data suggested better fit for Langmuir adsorption. Thus the synthesized nanocomposite ion exchanger was found to be a promising cation exchanger as well as an adsorbent for heavy metal ion and dye removal from textile industrial effluents. Graphical abstract: .

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Community Data Mining Approach for Surface Complexation Database Development

This paper presents a comprehensive data-to-model workflow, including a findable, accessible, interoperable, reusable (FAIR) community sorption database (newly developed LLNL Surface Complexation/Ion Exchange (L-SCIE) database) along with a data fitting workflow to efficiently optimize surface complexation reaction constants with multiple surface complexation model (SCM) constructs. This workflow serves as a universal framework to mine, compile, and analyze large numbers of published sorption data as well as to estimate reaction constants for parameterizing reactive transport models. Here the framework includes (1) data digitization from published papers, (2) data unification including unit conversions, and (3) data-model integration and reaction constant estimation using geochemical software PHREEQC coupled with the universal parameter estimation code PEST. We demonstrate our approach using an analysis of U(VI) sorption to quartz based on a first L-SCIE implementation, concluding that a multisite SCM construct with carbonate surface species yielded the best fit to community data. Surface complexation reaction constants extracted from this approach captured all available sorption data available in the literature and provided insight into previously published reaction constants and surface complexation model constructs. The L-SCIE sorption database presented herein allows for automating this approach across a wide range of metals and minerals and implementing novel machine learning approaches to reactive transport in the future.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Incorporating Cancer Stem Cells in Radiation Therapy Treatment Response Modeling and the Implication in Glioblastoma Multiforme Treatment Resistance

Purpose: To perform a preliminary exploration with a simplistic mathematical cancer stem cell (CSC) interaction model to determine whether the tumor-intrinsic heterogeneity and dynamic equilibrium between CSCs and differentiated cancer cells (DCCs) can better explain radiation therapy treatment response with a dual-compartment linear-quadratic (DLQ) model. Methods and Materials: The radiosensitivity parameters of CSCs and DCCs for cancer cell lines including glioblastoma multiforme (GBM), non–small cell lung cancer, melanoma, osteosarcoma, and prostate, cervical, and breast cancer were determined by performing robust least-square fitting using the DLQ model on published clonogenic survival data. Fitting performance was compared with the single-compartment LQ (SLQ) and universal survival curve models. The fitting results were then used in an ordinary differential equation describing the kinetics of DCCs and CSCs in response to 2- to 14.3-Gy fractionated treatments. The total dose to achieve tumor control and the fraction size that achieved the least normal biological equivalent dose were calculated. Results: Smaller cell survival fitting errors were observed using DLQ, with the exception of melanoma, which had a low α/β = 0.16 in SLQ. Ordinary differential equation simulation indicated lower normal tissue biological equivalent dose to achieve the same tumor control with a hypofractionated approach for 4 cell lines for the DLQ model, in contrast to SLQ, which favored 2 Gy per fraction for all cells except melanoma. The DLQ model indicated greater tumor radioresistance than SLQ, but the radioresistance was overcome by hypofractionation, other than the GBM cells, which responded poorly to all fractionations. Conclusion: The distinct radiosensitivity and dynamics between CSCs and DCCs in radiation therapy response could perhaps be one possible explanation for the heterogeneous intertumor response to hypofractionation and in some cases superior outcome from stereotactic ablative radiation therapy. The DLQ model also predicted the remarkable GBM radioresistance, a result that is highly consistent with clinical observations. The radioresistance putatively stemmed from accelerated DCC regrowth that rapidly restored compartmental equilibrium between CSCs and DCCs.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Optical response of columnar TiO{sub 2}:Co on Si substrate

Vertically arranged columnar TiO{sub 2}:Co structure has been deposited on Si (100) substrate using the metal organic chemical vapor deposition (MOCVD) technique. SEM images confirm the columnar structure of the TiO{sub 2}:Co layer on Si (100). Concentration of the Co dopant has been defined by means of the energy dispersive x-ray (EDX) spectra. Co concentrations are varied at 0.18, 0.33, and 0.68 at. % using pure TiO{sub 2} as the reference. X-ray diffraction (XRD) spectra show that the TiO{sub 2}:Co layer have a rutile phase. Further analysis of this structure has been performed by using room temperature spectroscopic ellipsometry (SE) measurements in the energy range of 0.5 - 6.5 eV. This optical response shows the unique feature of significantly high density fringe especially below the band gap energy region (lower than 3.3 eV) in contrast to the conventional TiO{sub 2} thin film structure. This unique feature is interpreted as a dominant contribution of the vertical interfaces at columnar boundaries towards reflected light from the sample. The evolution of Ψ and Δ will be discussed as the function of Co concentration and some optical properties will be extracted after reasonably well data fitting by multilayer modelling. Furthermore, by comparing the SE data from different measurement orientations, it has been shown that TiO{sub 2}:Co tends to be anisotropic. This study enables us to further evaluate structural properties through optical responses given by SE measurement.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Non-local gravity and comparison with observational datasets

We study the cosmological predictions of two recently proposed non-local modifications of General Relativity. Both models have the same number of parameters as ΛCDM, with a mass parameter m replacing the cosmological constant. We implement the cosmological perturbations of the non-local models into a modification of the CLASS Boltzmann code, and we make a full comparison to CMB, BAO and supernova data. We find that the non-local models fit these datasets very well, at the same level as ΛCDM. Among the vast literature on modified gravity models, this is, to our knowledge, the only example which fits data as well as ΛCDM without requiring any additional parameter. For both non-local models parameter estimation using Planck +JLA+BAO data gives a value of H{sub 0} slightly higher than in ΛCDM.

79 ASTRONOMY AND ASTROPHYSICS↗

Clustering with general photo- z uncertainties: application to Baryon Acoustic Oscillations

ABSTRACT Photometric data can be analysed using the 3D correlation function ξp to extract cosmological information via e.g. measurement of the Baryon Acoustic Oscillations (BAO). Previous studies modeled ξp assuming a Gaussian photo-z approximation. In this work we improve the modeling by incorporating realistic photo-z distribution. We show that the position of the BAO scale in ξp is determined by the photo-z distribution and the Jacobian of the transformation. The latter diverges at the transverse scale of the separation s⊥, and it explains why ξp traces the underlying correlation function at s⊥, rather than s, when the photo-z uncertainty σz/(1+ z) ≳ 0.02. We also obtain the Gaussian covariance for ξp. Due to photo-z mixing, the covariance of ξp shows strong off-diagonal elements. The high correlation of the data causes some issues to the data fitting. None the less, we find that either it can be solved by suppressing the largest eigenvalues of the covariance or it is not directly related to the BAO. We test our BAO fitting pipeline using a set of mock catalogs. The data set is dedicated for Dark Energy Survey Year 3 (DES Y3) BAO analyses and includes realistic photo-z distributions. The theory template is in good agreement with mock measurement. Based on the DES Y3 mocks, ξp statistic is forecast to constrain the BAO shift parameter α to be 1.001 ± 0.023, which is well consistent with the corresponding constraint derived from the angular correlation function measurements. Thus, ξp offers a competitive alternative for the photometric data analyses.

79 ASTRONOMY AND ASTROPHYSICS↗