Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sampling error”

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 973 records · Page 54

Are wide pre-main-sequence binaries coeval?

We have observed a sample of 39 wide (projected separations 400-6000 AU) pre-main-sequence binary pairs spectroscopically and with optical and near-infrared images. The observations enable us to place 26 of the pairs in an H-R diagram and to determine masses and ages of the primary and secondary according to three sets of pre-main-sequence evolutionary tracks. In two-thirds of the cases the primary and secondary lie along the same isochrone to within the observational errors. However, real age differences appear for about one-third of our sample pairs-there is no set of nonintersecting theoretical isochrones that can make the primary and secondary have the same age for all pairs in our sample. In the cases where there are significant age differences between the component stars, the less massive star is usually younger than the more massive star. There is no correlation of the age differences with the presence or absence of accretion disks around the young stars. Hence, while disk accretion may affect the evolutionary tracks of the pre-main-sequence stars H-R diagram, we see no clear evidence of this effect among the pairs in our sample. The age differences also do not depend systematically on the apparent separation, the mass ratio, or the ages of stars.

Hartigan, Patrick↗

A simple soil mass correction for a more accurate determination of soil carbon stock changes

Abstract Agricultural soils can act as a sink for large quantities of soil organic carbon (SOC) but can also be sources of carbon to the atmosphere. The international standard for assessing SOC stock and measuring stock change stipulates fixed depth sampling to at least 30 cm. The tendency of bulk density (BD) to decrease with decreasing disturbance and increasing SOC concentration and the assumption of constant SOC and BD within this depth profile promotes error in the estimates of SOC stock. A hypothetical but realistic change in BD from 1.5 to 1.1 g cm −3 from successive fixed depth sampling to 30 cm underestimates SOC stock change by 17%. Significant effort has been made to evaluate and reduce this fixed depth error by using the equivalent soil mass (ESM) approach, but with limited adoption. We evaluate the error in SOC stock assessment and change generated from fixed depth measurements over time relative to the ESM approach and propose a correction that can be readily adopted under current sampling and analytical methods. Our approach provides a more accurate estimate of SOC stock accumulation or loss that will help incentivize management practice changes that reduce the environmental impacts of agriculture and further legitimize the accounting practices used by the emerging carbon market and organizations that have pledged to reduce their supply chain greenhouse gas (GHG) footprints.

60 APPLIED LIFE SCIENCES↗

Optical Imaging and Radiometric Modeling and Simulation

OPTOOL software is a general-purpose optical systems analysis tool that was developed to offer a solution to problems associated with computational programs written for the James Webb Space Telescope optical system. It integrates existing routines into coherent processes, and provides a structure with reusable capabilities that allow additional processes to be quickly developed and integrated. It has an extensive graphical user interface, which makes the tool more intuitive and friendly. OPTOOL is implemented using MATLAB with a Fourier optics-based approach for point spread function (PSF) calculations. It features parametric and Monte Carlo simulation capabilities, and uses a direct integration calculation to permit high spatial sampling of the PSF. Exit pupil optical path difference (OPD) maps can be generated using combinations of Zernike polynomials or shaped power spectral densities. The graphical user interface allows rapid creation of arbitrary pupil geometries, and entry of all other modeling parameters to support basic imaging and radiometric analyses. OPTOOL provides the capability to generate wavefront-error (WFE) maps for arbitrary grid sizes. These maps are 2D arrays containing digital sampled versions of functions ranging from Zernike polynomials to combination of sinusoidal wave functions in 2D, to functions generated from a spatial frequency power spectral distribution (PSD). It also can generate optical transfer functions (OTFs), which are incorporated into the PSF calculation. The user can specify radiometrics for the target and sky background, and key performance parameters for the instrument s focal plane array (FPA). This radiometric and detector model setup is fairly extensive, and includes parameters such as zodiacal background, thermal emission noise, read noise, and dark current. The setup also includes target spectral energy distribution as a function of wavelength for polychromatic sources, detector pixel size, and the FPA s charge diffusion modulation transfer function (MTF).

Ha, Kong Q.↗

Systematic softening in universal machine learning interatomic potentials

