Search NASA⌕ Search

SEARCH · Search NASA

Results for “Function Estimation”

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,783 records · Page 99

Uncertainty analysis of signal deconvolution using a measured instrument response function

A common analysis procedure minimizes the ln-likelihood that a set of experimental observables matches a parameterized model of the observation. The model includes a description of the underlying physical process as well as the instrument response function (IRF). Here, we investigate the National Ignition Facility (NIF) neutron time-of-flight (nTOF) spectrometers, the IRF is constructed from measurements and models. IRF measurements have a finite precision that can make significant contributions to the uncertainty estimate of the physical model’s parameters. Finally, we apply a Bayesian analysis to properly account for IRF uncertainties in calculating the ln-likelihood function used to find the optimum physical parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Optimization-based scatter estimation using primary modulation for computed tomography

Purpose: Scatter reduces the image quality in computed tomography (CT), but scatter correction remains a challenge. A previously proposed primary modulation method simultaneously obtains the primary and scatter in a single scan. However, separating the scatter and primary in primary modulation is challenging because it is an underdetermined problem. In this study, an optimization-based scatter estimation (OSE) algorithm is proposed to estimate and correct scatter. Methods: In the concept of primary modulation, the primary is modulated, but the scatter remains smooth by inserting a modulator between the x-ray source and the object. In the proposed algorithm, an objective function is designed for separating the scatter and primary. Prior knowledge is incorporated in the optimization-based framework to improve the accuracy of the estimation: (1) the primary is always positive; (2) the primary is locally smooth and the scatter is smooth; (3) the location of penumbra can be determined; and (4) the scatter-contaminated data provide knowledge about which part is smooth. Results: The simulation study shows that the edge-preserving weighting in OSE improves the estimation accuracy near the object boundary. Simulation study also demonstrates that OSE outperforms the two existing primary modulation algorithms for most regions of interest in terms of the CT number accuracy and noise. The proposed method was tested on a clinical cone beam CT, demonstrating that OSE corrects the scatter even when the modulator is not accurately registered. Conclusions: The proposed OSE algorithm improves the robustness and accuracy in scatter estimation and correction. This method is promising for scatter correction of various kinds of x-ray imaging modalities, such as x-ray radiography, cone beam CT, and the fourth-generation CT.

60 APPLIED LIFE SCIENCES↗

Metagenomic analysis of the microbiota in the highly compartmented hindguts of six wood- or soil-feeding higher termites

Termites are important contributors to carbon and nitrogen cycling in tropical ecosystems. Higher termites digest lignocellulose in various stages of humification with the help of an entirely prokaryotic microbiota housed in their compartmented intestinal tract. Previous studies revealed fundamental differences in community structure between compartments, but the functional roles of individual lineages in symbiotic digestion are mostly unknown. Furthermore, we conducted a highly resolved analysis of the gut microbiota in six species of higher termites that feed on plant material at different levels of humification. Combining amplicon sequencing and metagenomics, we assessed similarities in community structure and functional potential between the major hindgut compartments (P1, P3, and P4). Cluster analysis of the relative abundances of orthologous gene clusters (COGs) revealed high similarities among woodand litter-feeding termites and strong differences to humivorous species. However, abundance estimates of bacterial phyla based on 16S rRNA genes greatly differed from those based on protein-coding genes. In conclusion, the community structure and functional potential of the microbiota in individual gut compartments are clearly driven by the digestive strategy of the host. The metagenomics libraries obtained in this study provide the basis for future studies that elucidate the fundamental differences in the symbiont-mediated breakdown of lignocellulose and humus by termites of different feeding groups. Furthermore, the high proportion of uncultured bacterial lineages in all samples calls for a reference-independent approach for the correct taxonomic assignment of protein-coding genes.

59 BASIC BIOLOGICAL SCIENCES↗

An Extended Group Additivity Method for Polycyclic Thermochemistry Estimation

