Search NASA⌕ Search

SEARCH · Search NASA

Results for “Population modeling”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Distinct critical behaviors from the same state in quantum spin and population dynamics perspectives

There is a deep connection between the ground states of transverse-field spin systems and the late-time distributions of evolving viral populations—within simple models, both are obtained from the principal eigenvector of the same matrix. However, that vector is the wave-function amplitude in the quantum spin model, whereas it is the probability itself in the population model. We show that this seemingly minor difference has significant consequences: Phase transitions that are discontinuous in the spin system become continuous when viewed through the population perspective, and transitions that are continuous become governed by new critical exponents. We introduce a more general class of models that encompasses both cases and that can be solved exactly in a mean-field limit. Numerical results are also presented for a number of one-dimensional chains with power-law interactions. We see that well-worn spin models of quantum statistical mechanics can contain unexpected new physics and insights when treated as population-dynamical models and beyond, motivating further studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Gas-Liquid Flow Modeling for Renewable Fuels Production

Aerobic/anaerobic and gas fermentation pathways have emerged as promising new technologies for the generation of renewable fuels/chemicals from biomass derived sugars, and mixtures of greenhouse/energy rich gas streams (CO2/CH4/H2/CO) via microbial action. Example pathways include sugars-to-ethanol conversion, biomethanation (CO2/H2 to CH4), biogas upgrading, CO fermentation and wet-waste conversion. Gas and liquid phase transport, mass-transfer, and mixing physics at large length scales can significantly affect microbial conversion rates, particularly when the microbial reaction requires a narrow set of conditions. These phenomena are difficult to study in small-scale bench-top reactors that are typically well-mixed. Predictive computational fluid dynamics (CFD) based simulations can therefore aid in the scale-up, design and optimization of these reactors. This work presents multiphase Euler-Euler CFD simulations of at-scale (~500 m3) bioreactors. Our mathematical model treats the gas and liquid as interpenetrating phases. This approach reduces the computational complexity of tracking individual gas bubbles that are several orders of magnitude smaller than reactor dimensions. We solve the Reynolds averaged Navier-Stokes (RANS) multiphase equations that account for phase and chemical species transport, interphase mass and momentum transfer and uses a phenomenological model for gas uptake by microbes. We use a customized solver derived from open-source CFD toolbox, OpenFOAM [1], to perform these simulations, which has been validated against small-scale reactors in our previous work [2]. There is currently a knowledge-gap regarding bubble-size distributions when using gas mixtures with vastly different properties, which can have a significant impact overall mass-transfer. For example, hydrogen bubbles are more buoyant compared to other relatively heavier gases (CO2/CH4/CO), resulting in a large distribution of residence times and bubble sizes. This work therefore develops a deeper understanding of bubble dynamics and interphase mass transfer in such heterogenous gas mixtures through well-resolved computational models. We use a population balance model (PBM) for bubble-size-distribution modeling that is validated against small-scale experiments in our solver with an uncertainty quantification study for bubble coalescence and break-up model parameters. Results pertaining to multiple simulations of gas-fermentation reactors are presented where gas mixtures with varying compositions of CO2/CH4/CO/H2 are imposed at the sparger boundaries. The spatio-temporal variations in bubble-size distribution and mass transfer coefficient are analyzed for varying superficial velocities and gas-compositions for varying sizes of bubble-column and airlift reactors. This work will also examine the performance of different reactor designs, viz. bubble column reactor, airlift reactor with an internal draft tube, and a stirred-tank reactor with Rushton impellers. Reactor mass-transfer coefficient, gas hold-up, and dissolved gas distribution are critically analyzed among reactors, and sensitivity studies pertaining to gas flow rates and reactor geometry will be presented. [1] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998. [2] Rahimi, M., Sitaraman, H., Humbird, D. and Stickel, J., Computational fluid dynamics study of full-scale aerobic bioreactors: Evaluation of gas-liquid mass transfer, oxygen uptake, and dynamic oxygen distribution, Chemical Engineering Research and Design, 139: 283-295.

BIOMASS FUELS↗
DOE OSTI2021

A Binomial Stochastic Framework for Efficiently Modeling Discrete Statistics of Convective Populations

