Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model 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 1,423 records · Page 79

Enhanced Jahn-Teller response induced by low-dose 10 MeV I{sup +} irradiation of La{sub 0.7}Ca{sub 0.3}MnO{sub 3−δ} films

The structural response of pulsed laser deposited La{sub 0.7}Ca{sub 0.3}Mn0{sub 3−δ} films to low-level 10 MeV I{sup +} irradiation is correlated with magnetic and electronic properties. A series of annealed pulsed laser deposited films were subjected to ion irradiation at doses of 0.5–2.0 × 10{sup 13} 10 MeV iodine ions/cm{sup 2}. X-ray diffraction measurements show no measurable increase in the full width at half maximum values of major reflections before and after irradiation. Extended x-ray absorption fine structure analyses (EXAFS), employing both a model independent fitting using expanded cumulants as well as a nonlinear least squares multiple scattering model based on FEFF generated theoretical EXAFS data, were employed to measure local properties relative to the Mn ions. Increasing irradiation doses lead to an unambiguous evolution of MnO{sub 6} octahedra from highly symmetric (i.e., undistorted) to a clear bimodal structure in which the MnO{sub 6} octahedra are measured to undergo a c-axis stretching bearing resemblance to a Jahn-Teller distortion. Concomitantly, the metal to insulator transition temperature decreases more than 270 K and the magnetoresistance increases from 10{sup 2}% to 10{sup 6}%.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Optical shielding of nickel nanoparticle by a bubble: Optical limiting gets limited

We have demonstrated that in a nickel nanoparticle colloid, the optical limiting action reduces if a vapor bubble forms around the nanoparticle. The energy-dependent transmission and z-scan measurements on nickel nanoparticles in toluene show the onset of an additional process. At high fluence excitation, the particle becomes less visible to the later part of the incoming pulse due to the heat generated bubble formed around it. We have proposed a simple “particle-in-bubble” model which fits the optical limiting and z-scan curves quite well. Using this model, we have also estimated that the bubble radius increases at a rate of 4.5 m/s.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

THE FERMI GBM GAMMA-RAY BURST SPECTRAL CATALOG: FOUR YEARS OF DATA

In this catalog we present the updated set of spectral analyses of gamma-ray bursts (GRBs) detected by the Fermi Gamma-Ray Burst Monitor during its first four years of operation. It contains two types of spectra, time-integrated spectral fits and spectral fits at the brightest time bin, from 943 triggered GRBs. Four different spectral models were fitted to the data, resulting in a compendium of more than 7500 spectra. The analysis was performed similarly but not identically to Goldstein et al. All 487 GRBs from the first two years have been re-fitted using the same methodology as that of the 456 GRBs in years three and four. We describe, in detail, our procedure and criteria for the analysis and present the results in the form of parameter distributions both for the observer-frame and rest-frame quantities. The data files containing the complete results are available from the High-Energy Astrophysics Science Archive Research Center.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of the theoretical bias in dark matter direct detection

Fitting the model ''A'' to dark matter direct detection data, when the model that underlies the data is ''B'', introduces a theoretical bias in the fit. We perform a quantitative study of the theoretical bias in dark matter direct detection, with a focus on assumptions regarding the dark matter interactions, and velocity distribution. We address this problem within the effective theory of isoscalar dark matter-nucleon interactions mediated by a heavy spin-1 or spin-0 particle. We analyze 24 benchmark points in the parameter space of the theory, using frequentist and Bayesian statistical methods. First, we simulate the data of future direct detection experiments assuming a momentum/velocity dependent dark matter-nucleon interaction, and an anisotropic dark matter velocity distribution. Then, we fit a constant scattering cross section, and an isotropic Maxwell-Boltzmann velocity distribution to the simulated data, thereby introducing a bias in the analysis. The best fit values of the dark matter particle mass differ from their benchmark values up to 2 standard deviations. The best fit values of the dark matter-nucleon coupling constant differ from their benchmark values up to several standard deviations. We conclude that common assumptions in dark matter direct detection are a source of potentially significant bias.

79 ASTRONOMY AND ASTROPHYSICS↗

Molecular radiotherapy: The NUKFIT software for calculating the time-integrated activity coefficient