Automatic kinetic mechanism generation, virtual high-throughput screening, and automatic transition state search are currently trending applications requiring exploration of a large molecule space. Large-scale search requires fast and accurate estimation of molecules’ properties of interest, such as thermochemistry. Existing approaches are not satisfactory for large polycyclic molecules: considering the number of molecules being screened, quantum chemistry (even cheap density functional theory methods) can be computationally expensive, and group additivity, though fast, is not sufficiently accurate. Here, this paper provides a fast and moderately accurate alternative by proposing a polycyclic thermochemistry estimation method that extends the group additivity method with two additional algorithms: similarity match and bicyclic decomposition. It significantly reduces H f (298 K) estimation error from over 60 kcal/mol (group additivity method) to around 5 kcal/mol, C p (298 K) error from 9 to 1 cal/mol/K,and S(298 K) error from 70 to 7 cal/mol/K. This method also works well for heteroatomic polycyclics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Calculation and interpretation of classical turning surfaces in solids

Abstract Classical turning surfaces of Kohn–Sham potentials separate classically allowed regions (CARs) from classically forbidden regions (CFRs). They are useful for understanding many chemical properties of molecules but need not exist in solids, where the density never decays to zero. At equilibrium geometries, we find that CFRs are absent in perfect metals, rare in covalent semiconductors at equilibrium, but common in ionic and molecular crystals. In all materials, CFRs appear or grow as the internuclear distances are uniformly expanded. They can also appear at a monovacancy in a metal. Calculations with several approximate density functionals and codes confirm these behaviors. A classical picture of conduction suggests that CARs should be connected in metals, and disconnected in wide-gap insulators, and is confirmed in the limits of extreme compression and expansion. Surprisingly, many semiconductors have no CFR at equilibrium, a key finding for density functional construction. Nonetheless, a strong correlation with insulating behavior can still be inferred. Moreover, equilibrium bond lengths for all cases can be estimated from the bond type and the sum of the classical turning radii of the free atoms or ions.

Chemistry↗

Physically based correction of systematic errors of Rotating Shadowband Irradiometers

Accurate measurements of direct normal, diffuse horizontal and global horizontal irradiance (DNI, DHI and GHI) are needed for meteorological studies and are essential for the solar resource assessment at potential solar power plant sites. Often, these potential sites are remote and hence require robust sensors that require minimal maintenance that are not affected strongly by soiling. Therefore, Rotating Shadowband Irradiometers (RSI) are widely used for resource assessment. To achieve the required accuracy, corrections for the raw values of RSIs depending on systematic temperature, incidence angle and spectral errors must be used, and a thorough calibration of the sensor head must be applied. The existing correction functions are derived from comparisons of RSIs to thermopile radiometers at selected sites and therefore empirical. Their accuracy is considered to be site dependent. In this work a new correction and calibration method is presented that removes the systematic errors using a physical approach. It is based on information of the sensor properties as well as measurements of its directional response, and incorporates the atmospheric conditions at the measurement site. In this case, no empiric relations obtained from a specific site are required. The method requires estimates of the current DHI and GHI spectra during each measurement of the RSI. Based on these spectra, a spectral correction, which includes a spectrum dependent temperature correction, can be made without employing empirical relationships. The new physical calibration and correction method is tested at three sites and reaches similar results compared to the empirical functions. This is already achieved with rudimentary estimations of the GHI and DHI spectra and we expect that these estimations can be improved in the future. The results indicate that the physical approach reduces the problematic location dependence of the current calibration and correction methods. The physical correction and calibration method show promise for a further improvement of the RSI accuracy.

47 OTHER INSTRUMENTATION↗

Atomistic Investigation of Plastic Deformation and Dislocation Motion in Uranium Mononitride