Abstract Understanding the coupling between convective clouds and the general circulation, as well as addressing the gray zone problem in convective parameterization, requires insight into the genesis and maintenance of spatial patterns in cumulus cloud populations. In this study, a simple toy model for recreating populations of interacting convective objects as distributed over a two‐dimensional Eulerian grid is formulated to this purpose. Key elements at the foundation of the model include i) a fully discrete formulation for capturing discrete behavior in convective properties at small population sample sizes, ii) object age‐dependence for representing life‐cycle effects, and iii) a prognostic number budget allowing for object interactions and co‐existence of multiple species. A primary goal is to optimize the computational efficiency of this system. To this purpose the object birth rate is represented stochastically through a spatially aware Bernoulli process. The same binomial stochastic operator is applied to horizontal advection of objects, conserving discreteness in object number. The applicability to atmospheric convection as well as behavior implied by the formulation is assessed. Various simple applications of the BiOMi model (Binomial Objects on Microgrids) are explored, suggesting that important convective behavior can be captured at low computational cost. This includes i) subsampling effects and associated powerlaw scaling in the convective gray zone, ii) stochastic predator‐prey behavior, iii) the downscale turbulent energy cascade, and iv) simple forms of spatial organization and convective memory. Consequences and opportunities for convective parameterization in next‐generation weather and climate models are discussed.

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

A Process-based Model with Temperature, Water, and Lab-derived Data Improves Predictions of Daily Culex pipiens/restuans Mosquito Density

While the number of human cases of mosquito-borne diseases has increased in North America in the last decade, accurate modeling of mosquito population density has remained a challenge. Longitudinal mosquito trap data over the many years needed for model calibration, and validation is relatively rare. In particular, capturing the relative changes in mosquito abundance across seasons is necessary for predicting the risk of disease spread as it varies from year to year. We developed a discrete, semi-stochastic, mechanistic process-based mosquito population model that captures life-cycle egg, larva, pupa, adult stages, and diapause for Culex pipiens (Diptera, Culicidae) and Culex restuans (Diptera, Culicidae) mosquito populations. This model combines known models for development and survival into a fully connected age-structured model that can reproduce mosquito population dynamics. Mosquito development through these stages is a function of time, temperature, daylight hours, and aquatic habitat availability. The time-dependent parameters are informed by both laboratory studies and mosquito trap data from the Greater Toronto Area. The model incorporates city-wide water-body gauge and precipitation data as a proxy for aquatic habitat. This approach accounts for the nonlinear interaction of temperature and aquatic habitat variability on the mosquito life stages. We demonstrate that the full model predicts the yearly variations in mosquito populations better than a statistical model using the same data sources. This improvement in modeling mosquito abundance can help guide interventions for reducing mosquito abundance in mitigating mosquito-borne diseases like West Nile virus.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Constraining the sources of ultra-high-energy cosmic rays across and above the ankle with the spectrum and composition data measured at the Pierre Auger Observatory

In this work we present the interpretation of the energy spectrum and mass composition data as measured by the Pierre Auger Collaboration above 6 × 10$^{17}$ eV. We use an astrophysical model with two extragalactic source populations to model the hardening of the cosmic-ray flux at around 5 × 10$^{18}$ eV (the so-called “ankle” feature) as a transition between these two components. We find our data to be well reproduced if sources above the ankle emit a mixed composition with a hard spectrum and a low rigidity cutoff. The component below the ankle is required to have a very soft spectrum and a mix of protons and intermediate-mass nuclei. The origin of this intermediate-mass component is not well constrained and it could originate from either Galactic or extragalactic sources.To the aim of evaluating our capability to constrain astrophysical models, we discuss the impact on the fit results of the main experimental systematic uncertainties and of the assumptions about quantities affecting the air shower development as well as the propagation and redshift distribution of injected ultra-high-energy cosmic rays (UHECRs).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Improved biomass feedstock materials handling and feeding engineering data sets, design methods, and modeling/simulation tools

Forest Concepts led a project to develop a set of tools which enable feedstock handing equipment designers to better understand and model the flowability of bulk particulate biomass materials. The work products included improved flowability mathematical models, data sets used to populate the models, and two new laboratory devices – a biomass-scale true cubical triaxial tester and a biomass-scale gas pycnometer. This work was funded in part by the US Department of Energy under contact DE-EE0008254.

09 BIOMASS FUELS↗
DOE OSTI2024

A strong blend in the morning: studying the circumgalactic medium before cosmic noon with strong, blended Ly α forest systems

