Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sampling variability”

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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The Median Isn’t the Message: soil nutrient hot spots have a disproportionate influence on biogeochemical structure across years, seasons, and depths

Abstract Soil nutrient distribution is heterogeneous in space and time, potentially altering nutrient acquisition by trees and microorganisms. Ecologists have distinguished “hot spots” (HSs) as areas with enhanced and sustained rates of nutrient fluxes relative to the surrounding soil matrix. We evaluated the spatial and temporal patterns in nutrient flux HSs in two mixed-conifer forest soils by repeatedly sampling the soil solution at the same spatial locations (horizontally and vertically) over multiple seasons and years using ion exchange resins incubated in situ. The climate of these forests is Mediterranean, with intense fall rains occurring following summers with little precipitation, and highly variable winter snowfall. Hot spots formed most often for NO 3 − and Na + . Although nutrient HSs often occurred in the same spatial location multiple times, HSs persisted more often for PO 4 3− NH 4 + , and NO 3 − , and were more transient for Ca 2+ , Mg 2+ , and Na + . Sampling year (annual precipitation ranged from 558 to 1223 mm) impacted the occurrence of HSs for most nutrients, but season was only significant for PO 4 3− , NH 4 + , NO 3 − , and Na + , with HSs forming more often after fall rains than after spring snowmelt. The frequency of HSs significantly decreased with soil depth for all nutrients, forming most commonly immediately below the surficial organic horizon. Although HSs accounted for less than 17% of the sampling volume, they were responsible for 56–88% of PO 4 3− , NH 4 + , and NO 3 − resin fluxes. Our results suggest that macronutrient HSs have a disproportional contribution to soil biogeochemical structure, with implications for vegetation nutrient acquisition strategies and biogeochemical models. Graphical abstract

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

The SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Quasar target selection

As part of the Sloan Digital Sky Survey (SDSS) IV the extended Baryon Oscillation Spectroscopic Survey (eBOSS) will improve measurements of the cosmological distance scale by applying the Baryon Acoustic Oscillation (BAO) method to quasar samples. eBOSS will adopt two approaches to target quasars over 7500 deg{sup 2}. First, a “CORE” quasar sample will combine the optical selection in ugriz using a likelihood-based routine called XDQSOz, with a mid-IR-optical color cut. eBOSS CORE selection (to g < 22 or r < 22) should return ∼70 deg{sup −2} quasars at redshifts 0.9 < z < 2.2 and ∼7 deg{sup −2}z > 2.1 quasars. Second, a selection based on variability in multi-epoch imaging from the Palomar Transient Factory should recover an additional ∼3–4 deg{sup −2}z > 2.1 quasars to g < 22.5. A linear model of how imaging systematics affect target density recovers the angular distribution of eBOSS CORE quasars over 96.7% (76.7%) of the SDSS north (south) Galactic Cap area. The eBOSS CORE quasar sample should thus be sufficiently dense and homogeneous over 0.9 < z < 2.2 to yield the first few-percent-level BAO constraint near z-bar ∼1.5.eBOSS quasars at z > 2.1 will be used to improve BAO measurements in the Lyα Forest. Beyond its key cosmological goals, eBOSS should be the next-generation quasar survey, comprising >500,000 new quasars and >500,000 uniformly selected spectroscopically confirmed 0.9 < z < 2.2 quasars. At the conclusion of eBOSS, the SDSS will have provided unique spectra for more than 800,000 quasars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Combined Machine Learning and High-Energy X-ray Diffraction Approach to Understanding Liquid and Amorphous Metal Oxides