Purpose: Calculation of the time-integrated activity coefficient (residence time) is a crucial step in dosimetry for molecular radiotherapy. However, available software is deficient in that it is either not tailored for the use in molecular radiotherapy and/or does not include all required estimation methods. The aim of this work was therefore the development and programming of an algorithm which allows for an objective and reproducible determination of the time-integrated activity coefficient and its standard error.Methods: The algorithm includes the selection of a set of fitting functions from predefined sums of exponentials and the choice of an error model for the used data. To estimate the values of the adjustable parameters an objective function, depending on the data, the parameters of the error model, the fitting function and (if required and available) Bayesian information, is minimized. To increase reproducibility and user-friendliness the starting values are automatically determined using a combination of curve stripping and random search. Visual inspection, the coefficient of determination, the standard error of the fitted parameters, and the correlation matrix are provided to evaluate the quality of the fit. The functions which are most supported by the data are determined using the corrected Akaike information criterion. The time-integrated activity coefficient is estimated by analytically integrating the fitted functions. Its standard error is determined assuming Gaussian error propagation. The software was implemented using MATLAB.Results: To validate the proper implementation of the objective function and the fit functions, the results of NUKFIT and SAAM numerical, a commercially available software tool, were compared. The automatic search for starting values was successfully tested for reproducibility. The quality criteria applied in conjunction with the Akaike information criterion allowed the selection of suitable functions. Function fit parameters and their standard error estimated by using SAAM numerical and NUKFIT showed differences of <1%. The differences for the time-integrated activity coefficients were also <1% (standard error between 0.4% and 3%). In general, the application of the software is user-friendly and the results are mathematically correct and reproducible. An application of NUKFIT is presented for three different clinical examples.Conclusions: The software tool with its underlying methodology can be employed to objectively and reproducibly estimate the time integrated activity coefficient and its standard error for most time activity data in molecular radiotherapy.

61 RADIATION PROTECTION AND DOSIMETRY↗

WE-D-BRF-05: Quantitative Dual-Energy CT Imaging for Proton Stopping Power Computation

Purpose: To extend the two-parameter separable basis-vector model (BVM) to estimation of proton stopping power from dual-energy CT (DECT) imaging. Methods: BVM assumes that the photon cross sections of any unknown material can be represented as a linear combination of the corresponding quantities for two bracketing basis materials. We show that both the electron density (ρe) and mean excitation energy (Iex) can be modeled by BVM, enabling stopping power to be estimated from the Bethe-Bloch equation. We have implemented an idealized post-processing dual energy imaging (pDECT) simulation consisting of monogenetic 45 keV and 80 keV scanning beams with polystyrene-water and water-CaCl2 solution basis pairs for soft tissues and bony tissues, respectively. The coefficients of 24 standard ICRU tissue compositions were estimated by pDECT. The corresponding ρe, Iex, and stopping power tables were evaluated via BVM and compared to tabulated ICRU 44 reference values. Results: BVM-based pDECT was found to estimate ρe and Iex with average and maximum errors of 0.5% and 2%, respectively, for the 24 tissues. Proton stopping power values at 175 MeV, show average/maximum errors of 0.8%/1.4%. For adipose, muscle and bone, these errors result range prediction accuracies less than 1%. Conclusion: A new two-parameter separable DECT model (BVM) for estimating proton stopping power was developed. Compared to competing parametric fit DECT models, BVM has the comparable prediction accuracy without necessitating iterative solution of nonlinear equations or a sample-dependent empirical relationship between effective atomic number and Iex. Based on the proton BVM, an efficient iterative statistical DECT reconstruction model is under development.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

CD8 + lymphocyte control of SIV infection during antiretroviral therapy

CD8 + lymphocytes play an important role in suppressing in vivo viral replication in HIV infection. However, both the extent to which and the mechanisms by which CD8 + lymphocytes contribute to viral control are not completely understood. A recent experiment depleted CD8 + lymphocytes in simian immunodeficiency virus (SIV)-infected rhesus macaques (RMs) on antiretroviral treatment (ART) to study the role of CD8 + lymphocytes. CD8 + lymphocytes depletion resulted in temporary plasma viremia in all studied RMs. Viral control was restored when CD8 + lymphocytes repopulated. We developed a viral dynamic model to fit the viral load (VL) data from the CD8 depletion experiment. We explicitly modeled the dynamics of the latent reservoir and the SIV-specific effector cell population including their exhaustion and their potential cytolytic and noncytolytic functions. We found that the latent reservoir significantly contributes to the size of the peak VL after CD8 depletion, while drug efficacy plays a lesser role. Our model suggests that the overall CD8 + lymphocyte cytolytic killing rate is dynamically changing depending on the levels of antigen-induced effector cell activation and exhaustion. Based on estimated parameters, our model suggests that before ART or without ART the overall CD8 cytolytic killing rate is small due to exhaustion. However, after the start of ART, the overall CD8 cytolytic killing rate increases due to an expansion of SIV-specific CD8 effector cells. Further, we estimate that the cytolytic killing rate can be significantly larger than the cytopathic death rate in some animals during the second phase of ART-induced viral decay. Lastly, our model provides a new explanation for the puzzling findings by Klatt et al. and Wong et al. that CD8 depletion done immediately before ART has no noticeable effect on the first phase viral decay slope seen after ART initiation Overall, by incorporating effector cells and their exhaustion, our model can explain the effects of CD8 depletion on VL during ART, reveals a detailed dynamic role of CD8 + lymphocytes in controlling viral infection, and provides a unified explanation for CD8 depletion experimental data.