Uranium mononitride (UN) is a promising advanced nuclear fuel due to its high thermal conductivity and high fissile density. However, many aspects of its mechanical behavior, particularly at reactor-relevant conditions, remain unclear. In this study, molecular dynamics (MD) simulations were employed to investigate the deformation behavior and dislocation motion in UN. We found that the Kocevski potential predicts the principal slip system as $\frac{1}{2}$ $\langle110\rangle${110}, aligning with experimental data. On the other hand, the Tseplyaev potential predicts slip to primarily occur on $\frac{1}{2}$ $\langle110\rangle${111}. MD simulations of stress–strain behavior were used to estimate the nanoindentation hardness, revealing that the Kocevski potential accurately predicts hardness even though it fails to model dynamic plasticity. Complete dislocation mobility functions have been fitted for the edge and screw dislocations in both the thermally activated and phonon-drag regimes. The 300 K linear mobility of the edge dislocation using the Tseplyaev potential was found to be 817 Pa -1 ·s -1 , whereas that of the screw dislocation using the Kocevski potential was found to be 4546 Pa -1 ·s -1 . At intermediate stresses, we observed that the subsonic steady-state motion of the edge dislocation in UN is intermittently interrupted by velocity jumps, reaching the average sound velocity. Finally, the threshold Schmid stress is calculated as 179–197 MPa, which gives an upper-limit estimate of the uniaxial yield stress of polycrystalline UN of 548–603 MPa. These findings, including the fitted dislocation mobility function, provide essential input for future plasticity and dislocation dynamics models of nuclear fuels.

36 MATERIALS SCIENCE↗

Predictions for Sivers single spin asymmetries in one- and two-hadron electroproduction at CLAS12 and EIC

The study of the Sivers effect, describing correlations between the transverse polarization of the nucleon and its constituent (unpolarized) parton's transverse momentum, has been the topic of a great deal of experimental, phenomenological and theoretical effort in recent years. Semi-Inclusive Deep Inelastic Scattering measurements of the corresponding single spin asymmetries (SSA) at the upcoming CLAS12 experiment at JLab and the proposed Electron-Ion Collider will help to pinpoint the flavor structure and the momentum dependence of the Sivers parton distribution function describing this effect. Here we describe a modified version of the $\tt{PYTHIA}$ Monte Carlo event generator that includes the Sivers effect. Then we use it to estimate the size of these SSAs, in the kinematics of these experiments, for both one and two hadron final states of pions and kaons. For this purpose we utilize the existing Sivers parton distribution function (PDF) parametrization extracted from HERMES and COMPASS experiments. Furthermore, we also show that the the leading order approximation commonly used in such extractions provides significantly underestimated values of Sivers PDFs, as the omitted parton showers and non-DIS processes play an important role in these SSAs at lower light-cone momentum fraction, for example in the COMPASS kinematics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Neutron light output function and resolution investigation of the deuterated organic liquid scintillator EJ-315

Here, we investigate the light response to fast neutrons and estimate the pulse height resolution of a deuterated liquid organic scintillator, EJ-315, considering the detector's non-linear light response to gamma-rays. Initially, collision data and a neutron beam trigger are recorded in coincidence mode, and incident neutron energy is calculated with a time-of-flight technique. Fast neutrons are further discriminated from gamma-ray background based on the scintillation material decay patterns using a pulse shape discrimination algorithm. We also derived a light response matrix composed of multiple neutron energy and their corresponding light outputs. The pulse height resolution property of the EJ-315 is characterized utilizing the derivatives of the pulse height distributions with corrections of the mea-surements setup uncertainties. In addition, the EJ-315's pulse height resolution is also characterized by comparing the smoothed derivatives of quasi-monoenergetic neutron pulse height distributions, given by the Peierls-formula-based analytic model, to match the measurement data. Results show rather consistent 10e13% pulse resolution for mono-energetic neutrons with kinetic energy above 2 MeV. The resolution decreases slightly with an increase in neutron energy indicating the improved resolution performance of EJ-315 in the higher energy events.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Effects of Improved 17 O Correction on Interlaboratory Agreement in Clumped Isotope Calibrations, Estimates of Mineral-Specific Offsets, and Temperature Dependence of Acid Digestion Fractionation

