Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model fitting”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Carbon exchange in an Amazon forest: from hours to years

Abstract. In Amazon forests, the relative contributions of climate, phenology, and disturbance to net ecosystem exchange of carbon (NEE) are not well understood. To partition influences across various timescales, we use a statistical model to represent eddy-covariance-derived NEE in an evergreen eastern Amazon forest as a constant response to changing meteorology and phenology throughout a decade. Our best fit model represented hourly NEE variations as changes due to sunlight, while seasonal variations arose from phenology influencing photosynthesis and from rainfall influencing ecosystem respiration, where phenology was asynchronous with dry-season onset. We compared annual model residuals with biometric forest surveys to estimate impacts of drought disturbance. We found that our simple model represented hourly and monthly variations in NEE well ( R 2 = 0.81 and 0.59, respectively). Modeled phenology explained 1 % of hourly and 26 % of monthly variations in observed NEE, whereas the remaining modeled variability was due to changes in meteorology. We did not find evidence to support the common assumption that the forest phenology was seasonally light- or water-triggered. Our model simulated annual NEE well, with the exception of 2002, the first year of our data record, which contained 1.2 MgC ha −1 of residual net emissions, because photosynthesis was anomalously low. Because a severe drought occurred in 1998, we hypothesized that this drought caused a persistent, multi-year depression of photosynthesis. Our results suggest drought can have lasting impacts on photosynthesis, possibly via partial damage to still-living trees.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

XRISM Spectroscopy of the Fe Kα Emission Line in the Seyfert Active Galactic Nucleus NGC 4151 Reveals the Disk, Broad-line Region, and Torus

Abstract We present an analysis of the first two XRISM/Resolve spectra of the well-known Seyfert-1.5 active galactic nucleus (AGN) in NGC 4151, obtained in 2023 December. Our work focuses on the nature of the narrow Fe K α emission line at 6.4 keV, the strongest and most common X-ray line observed in AGN. The total line is found to consist of three components. Even the narrowest component of the line is resolved with evident Fe K α ,1 (6.404 keV) and K α ,2 (6.391 keV) contributions in a 2:1 flux ratio, fully consistent with neutral gas with negligible bulk velocity. Subject to the limitations of our models, the narrowest and intermediate-width components are consistent with emission from optically thin gas, suggesting that they arise in a disk atmosphere and/or wind. Modeling the three line components in terms of Keplerian broadening, they are readily associated with (1) the inner wall of the “torus,” (2) the innermost optical “broad-line region” (or “X-ray BLR”), and (3) a region with a radius of r ≃ 100 GM / c 2 that may signal a warp in the accretion disk. Viable alternative explanations of the broadest component include a fast-wind component and/or scattering; however, we find evidence of variability in the narrow Fe K α line complex on timescales consistent with small radii. The best-fit models are statistically superior to simple Voigt functions, but when fit with Voigt profiles the time-averaged lines are consistent with a projected velocity broadening of FWHM = 1600 − 200 + 400 km s − 1 . Overall, the resolution and sensitivity of XRISM show that the narrow Fe K line in AGN is an effective probe of all key parts of the accretion flow, as it is currently understood. We discuss the implications of these findings for our understanding of AGN accretion, future studies with XRISM, and X-ray-based black hole mass measurements.

Audard, Marc (ORCID:000000034721034X)↗
DOE OSTI2026

Low energy neutron light output characterization of EJ301D and deuterated stilbene with a comparison of light output characterization methods