ABSTRACT We study of the properties of a new class of circumgalactic medium absorbers identified in the Ly α forest: ‘Strong, Blended Lyman-α’ (or SBLA) absorption systems. We study SBLAs at 2.4 < z < 3.1 in SDSS-IV/eBOSS spectra by their strong extended Ly α absorption complexes covering 138 $\, \, {\rm km}\, {\rm s}^{-1}$ with an integrated $\log (N_{\rm H\, {\small I}}/\mathrm{cm}^{-2}) =16.04$$\substack{+0.05 \\ -0.06}$ and Doppler parameter b = 18.1$\substack{+0.7 \\ -0.4}$$\, \, {\rm km}\, {\rm s}^{-1}$. Clustering with the Ly α forest provides a large-scale structure bias of b = 2.34 ± 0.06 and halo mass estimate of $M_h \approx 10^{12}\, h^{-1}\, {\rm M_{\odot }}$ for our SBLA sample. We measure the ensemble mean column densities of 22 metal features in the SBLA composite spectrum and find that no single-population multiphase model for them is viable. We therefore explore the underlying SBLA population by forward modelling the SBLA absorption distribution. Based on covariance measurements and favoured populations we find that ≈25 per cent of our SBLAs have stronger metals. Using silicon only we find that our strong metal SBLAs trace gas with a log (nH/cm−3) > −2.40 for T = 103.5 K and show gas clumping on <210 parsec scales. We fit multiphase models to this strong subpopulation and find a low ionization phase with nH = 1 cm−3, T = 103.5 K, and [X/H] = 0.8, an intermediate ionization phase with log (nH/cm−3) = −3.05, T = 103.5 K and [X/H] = −0.8, and a poorly constrained higher ionization phase. We find that the low ionization phase favours cold, dense super-solar metallicity gas with a clumping scale of just 0.009 parsecs.

Morrison, Sean (ORCID:0000000267702627)↗
DOE OSTI2024

LandScan Global 2023: Silver Edition

For a quarter of a century, the LandScan Global (LSG) project has annually released a global, high-resolution gridded population dataset representing the ambient or unwarned population at a 30 arcsecond resolution. LSG supports a range of applications such as emergency management, disaster response, and human health and security for understanding populations at risk. The 2023 release of LSG, the LandScan Silver Edition, represents a major methodological leap forward while also leveraging previous knowledge—the previous year was the baseline for the current annual update carrying forward valuable knowledge of the built environment for the past quarter century—to train the machine learning models. Compared with annual releases over the past 24years, multiple advancements were made to different aspects of the methodology to achieve reproducibility, transparency, and consistent global propagation of solutions to modeling or population distribution issues identified during the review process. These novel changes include incorporation of the latest available geospatial inputs across the globe, machine learning models instead of manual modifications, population feature importance analysis, open-source solutions vs. proprietary software, generation of multiple global versions, analytic validations, and human-in-the-loop revisions to produce the final version. Additionally, algorithms—such as anomaly detection—were introduced to quickly identify areas of focus to develop a new and robust systematic review. Significant changes in modeled population distributions were observed between the 2022 and 2023 releases, largely attributable to improvements in data and methods and discussed thoroughly within this report. In summation, the LandScan Silver Edition leverages the best of the past quarter century of LSG legacy knowledge and continues a tradition of applying cutting-edge enhancements to serve as a new benchmark for accurate, actionable gridded population data

Lebakula, Viswadeep↗
DOE OSTI2024

Carbon- and Oxygen-rich stars in MaStar: identification and classification

ABSTRACT Carbon- and Oxygen-rich stars populating the Thermally Pulsing Asymptotic Giant Branch (TP-AGB) phase of stellar evolution are relevant contributors to the spectra of ∼1 Gyr old populations. Atmosphere models for these types are uncertain, due to complex molecules and mass-loss effects. Empirical spectra are then crucial, but samples are small due to the short (∼3 Myr) TP-AGB lifetime. Here we exploit the vastness of the MaNGA Stellar library MaStar (∼60 000 spectra) to identify C and O-rich type stars. We define an optical colour selection with cuts of (g − r) > 2 and (g − i) < 1.55(g − r) − 0.07, calibrated with known C- and O-rich spectra. This identifies C- and O-rich stars along clean, separated sequences. An analogue selection is found in V, R, I bands. Our equation identifies C-rich and O-rich spectra with predictive performance metric F1-scores of 0.72 and 0.74 (over 1), respectively. We finally identify 41 C- and 87 O-rich type AGB stars in MaStar, 5 and 49 of which do not have a SIMBAD counterpart. We also detect a sample of non-AGB, dwarf C-stars. We further design a fitting procedure to classify the spectra into broad spectral types, by using as fitting templates empirical C- and O-rich spectra. We find remarkably good fits for the majority of candidates and categorize them into C- and O-rich bins following existing classifications, which correlate to effective temperature. Our selection models can be applied to large photometric surveys (e.g. Euclid, Rubin). The classified spectra will facilitate future evolutionary population synthesis models.