Machine learning interatomic potentials (MLIPs) have introduced a new paradigm for atomic simulations. Recent advancements have led to universal MLIPs (uMLIPs) that are pre-trained on diverse datasets, providing opportunities for universal force fields and foundational machine learning models. However, their performance in extrapolating to out-of-distribution complex atomic environments remains unclear. In this study, we highlight a consistent potential energy surface (PES) softening effect in three uMLIPs: M3GNet, CHGNet, and MACE-MP-0, which is characterized by energy and force underprediction in atomic-modeling benchmarks including surfaces, defects, solid-solution energetics, ion migration barriers, phonon vibration modes, and general high-energy states. The PES softening behavior originates primarily from the systematically underpredicted PES curvature, which derives from the biased sampling of near-equilibrium atomic arrangements in uMLIP pre-training datasets. Our findings suggest that a considerable fraction of uMLIP errors are highly systematic, and can therefore be efficiently corrected. We argue for the importance of a comprehensive materials dataset with improved PES sampling for next-generation foundational MLIPs.

36 MATERIALS SCIENCE↗

Sampling on NISQ Devices: "Who’s the Fairest One of All?"

Modern NISQ devices are subject to a variety of biases and sources of noise that degrade the solution quality of computations carried out on these devices. A natural question that arises in the NISQ era, is how fairly do these devices sample ground state solutions. To this end, we run five fair sampling problems (each with at least three ground state solutions) that are based both on quantum annealing and, on the Grover Mixer, -QAOA algorithm for gate-based NISQ hardware. In particular, we use seven IBM Q devices, the Aspen-9 Rigetti device, the IonQ device, and three D-Wave quantum annealers. For each of the fair sampling problems, we measure the ground state probability, the relative fairness of the frequency of each ground state solution with respect to the other ground state solutions, and the aggregate error as given by each hardware provider. Overall, our results show that NISQ devices do not achieve fair sampling yet. Furthermore, we also observe differences in the software stack with a particular focus on compilation techniques that illustrate what work will still need to be done to achieve a seamless integration of frontend (i.e., quantum circuit description) and backend compilation.

Computer Science↗

Magnetic Shield Design Modeling and Validation on SWOT Spacecraft Applied to Mars Flux Pinning Orbiting Sample Design

A modeling methodology was developed for use on the Surface Water and Ocean Topography (SWOT) mission for the computational modeling and design of a magnetic shield for a 63 A-sq m source. Shield options were modeled and tested across various design parameters and validated with measurement. Measurement results fell within 10% of simulation in most cases of concern with sources of error well understood. These methods results informed a subsequent modeling activity for a future Mars Sample Return (MSR) mission concept with stringent magnetic cleanliness requirements to retain geologic integrity of the samples. One option for capturing the sample from orbit requires rare earth magnets in close proximity to the Martian samples. Magnetic modeling with a finite-element method solver estimated the magnetic environment and established the need for shielding. Further modeling then determined the shielding necessary to meet magnetic requirements, followed by the design of a mass-optimized solution.

Gonzales, Edward↗

Is the galactic bulge devoid of planets?

We consider a sample of 31 exoplanetary systems detected by gravitational microlensing and investigate whether or not the estimated distances to these systems conform to the Galactic distribution of planets expected from models. We derive the expected distribution of distances and relative proper motions from a simulated microlensing survey, correcting for the dominant selection effects that affect the sensitivity of planet detection as a function of distance, and compare it to the observed distribution using Anderson–Darling (AD) hypothesis testing. Taking the relative abundance of planets in the bulge to that in the disk, f{sub bulge}, as a model parameter, we find that our model is consistent with the observed distribution only for f{sub bulge}<0.54 (for a p-value threshold of 0.01) implying that the bulge may be devoid of planets relative to the disk. Allowing for a dependence of planet abundance on metallicity and host mass, or an additional dependence of planet sensitivity on event timescale, does not restore consistency for f{sub bulge}=1. We examine the distance estimates of some events in detail, and conclude that some parallax-based estimates could be significantly in error. Only by combining the removal of one problematic event from our sample and the inclusion of strong dependences of planet abundance or detection sensitivity on host mass, metallicity, and event timescale are we able to find consistency with the hypothesis that the bulge and disk have equal planet abundance.