The neutron-induced light yield of a 2.54 cm diameter by 2.54 cm long right circular cylinder of EJ301D and a (5.08 cm)3 custom made cube of deuterated trans-stilbene-d12 (d-stilbene) were measured over incident neutron energies from 300 keV to 2.2 MeV and 200 keV to 2.4 MeV, respectively. The measurements were performed using a time-of-flight experiment with a Cf-252 source and an approximately 1.5 m flight path. We compare three light output spectrum full energy deposition edge estimation methods: (1) simulating the neutron energy spectrum edge and fitting it to the light output spectrum, (2) using the inflection point of the light output spectrum edge (derivative method, a.k.a. Kornilov’s method), and (3) using an empirical model fit to the edge of the light output spectrum. Both the derivative and equation fit methods do not account for physical processes such as multiple neutron scattering in the detectors. They instead rely on assumptions about the linear shape continuum shape of the light output spectrum and the direct correlation between the location of the spectrum’s inflection point and maximum energy deposition. These assumptions were found to introduce bias into those methods when tested against simulated spectra with known edge locations. When tested against measured spectra the derivative method was found to differ from the simulation fit by greater than 30% at low energies with large discontinuities for adjacent data points above 800 keV incident neutron energy. The empirical equation fitting method was found to also exhibit bias of a similar magnitude, but with significantly more continuous behavior, especially with the lower count data of the smaller volumed EJ301D scintillator. Experimental light output yield for this neutron energy range is reported using the simulated spectrum fitting method because it includes physics neglected by the other methods, and did not exhibit the bias observed in the other methods

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Multiple emission components in the Cygnus cocoon detected from Fermi -LAT observations

Star-forming regions may play an important role in the life cycle of Galactic cosmic rays (CRs), notably as home to specific acceleration mechanisms and transport conditions. Gamma-ray observations of Cygnus X have revealed the presence of an excess of hard-spectrum gamma-ray emission, possibly related to a cocoon of freshly accelerated particles. We seek an improved description of the gamma-ray emission from the cocoon using ~13 yr of observations with the Fermi-Large Area Telescope (LAT) and use it to further constrain the processes and objects responsible for the young CR population. We developed an emission model for a large region of interest, including a description of interstellar emission from the background population of CRs and recent models for other gamma-ray sources in the field. Thus, we performed an improved spectro-morphological characterisation of the residual emission including the cocoon. The best-fit model for the cocoon includes two main emission components: an extended component FCES G78.74+1.56, described by a 2D Gaussian of extension r 68 = 4.4° ± 0.1° -0.1° +0.1° and a smooth broken power law spectrum with spectral indices 1.67 ± 0.05 -0.01 +0.02 and 2.12 ± 0.02 -0.01 +0.00 below and above 3.0 ± 0.6 -0.2 +0.0 GeV, respectively; and a central component FCES G80.00+0.50, traced by the distribution of ionised gas within the borders of the photo-dissociation regions and with a power law spectrum of index 2.19 ± 0.03 -0.01 +0.00 that is significantly different from the spectrum of FCES G78.74+1.56. An additional extended emission component FCES G78.83+3.57, located on the edge of the central cavities in Cygnus X and with a spectrum compatible with that of FCES G80.00+0.50, is likely related to the cocoon. For the two brightest components FCES G80.00+0.50 and FCES G78.74+1.56, spectra and radial-azimuthal profiles of the emission can be accounted for in a diffusion-loss framework involving one single population of non-thermal particles with a flat injection spectrum. Particles span the full extent of FCES G78.74+1.56 as a result of diffusion from a central source, and give rise to source FCES G80.00+0.50 by interacting with ionised gas in the innermost region. For this simple diffusion-loss model, viable setups can be very different in terms of energetics, transport conditions, and timescales involved, and both hadronic and leptonic scenarios are possible. The solutions range from long-lasting particle acceleration, possibly in prominent star clusters such as Cyg OB2 and NGC 6910, to a more recent and short-lived release of particles within the last 10–100 kyr, likely from a supernova remnant. The observables extracted from our analysis can be used to perform detailed comparisons with advanced models of particle acceleration and transport in star-forming regions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Functionalization of Single-walled Carbon Nanotubes with Thermo-reversible Block Copolymers and Characterization by Small-angle Neutron Scattering