The clumped isotopic composition of carbonate-CO 2 (denoted Δ 47 ) is a function of carbonate formation temperature and in natural samples can act as a recorder of paleoclimate, burial, or diagenetic conditions. The absolute abundance of heavy isotopes in the universal standards VPDB and VSMOW (defined by four parameters: R 13 VPDB , R 17 VSMOW , R 18 VSMOW , and λ ) impact calculated Δ 47 values. In this work, we investigate whether use of updated and more accurate values for these parameters can remove observed interlaboratory differences in the measured T-Δ 47 relationship. Using the updated parameters, we reprocess 14 published calibration data sets measured in 11 different laboratories, representing many mineralogies, bulk compositions, sample types, reaction temperatures, and sample preparation and analysis methods. Exploiting this large composite data set ( n = 1,253 sample replicates), we investigate the possibility for a “universal” clumped isotope calibration. We find that applying updated parameters improves the T-Δ 47 relationship (reduces residuals) within most labs and improves overall agreement but does not eliminate all interlaboratory differences. We reaffirm earlier findings that different mineralogies do not require different calibration equations and that cleaning procedures, method of pressure baseline correction, and mass spectrometer type do not affect interlaboratory agreement. We also present new estimates of the temperature dependence of the acid digestion fractionation for Δ 47 (Δ* 25-X ), based on combining reprocessed data from four studies, and new theoretical equilibrium values to be used in calculation of the empirical transfer function. Overall, we have ruled out a number of possible causes of interlaboratory disagreement in the T-Δ 47 relationship, but many more remain to be investigated.

58 GEOSCIENCES↗

Applications of the worldline Monte Carlo formalism in quantum mechanics

Highlights: Universal approach to the numerical computation of the Euclidean path integral. Inspired by recent work in relativistic quantum field theory. Here adapted to non-relativistic quantum mechanics. Worked out for the computation of propagators and ground-state energies. Special smoothing procedure for singular potentials. -- Abstract: In recent years efficient algorithms have been developed for the numerical computation of relativistic single-particle path integrals in quantum field theory. Here, we adapt this “worldline Monte Carlo” approach to the standard problem of the numerical approximation of the non-relativistic path integral, resulting in a formalism whose characteristic feature is the fast, non-recursive generation of an ensemble of trajectories that is independent of the potential, and is thus universally applicable. The numerical implementation discretises the trajectories with respect to their time parametrisation but maintains a continuous spatial domain. In the case of singular potentials, the discretised action gets adapted to the singularity through a “smoothing” procedure. We show for a variety of examples (the harmonic oscillator in various dimensions, the modified Pöschl–Teller potential, delta-function potentials, the Coulomb and Yukawa potentials) that the method allows one to obtain fast and reliable estimates for the Euclidean propagator and use them in a certain time window suitable for extracting the ground state energy. Usually Monte Carlo calculations in quantum mechanics are performed using algorithms that involve a more specific adaptation to the potential, and we briefly compare our method with the most widely used such algorithms, which are of Metropolis type. As an aside, we apply it for studying the classical limit where nearly classical trajectories are expected to dominate in the path integral. We expect the advances made here to be useful also in the relativistic case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

EVOLUTION OF STELLAR-TO-HALO MASS RATIO AT z = 0–7 IDENTIFIED BY CLUSTERING ANALYSIS WITH THE HUBBLE LEGACY IMAGING AND EARLY SUBARU/HYPER SUPRIME-CAM SURVEY DATA

We present clustering analysis results from 10,381 Lyman break galaxies (LBGs) at z ∼ 4–7, identified in the Hubble legacy deep imaging and new complimentary large-area Subaru/Hyper Suprime-Cam data. We measure the angular correlation functions of these LBGs at z ∼ 4, 5, 6, and 7 and fit these measurements using halo occupation distribution (HOD) models that provide an estimate of halo masses, M{sub h}∼(1−20)×10{sup 11} M{sub ⊙}. Our M{sub h} estimates agree with those obtained by previous clustering studies in a UV-magnitude versus M{sub h} plane and allow us to calculate stellar-to-halo mass ratios (SHMRs) of LBGs. By comparison with the z∼0 SHMR, we identify evolution of the SHMR from z∼0 to z∼4 and from z∼4 to z∼7 at the >98% confidence levels. The SHMR decreases by a factor of ∼2 from z∼0 to 4 and increases by a factor of ∼4 from z∼4 to 7 at the dark matter halo mass of M{sub h}∼10{sup 11} M{sub ⊙}. We compare our SHMRs with results of a hydrodynamic simulation and a semianalytic model and find that these theoretical studies do not predict the SHMR increase from z∼4 to 7. We obtain the baryon conversion efficiency (BCE) of LBGs at z∼4 and find that the BCE increases with increasing dark matter halo mass. Finally, we compare our clustering+HOD estimates with results from abundance matching techniques and conclude that the M{sub h} estimates of the clustering+HOD analyses agree with those of the simple abundance matching within a factor of 3, and that the agreement improves when using more sophisticated abundance matching techniques that include subhalos, incompleteness, and/or evolution in the star formation and stellar mass functions.

