Search NASA⌕ Search

SEARCH · Search NASA

Results for “spectral 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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The CXSFIT spectral fitting code: Past, present and future

Magnetically confined plasma experiments generate a wealth of spectroscopic data. The first step toward extracting physical parameters is to fit a spectral model to the often complex spectra. The CXSFIT (Charge eXchange Spectroscopy FITting) spectral fitting code was originally developed for fitting charge exchange spectra on JET from the late 1980s onward and has been further developed over decades to keep up with the needs of the users. The primary use is to efficiently fit a large number of spectra with many constrained Gaussian spectral lines of which the physical parameters can be coupled in a user-friendly manner. More recent additions to the code include time-dependent couplings between parameters, flexible background subtraction, and a non-linear coupling scheme between fit parameters. The latter was a pre-requisite for implementing Zeeman and motional Stark effect multiplets in the library of spectral features. The ability to save and replay “fit recipes,” even when multiple iterations are required, has ensured the traceability of the results and is one of the keys to the longevity and success of the code. The code is also in use on other tokamaks (AUG, ST-40) and to fit data from other spectroscopic diagnostics on JET. In this paper, we document the current capabilities and philosophy behind the structure of the code, including some of the algorithms used to calculate spectral features numerically efficiently. We also provide an outline of how CXSFIT could be transferred into a framework that would be able to meet the spectral fitting requirements of future devices, such as ITER.

Delabie, Ephrem G.↗
DOE OSTI2024

SCEC/USGS Community Stress-Drop Validation Study: How Spectral Fitting Approaches Influence Measured Source Parameters

Spectral source parameters used to estimate an earthquake’s stress drop (⁠Δσ⁠) can vary significantly across measurement approaches. The Statewide California Earthquake Center/U.S. Geological Survey Community Stress‐Drop Validation Study was initiated to compare source parameter estimates, focusing initially on a dataset from the 2019 Ridgecrest earthquake sequence. As part of that validation effort, here we focus on one potential source of uncertainty: whether spectral fitting approaches alone, applied to a common set of spectra from the 2019 Ridgecrest sequence result in different source parameter estimates. By using a common set of benchmark spectra analyzed across a consistent frequency band of 1–40 Hz, we eliminate many sources of variability. A subgroup of validation study participants volunteered to estimate the low‐frequency displacement (⁠Ω 0 ⁠) and corner frequency (⁠ƒ c ⁠) by fitting a smooth function to benchmark displacement spectra. Participants used linear‐ or log‐sampled spectra, assumed a Brune or Boatwright spectral model, and applied different misfit criteria. We compare 17 approaches used to estimate ⁠Ω 0 ⁠, ƒ c ⁠, and Δσ for 54 earthquake spectra. Our results reveal that 35% of events have Δσ estimates within a factor of two, whereas others exhibit variations exceeding an order of magnitude. The variability in and can largely be attributed to whether a spectrum is consistent with the smooth function of an idealized simple crack model. The trade‐off between Ω 0 and ƒ c may be more pronounced when using linearly sampled spectra, as higher frequency spectral bumps control the fits. As expected, methods that assumed a Boatwright model tended to have lower Ω 0 and somewhat higher ƒ c compared to those assuming a Brune model, although resulting Δσ estimates are similar. Finally, when compared to the overall validation study results, the fitting approach alone may account for between 5% and 90% (25% on average) of the total variability in spectral Δσ⁠.

58 GEOSCIENCES↗
NASA NTRS1991

XSPEC: An x ray spectral fitting package. Version 2 of the user's guide

The user guide for XSPEC, a command driven, interactive, X-ray spectral fitting program is presented. It is designed to be completely detector independent so it can be used for any X-ray spectral instrument. An overview of the program commands and a walk through of an XSPEC session is presented. Individual commands and descriptions of the spectral models are given. For the more experienced user details of some of the programs associated with XSPEC, the command parser, and the addition of models to XSPEC are included. The PLT plotting package used by XSPEC is described together with the maximum likelihood option for users possessing X-ray spectra with small numbers of counts per bin. Details on the various files used by XSPEC and on the general file structure are also given.

Shafer, R. A.↗
DOE OSTI2022

Semi-analytic spectral fitting: simultaneously modelling the mass accumulation and chemical evolution in MaNGA spiral galaxies

