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 865 records · Page 48

SPITZER IMAGING OF STRONGLY LENSED HERSCHEL-SELECTED DUSTY STAR-FORMING GALAXIES

We present the rest-frame optical spectral energy distribution (SED) and stellar masses of six Herschel-selected gravitationally lensed dusty, star-forming galaxies (DSFGs) at 1 < z < 3. These galaxies were first identified with Herschel/SPIRE imaging data from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). The targets were observed with Spitzer/IRAC at 3.6 and 4.5 μm. Due to the spatial resolution of the IRAC observations at the level of 2″, the lensing features of a background DSFG in the near-infrared are blended with the flux from the foreground lensing galaxy in the IRAC imaging data. We make use of higher resolution Hubble/WFC3 or Keck/NIRC2 Adaptive Optics imaging data to fit light profiles of the foreground lensing galaxy (or galaxies) as a way to model the foreground components, in order to successfully disentangle the foreground lens and background source flux densities in the IRAC images. The flux density measurements at 3.6 and 4.5 μm, once combined with Hubble/WFC3 and Keck/NIRC2 data, provide important constraints on the rest-frame optical SED of the Herschel-selected lensed DSFGs. We model the combined UV- to millimeter-wavelength SEDs to establish the stellar mass, dust mass, star formation rate, visual extinction, and other parameters for each of these Herschel-selected DSFGs. These systems have inferred stellar masses in the range 8 × 10{sup 10}–4 × 10{sup 11} M{sub ⊙} and star formation rates of around 100 M{sub ⊙} yr{sup −1}. This puts these lensed submillimeter systems well above the SFR-M* relation observed for normal star-forming galaxies at similar redshifts. The high values of SFR inferred for these systems are consistent with a major merger-driven scenario for star formation.

79 ASTRONOMY AND ASTROPHYSICS↗

Improved Fitting of Solution X-ray Scattering Data to Macromolecular Structures and Structural Ensembles by Explicit Water Modeling

A new procedure, AXES, is introduced for fitting small-angle X-ray scattering (SAXS) data to macromolecular structures and ensembles of structures. By using explicit water models to account for the effect of solvent, and by restricting the adjustable fitting parameters to those that dominate experimental uncertainties, including sample/buffer rescaling, detector dark current, and, within a narrow range, hydration layer density, superior fits between experimental high resolution structures and SAXS data are obtained. AXES results are found to be more discriminating than standard Crysol fitting of SAXS data when evaluating poorly or incorrectly modeled protein structures. AXES results for ensembles of structures previously generated for ubiquitin show improved fits over fitting of the individual members of these ensembles, indicating these ensembles capture the dynamic behavior of proteins in solution.

59 BASIC BIOLOGICAL SCIENCES↗

Constraints on Cosmological Parameters from the Angular Power Spectrum of a Combined 2500 deg 2 SPT-SZ and Planck Gravitational Lensing Map

Here, we report constraints on cosmological parameters from the angular power spectrum of a cosmic microwave background (CMB) gravitational lensing potential map created using temperature data from 2500 deg$^2$ of South Pole Telescope (SPT) data supplemented with data from Planck in the same sky region, with the statistical power in the combined map primarily from the SPT data. We fit the corresponding lensing angular power spectrum to a model including cold dark matter and a cosmological constant ($\Lambda$CDM), and to models with single-parameter extensions to $\Lambda$CDM. We find constraints that are comparable to and consistent with constraints found using the full-sky Planck CMB lensing data. Specifically, we find $\sigma_8 \Omega_{\rm m}^{0.25}=0.598 \pm 0.024$ from the lensing data alone with relatively weak priors placed on the other $\Lambda$CDM parameters. In combination with primary CMB data from Planck, we explore single-parameter extensions to the $\Lambda$CDM model. We find $\Omega_k = -0.012^{+0.021}_{-0.023}$ or $M_{\nu}< 0.70$eV both at 95% confidence, all in good agreement with results that include the lensing potential as measured by Planck over the full sky. We include two independent free parameters that scale the effect of lensing on the CMB: $A_{L}$, which scales the lensing power spectrum in both the lens reconstruction power and in the smearing of the acoustic peaks, and $A^{\phi \phi}$, which scales only the amplitude of the CMB lensing reconstruction power spectrum. We find $A^{\phi \phi} \times A_{L} =1.01 \pm 0.08$ for the lensing map made from combined SPT and Planck temperature data, indicating that the amount of lensing is in excellent agreement with what is expected from the observed CMB angular power spectrum when not including the information from smearing of the acoustic peaks.