We demonstrate a protocol for single-walled carbon nanotube functionalization using thermo-sensitive PEO-PPO-PEO triblock copolymers in an aqueous solution In a carbon nanotube/PEO105-PPO70-PEO105 (poloxamer 407) aqueous solution, the amphiphilic poloxamer 407 adsorbs onto the carbon nanotube surfaces and self-assembles into continuous layers, driven by intermolecular interactions between constituent molecules. The addition of 5-methylsalicylic acid changes the self-assembled structure from spherical-micellar to a cylindrical morphology. The fabricated poloxamer 407/carbon nanotube hybrid particles exhibit thermo-responsive structural features so that the density and thickness of poloxamer 407 layers are also reversibly controllable by varying temperature. The detailed structural properties of the poloxamer 407/carbon nanotube particles in suspension can be characterized by small-angle neutron scattering experiments and model fit analyses. The distinct curve shapes of the scattering intensities depending on temperature control or addition of aromatic additives are well described by a modified core-shell cylinder model consisting of a carbon nanotube core cylinder, a hydrophobic shell, and a hydrated polymer layer. In conclusion, this method can provide a simple but efficient way for the fabrication and in-situ characterization of carbon nanotube-based nano particles with a structure-tunable encapsulation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Complex Sources and Impacts on Remediation and Monitoring Approaches

Remedy selection, implementation, and coordination with tank farm activities for adjacent waste sites need to be facilitated by a technically defensible description of sources and their relative strength and temporal contributions to groundwater contamination. Detailed data describing contaminant source releases are typically not available, so source locations are highly uncertain. Furthermore, data on sources and model parameters are not readily available through direct measurements. Instead, source data are often derived by extrapolation from a few direct and/or indirect measurements, while model parameters are commonly estimated using inverse modeling or calibration methods. Because the conceptual model of sources will influence the mass and timing of the flux of contaminants to groundwater, uncertainty in the source terms should be considered. Even though a best-fit model can often be identified, it may not be the best representation of the sources given the sparse, uncertain nature of the data used for estimating source terms. To support future feasibility studies, a systematic approach is needed for identifying remediation and monitoring activities appropriate under conditions of uncertainty. In this study, an integrated assessment of alternative conceptual model elements has been examined using the B-Complex site in the Hanford Central Plateau to demonstrate how multiple sources in a waste management complex can be described to identify their relative contributions to groundwater. The heterogeneous structure of the system is modeled based on available field information on the lithological data for both the vadose zone and aquifer. The approach includes a simulation assessment and linkages to site data, with an uncertainty analysis applied to evaluate the effects of source term uncertainty on contaminant mass flux to groundwater. An integrated assessment of alternative conceptual model elements was evaluated to identify how parameter uncertainty impacts both temporal and spatial uncertainties associated with vadose zone source terms. Modeling scenarios identify important characteristics of remediation target zones to provide context for how surface barriers, in situ remedies, and water-table treatment zones can be designed to account for source dynamics. Ramifications for remediation and monitoring are discussed within the context of the source term uncertainty. The described approach facilitates systematically addressing alternative conceptual model questions in relation to remediation and monitoring and supports the remedial investigation/ feasibility study process and the remediation strategy development for the Hanford Central Plateau.

Rockhold, Mark L.↗
DOE OSTI2019

Complex Sources and Impacts on Remediation and Monitoring Approaches - 19058

Remedy selection, implementation, and coordination with tank farm activities for adjacent waste sites need to be facilitated by a technically defensible description of sources and their relative strength and temporal contributions to groundwater contamination. Detailed data describing contaminant source releases are typically not available, so source locations are highly uncertain. Furthermore, data on sources and model parameters are not readily available through direct measurements. Instead, source data are often derived by extrapolation from a few direct and/or indirect measurements, while model parameters are commonly estimated using inverse modeling or calibration methods. Because the conceptual model of sources will influence the mass and timing of the flux of contaminants to groundwater, uncertainty in the source terms should be considered. Even though a best-fit model can often be identified, it may not be the best representation of the sources given the sparse, uncertain nature of the data used for estimating source terms. To support future feasibility studies, a systematic approach is needed for identifying remediation and monitoring activities appropriate under conditions of uncertainty. In this study, an integrated assessment of alternative conceptual model elements has been examined using the B-Complex site in the Hanford Central Plateau to demonstrate how multiple sources in a waste management complex can be described to identify their relative contributions to groundwater. The heterogeneous structure of the system is modeled based on available field information on the lithological data for both the vadose zone and aquifer. The approach includes a simulation assessment and linkages to site data, with an uncertainty analysis applied to evaluate the effects of source term uncertainty on contaminant mass flux to groundwater. An integrated assessment of alternative conceptual model elements was evaluated to identify how parameter uncertainty impacts both temporal and spatial uncertainties associated with vadose zone source terms. Modeling scenarios identify important characteristics of remediation target zones to provide context for how surface barriers, in situ remedies, and water-table treatment zones can be designed to account for source dynamics. Ramifications for remediation and monitoring are discussed within the context of the source term uncertainty. The described approach facilitates systematically addressing alternative conceptual model questions in relation to remediation and monitoring and supports the remedial investigation/ feasibility study process and the remediation strategy development for the Hanford Central Plateau. (authors)

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