79 ASTRONOMY AND ASTROPHYSICS↗

Coding isotropic images

Rate distortion functions for two-dimensional homogeneous isotropic images are compared with the performance of 5 source encoders designed for such images. Both unweighted and frequency weighted mean square error distortion measures are considered. The coders considered are differential PCM (DPCM) using six previous samples in the prediction, herein called 6 pel (picutre element) DPCM; simple DPCM using single sample prediction; 6 pel DPCM followed by entropy coding; 8 x 8 discrete cosine transform coder, and 4 x 4 Hadamard transform coder. Other transform coders were studied and found to have about the same performance as the two transform coders above. With the mean square error distortion measure DPCM with entropy coding performed best. The relative performance of the coders changes slightly when the distortion measure is frequency weighted mean square error. The performance of all the coders was separated by only about 4 dB.

Oneal, J. B., Jr.↗

Linear discriminant analysis with misallocation in training samples

Linear discriminant analysis for a two-class case is studied in the presence of misallocation in training samples. A general appraoch to modeling of mislocation is formulated, and the mean vectors and covariance matrices of the mixture distributions are derived. The asymptotic distribution of the discriminant boundary is obtained and the asymptotic first two moments of the two types of error rate given. Certain numerical results for the error rates are presented by considering the random and two non-random misallocation models. It is shown that when the allocation procedure for training samples is objectively formulated, the effect of misallocation on the error rates of the Bayes linear discriminant rule can almost be eliminated. If, however, this is not possible, the use of Fisher rule may be preferred over the Bayes rule.

Chhikara, R.↗

Gamma-ray Burst Energetics an the Gamma-ray Burst Hubble Diagram: Promises and Limitations