We develop a novel semi-analytic spectral fitting approach to quantify the star formation histories (SFHs) and chemical enrichment histories (ChEHs) of individual galaxies. We construct simple yet general chemical evolution models that account for gas inflow and outflow processes as well as star formation, to investigate the evolution of merger-free star-forming systems. These models are fitted directly to galaxies’ absorption-line spectra, while their emission lines are used to constrain current gas phase metallicity and star formation rate. We apply this method to spiral galaxies selected from the SDSS-IV Mapping Nearby Galaxies at Apache Point Observatory survey. By fitting the co-added absorption-line spectra for each galaxy, and using the emission-line constraints on present-day metallicity and star formation, we reconstruct both the SFHs and the ChEHs for all objects in the sample. We can use these reconstructions to obtain archaeological measures of derived correlations such as the mass–metallicity relation at any redshift, which compare favourably with direct observations. We find that both the SFHs and ChEHs have strong mass dependence: massive galaxies accumulate their stellar masses and become enriched earlier. This mass dependence causes the observed flattening of the mass–metallicity relation at lower redshifts. The model also reproduces the observed gas-to-stellar mass ratio and its mass dependence. Moreover, we are able to determine that more massive galaxies have earlier gas infall times and shorter infall time-scales, and that the early chemical enrichment of low-mass galaxies is suppressed by strong outflows, while outflows are not very significant in massive galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2023

Use of Machine Learning and Principal Component Analysis to Retrieve Nitrogen Dioxide (NO 2 ) With Hyperspectral Imagers and Reduce Noise in Spectral Fitting

Nitrogen dioxide (NO 2 ) is an important trace-gas pollutant and climate agent whose presence also leads to spectral interference in ocean color retrievals. NO 2 column densities have been retrieved with satellite UV–Vis spectrometers such as the Ozone Monitoring Instrument (OMI) and the Tropospheric Monitoring Instrument (TROPOMI) that typically have spectral resolutions of the order of 0.5 nm or better and spatial footprints as small as 3.6 km × 5.6 km. These NO 2 observations are used to estimate emissions, monitor pollution trends, and study effects on human health. Here, we investigate whether it is possible to retrieve NO 2 amounts with lower-spectral-resolution hyperspectral imagers such as the Ocean Color Instrument (OCI) that will fly on the Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) satellite set for launch in early 2024. OCI will have a spectral resolution of 5 nm and a spatial resolution of ∼ 1 km with global coverage in 1–2 d. At this spectral resolution, small-scale spectral structure from NO 2 absorption is still present. We use real spectra from the OMI to simulate OCI spectra that are in turn used to estimate NO 2 slant column densities (SCDs) with an artificial neural network (NN) trained on target OMI retrievals. While we obtain good results with no noise added to the OCI simulated spectra, we find that the expected instrumental noise substantially degrades the OCI NO 2 retrievals. Nevertheless, the NO 2 information from OCI may be of value for ocean color retrievals. OCI retrievals can also be temporally averaged over timescales of the order of months to reduce noise and provide higher-spatial-resolution maps that may be useful for downscaling lower-spatial-resolution data provided by instruments such as OMI and TROPOMI; this downscaling could potentially enable higher-resolution emissions estimates and be useful for other applications. In addition, we show that NNs that use coefficients of leading modes of a principal component analysis of radiance spectra as inputs appear to enable noise reduction in NO 2 retrievals. Once trained, NNs can also substantially speed up NO 2 spectral fitting algorithms as applied to OMI, TROPOMI, and similar instruments that are flying or will soon fly in geostationary orbit.

NO2↗
NASA NTRS2019

FIGS: spectral fitting constraints on the star formation history of massive galaxies since the cosmic noon

We constrain the stellar population properties of a sample of 52 massive galaxies – with stellar mass log (M(s)/M(⊙)) ≳ 10.5 – over the redshift range 0.5 < z < 2 by use of observer-frame optical and near-infrared slitless spectra from Hubble Space Telescope’s ACS and WFC3 grisms. The deep exposures (∼100 ks) allow us to target individual spectra of massive galaxies to F160W = 22.5 AB. Our spectral fitting approach uses a set of six base models adapted to the redshift and spectral resolution of each observation, and fits the weights of the base models, including potential dust attenuation, via a Markov Chain Monte Carlo method. Our sample comprises a mixed distribution of quiescent (19) and star-forming galaxies (33). We quantify the width of the age distribution (Δt) that is found to dominate the variance of the retrieved parameters according to principal component analysis. The population parameters follow the expected trend towards older ages with increasing mass, and Δt appears to weakly anticorrelate with stellar mass, suggesting a more efficient star formation at the massive end. As expected, the redshift dependence of the relative stellar age (measured in units of the age of the Universe at the source) in the quiescent sample rejects the hypothesis of a single burst (aka monolithic collapse). Radial colour gradients within each galaxy are also explored, finding a wider scatter in the star-forming subsample, but no conclusive trend with respect to the population parameters.