59 BASIC BIOLOGICAL SCIENCES↗

A Model-Based Approach to Scintillator/Photomultiplier System Characterization

Scintillator/photomultiplier systems are prevalent in many instruments to measure a variety on nuclear events. When gamma energy from a nuclear event irradiates the scintillator, event radiation interacts with the scintillator material generating photons which are detected by the photomultiplier tube (PMT). The PMT photoelectrons or photocathode current is amplified in the PMT by a number of dynode stages producing a total charge out for a given radiation flux into scintillator. The resulting current is converted to a voltage when passed through a typical load impedance. This voltage is attenuated and provided as input to a digitizer through a variety of signal conditioners (filters, amplifiers, receivers) before analog-to-digital conversion. In this report, we concentrate on the development of analytic models for the scintillator/photomultiplier system and develop both a (1) statistical simulator; and (2) a model-based signal processor to extract the desired signal information while rejecting the noise and uncertainty. We develop the model analytically, transform it to the model-based framework (state-space), and then use actual experimental data to fit the model and apply it to our processing problem.

42 ENGINEERING↗

Deep galaxy counts in the K band with the Kech telescope

We present deep galaxy counts in the K (lambda 2.2 micrometer) band, obtained at the W. M. Kech 10 m telescope. The data reach limiting magnitudes K approximately 24 mag, about 5 times deeper than the deepest published K-band images to date. The counts are performed in three small (approximately 1 min), widely separated high-latitude fields. Extensive Monte Carlo tests were used to derive the comleteness corrections and minimize photometric biases. The counts continue to rise, with no sign of a turnover, down to the limits of our data, with the logarithmic slope of d log N/dm = 0.315 +/- 0.02 between K = 20 and 24 mag. This implies a cumulative surface density of approximately 5 x 10(exp 5) galaxies/sq deg, or approximately 2 x 10(exp 10) over the entire sky, down to K = 24 mag. Our counts are in good agreement with, although slightly lower than, those from the Hawaii Deep Survey by Cowie and collaborators; the discrepancies may be due to the small differences in the aperture corrections. We compare our counts with some of the available theoretical predictions. The data do not require models with a high value of Omega(sub 0), but can be well fitted by models with no (or little) evolution, and cosmologies with a low value of Omega(sub 0). Given the uncertainties in the models, it may be premature to put useful constrains on the value of Omega(sub 0) from the counts alone. Optical-to-IR colors are computed, using CCD data obtaind previously at Palomar. We find a few red galaxies with (r-K) approximately greater than 5 mag, or (i-K) approximately greater than 5 mag; these may be ellipticals at z approximately 1. While the redshift distribution of galaxies in our counts is still unknown, the flux limits reached would allow us to detect unobscured L(sub *) galaxies out to substantial redshifts (z greater than 3?).

Djorgovski, S.↗

Analytic Guided-Search Model of Human Performance Accuracy in Target- Localization Search Tasks

Current models of human visual search have extended the traditional serial/parallel search dichotomy. Two successful models for predicting human visual search are the Guided Search model and the Signal Detection Theory model. Although these models are inherently different, it has been difficult to compare them because the Guided Search model is designed to predict response time, while Signal Detection Theory models are designed to predict performance accuracy. Moreover, current implementations of the Guided Search model require the use of Monte-Carlo simulations, a method that makes fitting the model's performance quantitatively to human data more computationally time consuming. We have extended the Guided Search model to predict human accuracy in target-localization search tasks. We have also developed analytic expressions that simplify simulation of the model to the evaluation of a small set of equations using only three free parameters. This new implementation and extension of the Guided Search model will enable direct quantitative comparisons with human performance in target-localization search experiments and with the predictions of Signal Detection Theory and other search accuracy models.

Eckstein, Miguel P.↗

Shock-Stationary Application of Pseudoshock Models During High-Amplitude Combustion-Driven Unsteadiness