The Pantheon+ Analysis: Forward Modeling the Dust and Intrinsic Color Distributions of Type Ia Supernovae, and Quantifying Their Impact on Cosmological Inferences

Abstract Recent studies have shown that the observed color distributions of Type Ia supernovae (SNe Ia) can be well described by a combination of a dust distribution and an intrinsic color distribution. Using the Pantheon+ sample of 1701 SN Ia, we apply a new forward-modeling fitting method (Dust2Dust) to measure the parent dust and color distributions, including their dependence on host-galaxy mass. At each fit step, the SN Ia selection efficiency is determined from a large simulated sample that is reweighted to reflect the proposed distributions. We use five separate metrics to describe the goodness of fit: distribution of fitted light-curve color c , cosmological residual trends with c , cosmological residual scatter with c , fitted color–luminosity relationship β SALT2 , and intrinsic scatter σ int . We present the results and the uncertainty in 12-dimensional space. Furthermore, we measure that the uncertainty on this modeling propagates to an upper threshold uncertainty in the equation of state of dark energy w of 0.014(1) for the Pantheon+ cosmology analysis and contributes negligible uncertainty to the Hubble constant H 0 . The Dust2Dust code is made publicly available at https://github.com/djbrout/dustdriver .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

SEC ‐ SAXS / MC Ensemble Structural Studies of the Microtubule Binding Protein Cdt1 Show Monomeric, Folded‐Over Conformations

ABSTRACT Cdt1 is a mixed folded protein critical for DNA replication licensing and it also has a “moonlighting” role at the kinetochore via direct binding to microtubules and the Ndc80 complex. However, it is unknown how the structure and conformations of Cdt1 could allow it to participate in these multiple, unique sets of protein complexes. While robust methods exist to study entirely folded or unfolded proteins, structure–function studies of combined, mixed folded/disordered proteins remain challenging. In this work, we employ orthogonal biophysical and computational techniques to provide structural characterization of mitosis‐competent human Cdt1. Thermal stability analyses shows that both folded winged helix domains1 are unstable. CD and NMR show that the N‐terminal and linker regions are intrinsically disordered. DLS shows that Cdt1 is monomeric and polydisperse, while SEC‐MALS confirms that it is monomeric at high concentrations, but without any apparent inter‐molecular self‐association. SEC‐SAXS enabled computational modeling of the protein structures. Using the program SASSIE, we performed rigid body Monte Carlo simulations to generate a conformational ensemble of structures. We observe that neither fully extended nor extremely compact Cdt1 conformations are consistent with SAXS. The best‐fit models have the N‐terminal and linker disordered regions extended into the solution and the two folded domains close to each other in apparent “folded over” conformations. We hypothesize the best‐fit Cdt1 conformations could be consistent with a function as a scaffold protein that may be sterically blocked without binding partners. Our study also provides a template for combining experimental and computational techniques to study mixed‐folded proteins.

Cell Biology↗
DOE OSTI2025

Generating multi-scale Li-ion battery cathode particles with radial grain architectures using stereological generative adversarial networks

Abstract Understanding structure-property relationships of Li-ion battery cathodes is crucial for optimizing rate-performance and cycle-life resilience. However, correlating the morphology of cathode particles, such as in LiNi0.8Mn0.1Co0.1O2 (NMC811), and their inner grain architecture with electrode performance is challenging, particularly, due to the significant length-scale difference between grain and particle sizes. Experimentally, it is not feasible to image such a high number of particles with full granular detail. A second challenge is that sufficiently high-resolution 3D imaging techniques remain expensive and are sparsely available at research institutions. Here, we present a stereological generative adversarial network-based model fitting approach to tackle this, that generates representative 3D information from 2D data, enabling characterization of materials in 3D using cost-effective 2D data. Once calibrated, this multi-scale model can rapidly generate virtual cathode particles that are statistically similar to experimental data, and thus is suitable for virtual characterization and materials testing through numerical simulations. A large dataset of simulated particles with inner grain architecture has been made publicly available.