Ignacio Ferreras↗
NASA NTRS2015

Intensity Conserving Spectral Fitting

The detailed shapes of spectral line profiles provide valuable information about the emitting plasma, especially when the plasma contains an unresolved mixture of velocities, temperatures, and densities. As a result of finite spectral resolution, the intensity measured by a spectrometer is the average intensity across a wavelength bin of non-zero size. It is assigned to the wavelength position at the center of the bin. However, the actual intensity at that discrete position will be different if the profile is curved, as it invariably is. Standard fitting routines (spline, Gaussian, etc.) do not account for this difference, and this can result in significant errors when making sensitive measurements. Detection of asymmetries in solar coronal emission lines is one example. Removal of line blends is another. We have developed an iterative procedure that corrects for this effect. It can be used with any fitting function, but we employ a cubic spline in a new analysis routine called Intensity Conserving Spline Interpolation (ICSI). As the name implies, it conserves the observed intensity within each wavelength bin, which ordinary fits do not. Given the rapid convergence, speed of computation, and ease of use, we suggest that ICSI be made a standard component of the processing pipeline for spectroscopic data.

Klimchuk, J. A.↗
DOE OSTI2021

Particle transport constraints via Bayesian spectral fitting of multiple atomic lines

Optimized operation of fusion devices demands detailed understanding of plasma transport, a problem that must be addressed with advances in both measurement and data analysis techniques. In this work, we adopt Bayesian inference methods to determine experimental particle transport, leveraging opportunities from high-resolution He-like ion spectra in a tokamak plasma. The Bayesian spectral fitting code is used to analyze resonance (w), forbidden (z), intercombination (x, y), and satellite (k, j) lines of He-like Ca following laser blow-off injections on Alcator C-Mod. This offers powerful transport constraints since these lines depend differently on electron temperature and density, but also differ in their relation to Li-like, He-like, and H-like ion densities, often the dominant Ca charge states over most of the C-Mod plasma radius. Using synthetic diagnostics based on the AURORA package, we demonstrate improved effectiveness of impurity transport inferences when spectroscopic data from a progressively larger number of lines are included.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Particle transport constraints via Bayesian spectral fitting of multiple atomic lines

Optimized operation of fusion devices demands detailed understanding of plasma transport, a problem that must be addressed with advances in both measurement and data analysis techniques. In this work, we adopt Bayesian inference methods to determine experimental particle transport, leveraging opportunities from high-resolution He-like ion spectra in a tokamak plasma. The Bayesian spectral fitting code is used to analyze resonance (w), forbidden (z), intercombination (x, y), and satellite (k, j) lines of He-like Ca following laser blow-off injections on Alcator C-Mod. This offers powerful transport constraints since these lines depend differently on electron temperature and density, but also differ in their relation to Li-like, He-like, and H-like ion densities, often the dominant Ca charge states over most of the C-Mod plasma radius. Using synthetic diagnostics based on the AURORA package, we demonstrate improved effectiveness of impurity transport inferences when spectroscopic data from a progressively larger number of lines are included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2018

A Cloud Algorithm Based on the O2-O2 477 nm Absorption Band Featuring an Advanced Spectral Fitting Method and the Use of Surface Geometry-Dependent Lambertian-Equivalent Reflectivity