Determining the structure-property relations of liquid and amorphous metal oxides is challenging, due to their variable short-range order and polyhedral connectivity. To predict chemically realistic structures, we have developed a Machine Learned, Gaussian Approximation Potential (GAP) for HfO 2 , with a focus on enhanced sampling of the training database and accurate density functional theory calculations. By using training datasets for the GAP model at the level of Density Functional Theory-Strongly Constrained and Appropriately Normed (DFT-SCAN) level of theory, our results show that the topology of both the low viscosity liquid and the amorphous form are dominated by edge-shared chains and small corner-shared rings of polyhedra. This topology is shown to be consistent with the structure of other liquid and amorphous transition metal oxides of variable ion size, such as TiO 2 and ZrO 2 . Current limitations of the ML-GAP modeling method for obtaining glass structures and future perspectives are also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Modeling and Simulation of High-Frequency Solar Irradiance

As the use of solar power as a source of electricity is increasing, so is the interest in modeling radiation at high temporal resolutions. High-dimensional remote sensing data products depend on cloud cover variability, atmosphere aerosol levels, and other atmospheric parameters. Because of weather fronts and aerosols, it is difficult to quantify solar power variability based on distributed solar networks. The global horizontal irradiance (GHI) component in the National Solar Radiation Database (NSRDB) is available at a 30-min time resolution. The algorithm proposed in this paper produces 1-min-resolution GHI samples at locations where NSRDB data are available. The synthetic irradiance datasets are produced using log-additive non-Gaussian mixture models and resampling techniques. The model is trained over historical data, and predicted values are compared with in situ data. This approach allows for estimating the solar irradiance at subhourly temporal resolutions, while featuring variability for locations where measurements are otherwise not available.

14 SOLAR ENERGY↗
DOE OSTI2017

Installing Deep Groundwater Monitoring Wells in a Volcanic Dominated Rift Basin: Los Alamos National Laboratory, New Mexico - 17594

Los Alamos National Laboratory (LANL, or the Laboratory) is located in northern New Mexico near the west margin of the Espanola Basin. LANL is underlain by a complex basin-fill sequence of alluvial fan deposits, volcanic tuff, basaltic and dacitic lavas, and riverine deposits. On-going environmental remediation and characterization activities have evolved to improve the performance of groundwater monitoring wells, ensuring they yield representative groundwater and core samples. Groundwater occurs as canyon-fill alluvial water, perched intermediate-depth water associated with the basalts and other local-scale perching horizons, and as a regional-scale aquifer that is the main water supply for the area. Anthropogenic contaminants are found in all three groundwater zones. The regional aquifer below LANL is nominally 1000-1400 feet below ground surface (bgs). Monitoring wells at the Laboratory, are therefore, among the deepest routinely installed in the US. Logistical challenges at LANL include steep terrain, widespread cultural sites, and threatened and endangered species habitat restrictions. Additionally, the cost of installing monitoring wells is high because of groundwater depths and the variable drilling techniques that are required to characterize and install wells in a complex geologic environment. A major challenge to the characterization approach at LANL has been to meet sample quality requirements while minimizing drilling costs. In the late 1990's threaded drill-casing was advanced by rotary drilling methods using compressed air as the only circulation fluid. This resulted in relatively low impacts to the geochemical integrity of rock and groundwater samples, but led to frequent episodes of stuck casing. Additional fluids were introduced to the air-rotary program, improving penetration rates, but also resulting in non-representative groundwater samples due to long-term impacts to the aquifer in the vicinity of the well screen. Mud-rotary drilling was incorporated into the program, but it also had significant long-term impacts on groundwater quality, and often inhibited identification of perched groundwater zones and accurate identification of the water table in the regional aquifer. Current drilling practice makes limited use of fluid additives (e.g. foam) to supplement air circulation methods in order to advance drill casing through the vadose zone. Starting about 100 feet above the water table, drill casing is advanced to the target horizon in the regional aquifer using only air and municipal water (when needed) for circulation. Additionally, the threaded drill casing couplings have been replaced with welded connections, significantly reducing the occurrence of stuck casing. Careful observations of water production and water levels during drilling, coupled with the ability to retract casing for video and geophysical logs, now allows for robust characterization of perched groundwater systems and accurate definition of the top of regional saturation. Recent advances in the drilling program include combined dual-rotary/sonic drilling to depths up to 1150 feet and dual-rotary casing advance at angles up to 25 deg. from vertical. The former allows for core collection with minimal fluid addition, thereby preserving pore water chemistry. The latter method allows for well installation into portions of the aquifer otherwise inaccessible due to constraints at the surface (e.g., cultural sites or topography). Similar to the development of drilling techniques, both well design and groundwater sampling system design have evolved throughout the program. Current well design emphasizes a minimal annulus of filter pack, maximum screen slot size based on formation sieve analyses, and thorough well development to mitigate formation damage due to drilling. To sample groundwater at multiple levels in the same well, earlier sampling systems included Westbay{sup R} multi-zone systems and dual-screen single-pump and dual-screen dual-pump systems designed in collaboration with Baski, Inc. These multi-zone sampling systems proved difficult to develop and maintain given the depth to groundwater, and most current well installations are configured with a single screen and a submersible pump. The ability to purge from a well-developed single-screen well has substantially improved the collection of representative groundwater samples. A carefully crafted groundwater characterization program in target aquifers at depths ranging from 1,000 to 1,400 feet bgs within a complex geologic setting, while also retaining the ability to collect representative core and groundwater samples, may have application at a number of sites throughout the environmental industry. Optimizing each well to meet program objectives while reducing total project costs benefits all environmental investigations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2017

