Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Functional priors”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Zeta potential measurements of SiO 2 and TiO 2 particles in anionic and cationic surfactant solutions

Zeta potential measurements can provide valuable insights into colloidal stability and may reflect changes in the electrochemical double layer that influence electrocatalysis. In this study, we used zeta potential measurements to monitor and investigate the adsorption behavior of anionic sodium dodecyl sulfate (SDS), cationic dodecyl trimethyl ammonium bromide (DTAB), and cationic cetyl trimethyl ammonium bromide (CTAB) surfactants on SiO 2 and TiO 2 particles in DI water. The zeta potentials of SiO 2 and TiO 2 particles were measured as a function of particle concentration in 1.0 mM, 5.0 mM, and 10.0 mM surfactant solutions. Results indicate that there is no significant adsorption of DS - ions on SiO 2 particles across all SDS concentrations investigated in this study. The zeta potential of TiO 2 particles shifts from negative to positive as particle concentration increases at 1.0 mM SDS, indicating a decrease in DS - ions adsorption per TiO 2 particle. Adsorption plateaus at 5.0 mM SDS for all TiO 2 concentrations studied. In contrast, DTA + and CTA + ions showed clear adsorption on SiO 2 particles, but DTA + exhibited no apparent adsorption on TiO 2 particles. Interestingly, CTA + adsorption on TiO 2 was only evident at higher surfactant concentrations. This study presents a systematic investigation of zeta potential behavior as a function of particle and surfactant concentration, highlighting adsorption saturation and electrostatic reversal. Unlike prior work that was limited to single particle concentration measurements, our approach reveals how surfactant adsorption per particle decreases with increasing surface area. Despite the structural similarity between DTAB and CTAB, this study captured the different adsorption behaviors of both surfactants on SiO 2 and TiO 2 , emphasizing the role of the overall structure of the surfactant on surfactant adsorption behavior. These insights refine the interpretation of zeta potential data and advance understanding of surfactant adsorption on oxide particles. The insight provided in this study about surfactant adsorption on SiO 2 and TiO 2 particles could have important implications for the design of more effective colloidal and catalytic systems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Plasma-Initiated Graft Polymerization of Acrylic Acid onto Fluorine-Doped Tin Oxide as a Platform for Immobilization of Water-Oxidation Catalysts