The isolator pseudo-shock provides necessary compression within a dual-mode scramjet engine and buffers the engine system against unstart. Quasi-1D flux-conserved models are the state-of-the-art reduced-order model for optimization and online control of dual-mode scramjet engines. The stability and efficacy of this modeling approach is evaluated against data from a combustor-driven direct-connect experiment. The experiment exhibited strong combustor-driven unsteadiness that produced upstream propagating weak shocks into the isolator, interacting with the pseudo-shock. While this configuration resulted in unsteadiness that is atypical of standard operation, the experiment provided an opportunity to evaluate the modeling techniques in highly transient states. Such transients could occur during maneuvering or result from unexpected combustor events. A flexible quasi-1D formulation, the Fievet flux-conserved model is fit using Bayesian inference in the laboratory and shock-stationary reference frames. Model performance is analyzed using the Bayesian posteriors and model evaluations over the measured shock train speed range. It is concluded that to produce consistent isolator pressure profile estimates in this unsteady environment, the model must be implemented in a shock-stationary reference frame. Implementing this conclusion in model-based engine controllers may reduce needed unstart safety margins and increase maximum performance.

Bayesian↗

Measurements of the Sunyaev–Zel’Dovich effect in MACS J0647.7+7015 and MACS J1206.2–0847 at high angular resolution with MUSTANG

We present high resolution (9″) imaging of the Sunyaev–Zel’dovich Effect (SZE) toward two massive galaxy clusters, MACS J0647.7+7015 (z = 0.591) and MACS J1206.2–0847 (z = 0.439). We compare these 90 GHz measurements, taken with the Multiplexed Squid/TES Array at Ninety Gigahertz (MUSTANG ) receiver on the Green Bank Telescope, with generalized Navarro–Frenk–White (gNFW) models derived from Bolocam 140 GHz SZE data as well as maps of the thermal gas derived from Chandra X-ray observations. We adopt a serial-fitting approach, in which gNFW models are first fit to the Bolocam data and then compared to the MUSTANG data to determine an overall best-fit model. For MACS J0647.7+7015, we find a gNFW profile with core slope parameter γ = 0.9 fits the MUSTANG image with χ{sub red}{sup 2}=1.005 and probability to exceed (PTE) = 0.34. For MACS J1206.2–0847, we find γ=0.7, χ{sub red}{sup 2}=0.993, and PTE = 0.70. In addition, we find a significant (>3σ) residual SZE feature in MACS J1206.2–0847 coincident with a group of galaxies identified in Very Large Telescope data and filamentary structure found in a weak-lensing mass reconstruction. We suggest the detected sub-structure may be the SZE decrement from a low mass foreground group or an infalling group. Giant Metrewave Radio Telescope measurements at 610 MHz reveal diffuse extended radio emission to the west, which we posit is either an active galactic nucleus-driven radio lobe, a bubble expanding away from disturbed gas associated with the SZE signal, or a bubble detached and perhaps re-accelerated by sloshing within the cluster. Using the spectroscopic redshifts available, we find evidence for a foreground (z = 0.423) or infalling group, coincident with the residual SZE feature.

79 ASTRONOMY AND ASTROPHYSICS↗

The Atacama Cosmology Telescope: Likelihood for Small-Scale CMB Data

The Atacama Cosmology Telescope has measured the angular power spectra of microwave fluctuations to arcminute scales at frequencies of 148 and 218 GHz, from three seasons of data. At small scales the fluctuations in the primordial Cosmic Microwave Background (CMB) become increasingly obscured by extragalactic foregounds and secondary CMB signals. We present results from a nine-parameter model describing these secondary effects, including the thermal and kinematic Sunyaev-Zel'dovich (tSZ and kSZ) power; the clustered and Poisson-like power from Cosmic Infrared Background (CIB) sources, and their frequency scaling; the tSZ-CIB correlation coefficient; the extragalactic radio source power; and thermal dust emission from Galactic cirrus in two different regions of the sky. In order to extract cosmological parameters, we describe a likelihood function for the ACT data, fitting this model to the multi-frequency spectra in the multipole range 500 < l < 10000. We extend the likelihood to include spectra from the South Pole Telescope at frequencies of 95, 150, and 220 GHz. Accounting for different radio source levels and Galactic cirrus emission, the same model provides an excellent fit to both datasets simultaneously, with χ2/dof= 675/697 for ACT, and 96/107 for SPT. We then use the multi-frequency likelihood to estimate the CMB power spectrum from ACT in bandpowers, marginalizing over the secondary parameters. This provides a simplified 'CMB-only' likelihood in the range 500 < l < 3500 for use in cosmological parameter estimation