79 ASTRONOMY AND ASTROPHYSICS↗

Ring Current Development During Storm Main Phase

The development of the ring current ions in the inner magnetosphere during the main phase of a magnetic storm is studied. The temporal and spatial evolution of the ion phase space densities in a dipole field are calculated using a three dimensional ring current model, considering charge exchange and Coulomb losses along drift paths. The simulation starts with a quiet time distribution. The model is tested by comparing calculated ion fluxes with Active Magnetospheric Particle Tracer Explorers/CCE measurement during the storm main phase on May 2, 1986. Most of the calculated omnidirectional fluxes are in good agreement with the data except on the dayside inner edge (L less than 2.5) of the ring current, where the ion fluxes are underestimated. The model also reproduces the measured pitch angle distributions of ions with energies below 10 keV. At higher energy, an additional diffusion in pitch angle is necessary in order to fit the data. The role of the induced electric field on the ring current dynamics is also examined by simulating a series of substorm activities represented by stretching and collapsing the magnetic field lines. In response to the impulsively changing fields, the calculated ion energy content fluctuates about a mean value that grows steadily with the enhanced quiescent field.

Fok, Mei-Ching↗

PSR J0007+7303 in the CTA1 SNR: New Gamma-ray Results from Two Years of Fermi-LAT Observations

One of the main results of the Fermi Gamma-Ray Space Telescope is the discovery of gamma-ray selected pulsars. The high magnetic field pulsar, PSR J0007+7303 in CTA1, was the first ever to be discovered through its gamma-ray pulsations. Based on analysis of 2 years of LAT survey data, we report on the discovery of I-ray emission in the off-pulse phase interval at the approx. 6sigma level. The flux from this emission in the energy range E > or =::: 100 MeV is F(sub 100) = (1.73+/-0.40) x 10(exp -8) photons/sq cm/s and is best fitted by a power law with a photon index of Gamma = 2.54+/-0.14. The pulsed gamma-ray flux in the same energy range is F(sub 100) = (3.95+/-0.07) x 10(exp -7) photons/sq cm/s and is best fitted by an exponentially-cutoff power-law spectrum with a photon index of Gamma = 1.41+/-0.23 and a cutoff energy E(sub c) = 4.04+/-0.20 GeV. We find no flux variability neither at the 2009 May glitch nor in the long term behavior. We model the gamma-ray light curve with two high-altitude emission models, the outer gap and slot gap, and find that the model that best fits the data depends strongly on the assumed origin of the off-pulse emission. Both models favor a large angle between the magnetic axis and observer line of sight, consistent with the nondetection of radio emission being a geometrical effect. Finally we discuss how the LAT results bear on the understanding of the cooling of this neutron star.

Abdo, A. A.↗

GPS measurements of strain accumulation across the Imperial Valley, California: 1986-1989

The Global Positioning System (GPS) data collected in southern California from 1986 to 1989 indicate considerable strain accumulation across the Imperial Valley. Displacements are computed at 29 stations in and near the valley from 1986 to 1988, and at 11 sites from 1988 to 1989. The earlier measurements indicate 5.9 +/- 1.0 cm/yr right-lateral differential velocity across the valley, although the data are heavily influenced by the 1987 Superstition Hills earthquake sequence. Some measurements, especially the east-trending displacements, are suspects for large errors. The 1988 to 1989 GPS displacements are best modeled by 5.2 +/- 0.9 cm/yr of valley crossing deformation, but rates calculated from conventional geodetic measurements (3.4 to 4.3 cm/yr) fit the data nearly as well. There is evidence from GPS and Very Long Base Interferometry (VLBI) observations that the present slip rate along the southern San Andreas fault is smaller than the long-term geologic estimate, suggesting a lower earthquake potential than is currently assumed. Correspondingly, a higher earthquake potential is indicated for the San Jacinto fault. The Imperial Valley GPS sites form part of a 183 station network in southern California and northern Baja California, which spans a cross-section of the North American-Pacific plate boundary.