Western North Pacific Tropical Cyclone Model Tracks in Present and Future Climates

Western North Pacific tropical cyclone (TC) model tracks are analyzed in two large multimodel ensembles, spanning a large variety of models and multiple future climate scenarios. Two methodologies are used to synthesize the properties of TC tracks in this large data set: cluster analysis and mass moment ellipses. First, the models' TC tracks are compared to observed TC tracks' characteristics, and a subset of the models is chosen for analysis, based on the tracks' similarity to observations and sample size. Potential changes in track types in a warming climate are identified by comparing the kernel smoothed probability distributions of various track variables in historical and future scenarios using a Kolmogorov-Smirnov significance test. Two track changes are identified. The first is a statistically significant increase in the north-south expansion, which can also be viewed as a poleward shift, as TC tracks are prevented from expanding equatorward due to the weak Coriolis force near the equator. The second change is an eastward shift in the storm tracks that occur near the central Pacific in one of the multimodel ensembles, indicating a possible increase in the occurrence of storms near Hawaii in a warming climate. The dependence of the results on which model and future scenario are considered emphasizes the necessity of including multiple models and scenarios when considering future changes in TC characteristics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Protonation Dynamics of Confined Ethanol–Water Mixtures in H-ZSM-5 from Machine Learning-Driven Metadynamics

Zeolites are indispensable heterogeneous catalysts in industrial chemical processes, valued for their strong Brønsted acidity, well-defined microporous frameworks, and tunable pore structures. Their catalytic activity arises primarily from Brønsted acid sites (BAS), typically present as bridging hydroxyl groups (Si–OH–Al). Under aqueous reaction conditions, these protons interact dynamically with water and alcohol molecules, leading to complex solvation and protonation behavior within confined pores. In this study, we investigate the protonation equilibrium occurring between ethanol and water at the BAS of acidic zeolites under varying hydration levels, i.e., C2H5OH–(H2O)n, n=1–4. Local structure was analyzed through an adaptive-learning global optimization algorithm, while enhanced sampling molecular dynamics simulations with Well-Tempered Metadynamics (WMetaD) and machine learning interatomic potentials (MLPs) provide free-energy surfaces (FES) at variable hydration levels. The results reveal a strong dependence of proton localization on the degree of hydration. At low hydration (1 water molecule), the proton resides predominantly on ethanol; with 2 water molecules, it shifts toward water, and at higher hydration (3 or more water molecules), it becomes extensively delocalized over the water cluster. These findings underscore the critical role of solvation in modulating acid site behavior and suggest that a minimum of three water molecules is necessary to fully stabilize the proton on water within the zeolite framework. This solvation threshold has significant implications for catalytic processes, particularly in biomass conversion reactions where alcohol protonation is a key step in dehydration mechanisms.