79 ASTRONOMY AND ASTROPHYSICS↗

On the uncertainty of estimating photovoltaic soiling using nearby soiling data

The accumulation of dust on the surface of photovoltaic modules can reduce their performance and affect the cost competitiveness of this technology. This phenomenon is known as soiling and can be mitigated through appropriate corrective and/or preventive actions. In order to maximize its effectiveness, it is important to plan the soiling mitigation strategy even before the PV system is operational. This is typically done through a nearest neighbor approach, by estimating soiling using data from the nearest operational photovoltaic system. This work focuses on understanding the uncertainty related to this practice. For this purpose, the semi-variance function is used to study the dissimilarity between the soiling losses of two locations in California depending on their distance. The results show that, when the soiling loss at a nearby system is used to estimate soiling of a site, the uncertainty can be approximated to increase linearly at a rate of 0.08-0.10%/km up to 60 or 80 km. After this distance, the use of a nearest neighbor approach is no longer justified, as it produces an uncertainty as big as the average soiling loss of the sites in the dataset used in this study. In some conditions, uncertainties > 0% are found also for sites located within 25 km, meaning that even close-by systems might soil differently.

14 SOLAR ENERGY↗

z ∼ 1 Lyα emitters. I. The luminosity function , , ,

We construct a flux-limited sample of 135 candidate z ∼ 1 Lyα emitters (LAEs) from Galaxy Evolution Explorer (GALEX) grism data using a new data cube search method. These LAEs have luminosities comparable to those at high redshifts and lie within a 7 Gyr gap present in existing LAE samples. We use archival and newly obtained optical spectra to verify the UV redshifts of these LAEs. We use the combination of the GALEX UV spectra, optical spectra, and X-ray imaging data to estimate the active galactic nucleus (AGN) fraction and its dependence on Lyα luminosity. We remove the AGNs and compute the luminosity function (LF) from 60 z ∼ 1 LAE galaxies. We find that the best-fit LF implies a luminosity density increase by a factor of ∼1.5 from z ∼ 0.3 to z ∼ 1 and ∼20 from z ∼ 1 to z ∼ 2. We find a z ∼ 1 volumetric Lyα escape fraction of 0.7% ± 0.4%.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurements of the electron energy distribution function in partially magnetized low temperature plasmas

While Langmuir probes (LPs) are relatively simple and inexpensive plasma diagnostics for the electron density, temperature, and the electron energy distribution function (EEDF), the interpretation of the measured current–voltage (I–V) characteristic is complicated considerably by the presence of a magnetic field. In regimes where the electron gyroradius is comparable to the probe radius, the electron flux to the probe surface is retarded by reduced mobility across field lines and can no longer be described by a thermal model. Predicting the current collected by the probe in these regimes requires accurate estimates of the plasma diffusion coefficients, which are usually difficult to obtain. In this work, we measure electron energy distribution functions in E×B plasmas with magnetized electrons and non-magnetized ions in argon and krypton gases, using both a LP and laser Thomson scattering (LTS) at various magnetic fields. Using the LTS measurements to provide a robust benchmark for comparison, we compare existing theories describing the flux to the probe under magnetized conditions. We find that even when the electron gyroradius associated with the effective electron temperature is small compared to the probe radius, the EEDF computed using classical probe theory is still robust at energies higher than the energy at which the gyroradius becomes larger than the probe size. For plasmas that are not strongly non-Maxwellian, we formulate a method to extract robust density and temperature measurements using physics-informed fitting techniques to analyze EEDFs computed using classical theory.