Larsen, Shawn↗

Global Positioning System measurements of strain accumulation across the Imperial Valley, California - 1986-1989

The Global Positioning System (GPS) data collected in southern California from 1986 to 1989 indicate considerable strain accumulation across the Imperial Valley. Displacements are computed at 29 stations in and near the valley from 1986 to 1988, and at 11 sites from 1988 to 1989. The earlier measurements indicate 5.9 =/- 1.0 cm/yr right-lateral differential velocity across the valley, although the data are heavily influenced by the 1987 Superstition Hills earthquake sequence. Some measurements, especially the east-trending displacements, are suspects for large errors. The 1988 to 1989 GPS displacements are best modeled by 5.2 =/- 0.9 cm/yr of valley crossing deformation, but rates calculated from conventional geodetic measurements (3.4 to 4.3 cm/yr) fit the data nearly as well. There is evidence from GPS and Very Long Base Interferometry (VLBI) observations that the present slip rate along the southern San Andreas fault is smaller than the long-term geologic estimate, suggesting a lower earthquake potential than is currently assumed. Correspondingly, a higher earthquake potential is indicated for the San Jacinto fault. The Imperial Valley GPS sites form part of a 183 station network in southern California and northern Baja California, which spans a cross-section of the North American-Pacific plate boundary.

Larsen, Shawn↗

Predicting Human Cooperation

The Prisoner’s Dilemma has been a subject of extensive research due to its importance in understanding the ever-present tension between individual self-interest and social benefit. A strictly dominant strategy in a Prisoner’s Dilemma (defection), when played by both players, is mutually harmful. Repetition of the Prisoner’s Dilemma can give rise to cooperation as an equilibrium, but defection is as well, and this ambiguity is difficult to resolve. The numerous behavioral experiments investigating the Prisoner’s Dilemma highlight that players often cooperate, but the level of cooperation varies significantly with the specifics of the experimental predicament. We present the first computational model of human behavior in repeated Prisoner’s Dilemma games that unifies the diversity of experimental observations in a systematic and quantitatively reliable manner. Our model relies on data we integrated from many experiments, comprising 168,386 individual decisions. The model is composed of two pieces: the first predicts the first-period action using solely the structural game parameters, while the second predicts dynamic actions using both game parameters and history of play. Our model is successful not merely at fitting the data, but in predicting behavior at multiple scales in experimental designs not used for calibration, using only information about the game structure. As a result, we demonstrate the power of our approach through a simulation analysis revealing how to best promote human cooperation.

59 BASIC BIOLOGICAL SCIENCES↗

Characterizing Quantum Devices Using Model Selection

Quantum information processors are a new kind of computing device with the potential to solve certain problems very quickly. The quantum bits, or "ʺqubits,"ʺ that enable such processors are being developed at Sandia and around the world. In support of this effort, Sandia’s Center for Computing Research is designing methods and algorithms to characterize and debug the behavior of experimental qubits, such as trapped ions and electron spins in silicon. This is a challenge, because qubits are hard to characterize for exactly the same reason that they are powerful. Their behavior (and misbehavior) is far richer and more complicated than that of classical bits! The most reliable models of quantum errors have infinitely many parameters, and we never have enough data to fit all of them. Until now, researchers simply threw out most of these parameters resulting in a crude approximation that makes the problem tractable, but isn’t reliable or accurate enough for Sandia’s experimental work. Our new techniques use statistical model selection to identify “significant” error parameters on the fly, applying a data-­‐‑driven information criterion to choose a Hilbert space that describes the data well without overfitting. To find a reliable criterion, we needed to evaluate many candidates on an enormous corpus of test cases–for which we relied on Red Sky’s ability to run parallel Monte Carlo simulations. Our computations demonstrate how classic model selection methods have to be modified for quantum devices. These results will enable new techniques that identify the effective dimension of as-­‐‑built qubits on the fly, and then use this information to reliably pin down error behaviors and produce successively better and more useful qubits. This project highlights the close and productive relationship between Sandia and the University of New Mexico, where Travis Scholten is working towards a Ph.D. in physics. Both the intellectual environment of UNM and the high-­‐‑powered supercomputing available at Sandia were critical ingredients in achieving our research goals.