The discovery of new and versatile strategies for the immobilization of molecular water-oxidation catalysts (WOCs) is crucial for developing clean energy conversion devices [e.g., (photo)electrocatalytic cells for water splitting]. The traditional approach for surface attachment to transparent conductive oxides [e.g., fluorine doped tin oxide (FTO)] is via synthetic modification of the ligand architecture to incorporate functional groups such as carboxylic acids (-COOH) or phosphonates (-PO 3 H 2 ) prior to immobilization. However, challenges arising from desorption and the cumbersome derivatizations steps have limited the scope and applications of surface-bound WOCs. Herein, we report the successful immobilization of underivatized Ru(II)-based WOCs (Ru–Cat1 = [Ru(tpy) (bpy) (H 2 O)] 2+ (tpy = 2,2':6'2"–terpyridine and bpy = 2,2;-bipyridine) and Ru–Cat2 = [Ru(Mebimpy) (bpy) (H 2 O)] 2+ (Mebimpy = 2,6-bis(1-methylbenzimidazol-2-yl) pyridine)) and the Ru(II) polypyridyl chromophore Ru–C3 = [Ru(bpy) 3 ] 2+ onto a FTO plasma-grafted poly(acrylic acid) surface (PAA|FTO). Various characterization techniques such as attenuated total reflectance Fourier transform infrared spectroscopy, scanning electron microscopy, atomic force microscopy, and cyclic voltammetry measurements provide evidence for the plasma-induced grafted PAA|FTO film and immobilization. Surface stability and electrocatalytic properties of these new hybrid composite films upon cycling were investigated at different pH values. Immobilized Ru–Cat1 and Ru–Cat2 onto PAA|FTO displayed pH-dependent (RuIII/RuII) couples and onset potentials indicative of PCET (proton-coupled electron transfer) reactions. Based on cyclic voltammetry results and spectroscopic monitoring, the immobilized WOCs Ru–Cat1 and Ru–Cat2 exhibited a higher surface stability in neutral aqueous solutions relative to Ru–C3 upon electrochemical oxidation. Furthermore, we attribute the surface PCET and stability to the presence of a water ligand in the coordination sphere of immobilized Ru–Cat1 and Ru–Cat2 which can H-bond with negatively charged carboxylate groups of the cross-linked PAA brushes. Our findings demonstrate that the plasma-grafted polymeric network onto FTO offers a versatile platform to directly anchor unmodified homogeneous WOCs or chromophores for potential applications in solar-to-fuel energy conversion.

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

Constrained non-negative matrix factorization enabling real-time insights of in situ and high-throughput experiments

Non-negative matrix factorization (NMF) is an appealing class of methods for performing unsupervised learning on streaming spectral data, particularly in time-sensitive applications such as in situ characterization of materials. These methods seek to decompose a dataset into a small number of components and weights that can compactly represent the underlying signal while effectively reconstructing the observations with minimal error. However, canonical NMF methods have no underlying requirement that the reconstruction uses components or weights that are representative of the true physical processes. In this work, we demonstrate how constraining a subset of the NMF weights or components as rigid priors, provided as known or assumed values, can provide significant improvement in revealing true underlying phenomena. We present a PyTorch-based method for efficiently applying constrained NMF and demonstrate its application to several synthetic examples. Our implementation allows an expert researcher-in-the-loop to provide and dynamically adjust the constraints during a live experiment involving streaming spectral data. Such interactive priors allow researchers to specify known or identified independent components, as well as functional expectations about the mixing or transitions between the components. We further demonstrate the application of this method to measured synchrotron x-ray total scattering data from in situ beamline experiments. In such a context, constrained NMF can result in a more interpretive and scientifically relevant decomposition than canonical NMF or other decomposition techniques. As a result, the details of the method are provided, along with general guidance for employing constrained NMF in the extraction of critical information and insights during time-sensitive experimental applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Thiol-based pathways in the thylakoid lumen and their role in photoprotection

The long-term goal of this research is to elaborate the catalysis of thiol-disulfide transactions in the thylakoid lumen, a compartment required for transducing energy via photophosphorylation. The molecular identity of the redox players catalyzing thiol-disulfide exchanges, their relevant targets of action in vivo and how thiol-disulfide chemistry in general controls photosynthesis are the experimental questions this research aims to address. Through previous work, we documented the requirement of catalyzed disulfide formation and reduction for the biogenesis of two photosynthetic enzymes, namely Photosystem II and the cytochrome b6f complex. We established that dedicated trans-thylakoid pathways operate in these processes by delivering reducing and oxidizing power to cysteine-containing subunits of the photosynthetic complexes, which are localized in the lumen compartment. Here we focused on CCDA, CCS4 and CCS5/HCF164, the components of the disulfide-reducing pathway which maintains the heme binding cysteines of apoforms of plastid cytochrome c in the reduced form prior to covalent heme ligation. We showed that CCS4 and CCS5 operate in functionally redundant pathways for the supply of reducing power and postulate a role for CCS4 in recruiting the source of reductants to CCDA or controlling the CCDA-dependent transduction of reducing equivalents. We also demonstrate the disulfide reducing pathway operates to counter thiol oxidation of the heme-binding cysteines by disulfide forming enzyme LTO1 in the thylakoid lumen.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Functional roles of the [2Fe-2S] clusters in Synechocystis PCC 6803 Hox [NiFe]-hydrogenase reactivity with ferredoxins

The HoxEFUYH complex of Synechocystis PCC 6803 (S. 6803) consists of a HoxEFU ferredoxin:NAD(P)H oxidoreductase subcomplex and a HoxYH [NiFe]-hydrogenase subcomplex that catalyzes reversible H2 oxidation. Prior studies have suggested that the presence of HoxE is required for reactivity with ferredoxin; however, it is unknown how HoxE is functionally integrated into the electron transfer network of the HoxEFU:ferredoxin complex. Deciphering electron transfer pathways is challenged by the rich iron-sulfur cluster content of HoxEFU, which includes a [2Fe-2S] cluster in each subunit, along with multiple [4Fe-4S] clusters and a flavin cofactor. To resolve the role of HoxE, we determined the biophysical and thermodynamic properties of each [2Fe-2S] cluster in HoxEFU using steady-state and potentiometric EPR analysis in combination with square wave voltammetry (SWV). The temperature-dependence of the EPR signal for HoxE confirmed the coordination of a single [2Fe-2S] cluster that was shown by SWV to have an E m = -424 mV (versus SHE). Strikingly, when the E m of the HoxE [2Fe-2S] cluster was analyzed in HoxEFU titrations, it was shifted by >100 mV to an E m < -525 mV (versus SHE). EPR titrations of HoxEFU gave an E m value for the [2Fe-2S] cluster of HoxF, E m = -419 mV and HoxU, E m = -349 mV. These values were used to re-analyze the diaphorase kinetics in reactions performed with ferredoxins with varying E m 's. The results are formulated into a model of HoxEFU:ferredoxin reactivity and the role of HoxE in mediating electron transfer within the HoxEFU:ferredoxin complex.

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

Analysis and optimization of seismic monitoring networks with Bayesian optimal experimental design

SUMMARY Monitoring networks increasingly aim to assimilate data from a large number of diverse sensors covering many sensing modalities. Bayesian optimal experimental design (OED) seeks to identify data, sensor configurations or experiments which can optimally reduce uncertainty and hence increase the performance of a monitoring network. Information theory guides OED by formulating the choice of experiment or sensor placement as an optimization problem that maximizes the expected information gain (EIG) about quantities of interest given prior knowledge and models of expected observation data. Therefore, within the context of seismo-acoustic monitoring, we can use Bayesian OED to configure sensor networks by choosing sensor locations, types and fidelity in order to improve our ability to identify and locate seismic sources. In this work, we develop the framework necessary to use Bayesian OED to optimize a sensor network’s ability to locate seismic events from arrival time data of detected seismic phases at the regional-scale. This framework requires five elements: (i) A likelihood function that describes the distribution of detection and traveltime data from the sensor network, (ii) A prior distribution that describes a priori belief about seismic events, (iii) A Bayesian solver that uses a prior and likelihood to identify the posterior distribution of seismic events given the data, (iv) An algorithm to compute EIG about seismic events over a data set of hypothetical prior events, (v) An optimizer that finds a sensor network which maximizes EIG. Once we have developed this framework, we explore many relevant questions to monitoring such as: how to trade off sensor fidelity and earth model uncertainty; how sensor types, number and locations influence uncertainty; and how prior models and constraints influence sensor placement.

58 GEOSCIENCES↗
DOE OSTI2022

Physics makes the difference: Bayesian optimization and active learning via augmented Gaussian process

Abstract Both experimental and computational methods for the exploration of structure, functionality, and properties of materials often necessitate the search across broad parameter spaces to discover optimal experimental conditions and regions of interest in the image space or parameter space of computational models. The direct grid search of the parameter space tends to be extremely time-consuming, leading to the development of strategies balancing exploration of unknown parameter spaces and exploitation towards required performance metrics. However, classical Bayesian optimization (BO) strategies based on the Gaussian process (GP) do not readily allow for the incorporation of the known physical behaviors or past knowledge. Here we explore a hybrid optimization/exploration algorithm created by augmenting the standard GP with a structured probabilistic model of the expected system’s behavior. This approach balances the flexibility of the non-parametric GP approach with a rigid structure of physical knowledge encoded into the parametric model. The fully Bayesian treatment of the latter allows additional control over the optimization via the selection of priors for the model parameters. The method is demonstrated for a noisy version of a standard univariate test function used to evaluate optimization algorithms and further extended to physical lattice models. This methodology is expected to be universally suitable for injecting prior knowledge in the form of physical models and past data in the BO framework.

42 ENGINEERING↗
DOE OSTI2022

Qualitative trend analysis based on a mixed-integer representation

Shape constrained spline fitting is a useful method to impose prior knowledge onto flexible semi-parametric models during parameter estimation. Most typically, the function shape is imposed through order restrictions on the regression coefficients. The intended shape is considered known or selected based on heuristic rules. In this study, we present a method to estimate the optimal set of order restrictions to segment a univariate data series into episodes with distinct shapes. This is also known as the qualitative trend analysis (QTA) problem. The obtained solution uses a trade-off between lack-of-fit and model complexity. Further, our practical implementation takes inspiration from the generalized order restricted information criterion (GORIC) for inequality-constrained model selection. From this, one learns (a) that QTA can be formulated as a mixed-integer quadratic program (MIQP) and (b) that the newly proposed mixed order restricted information criterion (MORIC) enables optimal segmentation. This is illustrated through didactic case studies.

42 ENGINEERING↗
DOE OSTI2025

Causal discovery from data assisted by large language models

Knowledge-driven discovery of novel materials necessitates the development of causal models for property emergence. While in the classical physical paradigm, the causal relationships are deduced based on physical principles or via experiment, the rapid accumulation of observational data necessitates learning causal relationships between dissimilar aspects of material structure and functionalities based on observations. For this, it is essential to integrate experimental data with prior domain knowledge. Here, we demonstrate this approach by combining high-resolution scanning transmission electron microscopy data with insights derived from large language models (LLMs). By applying ChatGPT to domain-specific literature, such as arXiv papers on ferroelectrics, and combining the obtained information with data-driven causal discovery, we construct adjacency matrices for directed acyclic graphs that map the causal relationships between structural, chemical, and polarization degrees of freedom in Sm-doped BiFeO 3 . This approach enables us to hypothesize how synthesis conditions influence material properties and guides experimental validation. Furthermore, the ultimate objective of this work is to develop a unified framework that integrates LLM-driven literature analysis with data-driven discovery, facilitating the precise engineering of ferroelectric materials by establishing clear connections between synthesis conditions and their resulting material properties.

Causal inference↗
DOE OSTI2022

Scientific Computational Imaging Code (SCICO)

Scientific Computational Imaging Code (SCICO) is a Python package for solving the inverse problems that arise in scientific imaging applications. Its primary focus is providing methods for solving ill-posed inverse problems by using an appropriate prior model of the reconstruction space. SCICO includes a growing suite of operators, cost functionals, regularizers, and optimization routines that may be combined to solve a wide range of problems, and is designed so that it is easy to add new building blocks. SCICO is built on top of JAX rather than NumPy, enabling GPU/TPU acceleration, just-in-time compilation, and automatic gradient functionality, which is used to automatically compute the adjoints of linear operators. An example of how to solve a multi-channel tomography problem with SCICO is shown in Figure 1. The SCICO source code is available from GitHub, and pre-built packages are available from PyPI. It has extensive online documentation, including API documentation and usage examples, which can be run online at Google Colab and binder.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Time-dependent phenomena in correlated materials

Understanding time-dependent processes and light-matter interaction in strongly correlated materials, and the interplay between electronic, orbital, vibrational, and spin degrees of freedom, is a cornerstone of condensed matter. These mechanisms can be proven by measuring the response of the systems to time-dependent perturbations. The corresponding time scales are dictated by the way light couples to the different excitations, and how these excitations exchange energy and momentum. Our research advances our understanding of these processes, and the interpretation of different equilibrium and time-resolved spectroscopies. Our project encompasses two main themes: (i) developing and refining computational techniques to study non-equilibrium spectroscopies including non-perturbative effects and (ii) applications to non-equilibrium phenomena. We have developed a new computational approach that works directly in the time domain: by including all the degrees of freedom involved in the scattering process (e.g. electrons, photons, neutrons), we solve the time dependent problem: a faithful numerical simulation of the experiment. By measuring the energy and momentum of the outgoing particles, we can extract information about the energy and momentum absorbed by the system. Prior to our work, people attempting to model and calculate non-equilibrium spectral functions relied on a description of the scattering cross section based on a formulation in the frequency domain, a treatment that is extremely cumbersome and complex. Our technique works in and out of equilibrium and can reproduce spectra by several spectroscopic techniques, such as time-resolved photoemission, neutron scattering, Raman, X-ray spectroscopies (RIXS, Auger, XAS, XMCD), and, by not relying on analytical approximations, yields results that reveal novel overlooked transient mechanisms. These tools provide sorely needed intuition for understanding the phenomenology of strongly correlated materials and will help experimentalists in identifying signatures of relevant excitations in pump-probe experiments, such as those conducted in DOE supported facilities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

MOSAIC-CONUS: A Multimodal, Multi-Temporally Paired Dataset for Earth Sciences

Earth embeddings—vector representations of geographic locations indexed in space and time—are emerging as a unifying interface for geospatial AI. However, their quality depends not only on model design, but on how multimodal Earth observation (EO) data are spatially indexed, temporally aligned, and cross-modally associated during pretraining. We introduce MOSAIC-CONUS (Multimodal Observations with Spatially Aligned Imagery, Urban Points of Interest, In-Situ Measurements and Text Captions), a large-scale EO dataset over the contiguous United States, organized around 250,000 stratified point indices that serve as stable spatial keys across seven modalities: active radar, passive optical imagery, lidar-derived elevation, land cover, functional context, hydrometeorological measurements, and textual summaries. Unlike existing EO datasets, MOSAIC-CONUS introduces four contributions not jointly addressed in prior work: 1. an open-source, large-scale multimodal EO corpus structured around point-indexed data designed to support Earth embedding learning; 2. explicit radar-optical pairing tables spanning twelve temporal alignment regimes, formalizing cross-sensor alignment as a controllable variable for analyzing how temporal mismatch across modalities influences learned embeddings quality; 3. a benchmark suite spanning cross-modal retrieval, annual nightlights regression, and basin-held-out streamflow prediction, positioning MOSAIC-CONUS as a benchmark-ready resource for multimodal AI systems; and 4. a language-based embedding layer through co-registered textual summaries, enabling Earth embeddings to function as a queryable interface for agentic AI systems. The dataset and pairing protocols are publicly released.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Tale of Two Catchments: Causality Analysis and Isotope Systematics Reveal Mountainous Watershed Traits That Regulate the Retention and Release of Nitrogen

Abstract Mountainous watersheds are characterized by variability in functional traits, including vegetation, topography, geology, and geomorphology, which determine nitrogen (N) retention, and release. Coal Creek and East River are two contrasting catchments within the Upper Colorado River Basin that differ markedly in total nitrate (NO 3 − ) export. The East River has a diverse vegetation cover, and sinuous floodplains, and is underlain by N‐rich marine shale. At 0.21 ± 0.14 kg ha −1 yr −1 , the East River exports ∼3.5 times more NO 3 − relative to the conifer‐dominated Coal Creek (0.06 ± 0.02 kg ha −1 yr −1 ). While this can partly be explained by the larger size of the East River, the distinct watershed traits of these two catchments imply different mechanisms controlling the aggregate N‐export signal. A causality analysis shows physical and biogenic processes were critical in determining NO 3 − export from the East River catchment. Stable isotope ratios of NO 3 − (δ 15 N NO3 and δ 18 O NO3 ) show the East River catchment is a strong hotspot for biogeochemical processing of NO 3 − at the hillslope soil‐saprolite. By contrast, the conifer‐dominated Coal Creek retained nearly all atmospherically deposited NO 3 − , and its export was controlled by catchment hydrological traits (i.e., snowmelt periods and water table depth). The conservative N‐cycle within Coal Creek is likely due to the abundance of conifer trees, and smaller riparian regions, retaining more NO 3 − overall and reduced processing prior to export. This study highlights the value of integrating isotope systematics to link watershed functional traits to mechanisms of watershed element retention and release.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Analysis and Optimization of Seismo-Acoustic Monitoring Networks with Bayesian Optimal Experimental Design

The Bayesian optimal experimental design (OED) problem seeks to identify data, sensor configurations, or experiments which can optimally reduce uncertainty. The goal of OED is to find an experiment that maximizes the expected information gain (EIG) about quantities of interest given prior knowledge about expected data. Therefore, within the context of seismic monitoring, we can use Bayesian OED to configure sensor networks by choosing sensor locations, types, and fidelity in order to improve our ability to identify and locate seismic sources. In this work, we develop the framework necessary to use Bayesian OED to optimize the ability to locate seismic events from arrival time data of detected seismic phases. In order to do utilize Bayesian OED we must develop four elements:1. A likelihood function that describes the uncertainty of detection and travel times; 2. A Bayesian solver that takes a prior and likelihood to identify the posterior; 3. An algorithm to compute EIG; and, 4. An optimizer that finds a sensor network which maximizes EIG. Once we have developed this framework, we can explore many relevant questions to monitoring such as: how and what multiphenomenology data can be used to optimally reduce uncertainty, how to trade off sensor fidelity and earth model uncertainty, and how sensor types, number, and locations influence uncertainty

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗
DOE OSTI2025

Analysis of the Trusted Inertial Terrain-Aided Navigation Measurement Function

The trusted inertial terrain-aided navigation (TITAN) algorithm leverages an airborne vertical synthetic aperture radar to measure the range to the closest ground points along several prescribed iso-Doppler contours. These TITAN minimum-range, prescribed-Doppler measurements are the result of a constrained nonlinear optimization problem whose optimization function and constraints both depend on the radar position and velocity. Owing to the complexity of this measurement definition, analysis of the TITAN algorithm is lacking in prior work. This publication offers such an analysis, making the following three contributions: (1) an analytical solution to the TITAN constrained optimization measurement problem, (2) a derivation of the TITAN measurement function Jacobian, and (3) a derivation of the Cramér-Rao lower bound on the estimated position and velocity error covariance. These three contributions are verified via Monte Carlo simulations over synthetic terrain, which further reveal two remarkable properties of the TITAN algorithm: (1) the along-track positioning errors tend to be smaller than the cross-track positioning errors, and (2) the cross-track positioning errors are independent of the terrain roughness.

TITAN↗
DOE OSTI2022

How to Obtain the Redshift Distribution from Probabilistic Redshift Estimates

Abstract A reliable estimate of the redshift distribution n ( z ) is crucial for using weak gravitational lensing and large-scale structures of galaxy catalogs to study cosmology. Spectroscopic redshifts for the dim and numerous galaxies of next-generation weak-lensing surveys are expected to be unavailable, making photometric redshift (photo- z ) probability density functions (PDFs) the next best alternative for comprehensively encapsulating the nontrivial systematics affecting photo- z point estimation. The established stacked estimator of n ( z ) avoids reducing photo- z PDFs to point estimates but yields a systematically biased estimate of n ( z ) that worsens with a decreasing signal-to-noise ratio, the very regime where photo- z PDFs are most necessary. We introduce Cosmological Hierarchical Inference with Probabilistic Photometric Redshifts ( CHIPPR ), a statistically rigorous probabilistic graphical model of redshift-dependent photometry that correctly propagates the redshift uncertainty information beyond the best-fit estimator of n ( z ) produced by traditional procedures and is provably the only self-consistent way to recover n ( z ) from photo- z PDFs. We present the chippr prototype code, noting that the mathematically justifiable approach incurs computational cost. The CHIPPR approach is applicable to any one-point statistic of any random variable, provided the prior probability density used to produce the posteriors is explicitly known; if the prior is implicit, as may be the case for popular photo- z techniques, then the resulting posterior PDFs cannot be used for scientific inference. We therefore recommend that the photo- z community focus on developing methodologies that enable the recovery of photo- z likelihoods with support over all redshifts, either directly or via a known prior probability density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

NREM sleep as a novel protective cognitive reserve factor in the face of Alzheimer's disease pathology

Alzheimer’s disease (AD) pathology impairs cognitive function. Yet some individuals with high amounts of AD pathology suffer marked memory impairment, while others with the same degree of pathology burden show little impairment. Why is this? One proposed explanation is cognitive reserve i.e., factors that confer resilience against, or compensation for the effects of AD pathology. Deep NREM slow wave sleep (SWS) is recognized to enhance functions of learning and memory in healthy older adults. However, that the quality of NREM SWS (NREM slow wave activity, SWA) represents a novel cognitive reserve factor in older adults with AD pathology, thereby providing compensation against memory dysfunction otherwise caused by high AD pathology burden, remains unknown. Here, we tested this hypothesis in cognitively normal older adults (N = 62) by combining 11 C-PiB (Pittsburgh compound B) positron emission tomography (PET) scanning for the quantification of β-amyloid (Aβ) with sleep electroencephalography (EEG) recordings to quantify NREM SWA and a hippocampal-dependent face-name learning task. We demonstrated that NREM SWA significantly moderates the effect of Aβ status on memory function. Specifically, NREM SWA selectively supported superior memory function in individuals suffering high Aβ burden, i.e., those most in need of cognitive reserve (B = 2.694, p = 0.019). In contrast, those without significant Aβ pathological burden, and thus without the same need for cognitive reserve, did not similarly benefit from the presence of NREM SWA (B = -0.115, p = 0.876). This interaction between NREM SWA and Aβ status predicting memory function was significant after correcting for age, sex, Body Mass Index, gray matter atrophy, and previously identified cognitive reserve factors, such as education and physical activity (p = 0.042). These findings indicate that NREM SWA is a novel cognitive reserve factor providing resilience against the memory impairment otherwise caused by high AD pathology burden. Furthermore, this cognitive reserve function of NREM SWA remained significant when accounting both for covariates, and factors previously linked to resilience, suggesting that sleep might be an independent cognitive reserve resource. Beyond such mechanistic insights are potential therapeutic implications. Unlike many other cognitive reserve factors (e.g., years of education, prior job complexity), sleep is a modifiable factor. As such, it represents an intervention possibility that may aid the preservation of cognitive function in the face of AD pathology, both present moment and longitudinally.

60 APPLIED LIFE SCIENCES↗