Cosmology↗

Does the first chaotic inflation model in supergravity provide the best fit to the Planck data?

I describe the first model of chaotic inflation in supergravity, which was proposed by Goncharov and the present author in 1983. The inflaton potential of this model has a plateau-type behavior V{sub 0}(1−(8/3) e{sup −√6|ϕ|}) at large values of the inflaton field. This model predicts n{sub s}=1−(2/N)≈0.967 and r=(4/(3N{sup 2}))≈4×10{sup −4}, in good agreement with the Planck data. I propose a slight generalization of this model, which allows to describe not only inflation but also dark energy and supersymmetry breaking.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A viscous blast-wave model for high energy heavy-ion collisions

Employing a viscosity-based survival scale for initial geometrical perturbations formed in relativistic heavy-ion collisions, we model the radial flow velocity at freeze-out. Subsequently, we use the Cooper-Frye freeze-out prescription, with viscous corrections to the distribution function, to extract the transverse momentum dependence of particle yields and flow harmonics. We fit the model parameters for central collisions, by fitting the spectra of identified particles at the Large Hadron Collider (LHC), and estimate them for other centralities using simple hydrodynamic relations. We use the results of Monte Carlo Glauber model for initial eccentricities. We demonstrate that this improved viscous blast-wave model leads to good agreement with transverse momentum distribution of elliptic and triangular flow for all centralities and estimate the shear viscosity to entropy density ratio η/s ≃ 0.24 at the LHC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Full Modeling and parameter compression methods in configuration space for DESI 2024 and beyond

In the contemporary era of high-precision spectroscopic surveys, led by projects like DESI, there is an increasing demand for optimizing the extraction of cosmological information from clustering data. This work conducts a thorough comparison of various methodologies for modeling the full shape of the two-point statistics in configuration space. We investigate the performance of both direct fits (Full Modeling) and the parameter compression approaches (ShapeFit and Standard). We utilize the ABACUS-SUMMIT simulations, tailored to exceed DESI's precision requirements. Particularly, we fit the two-point statistics of three distinct tracers (LRG, ELG, and QSO), by employing a Gaussian Streaming Model in tandem with Convolution Lagrangian Perturbation Theory and Effective Field Theory. We explore methodological setup variations, including the range of scales, the set of galaxy bias parameters, the inclusion of the hexadecapole, as well as model extensions encompassing varying ns and allowing for w 0 w a CDM dark energy model. Throughout these varied explorations, while precision levels fluctuate and certain configurations exhibit tighter parameter constraints, our pipeline consistently recovers the parameter values of the mocks within 1σ in all cases for a 1-year DESI volume. Additionally, we compare the performance of configuration space analysis with its Fourier space counterpart using three models: PyBird, FOLPS and velocileptors, presented in companion papers. We find good agreement with the results from all these models.

79 ASTRONOMY AND ASTROPHYSICS↗

Observations of ultraviolet-bright stars in globular clusters with the Hopkins Ultraviolet Telescope

Two UV-bright stars in globular clusters, UV5 in NGC 1851 and vZ 1128 in M3, were observed with the Hopkins Ultraviolet Telescope (HUT) during the Astro-1 space shuttle mission in 1990 December. The stars' spectra show weak absorption features and no significant emission features other than well known geocoronal lines. Detailed fitting of Kurucz (1991) stellar atmosphere models using a chi(sup 2) minimization technique indicates T(sub eff) = 16 000 K, log g = 2.5, and abundance (-1.0) (logarithm of abundance of elements heavier than helium relative to solar) for UV5, and T(sub eff) = 35 000 K, log g = 4.0, and abundance (-3.5) for vZ 1128. The Kurucz model which best fits vZ 1128 overpredicts the flux in the region below approximately 1000 A, an effect seen in previous models of O-type stars. Our results are robust (to within approximately 1000 K) with respect to uncertainties in interstellar reddening and atomic and molecular hydrogen column densities. We do not see significant molecular hydrogen absorption, which might have indicated material in a circumstellar shell, in either star's spectrum. We estimate the stellar luminosities to be log L/solar luminosity = 3.33 +/- 0.15 for UV5 and log L/solar luminosity = 3.21 +/- 0.12 for vZ 1128. These atmospheric parameters place both stars on the Schoenberner (post-AGB) tracks, though the stellar masses derived from the best-fitting Kurucz models are somewhat less than those predicted by the Schoenberner models. Examination of individual absorption line strengths reveals no significant abundance anomalies in either star.

Van Dyke Dixon, W.↗