97 MATHEMATICS AND COMPUTING↗

Inverting infrasound data for the seismoacoustic source time functions and surface spall at the Source Physics Experiments Phase II: Dry Alluvium Geology

This report presents the infrasound data recorded as part of the Source Physics Experiment - Phase 2, Dry Alluvium Geology. This experiment, also known colloquially as DAG, consisted of four underground chemical explosions at the Nevada National Security Site. We focus our analysis on only the fourth explosion (DAG-4) as we determined that this was the only event that produced clear source-generated infrasound energy as recorded by the DAG sensors. We analyze the data using two inversion methods. The first method is designed to estimate the point-source seismoacoustic source time functions, and the second inversion method is designed to estimate the first-order characteristics (e.g. horizontal dimensions and maximum amplitude) of the actual spall surface. For both analysis methods, we are able to fit the data reasonably well, with various assumptions of the source model. The estimated seismoacoustic source appears to be a combination of a buried, isotropic explosion with a maximum amplitude of ~2 x 10 9 Nm and a vertically oriented force, applied to the Earth's surface with a maximum amplitude of 4 x 10 7 N. We use the vertically oriented force to simulate surface spall. The estimated spall surface has an approximate radius of ~40 m with a maximum acceleration magnitude in the range of 0.8 to 1.5 m/s/s. These estimates are approximately similar to the measured surface acceleration at the site.

58 GEOSCIENCES↗

Computer codes for the evaluation of thermodynamic properties, transport properties, and equilibrium constants of an 11-species air model

The computer codes developed provide data to 30000 K for the thermodynamic and transport properties of individual species and reaction rates for the prominent reactions occurring in an 11-species nonequilibrium air model. These properties and the reaction-rate data are computed through the use of curve-fit relations which are functions of temperature (and number density for the equilibrium constant). The curve fits were made using the most accurate data believed available. A detailed review and discussion of the sources and accuracy of the curve-fitted data used herein are given in NASA RP 1232.

Thompson, Richard A.↗

MB-Fit: Software infrastructure for data-driven many-body potential energy functions

Many-body potential energy functions (MB-PEFs), which integrate data-driven representations of many-body short-range quantum mechanical interactions with physics-based representations of many-body polarization and long-range interactions, have recently been shown to provide high accuracy in the description of molecular interactions from the gas to the condensed phase. Here, we present MB-Fit, a software infrastructure for the automated development of MB-PEFs for generic molecules within the TTM-nrg (Thole-type model energy) and MB-nrg (many-body energy) theoretical frameworks. Besides providing all the necessary computational tools for generating TTM-nrg and MB-nrg PEFs, MB-Fit provides a seamless interface with the MBX software, a many-body energy and force calculator for computer simulations. Given the demonstrated accuracy of the MB-PEFs, particularly within the MB-nrg framework, we believe that MB-Fit will enable routine predictive computer simulations of generic (small) molecules in the gas, liquid, and solid phases, including, but not limited to, the modeling of quantum isomeric equilibria in molecular clusters, solvation processes, molecular crystals, and phase diagrams.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

HERSCHEL PACS OBSERVATIONS AND MODELING OF DEBRIS DISKS IN THE TUCANA-HOROLOGIUM ASSOCIATION