Hill, Lewis (ORCID:0000000180528080)↗
DOE OSTI2023

Simulation of radioactive plume transport in the atmosphere including dynamics of particle aggregation and breakup

Accurate prediction of the atmospheric transport of debris particles relies heavily on our knowledge of the size distribution of the particles within a debris cloud. Assuming a fixed particle size during simulations is not always viable since the size distribution of the debris can change during transport. Various microphysical processes, such as aggregation and breakup, influence debris particles and dictate any changes to the size distribution. To track those changes that can occur, a population balance model can be adopted and instituted within a model framework. Nonetheless, many of the models that simulate the transport of radioactive debris following a device-driven fission incident have historically neglected to consider these processes. As such, this work describes our effort to develop a modeling framework capable of simulating the transport and deposition of a radioactive plume generated from a fission incident with a dynamic population balance including particle aggregation and breakup. Herein the impact of aggregation and breakup, individually and collectively, on the particle size distribution is explored using the developed framework. When simulating aggregation, for example, six mechanisms, including Brownian coagulation, the convective enhancement to Brownian coagulation, van der Waals-viscous force correction for Brownian coagulation, gravitational collection, turbulent inertial motion, and turbulent shear, are considered. Brownian coagulation and its corrections have, as one would expect, a large impact on relatively small aggregates. Aggregates with a diameter that is less than or equal to 1.0 μm, for instance, comprise 50.6 vol % of all aggregates in the absence of aggregation and 31.2 vol % when Brownian coagulation and its corrections are considered. Gravitational collection and, to a much lesser extent, turbulent shear and turbulent inertial motion are, conversely, of great importance to relatively large aggregates (i.e., diameter greater than 3.0 μm). Additionally, the individual effects of atmospheric and particle parameters, such as wind speed and particle density, are examined. Of the parameters examined, turbulent energy dissipation and aggregate fractal dimension (i.e., aggregate shape with lower values representing more irregular particles) were of substantial importance since both terms directly impact aggregate stability and, by extension, the breakup rate. Large-scale transport and deposition simulations in a dry atmosphere are also presented and discussed as a proof of concept.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2020

National population mapping from sparse survey data: A hierarchical Bayesian modeling framework to account for uncertainty

Population estimates are critical for government services, development projects, and public health campaigns. Such data are typically obtained through a national population and housing census. However, population estimates can quickly become inaccurate in localized areas, particularly where migration or displacement has occurred. Some conflict-affected and resource-poor countries have not conducted a census in over 10 y. We developed a hierarchical Bayesian model to estimate population numbers in small areas based on enumeration data from sample areas and nationwide information about administrative boundaries, building locations, settlement types, and other factors related to population density. We demonstrated this model by estimating population sizes in every 10- m grid cell in Nigeria with national coverage. These gridded population estimates and areal population totals derived from them are accompanied by estimates of uncertainty based on Bayesian posterior probabilities. The model had an overall error rate of 67 people per hectare (mean of absolute residuals) or 43% (using scaled residuals) for predictions in out-of-sample survey areas (approximately 3 ha each), with increased precision expected for aggregated population totals in larger areas. This statistical approach represents a significant step toward estimating populations at high resolution with national coverage in the absence of a complete and recent census, while also providing reliable estimates of uncertainty to support informed decision making.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Developing an Evolutionary Baseline Model for Humans: Jointly Inferring Purifying Selection with Population History