machine learning↗
DOE OSTI2023

Automatic detection of cataclysmic variables from SDSS images

Abstract Investigating rare and new objects have always been an important direction in astronomy. Cataclysmic variables (CVs) are ideal and natural celestial bodies for studying the accretion process of semi-detached binaries with accretion processes. However, the sample size of CVs must increase because a lager gap exists between the observational and the theoretical expanding CVs. Astronomy has entered the big data era and can provide massive images containing CV candidates. CVs as a type of faint celestial objects, are highly challenging to be identified directly from images using automatic manners. Deep learning has rapidly developed in intelligent image processing and has been widely applied in some astronomical fields with excellent detection results. YOLOX, as the latest YOLO framework, is advantageous in detecting small and dark targets. This work proposes an improved YOLOX-based framework according to the characteristics of CVs and Sloan Digital Sky Survey (SDSS) photometric images to train and verify the model to realise CV detection. We use the Convolutional Block Attention Module to increase the number of output features with the feature extraction network and adjust the feature fusion network to obtain fused features. Accordingly, the loss function is modified. Experimental results demonstrate that the improved model produces satisfactory results, with average accuracy (mean average Precision at 0.5) of 92.0%, Precision of 92.9%, Recall of 94.3%, and $F1-score$ of 93.6% on the test set. The proposed method can efficiently achieve the identification of CVs in test samples and search for CV candidates in unlabeled images. The image data vastly outnumber the spectra in the SDSS-released data. With supplementary follow-up observations or spectra, the proposed model can help astronomers in seeking and detecting CVs in a new manner to ensure that a more extensive CV catalog can be built. The proposed model may also be applied to the detection of other kinds of celestial objects.

Astronomy & Astrophysics↗
DOE OSTI2019

Radiation damage effects in chlorite investigated using microfocus synchrotron techniques

A detailed understanding of the mechanisms and effects of radiation damage in phyllosilicate minerals is a necessary component of the evaluation of the safety case for a deep geological disposal facility (GDF) for radioactive waste. Structural and chemical changes induced by a-particle damage will affect these mineral’s performance as a reactive barrier material (both in the near and far-field) over timescales relevant to GDF integrity. In this study, two examples of chlorite group minerals have been irradiated at a-particle doses comparable to those predicted to be experienced by the clay buffer material surrounding high level radioactive waste canisters. Crystallographic aberrations induced by the focused 4He2+ ion beam are revealed via high-resolution, microfocus X-ray diffraction mapping. Interlayer collapse by up to 0.5 Å is prevalent across both macro- and micro-crystalline samples, with the macro-crystalline specimen displaying a breakdown of the phyllosilicate structure into loosely-connected, multi-oriented crystallites displaying variable lattice parameters. The damaged lattice parameters suggest a localised breakdown and collapse of the OH- rich, ‘brucite-like’ interlayer. Microfocus Fe K-edge X-ray absorption spectroscopy illustrates this defect accumulation, manifest as a severe damping of the X-ray absorption edge. Subtle Fe2+/Fe3+ speciation changes are apparent across the damaged structures. A trend towards Fe reduction is evident at depth in the damaged structures at certain doses (8.76 x 1015 alpha particles/ cm2). Interestingly, this reductive trend does not increase with radiation dose, indeed at the maximum dose (1.26 x 1016 alpha particles/ cm2) administered in this study, there is evidence for a slight increase in Fe binding energy, suggesting the development of a depth-dependant redox gradient concurrent with light ion damage. At the doses examined here, these damaged structures are likely highly reactive, as sorption capacity will, to an extent, be largely enhanced by lattice disruption and an increase in available ‘edge’ sites.

Bower, William R.↗
DOE OSTI2024

Metadynamics investigation of lanthanide solvation free energy landscapes and insights into separations energetics