25 ENERGY STORAGE↗
DOE OSTI2015

Hubble Space Telescope near-ultraviolet spectroscopy of bright CEMP-s stars

We present an elemental-abundance analysis, in the near-ultraviolet (NUV) spectral range, for the bright carbon-enhanced metal-poor (CEMP) stars HD 196944 (V=8.40, [Fe/H] = −2.41) and HD 201626 (V=8.16, [Fe/H] = −1.51), based on data acquired with the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope. Both of these stars belong to the sub-class CEMP-s, and exhibit clear over-abundances of heavy elements associated with production by the slow neutron-capture process. HD 196944 has been well-studied in the optical region, but we add abundance results for six species (Ge, Nb, Mo, Lu, Pt, and Au) that are only accessible in the NUV. In addition, we provide the first determination of its orbital period, P = 1325 days. HD 201626 has only a limited number of abundance results based on previous optical work—here we add five new species from the NUV, including Pb. We compare these results with models of binary-system evolution and s-process element production in stars on the asymptotic giant branch, with the goal of explaining their origin and evolution. Our best-fitting models for HD 196944 (M{sub 1,i}=0.9M{sub ⊙}, M{sub 2,i}=0.86M{sub ⊙}, for [Fe/H] = −2.2), and HD 201626 (M{sub 1,i}=0.9M{sub ⊙}, M{sub 2,i}=0.76M{sub ⊙}, for [Fe/H] = −2.2; M{sub 1,i}=1.6M{sub ⊙}, M{sub 2,i}=0.59M{sub ⊙}, for [Fe/H] = −1.5) are consistent with the current accepted scenario for the formation of CEMP-s stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Review of PBX 9502 Pop Plot data

Pop plot data (distance of run and time to detonation as a function of shock pressure) is a simple characterization of shock initiation of a high explosive. It is used, along with other data, to calibrate reactive burn models. For the insensitive TATB based explosive PBX 9502, there are data at 5 initial temperatures (77 K, -55 C, 23 C, 76 C, 130 C) and for several lots. Here we review the dependence of the Pop plot on temperature and the variation with lot. There is a significant increase in shock sensitivity with temperature.Since the Pop plot is important for modeling, fit parameters are interpolated to cover the temperature range –55 < T < 76 C. The scatter in the data and variation with lot is important for assessing the accuracy to which a model can be calibrated.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2016

Investigation of SiPM physics parameters down to cryogenic temperatures and for a bio-medical application