Building evolutionarily appropriate baseline models for natural populations is not only important for answering fundamental questions in population genetics—including quantifying the relative contributions of adaptive versus nonadaptive processes—but also essential for identifying candidate loci experiencing relatively rare and episodic forms of selection (e.g., positive or balancing selection). Here, a baseline model was developed for a human population of West African ancestry, the Yoruba, comprising processes constantly operating on the genome (i.e., purifying and background selection, population size changes, recombination rate heterogeneity, and gene conversion). Specifically, to perform joint inference of selective effects with demography, an approximate Bayesian approach was employed that utilizes the decay of background selection effects around functional elements, taking into account genomic architecture. This approach inferred a recent 6-fold population growth together with a distribution of fitness effects that is skewed towards effectively neutral mutations. Importantly, these results further suggest that, although strong and/or frequent recurrent positive selection is inconsistent with observed data, weak to moderate positive selection is consistent but unidentifiable if rare.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Surrogate modelling the Baryonic Universe II: On forward modelling the colours of individual and populations of galaxies

ABSTRACT Among the properties shaping the light of a galaxy, the star formation history (SFH) is one of the most challenging to model due to the variety of correlated physical processes regulating star formation. In this work, we leverage the stellar population synthesis model fsps, together with SFHs predicted by the hydrodynamical simulation IllustrisTNG and the empirical model universemachine, to study the impact of star formation variability on galaxy colours. We start by introducing a model-independent metric to quantify the burstiness of a galaxy formation model, and we use this metric to demonstrate that universemachine predicts SFHs with more burstiness relative to IllustrisTNG. Using this metric and principal component analysis, we construct families of SFH models with adjustable variability, and we show that the precision of broad-band optical and near-infrared colours degrades as the level of unresolved short-term variability increases. We use the same technique to demonstrate that variability in metallicity and dust attenuation presents a practically negligible impact on colours relative to star formation variability. We additionally provide a model-independent fitting function capturing how the level of unresolved star formation variability translates into imprecision in predictions for galaxy colours; our fitting function can be used to determine the minimal SFH model that reproduces colours with some target precision. Finally, we show that modelling the colours of individual galaxies with per cent-level precision demands resorting to complex SFH models, while producing precise colours for galaxy populations can be achieved using models with just a few degrees of freedom.

Chaves-Montero, Jonás (ORCID:0000000295534261)↗
DOE OSTI2025

Bayesian calibration of bubble size dynamics applied to CO 2 gas fermenters

To accelerate the scale-up of gaseous CO 2 fermentation reactors, computational models need to predict gas-to-liquid mass transfer which requires capturing the bubble size dynamics, i.e. bubble breakup and coalescence. However, the applicability of existing models beyond air–water mixtures remains to be established. Here, an inverse modeling approach, accelerated with a neural network surrogate, calibrates the breakup and coalescence closure models, that are used in class methods for population balance modeling (PBM). The calibration is performed based on experimental results obtained in a CO 2 -air–water-coflowing bubble column reactor. Bayesian inference is used to account for noise in the experimental dataset and bias in the simulation results. To accurately capture gas holdup and interphase mass transfer, the results show that the breakage rate needs to be increased by one order of magnitude. In conclusion, the inferred model parameters are then used on a separate configuration and shown to also improve bubble size distribution predictions.

09 BIOMASS FUELS↗
DOE OSTI2022

Behavioral and Population Data-Driven Distribution System Load Modeling

Distribution system residential load modeling and analysis for different geographic areas within a utility or an independent system operator territory are critical for enabling small-scale, aggregated distributed energy resources to participate in grid services under Federal Energy Regulatory Commission Order No. 2222 [1]. In this study, we develop a methodology of modeling residential load profiles in different geographic areas with a focus on human behavior impact. First, we construct a behavior-based load profile model leveraging state-of-the-art appliance models. We simulate human activity and occupancy using Markov chain Monte Carlo methods calibrated with the American Time Use Survey data set. Second, we link our model with cleaned Current Population Survey data from the U.S. Census Bureau. Finally, we populate two sets of 500 households using California and Texas census data, respectively, to perform an initial analysis of the load in different geographic areas with various group features (e.g., different income levels). To distinguish the effect of population behavior differences on aggregated load, we simulate load profiles for both sets assuming fixed physical household parameters and weather data. Analysis shows that average daily load profiles vary significantly by income and income dependency varies by locality.

American Time Use Survey↗
DOE OSTI2020

Potentially adaptive SARS-CoV-2 mutations discovered with novel spatiotemporal and explainable AI models