Lanthanide ion solvation chemistry in nonaqueous phases is key to understanding and developing effective separation processes for these critical materials. Due to the complexity and inherent disorder of the solution phase, a comprehensive picture of the solvated metal ion is often difficult to generate solely from conventional spectroscopic approaches and electronic structure calculations, particularly in the extractant phase. In this work, we use classical molecular dynamics (MD) simulation with an advanced sampling technique, metadynamics, supplemented by experimental spectroscopy and speciation analysis, to measure lanthanide solvation free energy landscapes. We define coordination-based collective variables to probe the entire range of solvation configurations in the organic phase of lanthanum (La), europium (Eu), and lutetium (Lu) nitrate salts bound with a commonly used extractant, N,N′-dimethyl, N,N′-dioctylhexylethoxymalonamide (DMDOHEMA). The known lanthanide extraction trend of La ≈ Eu > Lu is readily explained by the measured free energy surfaces, which show consistent DMDOHEMA coordination from La to Eu, followed by loss of DMDOHEMA coordination from Eu to Lu. These simulations suggest how ligand crowding at the metal center can control selectivity, in this case resulting in the opposite extraction trend as observed with other conventional extractants, where the enthalpic contribution from increasing lanthanide charge density across the series dominates the extraction energetics. We also find that the presence of inner-sphere water, verified by time-resolved fluorescence, diversifies the accessible solvation structures. As a result, understanding solvation requires consideration of an entire thermodynamic ensemble, rather than the single dominant lowest-energy structure, as is often considered out of necessity in interpretation of spectroscopic data or in electronic structure-based ligand design approaches. In general, we demonstrate how metadynamics uniquely enables investigation of complex, multidimensional solvation energetic landscapes, and how it can explain selectivity trends where extraction is controlled by more complex mechanisms than simple charge density-based selectivity.

Wang, Xiaoyu↗
DOE OSTI2013

Making use of x-ray optical effects in photoelectron-, Auger electron-, and x-ray emission spectroscopies: Total reflection, standing-wave excitation, and resonant effects

We present a general theoretical methodology and related open-access computer program for carrying out the calculation of photoelectron, Auger electron, and x-ray emission intensities in the presence of several x-ray optical effects, including total reflection at grazing incidence, excitation with standing-waves produced by reflection from synthetic multilayers and at core-level resonance conditions, and the use of variable polarization to produce magnetic circular dichroism. Calculations illustrating all of these effects are presented, including in some cases comparisons to experimental results. Sample types include both semi-infinite flat surfaces and arbitrary multilayer configurations, with interdiffusion/roughness at their interfaces. These x-ray optical effects can significantly alter observed photoelectron, Auger, and x-ray intensities, and in fact lead to several generally useful techniques for enhancing surface and buried-layer sensitivity, including layer-resolved densities of states and depth profiles of element-specific magnetization. The computer program used in this study should thus be useful for a broad range of studies in which x-ray optical effects are involved or are to be exploited in next-generation surface and interface studies of nanoscale systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

ENIGMAS FROM THE SLOAN DIGITAL SKY SURVEY DR7 KLEINMAN WHITE DWARF CATALOG