We present Herschel PACS photometry of 17 B- to M-type stars in the 30 Myr old Tucana-Horologium Association. This work is part of the Herschel Open Time Key Programme 'Gas in Protoplanetary Systems'. 6 of the 17 targets were found to have infrared excesses significantly greater than the expected stellar IR fluxes, including a previously unknown disk around HD30051. These six debris disks were fitted with single-temperature blackbody models to estimate the temperatures and abundances of the dust in the systems. For the five stars that show excess emission in the Herschel PACS photometry and also have Spitzer IRS spectra, we fit the data with models of optically thin debris disks with realistic grain properties in order to better estimate the disk parameters. The model is determined by a set of six parameters: surface density index, grain size distribution index, minimum and maximum grain sizes, and the inner and outer radii of the disk. The best-fitting parameters give us constraints on the geometry of the dust in these systems, as well as lower limits to the total dust masses. The HD105 disk was further constrained by fitting marginally resolved PACS 70 {mu}m imaging.

79 ASTRONOMY AND ASTROPHYSICS↗

Reanalysis of x ray emission from M87 II: Multiphase models

In a previous paper by the authors (hereafter TB), it was shown that the hot intracluster gas around M87 could not be adequately explained in terms of a spherically symmetric, single phase model. It was found that although data from the Einstein satellite High Resolution Imager (HRI), Imaging Proportional Counter (IPC), and Focal Point Crystal Spectrometer (FPCS) could indeed be simultaneously explained by a single phase model, data from the Solid State Spectrometer (SSS) and optically determined mass estimates could not be similarly explained. A qualitative discussion was then given to indicate how the adoption of a specific multiphase model could plausibly improve upon the single phase model, but a detailed discussion was delayed until now. The x ray data is reanalyzed for the gas around M87 in the same spirit as the work in TB. That is, by assuming several cooling flow models, the shortcomings are improved upon. The procedure of TB was adopted and the data reanalyzed in the following manner: a multiphase model is assumed which can compute the surface brightness as seen by HRI and IPC. Also, the line fluxes can be computed of the lines seen by the FPCS and the mass profile after assuming hydrostatic equilibrium. The parameters of the model are then adjusted to fit the data. A check is then made to see whether the resulting model is consistent with the SSS spectrum and the equivalent width of the 7 keV complex of Fe lines as seen by large field of view instruments.

Tsai, John C.↗

MAGIC gamma-ray and multi-frequency observations of flat spectrum radio quasar PKS 1510-089 in early 2012

Aims. We report amongst more than fifty blazars detected in very high energy (VHE, E> 100 GeV) γ rays, only three belong to the subclass of flat spectrum radio quasars (FSRQs). The detection of FSRQs in the VHE range is challenging, mainly because of their soft spectra in the GeV-TeV regime. MAGIC observed PKS 1510-089 (z = 0.36) starting 2012 February 3 until April 3 during a high activity state in the high energy (HE, E> 100 MeV) γ-ray band observed by AGILE and Fermi. MAGIC observations result in the detection of a source with significance of 6.0 standard deviations (σ). We study the multi-frequency behaviour of the source at the epoch of MAGIC observation, collecting quasi-simultaneous data at radio and optical (GASP-WEBT and F-Gamma collaborations, REM, Steward, Perkins, Liverpool, OVRO, and VLBA telescopes), X-ray (Swift satellite), and HE γ-ray frequencies. Methods. We study the VHE γ-ray emission, together with the multi-frequency light curves, 43 GHz radio maps, and spectral energy distribution (SED) of the source. The quasi-simultaneous multi-frequency SED from the millimetre radio band to VHE γ rays is modelled with a one-zone inverse Compton model. We study two different origins of the seed photons for the inverse Compton scattering, namely the infrared torus and a slow sheath surrounding the jet around the Very Long Baseline Array (VLBA) core. Results. We find that the VHE γ-ray emission detected from PKS 1510-089 in 2012 February-April agrees with the previous VHE observations of the source from 2009 March-April. We find no statistically significant variability during the MAGIC observations on daily, weekly, or monthly time scales, while the other two known VHE FSRQs (3C 279 and PKS 1222+216) have shown daily scale to sub-hour variability. The γ-ray SED combining AGILE, Fermi and MAGIC data joins smoothly and shows no hint of a break. The multi-frequency light curves suggest a common origin for the millimetre radio and HE γ-ray emission, and the HE γ-ray flaring starts when the new component is ejected from the 43 GHz VLBA core and the studied SED models fit the data well. However, the fast HE γ-ray variability requires that within the modelled large emitting region, more compact regions must exist. Lastly, we suggest that these observed signatures would be most naturally explained by a turbulent plasma flowing at a relativistic speed down the jet and crossing a standing conical shock.