Silicon PhotoMultiplier (SiPM) detector has become a suitable visible light/photon detector for many applications like high energy physics and neutrino experiments, fluorescence detection, bio-photonics and medical imaging. The first part of my thesis was oriented to the studies of SiPM physics parameters as a function of temperature. Particularly, recent KETEK devices (year 2015) with different technological characteristics like p/n and n/p junctions, with and without trench technology, and different widths of epitaxial layer were studied in the temperature range from 308.15 K (+35°C) down to 238.15 K (-35°C). In addition, the Hamamatsu devices from 2011 production run as well as new devices from 2015 year, with improved technological characteristics inducing a reduced noise, were investigated in a wider temperature range from 318.15 K (+45°C) down to 98.15 K (-175°C). For these purposes, I participated to the design, installation, commissioning and calibration of a cryogenic experimental setup dedicated to electrical, optical and temperature studies of SiPM devices. Also, I have developed an automatic analysis procedure able to handle in a short time an impressive quantity of experimental data (i.e. tens of Gb/device) and to give a precise and fast information on main SiPM parameters and their temperature dependence. I have also developed a physical modeldescribing the DC I-V curves of SiPM detectors at different temperatures. The proposed model fits well the shape of IV curve in a very large currents range from 10⁻¹² A up to 10⁻⁵ A over the full working range of various devices. Consequently, the IV model can be used as a simple and fast method for determination of SiPM parameters like breakdown voltage VBD, the shape of Geiger triggering probability PGeiger as a function of Vbias as well as the Vbias working range. The comparison of these parameters with those calculated from AC measurements and analyzed by the automatic procedure showed a good agreement. The second part of my thesis was oriented to the study of SiPM devices and their physical parameters required to build a prototype of betasensitive intracerebral probe. Such probe is dedicated to measure the local concentration of radiolabeled molecules on awake and freely moving animal and to study new animal models of human disorders (neurodegenerative diseases, tumor growth, and neuropsychiatric disorders). It is composed of small size, low-noise SiPM device coupled to a scintillating fiber and readout by a dedicated miniaturized low-power consumption electronics. Three SiPM devices have been chosen as the most adapted for our application: two small KETEK devices of 0.5×0.5 mm² size (with and without optical trenches, specially developed by KETEK to fulfill our requirements) and a standard Hamamatsu device of 1.3×1.3 mm² size, all devices having 50 × 50 μm² μcell size. For each SiPM the gain G, dark count rate DCR and beta sensitivity were measured as a function of Vbias and temperature. The obtained results showed that the small field of view and newly developed structure of the KETEK devices allow a large decrease of the dark count rate DCR. However, this small field of view also leads to a reduced light collection due to the thickness of the epoxy protection resin on top of the SiPM and the acceptance angle of the fiber. Since the beta sensitivity represents a tradeoff between photon detection efficiency PDE and dark count rate DCR, KETEK SiPMs exhibit similar performances in comparison with the Hamamatsu device. Preliminary results demonstrate that the beta sensitivity of KETEK devices can be significantly improved by using focusing lens between the scintillating fiber and the SiPM or by reducing the thickness of its epoxy protection resin.

Nagai, Andrii↗
DOE OSTI2018

An Indication of Anisotropy in Arrival Directions of Ultra-high-energy Cosmic Rays through Comparison to the Flux Pattern of Extragalactic Gamma-Ray Sources

A new analysis of the dataset from the Pierre Auger Observatory provides evidence for anisotropy in the arrival directions of ultra-high-energy cosmic rays on an intermediate angular scale, which is indicative of excess arrivals from strong, nearby sources. The data consist of 5514 events above 20 EeV with zenith angles up to 80 deg recorded before 2017 April 30. Sky models have been created for two distinct populations of extragalactic gamma-ray emitters: active galactic nuclei from the second catalog of hard Fermi-LAT sources (2FHL) and starburst galaxies from a sample that was examined with Fermi-LAT. Flux-limited samples, which include all types of galaxies from the Swift-BAT and 2MASS surveys, have been investigated for comparison. The sky model of cosmic-ray density constructed using each catalog has two free parameters, the fraction of events correlating with astrophysical objects and an angular scale characterizing the clustering of cosmic rays around extragalactic sources. A maximum-likelihood ratio test is used to evaluate the best values of these parameters and to quantify the strength of each model by contrast with isotropy. It is found that the starburst model fits the data better than the hypothesis of isotropy with a statistical significance of 4.0 sigma, the highest value of the test statistic being for energies above 39 EeV. The three alternative models are favored against isotropy with 2.7-3.2 sigma significance. The origin of the indicated deviation from isotropy is examined and prospects for more sensitive future studies are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Quantitative assessment of soil CO 2 concentration and stable carbon isotope for leakage detection at geological carbon sequestration sites

This paper presents a quantitative method to evaluate CO 2 concentration ([CO 2 ]) and stable carbon isotope ratio (δ 13 C) as indicators for leakage detection at a geological carbon sequestration site by combining use of field release tests and a numerical modeling approach. A numerical model was developed to simulate CO 2 dynamics by considering diffusion, dissolution in soil water, and soil respiration. The numerical model fits the background dynamics of [CO 2 ] (360 to 550 ppm) and δ 13 C (-16‰ to -6‰) well and reproduces fairly the overall trend observed during the CO 2 release test. The model was further applied to assess detection probability (DP) of [CO 2 ] and δ 13 C for leakage detection in terms of various factors, such as CO 2 leakage rate, background variations, δ 13 C of the CO 2 leaked, and the threshold value of signal-to-noise ratio. Modeling results suggest that δ 13 C may have a higher DP than [CO 2 ]. This study also shows that DP of δ 13 C for the IEA Weyburn project is close to 0, implying δ 13 C is inappropriate to be an indicator for CO 2 leakage at the site. Finally, the quantitative method developed can also be used to design a monitoring plan or strategy in a near-surface environment for geological carbon sequestration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Barium stars as tracers of s -process nucleosynthesis in AGB stars: II. Using machine learning techniques on 169 stars