We present a complete sample of 29 gamma-ray bursts (GRBs) for which it has been possible to determine temporal breaks (or limits) from their afterglow light curves. We interpret these breaks within the framework of the uniform conical jet model, incorporating realistic estimates of the ambient density and propagating error estimates on the measured quantities. In agreement with our previous analysis of a smaller sample, the derived jet opening angles of those 16 bursts with redshifts result in a narrow clustering of geometrically corrected gamma-ray energies about 1.33 x 10(exp 51) ergs; the burst-to-burst variance about this value is 0.35 dex, a factor of 2.2. Despite this rather small scatter, we demonstrate in a series of GRB Hubble diagrams that the current sample cannot place meaningful constraints upon the fundamental parameters of the universe. Indeed, for GRBs to ever be useful in cosmographic measurements, we argue the necessity of two directions. First, GRB Hubble diagrams should be based upon fundamental physical quantities such as energy, rather than empirically derived and physically ill-understood distance indicators (such as those based upon prompt burst time-profiles and spectra). Second, a more homogeneous set should be constructed by culling subclasses from the larger sample. These subclasses, although now first recognizable by deviant energies, ultimately must be identifiable by properties other than those directly related to energy. We identify a new subclass of GRBs (" f-GRBs ") that appear both underluminous by factors of at least 10 and exhibit a rapid fading (f(sub nu is proportional to t(sup -2) at early times (t < or = 0.5 day). About 10%-20% of observed long-duration bursts appear to be f-GRBs.

Bloom, J. S.↗

CPC data from TAMU TRACER campaign in the Houston, TX region from July to September 2022

During TRACER, the Texas A&M Rapid Onsite Atmospheric Measurements Van (ROAM-V) was deployed to capture airmasses behind (maritime) and ahead (continental) of the passage of the sea-breeze front through Houston. On select sampling days, ROAM-V sampled in the morning/mid-day on the coast and then transited to a second inland site for the afternoon/evening. The suite of instruments deployed on ROAM-V included a Condensation Particle Counter (CPC; GRIMM Model 5.403 CPC), Scanning Mobility Particle Sizer (SMPS; TSI 3750 detector, TSI 3082 classifier, TSI 3088 neutralizer, TSI 3081A Differential Mobility Analyzer), Cloud Condensation Nuclei counter (Droplet Measurement Technologies CCN Counter), micro pulse lidar (Droplet Measurement Technologies Micro Pulse LiDAR (miniMPL)), and a Davis Rotating Uniform size-cut Monitor (DRUM; DRUMAir 4-DRUM). Before sampling at each location, the latitude and longitude were recorded using the GPS on the phone application “My Altitude”. Onboard the ROAM-V, aerosol samples are drawn through a shared isokinetic inlet at a flow rate ranging from 3.5 to 7.0 LPM. A portion of this flow is directed through a cyclone impactor (Brechtel, Inc. Model SCC 0.732) and 0.3 LPM is directed to the CPC. To correct the data for particle losses, we used a two-step method. First, the measured SMPS size distributions were used to calculate total particle loss through the inlet for every SMPS scan at a single deployment location. Second, the CPC data was corrected for particle losses using the average of the total losses per scan. This calculation was done separately for each deployment location due to changes in the measured size distributions between locations. Particle loss from diffusion (based on Kesten, 1991 and Gormley, 1949), inertial impaction in 90-degree bends (based on Aerosol Measurement, 2011 and Crane, 1977), and cyclone impactor efficiency (based on Dirgo, 1985) were included in the loss calculation. When the SMPS was not sampling at a location (in the case of an instrument malfunction or operator error), the reported CPC data was corrected with an average of the total losses for the entire campaign at the specified deployment location (e.g., if we needed to correct Galveston data, then the average of all calculated losses at Galveston was taken). These flatline corrections were used for all data on 22/07/13, 22/07/20, 22/07/22, and the data from Galveston on 22/08/09. An additional flatline correction of 14% was applied to all CPC data based on laboratory calibration prior to and after the campaign. In laboratory calibrations, we identified that for the same sample of air the field CPC (GRIMM Model 5.403 CPC) undercounted the total concentration of particles by 14% compared to the ground truth laboratory CPC (TSI Model 3750). This data was collected for ARM Field Campaign AFC07055 and supported by DOE ASR grant DE-SC0021047. For any further questions, please feel free to contact the instrument PI, Sarah D. Brooks, sbrooks@tamu.edu. Kesten et. al. Calibration of a TSI Model 3025 Ultrafine Condensation Particle Counter. Aerosol Science and Technology, 15:2, 107-111, 1991. Gormley et. al. Diffusion from a Stream Flowing through a Cylindrical Tube. Proceedings of the Royal Irish Academy, Vol 52, 163-169, 1948. Aerosol Measurement: Principles, Techniques, and Applications, Third Edition. John Wiley & Sons, Inc, 2011. Crane et. al. Inertial Deposition of Particles in a Bent Pipe. Journal of Aerosol Science, Vol 8, 161-170, 1977. Dirgo et. al. Cyclone Collection Efficiency: Comparison of Experimental Results with Theoretical Predictions. Aerosol Science and Technology, 4:4, 401-415, 1985.

54 ENVIRONMENTAL SCIENCES↗

Northern dwarf and low surface brightness galaxies. II - The Green Bank neutral hydrogen survey

The paper reports neutral hydrogen observations of a large sample of dwarf and other low surface brightness galaxies. A detailed discussion and error analysis of the observations are presented, and spectra are displayed for 329 galaxies detected for the first time, or detected with substantially better signal-to-noise ratios than achieved previously. The positions on the sky of 667 galaxies meeting the present selection criteria north of delta = 38 deg are shown. The distribution of the redshifts of galaxies detected at Green Bank is illustrated. The Green Bank detections tapered off strongly below the median H I flux of 3.7 Jy km/s detected at Arecibo: only 12 percent of the Green Bank sample was detected with smaller fluxes.

Schneider, Stephen E.↗

Predicting critical heat flux with uncertainty quantification and domain generalization using conditional variational autoencoders and deep neural networks

Deep generative models (DGMs) can generate synthetic data samples that closely resemble the original dataset, addressing data scarcity. In this work, we developed a conditional variational autoencoder (CVAE) to augment critical heat flux (CHF) data used for the 2006 Groeneveld lookup table. To compare with traditional methods, a fine-tuned deep neural network (DNN) regression model was evaluated on the same dataset. Both models achieved small mean absolute relative errors, with the CVAE showing more favorable results. Uncertainty quantification (UQ) was performed using repeated CVAE sampling and DNN ensembling. The DNN ensemble improved performance over the baseline, while the CVAE maintained consistent results with less variability and higher confidence. Both models achieved small errors inside and outside the training domain, with slightly larger errors outside. Altogether, the CVAE performed better than the DNN in predicting CHF and exhibited better uncertainty behavior.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The velocity distribution function of galaxy clusters as a cosmological probe

We present a new approach for quantifying the abundance of galaxy clusters and constraining cosmological parameters using dynamical measurements. In the standard method, galaxy line-of-sight velocities, v, or velocity dispersions are used to infer cluster masses, M, to quantify the halo mass function (HMF), dn(M)/dlog(M), which is strongly affected by mass measurement errors. In our new method, the probability distributions of velocities for each cluster in the sample are summed to create a new statistic called the velocity distribution function (VDF), dn(v)/dv. The VDF can be measured more directly and precisely than the HMF and can be robustly predicted with cosmological simulations that capture the dynamics of subhalos or galaxies. We apply these two methods to realistic (ideal) mock cluster catalogs with (without) interlopers and forecast the bias and constraints on the matter density parameter Ω{sub m} and the amplitude of matter fluctuations σ{sub 8} in flat ΛCDM cosmologies. For an example observation of 200 massive clusters, the VDF with (without) interloping galaxies constrains the parameter combination σ{sub 8} Ω{sub m}{sup 0.29(0.29)}=0.589±0.014 (0.584±0.011) and shows only minor bias. However, the HMF with interlopers is biased to low Ω{sub m} and high σ{sub 8} and the fiducial model lies well outside of the forecast constraints, prior to accounting for Eddington bias. When the VDF is combined with constraints from the cosmic microwave background, the degeneracy between cosmological parameters can be significantly reduced. Upcoming spectroscopic surveys that probe larger volumes and fainter magnitudes will provide clusters for applying the VDF as a cosmological probe.

79 ASTRONOMY AND ASTROPHYSICS↗

Algorithm and program for precise determination of unit-cell parameters of single crystal taking into account the sample eccentricity

A technique has been developed to refine the unit-cell parameters of single crystals with minimization of the influence of instrumental errors on the result. The corresponding computational procedure HuberUB is added to the software package of Huber-5042 diffractometer with a point detector and closedcycle helium cryostat Displex DE-202. The parameters of unit cell, its orientation, the goniometer zero angles, the sample eccentricity, the distances in the goniometer, and the radiation wavelength were refined by the nonlinear least-squares method, which allows imposition of constraints on the unit-cell parameters, depending on the crystal symmetry. The technique is approved on a LuB{sub 12} single crystal. The unit-cell parameters are determined in a temperature range of 20–295 K, with an absolute error not larger than 0.0004 Å (the relative error is of ~5 × 10{sup –5}). The estimates of the unit-cell parameters obtained by the proposed method are evidenced to be unbiased. Some specific features of the behavior of parameters in the ranges of 120–140 and 20–50 K are revealed, which correlate with the anomalies of the physical properties of the crystal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On the Ability of Space- Based Passive and Active Remote Sensing Observations of CO2 to Detect Flux Perturbations to the Carbon Cycle

Space-borne observations of CO2 are vital to gaining understanding of the carbon cycle in regions of the world that are difficult to measure directly, such as the tropical terrestrial biosphere, the high northern and southern latitudes, and in developing nations such as China. Measurements from passive instruments such as GOSAT (Greenhouse Gases Observing Satellite) and OCO-2 (Orbiting Carbon Observatory 2), however, are constrained by solar zenith angle limitations as well as sensitivity to the presence of clouds and aerosols. Active measurements such as those in development for the Active Sensing of CO2 Emissions over Nights, Days and Seasons (ASCENDS) mission show strong potential for making measurements in the high-latitude winter and in cloudy regions. In this work we examine the enhanced flux constraint provided by the improved coverage from an active measurement such as ASCENDS. The simulation studies presented here show that with sufficient precision, ASCENDS will detect permafrost thaw and fossil fuel emissions shifts at annual and seasonal time scales, even in the presence of transport errors, representativeness errors, and biogenic flux errors. While OCO-2 can detect some of these perturbations at the annual scale, the seasonal sampling provided by ASCENDS provides the stronger constraint. Plain Language Summary: Active and passive remote sensors show the potential to provide unprecedented information on the carbon cycle. With the all-season sampling, active remote sensors are more capable of constraining high-latitude emissions. The reduced sensitivity to cloud and aerosol also makes active sensors more capable of providing information in cloudy and polluted scenes with sufficient accuracy. These experiments account for errors that are fundamental to the top-down approach for constraining emissions, and even including these sources of error, we show that satellite remote sensors are critical for understanding the carbon cycle.

tropical terrestrial biosphere↗

Taylor approximation variance reduction for approximation errors in PDE-constrained Bayesian inverse problems

In numerous applications, surrogate models are used as a replacement for accurate parameter-to-observable mappings when solving large-scale inverse problems governed by partial differential equations (PDEs). The surrogate model may be a computationally cheaper alternative to the accurate parameter-to-observable mappings and/or may ignore additional unknowns or sources of uncertainty. The Bayesian approximation error (BAE) approach provides a means to account for the induced uncertainties and approximation errors, i.e. the errors between the accurate parameter-to-observable mapping and the surrogate. The statistics of these errors are, however, in general unknown a priori, and are thus calculated using Monte Carlo sampling. Although the sampling is typically carried out offline, i.e. before considering the data, the process can still represent a computational bottleneck. In this work, we develop a scalable computational approach for reducing the costs associated with the sampling stage of the BAE approach. Specifically, we consider the Taylor expansion of the accurate and surrogate forward models with respect to the uncertain parameter fields either as a control variate for variance reduction or as a means to directly and efficiently approximate the mean and covariance of the approximation errors. We propose efficient methods for evaluating the expressions for the mean and covariance of the Taylor approximations based on linear(-ized) PDE solves. Furthermore, the proposed approach is independent of the dimension of the uncertain parameter, depending instead on the intrinsic dimension of the data, ensuring scalability to high-dimensional problems. The potential benefits of the proposed approach are demonstrated for two high-dimensional inverse problems governed by PDE examples, namely for the estimation of a distributed Robin boundary coefficient in a linear diffusion problem, and for a coefficient estimation problem governed by a nonlinear diffusion problem.

Bayesian approximation error↗

Non-local Thermodynamic Equilibrium Stellar Spectroscopy with 1D and 〈3D〉 Models. I. Methods and Application to Magnesium Abundances in Standard Stars

We determine Mg abundances in six Gaia benchmark stars using theoretical one-dimensional (1D) hydrostatic model atmospheres, as well as temporally and spatially averaged three-dimensional (〈3D〉) model atmospheres. The stars cover a range of T{sub eff} from 4700 to 6500 K, logg from 1.6 to 4.4 dex, and [Fe/H] from −3.0 dex to solar. Spectrum synthesis calculations are performed in local thermodynamic equilibrium (LTE) and in non-LTE (NLTE) using the oscillator strengths recently published by Pehlivan Rhodin et al. We find that: (a) Mg abundances determined from the infrared spectra are as accurate as the optical diagnostics, (b) the NLTE effects on Mg i line strengths and abundances in this sample of stars are minor (although for a few Mg i lines the NLTE effects on abundance exceed 0.6 dex in 〈3D〉 and 0.1 dex in 1D, (c) the solar Mg abundance is 7.56±0.05 dex (total error), in excellent agreement with the Mg abundance measured in CI chondritic meteorites, (d) the 1D NLTE and 〈3D〉 NLTE approaches can be used with confidence to analyze optical Mg i lines in spectra of dwarfs and sub-giants, but for red giants the Mg i 5711 Å line should be preferred, (e) low-excitation Mg i lines are sensitive to the atmospheric structure; for these lines, LTE calculations with 〈3D〉 models lead to significant systematic abundance errors. The methods developed in this work will be used to study Mg abundances of a large sample of stars in the next paper in the series.

79 ASTRONOMY AND ASTROPHYSICS↗