79 ASTRONOMY AND ASTROPHYSICS↗

Infrared Communications for Small Spacecraft: From a Wireless Bus to Cluster Concepts

Nanosatellites operating singly or in clusters are anticipated for future space science missions. To implement this new communications paradigm, we are approaching cluster communications by first developing an infrared (IR) intra-craft wireless bus capability, following initially the MIL-STD-1553B protocol. Benefits of an IR wireless bus are low mass, size, power, and cost, simplicity of implementation, ease of use, minimum EMI, and efficient and reliable data transfer. Our goals are to maximize the reliable link margin in order to afford greater flexibility in receiver placement, which will ease technology insertion. We have developed a concept demonstration using a high-speed visible-band silicon PIN photodiode and a high-efficiency visible LED operating at a data rate up to 4 Mb/sec. In designing an internal IR wireless bus, we have characterized various candidate materials, emitters, and geometries, assuming a single reflection. Thus, we have measured the bidirectional reflectance distribution function (BRDF) for five different materials characteristic of typical spacecraft structures, which range from nearly Lambertian to highly specular. We have fit our data to empirical BRDF functions and modeled the detected irradiance anywhere in the plane of incidence for a divergent (LED) emitter. We have also determined the angular limits on the link geometry to remain within the required bit error rate by determining the received signal-to-noise ratio (SNR) for minimum values of irradiance received at the detector.

Webb, Suzanne C.↗

Neutrino masses and cosmology with Lyman-alpha forest power spectrum

We present constraints on neutrino masses, the primordial fluctuation spectrum from inflation, and other parameters of the ΛCDM model, using the one-dimensional Lyα-forest power spectrum measured by [1] from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey (SDSS-III), complemented by Planck 2015 cosmic microwave background (CMB) data and other cosmological probes. This paper improves on the previous analysis by [2] by using a more powerful set of calibrating hydrodynamical simulations that reduces uncertainties associated with resolution and box size, by adopting a more flexible set of nuisance parameters for describing the evolution of the intergalactic medium, by including additional freedom to account for systematic uncertainties, and by using Planck 2015 constraints in place of Planck 2013. Fitting Lyα data alone leads to cosmological parameters in excellent agreement with the values derived independently from CMB data, except for a weak tension on the scalar index n{sub s}. Combining BOSS Lyα with Planck CMB constrains the sum of neutrino masses to ∑ m{sub ν} < 0.12 eV (95% C.L.) including all identified systematic uncertainties, tighter than our previous limit (0.15 eV) and more robust. Adding Lyα data to CMB data reduces the uncertainties on the optical depth to reionization τ, through the correlation of τ with σ{sub 8}. Similarly, correlations between cosmological parameters help in constraining the tensor-to-scalar ratio of primordial fluctuations r. The tension on n{sub s} can be accommodated by allowing for a running dn{sub s}/d ln k. Allowing running as a free parameter in the fits does not change the limit on ∑ m{sub ν}. We discuss possible interpretations of these results in the context of slow-roll inflation.

79 ASTRONOMY AND ASTROPHYSICS↗

Constrained curve fitting for semi-parametric models with radial basis function networks

Common to many analysis pipelines in lattice gauge theory and the broader scientific discipline is the need to fit a semi-parametric model to data. We propose a fit method that utilizes a radial basis function network to approximate the non-parametric component of such models. The approximate parametric model is fit to data using the basin hopping global optimization algorithm. Parameter constraints are enforced through Gaussian priors. The viability of our method is tested by examining its use in a finite-size scaling analysis of the $q$-state Potts model and $p$-state clock model with $q=2,3$ and $p=4,\infty$.

Peterson, Curtis T.↗