Barium (Ba) stars are characterised by an abundance of heavy elements made by the slow neutron capture process (s-process). This peculiar observed signature is due to the mass transfer from a stellar companion, bound in a binary stellar system, to the Ba star observed today. The signature is created when the stellar companion is an asymptotic giant branch (AGB) star. We aim to analyse the abundance pattern of 169 Ba stars using machine learning techniques and the AGB final surface abundances predicted by the FRUITY and Monash stellar models. We developed machine learning algorithms that use the abundance pattern of Ba stars as input to classify the initial mass and metallicity of each Ba star’s companion star using stellar model predictions. We used two algorithms. The first exploits neural networks to recognise patterns, and the second is a nearest-neighbour algorithm that focuses on finding the AGB model that predicts the final surface abundances closest to the observed Ba star values. In the second algorithm, we included the error bars and observational uncertainties in order to find the best-fit model. The classification process was based on the abundances of Fe, Rb, Sr, Zr, Ru, Nd, Ce, Sm, and Eu. We selected these elements by systematically removing s-process elements from our AGB model abundance distributions and identifying the elements whose removal had the biggest positive effect on the classification. We excluded Nb, Y, Mo, and La. Our final classification combined the output of both algorithms to identify an initial mass and metallicity range for each Ba star companion. With our analysis tools, we identified the main properties for 166 of the 169 Ba stars in the stellar sample. The classifications based on both stellar sets of AGB final abundances show similar distributions, with an average initial mass of M = 2.23 M ⊙ and 2.34 M ⊙ and an average [Fe/H] = –0.21 and –0.11, respectively. We investigated why the removal of Nb, Y, Mo, and La improves our classification and identified 43 stars for which the exclusion had the biggest effect. We found that these stars have statistically significant and different abundances for these elements compared to the other Ba stars in our sample. We discuss the possible reasons for these differences in the abundance patterns.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Neptune's Latitudinal Variations as Viewed with ALMA

We present spatially resolved millimeter maps of Neptune between 95 and 242 GHz taken with the Atacama Large Millimeter/submillimeter Array (ALMA) in 2016–2017. The millimeter weighting functions peak between 1 and 10 bar on Neptune, lying in between the altitudes probed at visible/infrared and centimeter wavelengths. Thus, these observations provide important constraints on the atmospheric structure and dynamics of Neptune. We identify seven well-resolved latitudinal bands of discrete brightness temperature variations, on the order of 0.5–3 K in all three observed ALMA spectral bands. We model Neptune’s brightness temperature using the radiative-transfer code Radio-BEAR and compare how various H{sub 2}S, CH{sub 4}, and ortho-/para-H{sub 2} abundance profiles can fit the observed temperature variations across the disk. We find that observed variations in brightness temperature with latitude can be explained by variations in the H{sub 2}S profile that range from sub- to supersaturations at altitudes above the 10 bar pressure level, while variations in CH{sub 4} improve the quality of fit near the equator. At the south polar cap, our best-fit model has a depleted deep atmospheric abundance of H{sub 2}S from 30 to only 1.5 times the protosolar value, while simultaneously depleting the CH{sub 4} abundance. This pattern of enhancement and depletion of condensible species is consistent with a global circulation structure where enriched air rises at the midlatitudes (32°–12°S) and north of the equator (2°–20°N), and dry air descends at the poles (90°–66°S) and just south of the equator (12°S–2°N). Our analysis finds more complex structure near the equator than accounted for in previous circulation models.

79 ASTRONOMY AND ASTROPHYSICS↗