Abstract Background A mechanistic understanding of the spread of SARS-CoV-2 and diligent tracking of ongoing mutagenesis are of key importance to plan robust strategies for confining its transmission. Large numbers of available sequences and their dates of transmission provide an unprecedented opportunity to analyze evolutionary adaptation in novel ways. Addition of high-resolution structural information can reveal the functional basis of these processes at the molecular level. Integrated systems biology-directed analyses of these data layers afford valuable insights to build a global understanding of the COVID-19 pandemic. Results Here we identify globally distributed haplotypes from 15,789 SARS-CoV-2 genomes and model their success based on their duration, dispersal, and frequency in the host population. Our models identify mutations that are likely compensatory adaptive changes that allowed for rapid expansion of the virus. Functional predictions from structural analyses indicate that, contrary to previous reports, the Asp 614 Gly mutation in the spike glycoprotein (S) likely reduced transmission and the subsequent Pro 323 Leu mutation in the RNA-dependent RNA polymerase led to the precipitous spread of the virus. Our model also suggests that two mutations in the nsp13 helicase allowed for the adaptation of the virus to the Pacific Northwest of the USA. Finally, our explainable artificial intelligence algorithm identified a mutational hotspot in the sequence of S that also displays a signature of positive selection and may have implications for tissue or cell-specific expression of the virus. Conclusions These results provide valuable insights for the development of drugs and surveillance strategies to combat the current and future pandemics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Chemical Cartography with APOGEE: Mapping Disk Populations with a 2-process Model and Residual Abundances

We apply a novel statistical analysis to measurements of 16 elemental abundances in 34,410 Milky Way disk stars from the final data release (DR17) of APOGEE-2. Building on recent work, we fit median abundance ratio trends [X/Mg] versus [Mg/H] with a 2-process model, which decomposes abundance patterns into a "prompt" component tracing core-collapse supernovae and a "delayed" component tracing Type Ia supernovae. For each sample star, we fit the amplitudes of these two components, then compute the residuals Δ[X/H] from this two-parameter fit. The rms residuals range from ~0.01–0.03 dex for the most precisely measured APOGEE abundances to ~0.1 dex for Na, V, and Ce. The correlations of residuals reveal a complex underlying structure, including a correlated element group comprised of Ca, Na, Al, K, Cr, and Ce and a separate group comprised of Ni, V, Mn, and Co. Selecting stars poorly fit by the 2-process model reveals a rich variety of physical outliers and sometimes subtle measurement errors. Residual abundances allow for the comparison of populations controlled for differences in metallicity and [α/Fe]. Relative to the main disk (R = 3–13 kpc), we find nearly identical abundance patterns in the outer disk (R = 15–17 kpc), 0.05–0.2 dex depressions of multiple elements in LMC and Gaia Sausage/Enceladus stars, and wild deviations (0.4–1 dex) of multiple elements in ω Cen. The residual abundance analysis opens new opportunities for discovering chemically distinctive stars and stellar populations, for empirically constraining nucleosynthetic yields, and for testing chemical evolution models that include stochasticity in the production and redistribution of elements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

DAmodel: hierarchical Bayesian modelling of DA white dwarfs for spectrophotometric calibration

We use hierarchical Bayesian modelling to calibrate a network of 32 all-sky faint DA white dwarf (DA WD) spectrophotometric standards (⁠16.5 < V , 19.5⁠) alongside three CALSPEC standards, from 912 Å to 32 μm. The framework is the first of its kind to jointly infer photometric zero points and WD parameters (surface gravity log g⁠, effective temperature T eff ⁠, extinction A V ⁠, dust relation parameter R V ) by simultaneously modelling both photometric and spectroscopic data. We model panchromatic Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) UVIS and IR photometry, HST/STIS UV spectroscopy, and ground-based optical spectroscopy to sub-per cent precision. Photometric residuals for the sample are the lowest yet yielding < 0.004 mag RMS on average from the UV to the NIR, achieved by jointly inferring time-dependent changes in system sensitivity and WFC3/IR count-rate nonlinearity. Our GPU-accelerated implementation enables efficient sampling via Hamiltonian Monte Carlo, critical for exploring the high-dimensional posterior space. The hierarchical nature of the model enables population analysis of intrinsic WD and dust parameters. Inferred spectral energy distributions from this model will be essential for calibrating the James Webb Space Telescope as well as next-generation surveys, including Vera Rubin Observatory’s Legacy Survey of Space and Time and the Nancy Grace Roman Space Telescope.

methods: statistical↗