We report results from a continuation of our searches for high field magnetic white dwarfs (WDs) paired in a detached binary with non-degenerate companions. We made use of the Sloan Digital Sky Survey DR7 catalog of Kleinman et al. (2013) with 19,712 spectroscopically identified WDs. These include 1735 WD plus M dwarf detached pairs (almost 10% of the Kleinman et al.'s list). No new pairs were found, although we did recover the polar (AM Herculis system) ST LMi in a low state of accretion. With the larger sample the original situation reported 10 yr ago remains intact now at a much higher level of statistical significance: in the selected SDSS sample, high field magnetic WDs are not found in an apparently detached pairing with an M dwarf, unless they are a magnetic cataclysmic variable (CV) in a low state of accretion. This finding strengthens the case that the fields in the isolated high field magnetic WDs are generated by stellar mergers but also raises questions on the nature of the progenitors of the magnetic CVs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A Novel Front-End ASIC With Post Digital Filtering and Calibration for CZT-Based PET Detector

This paper presents a novel front-end electronics based on a front-end ASIC with post digital filtering and calibration dedicated to CZT detectors for PET imaging. A cascade amplifier based on split-leg topology is selected to realize the charge-sensitive amplifier (CSA) for the sake of low noise performances and the simple scheme of the power supplies. The output of the CSA is connected to a variable-gain amplifier to generate the compatible signals for the A/D conversion. A multi-channel single-slope ADC is designed to sample multiple points for the digital filtering and shaping. The digital signal processing algorithms are implemented by a FPGA. To verify the proposed scheme, a front-end readout prototype ASIC is designed and implemented in 0.35 μm CMOS process. In a single readout channel, a CSA, a VGA, a 10-bit ADC and registers are integrated. Two dummy channels, bias circuits, and time controller are also integrated. The die size is 2.0 mm x 2.1 mm. The input range of the ASIC is from 2000 e{sup -} to 100000 e{sup -}, which is suitable for the detection of the X-and gamma ray from 11.2 keV to 550 keV. The linearity of the output voltage is less than 1 %. The gain of the readout channel is 40.2 V/pC. The static power dissipation is about 10 mW/channel. The above tested results show that the electrical performances of the ASIC can well satisfy PET imaging applications. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Mechanical Response of Castlegate Sandstone under Hydrostatic Cyclic Loading

The stress history of rocks in the subsurface affects their mechanical and petrophysical properties. Rocks can often experience repeated cycles of loading and unloading due to fluid pressure fluctuations, which will lead to different mechanical behavior from static conditions. This is of importance for several geophysical and industrial applications, for example, wastewater injection and reservoir storage wells, which generate repeated stress perturbations. Laboratory experiments were conducted with Castlegate sandstone to observe the effects of different cyclic pressure loading conditions on a common reservoir analogue. Each sample was hydrostatically loaded in a triaxial cell to a low effective confining pressure, and either pore pressure or confining pressure was cycled at different rates over the course of a few weeks. Fluid permeability was measured during initial loading and periodically between stress cycles. Samples that undergo cyclic loading experience significantly more inelastic (nonrecoverable) strain compared to samples tested without cyclic hydrostatic loading. Permeability decreases rapidly for all tests during the first few days of testing, but the decrease and variability of permeability after this depend upon the loading conditions of each test. Cycling conditions do affect the mechanical behavior; the elastic moduli decrease with the increasing loading rate and stress cycling. The degree of volumetric strain induced by stress cycles is the major control on permeability change in the sandstones, with less compaction leading to more variation from measurement to measurement. The data indicate that cyclic loading degrades permeability and porosity more than static conditions over a similar period, but the petrophysical properties are dictated more by the hydrostatic loading rate rather than the total length of time stress cycling is imposed.

58 GEOSCIENCES↗
DOE OSTI2021

Anomaly Attribution with Likelihood Compensation

This paper addresses the task of explaining anomalous predictions of a black-box regression model. When using a black-box model, such as one to predict building energy consumption from many sensor measurements, we often have a situation where some observed samples may significantly deviate from their prediction. It may be due to a sub-optimal black-box model, or simply because those samples are outliers. In either case, one would ideally want to compute a responsibility score indicative of the extent to which an input variable is responsible for the anomalous output. In this work, we formalize this task as a statistical inverse problem: Given model deviation from the expected value, infer the responsibility score of each of the input variables. We propose a new method called likelihood compensation (LC), which is founded on the likelihood principle and computes a correction to each input variable. To the best of our knowledge, this is the first principled framework that computes a responsibility score for real valued anomalous model deviations. We apply our approach to a real-world building energy prediction task and confirm its utility based on expert feedback.

Idé, Tsuyoshi↗
DOE OSTI2017

Stan : A Probabilistic Programming Language

Stan is a probabilistic programming language for specifying statistical models. A Stan program imperatively defines a log probability function over parameters conditioned on specified data and constants. As of version 2.14.0, Stan provides full Bayesian inference for continuous-variable models through Markov chain Monte Carlo methods such as the No-U-Turn sampler, an adaptive form of Hamiltonian Monte Carlo sampling. Penalized maximum likelihood estimates are calculated using optimization methods such as the limited memory Broyden-Fletcher-Goldfarb-Shanno algorithm. Stan is also a platform for computing log densities and their gradients and Hessians, which can be used in alternative algorithms such as variational Bayes, expectation propagation, and marginal inference using approximate integration. To this end, Stan is set up so that the densities, gradients, and Hessians, along with intermediate quantities of the algorithm such as acceptance probabilities, are easily accessible. Stan can also be called from the command line using the cmdstan package, through R using the rstan package, and through Python using the pystan package. All three interfaces support sampling and optimization-based inference with diagnostics and posterior analysis. rstan and pystan also provide access to log probabilities, gradients, Hessians, parameter transforms, and specialized plotting.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Microbiomes Associated With Foods From Plant and Animal Sources

Food microbiome composition impacts food safety and quality. The resident microbiota of many food products is influenced throughout the farm to fork continuum by farming practices, environmental factors, and food manufacturing and processing procedures. Currently, most food microbiology studies rely on culture-dependent methods to identify bacteria. However, advances in high-throughput DNA sequencing technologies have enabled the use of targeted 16S rRNA gene sequencing to profile complex microbial communities including non-culturable members. In this study we used 16S rRNA gene sequencing to assess the microbiome profiles of plant and animal derived foods collected at two points in the manufacturing process; post-harvest/pre-retail (cilantro) and retail (cilantro, masala spice mixes, cucumbers, mung bean sprouts, and smoked salmon). Our findings revealed microbiome profiles, unique to each food, that were influenced by the moisture content (dry spices, fresh produce), packaging methods, such as modified atmospheric packaging (mung bean sprouts and smoked salmon), and manufacturing stage (cilantro prior to retail and at retail). The masala spice mixes and cucumbers were comprised mainly of Proteobacteria, Firmicutes, and Actinobacteria. Cilantro microbiome profiles consisted mainly of Proteobacteria, followed by Bacteroidetes, and low levels of Firmicutes and Actinobacteria. The two brands of mung bean sprouts and the three smoked salmon samples differed from one another in their microbiome composition, each predominated by either by Firmicutes or Proteobacteria. These data demonstrate diverse and highly variable resident microbial communities across food products, which is informative in the context of food safety, and spoilage where indigenous bacteria could hamper pathogen detection, and limit shelf life.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

A SEARCH FOR SHORT-PERIOD ROCKY PLANETS AROUND WDs WITH THE COSMIC ORIGINS SPECTROGRAPH (COS)

The search for transiting habitable exoplanets has broadened to include several types of stars that are smaller than the Sun in an attempt to increase the observed transit depth and hence the atmospheric signal of the planet. Of all spectral types, white dwarfs (WDs) are the most favorable for this type of investigation. The fraction of WDs that possess close-in rocky planets is unknown, but several large angle stellar surveys have the photometric precision and cadence to discover at least one if they are common. Ultraviolet observations of WDs may allow for detection of molecular oxygen or ozone in the atmosphere of a terrestrial planet. We use archival Hubble Space Telescope data from the Cosmic Origins Spectrograph to search for transiting rocky planets around UV-bright WDs. In the process, we discovered unusual variability in the pulsating WD GD 133, which shows slow sinusoidal variations in the UV. While we detect no planets around our small sample of targets, we do place stringent limits on the possibility of transiting planets, down to sub-lunar radii. We also point out that non-transiting small planets in thermal equilibrium are detectable around hotter WDs through infrared excesses, and identify two candidates.

79 ASTRONOMY AND ASTROPHYSICS↗