Devin, E. G. [Princeton Plasma Physics Laboratory ↗

Over 150 years of long-term fertilization alters spatial scaling of microbial biodiversity

Spatial scaling is a critical issue in ecology, but how anthropogenic activities like fertilization affect spatial scaling is poorly understood, especially for microbial communities. Here, we determined the effects of long-term fertilization on the spatial scaling of microbial functional diversity and its relationships to plant diversity in the 150-year-old Park Grass Experiment, the oldest continuous grassland experiment in the world. Nested samples were taken from plots with contrasting inorganic fertilization regimes, and community DNAs were analyzed using the GeoChip-based functional gene array. The slopes of microbial gene-area relationships (GARs) and plant species-area relationships (SARs) were estimated in a plot receiving nitrogen (N), phosphorus (P), and potassium (K) and a control plot without fertilization. Our results indicated that long-term inorganic fertilization significantly increased both microbial GARs and plant SARs. Microbial spatial turnover rates (i.e., z values) were less than 0.1 and were significantly higher in the fertilized plot (0.0583) than in the control plot (0.0449) (P < 0.0001). The z values also varied significantly with different functional genes involved in carbon (C), N, P, and sulfur (S) cycling and with various phylogenetic groups (archaea, bacteria, and fungi). Similarly, the plant SARs increased significantly (P < 0.0001), from 0.225 in the control plot to 0.419 in the fertilized plot. Soil fertilization, plant diversity, and spatial distance had roughly equal contributions in shaping the microbial functional community structure, while soil geochemical variables contributed less. Results indicated that long-term agricultural practice could alter the spatial scaling of microbial biodiversity. Determining the spatial scaling of microbial biodiversity and its response to human activities is important but challenging in microbial ecology. Most studies to date are based on different sites that may not be truly comparable or on short-term perturbations, and hence, the results observed could represent transient responses. This study examined the spatial patterns of microbial communities in response to different fertilization regimes at the Rothamsted Research Experimental Station, which has become an invaluable resource for ecologists, environmentalists, and soil scientists. The current study is the first showing that long-term fertilization has dramatic impacts on the spatial scaling of microbial communities. In addition, by identifying the spatial patterns in response to long-term fertilization and their underlying mechanisms, the study makes fundamental contributions to predictive understanding of microbial biogeography.

54 ENVIRONMENTAL SCIENCES↗

CaCO 3 Crystallite Evolution during CaO Carbonation: Critical Crystallite Size and Rate Constant Measurement by In-Situ Synchrotron Radiation X-ray Powder Diffraction

In this work, the evolution of the CaCO 3 crystalline phase during the CaO–CO 2 reaction was investigated by means of in situ synchrotron radiation X-ray powder diffraction performed at the Advanced Photon Source (APS) facilities of the Argonne National Laboratory. CO 2 absorption experiments were carried out in a high temperature reaction capillary with a controlled atmosphere of pure carbon dioxide (CO 2 partial pressure of 1 bar) and in the temperature range between 450 and 750 °C, using CaO-based sorbents obtained by calcination of commercial calcium carbonate. The Rietveld refinement method was applied to estimate the average size of the CaCO 3 crystallites formed during a carbonation time of 20 min, as a function of the carbonation temperature and of the initial calcination conditions. Local maxima were observed in the CaCO 3 crystallite size versus time curves and were identified as the critical CaCO 3 crystallite sizes, marking the transition between the first fast carbonation stage and the second reaction stage controlled by product-layer diffusion. A relationship between this parameter and the reaction temperature, as well as with the initial (at the beginning of carbonation) CaO crystallite size, were found. The CaCO 3 critical crystallite sizes were used to estimate the active surface areas of the CaO sorbent particles where CaCO 3 crystals form and grow. Here, the computed active surface areas were utilized to calculate the kinetic parameters of the surface carbonation reaction: a reaction rate constant of 4.41 × 10 –4 mol/m 2 s, with zero-activation energy, was obtained.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