We discuss a new cloud algorithm that retrieves an effective cloud pressure, also known as cloud optical centroid pressure (OCP), from oxygen dimer (O2-O2) absorption at 477nm after determining an effective cloud fraction (ECF) at 466nm, a wavelength not significantly affected by trace-gas absorption and rotational Raman scattering. The retrieved cloud products are intended for use as inputs to the operational nitrogen dioxide (NO2) retrieval algorithm for the Ozone Monitoring Instrument (OMI) flying on the Aura satellite. The cloud algorithm uses temperature-dependent O2-O2 cross sections and incorporates flexible spectral fitting techniques that account for specifics of the surface reflectivity. The fitting procedure derives O2-O2 slant column densities (SCDs) from radiances after O3, NO2, and H2O absorption features have been removed based on estimates of the amounts of these species from independent OMI algorithms. The cloud algorithm is based on the frequently used mixed Lambertian-equivalent reflectivity (MLER) concept. A geometry-dependent Lambertian-equivalent reflectivity (GLER), which is a proxy of surface bidirectional reflectance, is used for the ground reflectivity in our implementation of the MLER approach. The OCP is derived from a match of the measured O2-O2 SCD to that calculated with the MLER method. Temperature profiles needed for computation of vertical column densities are taken from the Global Modeling Initiative (GMI) model. We investigate the effect of using GLER instead of climatological LER on the retrieved ECF and OCP. For evaluation purposes, the retrieved ECFs and OCPs are compared with those from the operational OMI cloud product, which is also based on the same O2-O2 absorption band. Impacts of the application of the newly developed cloud algorithm to the OMI NO2 retrieval are discussed.

cloud fraction and pressure↗
NASA NTRS1992

The X-ray spectra of galaxies. I - Spectral fits of individual galaxies and X-ray colors

The X-ray spectra of normal galaxies is investigated systematically by using the Imaging Proportional Counter (IPC) data in the Einstein data bank. In addition to the standard model fitting technique, X-ray colors are introduced in order to extract spectral information from the fainter X-ray sources. Spectral parameters for 43 galaxies and X-ray colors for 127 galaxies are presented.

Kim, D.-W.↗
DOE OSTI2022

Hybrid time–frequency domain dual-probe coherent anti-Stokes Raman scattering for simultaneous temperature and pressure measurements in compressible flows via spectral fitting

We demonstrate a hybrid time–frequency spectroscopic method for simultaneous temperature/pressure measurements in nonreacting compressible flows with known gas composition. Hybrid femtosecond–picosecond, pure-rotational coherent anti-Stokes Raman scattering (CARS), with two independent, time-delayed probe pulses, is deployed for single-laser-shot measurements of temperature and pressure profiles along an ∼ <#comment/> 5 -mm line. The theory of dual-probe CARS is presented, along with a discussion of the iterative fitting of experimental spectra. Temperature is obtained from spectra acquired with an early, near-collision-free probe time delay ( τ <#comment/> 1 = 0 p s ) and pressure from spectra obtained at probe delays of τ <#comment/> 2 = 150 − <#comment/> 1000 p s , where collisions significantly impact the spectral profile. Unique solutions for temperature and pressure are obtained by iteratively fitting the two spectra to account for small collisional effects observed for the near zero probe delay spectrum. A dual-probe pure-rotational CARS system, in a 1D line-imaging configuration, is developed to demonstrate effectively the simultaneous temperature and pressure profiles recorded along the axial centerline of a highly underexpanded jet. The underexpanded air jet permits evaluation of this hybrid time–frequency domain approach for temperature and pressure measurements across a wide range of low-temperature–low-pressure conditions of interest in supersonic ground-test facilities. Single-laser-shot measurement precisions in both quantities and pressure measurement accuracy are systematically evaluated in the quiet zone upstream of the Mach disk. Precise thermometry approaching 1%−2% is observed in regions of high CARS signal-to-noise ratios. Pressure measurements are optimized at probe time delays where the ratio of the late probe delay to the Raman lifetime exceeds four ( τ <#comment/> 2 / τ <#comment/> R > <#comment/> 4 ). The impact of low-temperature Raman linewidths on CARS pressure measurements is evaluated, and comparisons of CARS pressures obtained with our recent low-temperature pure-rotational Raman linewidth data and extrapolated high-temperature Q -branch linewidths are presented. Considering all measurements with τ <#comment/> 2 / τ <#comment/> R ≥ <#comment/> 4.0 , measured pressures were on average 7.9% of the computed isentropic values with average shot-to-shot deviations representing a combination of instrument noise and fluid fluctuations of 5.0%.

Retter, Jonathan E.↗
NASA NTRS1991

Spectral fitting, shock layer modeling, and production of nitrogen oxides and excited nitrogen

An analysis was made of N2 emission from 8.72 MJ/kg shock layer at 2.54, 1.91, and 1.27 cm positions and vibrational state distributions, temperatures, and relative electronic state populations was obtained from data sets. Other recorded arc jet N2 and air spectral data were reviewed and NO emission characteristics were studied. A review of operational procedures of the DSMC code was made. Information on other appropriate codes and modifications, including ionization, were made as well as a determination of the applicability of codes reviewed to task requirement. A review was also made of computational procedures used in CFD codes of Li and other codes on JSC computers. An analysis was made of problems associated with integration of specific chemical kinetics applicable to task into CFD codes.

Blackwell, H. E.↗
NASA NTRS1994

Soft X-ray spectral fits of Geminga with model neutron star atmospheres

The spectrum of the soft X-ray pulsar Geminga consists of two components, a softer one which can be interpreted as thermal-like radiation from the surface of the neutron star, and a harder one interpreted as radiation from a polar cap heated by relativistic particles. We have fitted the soft spectrum using a detailed magnetized hydrogen atmosphere model. The fitting parameters are the hydrogen column density, the effective temperature T(sub eff), the gravitational redshift z, and the distance to radius ratio, for different values of the magnetic field B. The best fits for this model are obtained when B less than or approximately 1 x 10(exp 12) G and z lies on the upper boundary of the explored range (z = 0.45). The values of T(sub eff) approximately = (2-3) x 10(exp 5) K are a factor of 2-3 times lower than the value of T(sub eff) obtained for blackbody fits with the same z. The lower T(sub eff) increases the compatibility with some proposed schemes for fast neutrino cooling of neutron stars (NSs) by the direct Urca process or by exotic matter, but conventional cooling cannot be excluded. The hydrogen atmosphere fits also imply a smaller distance to Geminga than that inferred from a blackbody fit. An accurate evaluation of the distance would require a better knowledge of the ROSAT Position Sensitive Proportional Counter (PSPC) response to the low-energy region of the incident spectrum. Our modeling of the soft component with a cooler magnetized atmosphere also implies that the hard-component fit requires a characteristic temperature which is higher (by a factor of approximately 2-3) and a surface area which is smaller (by a factor of 10(exp 3), compared to previous blackbody fits.

Meyer, R. D.↗
NASA NTRS2010

CARS Spectral Fitting with Multiple Resonant Species using Sparse Libraries

The dual pump CARS technique is often used in the study of turbulent flames. Fast and accurate algorithms are needed for fitting dual-pump CARS spectra for temperature and multiple chemical species. This paper describes the development of such an algorithm. The algorithm employs sparse libraries, whose size grows much more slowly with number of species than a conventional library. The method was demonstrated by fitting synthetic "experimental" spectra containing 4 resonant species (N2, O2, H2 and CO2), both with noise and without it, and by fitting experimental spectra from a H2-air flame produced by a Hencken burner. In both studies, weighted least squares fitting of signal, as opposed to least squares fitting signal or square-root signal, was shown to produce the least random error and minimize bias error in the fitted parameters.

Cutler, Andrew D.↗
NASA NTRS2013

An Archival Chandra and XMM-Newton Survey of Type 2 Quasars

In order to investigate obscuration in high-luminosity type 2 active galactic nuclei (AGNs), we analyzed Chandra and XMM-Newton archival observations for 71 type 2 quasars detected at 0.05 < z < 0.73, which were selected based on their [O III] lambda5007 emission lines. For 54 objects with good spectral fits, the observed hard X-ray luminosity ranges from 2 × 10(exp 41) to 5.3 × 10(exp 44) erg s(exp −1), with a median of 1.1 × 10(exp 43) erg s(exp −1). We find that the means of the column density and photon index of our sample are log N(sub H) = 22.9 cm(exp −2) and gamma = 1.87, respectively. From simulations using a more physically realistic model, we find that the absorbing column density estimates based on simple power-law models significantly underestimate the actual absorption in approximately half of the sources. Eleven sources show a prominent Fe K alpha emission line (EW>100 eV in the rest frame) and we detect this line in the other sources through a joint fit (spectral stacking). The correlation between the Fe K alpha and [O III] fluxes and the inverse correlation of the equivalent width of the Fe Kα line with the ratio of hard X-ray and [O III] fluxes is consistent with previous results for lower luminosity Seyfert 2 galaxies. We conclude that obscuration is the cause of the weak hard X-ray emission rather than intrinsically low X-ray luminosities. We find that about half of the population of optically selected type 2 quasars are likely to be Compton thick. We also find no evidence that the amount of X-ray obscuration depends on the AGN luminosity (over a range of more than three orders of magnitude in luminosity).

quasars: general↗