Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Inference mechanisms”

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

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

At least 1,063 records · Page 59

MONTI: A Multi-Omics Non-negative Tensor Decomposition Framework for Gene-Level Integrative Analysis

Multi-omics data is frequently measured to enrich the comprehension of biological mechanisms underlying certain phenotypes. However, due to the complex relations and high dimension of multi-omics data, it is difficult to associate omics features to certain biological traits of interest. For example, the clinically valuable breast cancer subtypes are well-defined at the molecular level, but are poorly classified using gene expression data. Here, we propose a multi-omics analysis method called MONTI (Multi-Omics Non-negative Tensor decomposition for Integrative analysis), which goal is to select multi-omics features that are able to represent trait specific characteristics. Here, we demonstrate the strength of multi-omics integrated analysis in terms of cancer subtyping. The multi-omics data are first integrated in a biologically meaningful manner to form a three dimensional tensor, which is then decomposed using a non-negative tensor decomposition method. From the result, MONTI selects highly informative subtype specific multi-omics features. MONTI was applied to three case studies of 597 breast cancer, 314 colon cancer, and 305 stomach cancer cohorts. For all the case studies, we found that the subtype classification accuracy significantly improved when utilizing all available multi-omics data. MONTI was able to detect subtype specific gene sets that showed to be strongly regulated by certain omics, from which correlation between omics types could be inferred. Furthermore, various clinical attributes of nine cancer types were analyzed using MONTI, which showed that some clinical attributes could be well explained using multi-omics data. We demonstrated that integrating multi-omics data in a gene centric manner improves detecting cancer subtype specific features and other clinical features, which may be used to further understand the molecular characteristics of interest. The software and data used in this study are available at: https://github.com/inukj/MONTI.

59 BASIC BIOLOGICAL SCIENCES↗

Wildfire Segmentation From Remotely Sensed Data Using Quantum-Compatible Conditional Vector Quantized-Variational Autoencoders

Wildfires represent a critical environmental hazard with multifaceted implications for ecosystems, communities, and public health [1]. The escalating frequency and intensity of wildfires globally have intensified the urgency for robust segmentation methodologies to facilitate effective mitigation, response, and recovery strategies [2]. Accurate wildfire segmentation is pivotal for delineating fire boundaries, assessing progression patterns, and prioritizing resource allocation during emergency scenarios. Furthermore, precise segmentation enables stakeholders, including policymakers, environmental scientists, and emergency responders, to formulate evidence-based strategies, thereby minimizing socio-economic disruptions and ecological degradation. Consequently, advancing wildfire segmentation techniques through innovative technological interventions remains a paramount research imperative. Although foundational in wildfire segmentation, traditional deterministic models exhibit inherent limitations that compromise their efficacy in dynamic and uncertain environments. These models often operate on rigid algorithms prioritizing deterministic classifications, thereby overlooking the inherent complexities and uncertainties associated with wildfire behavior and satellite data variability. Such deterministic frameworks tend to produce oversimplified representations that fail to capture the intricate nuances of evolving fire dynamics, spatial heterogeneity, and environmental interactions [1]. Consequently, the deterministic approach’s propensity for uncertainty collapsing [1, 3] hampers the accuracy, reliability, and applicability of segmentation outcomes in real-world scenarios. Contrastingly, stochastic models offer a more nuanced and adaptable framework for wildfire segmentation. By integrating probabilistic elements into the modeling paradigm, stochastic approaches, particularly probabilistic approaches such as variational auto encoders (VAEs) [4], facilitate comprehensive uncertainty assessment, enabling researchers to quantify and incorporate uncertainties into segmentation outcomes effectively. This probabilistic nature empowers stochastic models to encapsulate variability, account for data inconsistencies, and adapt to evolving environmental conditions, enhancing segmentation accuracy, reliability, and robustness. Embracing stochastic methodologies thus catalyzes advancements in wildfire science by fostering a more holistic, adaptive, and resilient segmentation framework. Despite VAEs demonstrating significant promise in various applications, they come with inherent limitations that have garnered attention within the machine learning community. One of the primary drawbacks lies in their reliance on static priors, which essentially assume a fixed distribution for latent variables, thereby limiting the model’s flexibility to capture complex data structures effectively [5]. This static nature leads to suboptimal representations, especially when dealing with complex and high-dimensional data. Additionally, VAEs often struggle with generating sharp and realistic samples, a phenomenon commonly referred to as mode collapse [5, 7, 6]. Furthermore, the optimization process in VAEs, which involves balancing the reconstruction loss and the regularization term, can sometimes be challenging to fine-tune [7]. In recent efforts to address these shortcomings, alternative approaches like Vector Quantized Variational Auto encoders(VQ-VAEs) [7], address the challenges by incorporating discrete latent variables and leveraging techniques that enhance the quality and diversity of generated samples while maintaining efficient training dynamics. VQ-VAEs propose a dynamic prior distribution generation mechanism that diverges from the static priors commonly associated with traditional VAEs. This dynamic approach allows for more adaptive and context-aware latent variable representations, thereby potentially capturing complex data structures more effectively. Unlike autoregressive prior models such as PixelCNN, which, despite their ability to model dependencies across data dimensions, suffer from significant computational inefficiencies and lack flexibility in handling diverse datasets. In our work, we propose to use a generative quantum-compatible approach to help alleviate the shortcomings of autoregressive prior model in VQ-VAEs. Restricted Boltzmann Machines (RBMs) are a viable alternative prior model that can learn prior distributions in a faster and more flexible manner. In this research endeavor, we meticulously curate a state-of-the-art dataset leveraging satellite MODIS data in conjunction with VIIRS fire masks, derived from Fire Radiative Power (FRP), thereby encapsulating diverse wildfire scenarios and environmental contexts. We developed a conditional VQ-VAE architecture with the RBM prior model that is trained in a supervised manner for segmenting wildfire masks. This innovative approach synergistically harnesses deep learning capabilities, enabling the generation of segmentation maps characterized by heightened precision, granularity, and contextual relevance. Furthermore, replacing the autoregressive prior learning method proposed by the original VQ-VAE with a prior density approximation via quantum-compatible RBM facilitates expedited inference processes, augments flexibility in prior sampling, optimizes computational efficiency and establishes a groundbreaking benchmark in wildfire segmentation methodologies.

quantum machine learning↗

A seamless multiscale operator neural network for inferring bubble dynamics

Modelling multiscale systems from nanoscale to macroscale requires the use of atomistic and continuum methods and, correspondingly, different computer codes. Here, we develop a seamless method based on DeepONet, which is a composite deep neural network (a branch and a trunk network) for regressing operators. In particular, we consider bubble growth dynamics, and we model tiny bubbles of initial size from 100 nm to 10 $\mathrm {\mu }\textrm {m}$ , modelled by the Rayleigh–Plesset equation in the continuum regime above 1 $\mathrm {\mu }\textrm {m}$ and the dissipative particle dynamics method for bubbles below 1 $\mathrm {\mu }\textrm {m}$ in the atomistic regime. After an offline training based on data from both regimes, DeepONet can make accurate predictions of bubble growth on-the-fly (within a fraction of a second) across four orders of magnitude difference in spatial scales and two orders of magnitude in temporal scales. The framework of DeepONet is general and can be used for unifying physical models of different scales in diverse multiscale applications.

Mechanics↗

Implementation of a Binary Neural Network on a Passive Array of Magnetic Tunnel Junctions

The increasing scale of neural networks and their growing application space have produced demand for more energy- and memory-efficient artificial-intelligence-specific hardware. Avenues to mitigate the main issue, the von Neumann bottleneck, include in-memory and near-memory architectures, as well as algorithmic approaches. In this report we leverage the low-power and the inherently binary operation of magnetic tunnel junctions (MTJs) to demonstrate neural network hardware inference based on passive arrays of MTJs. In general, transferring a trained network model to hardware for inference is confronted by degradation in performance due to device-to-device variations, write errors, parasitic resistance, and nonidealities in the substrate. To quantify the effect of these hardware realities, we benchmark 300 unique weight matrix solutions of a two-layer perceptron to classify the Wine dataset for both classification accuracy and write fidelity. Despite device imperfections, we achieve software-equivalent accuracy of up to 95.3% with proper tuning of network parameters in 15 x 15 MTJ arrays having a range of device sizes. The success of this tuning process shows that new metrics are needed to characterize the performance and quality of networks reproduced in mixed signal hardware.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Paleoclimatic implications of glacial fluctuations in the Sierra Nevada del Cocuy, northern Andes, Colombia, during the Lateglacial and Holocene

The reconstruction of former mountain glaciers from geomorphic mapping and cosmogenic-nuclide surface-exposure dating provides a unique opportunity to infer patterns of past terrestrial climate variability. Tropical mountain glaciers are particularly valuable as there are comparatively few terrestrial climate proxies at equatorial latitudes relative to higher latitudes. As the single largest climate zone on Earth, the tropics play an outsized role in mediating global climate via the ocean-atmosphere transfer of latent heat and water vapor. Nonetheless, there remains a persistent gap in our understanding of how the tropics influenced – or were influenced by – the high-magnitude climate shifts of the Late Pleistocene, and whether this high-energy region simply responded to extratropical forcing or was itself a driver of global climatic change. To help address this knowledge gap, we analyzed geologic evidence for past glacial fluctuations in three adjacent valleys in the Sierra Nevada del Cocuy, the highest subrange of the Eastern Cordillera in the Colombian Andes, to provide a terrestrial record of atmospheric temperature during the latter part of Termination 1. Coupled with geomorphic mapping and paleo-snowline reconstructions, our beryllium-10 glacial chronology indicates that glaciers in the humid inner tropics underwent pronounced growth and gradual decay during the Antarctic Cold Reversal (14.5–12.8 ka) and Younger Dryas (12.8–11.7 ka) periods, respectively, following a trend that, according to directly dated moraine records from throughout both polar hemispheres, appears to have been global. While the specific mechanism(s) behind this large-scale behavior remains to be corroborated, we revisit the hypothesis that ocean-atmosphere heat transfer and water vapor flux are key drivers of abrupt Lateglacial temperature fluctuations. Subsequent to the Lateglacial, deglaciation of the Sierra Nevada del Cocuy accelerated during the Early Holocene, a pattern also observed in other tropical glacier records. More recently, the magnitude of snowline rise and glacier retreat over the last two centuries supports the view that modern tropospheric warming is anomalously strong at least relative to the last ∼16,000 years.

Andes↗

Field intercomparison of ice nucleation measurements: the Fifth International Workshop on Ice Nucleation Phase 3 (FIN-03)

Abstract. The third phase of the Fifth International Ice Nucleation Workshop (FIN-03) was conducted at the Storm Peak Laboratory in Steamboat Springs, Colorado, in September 2015 to facilitate the intercomparison of instruments measuring ice-nucleating particles (INPs) in the field. Instruments included two online and four offline measurement systems for INPs, which are a subset of those utilized in the laboratory study that comprised the second phase of FIN (FIN-02). The composition of the total aerosols was characterized using the Particle Analysis by Laser Mass Spectrometry (PALMS) and Wideband Integrated Bioaerosol Sensor (WIBS) instruments, and aerosol size distributions were measured by a laser aerosol spectrometer (LAS). The dominant total particle compositions present during FIN-03 were composed of sulfates, organic compounds, and nitrates, as well as particles derived from biomass burning. Mineral-dust-containing particles were ubiquitous throughout and represented 67 % of supermicron particles. Total WIBS fluorescing particle concentrations for particles with diameters of > 0.5 µm were 0.04 ± 0.02 cm−3 (0.1 cm−3 highest; 0.02 cm−3 lowest), typical of the warm season in this region and representing ≈ 9 % of all particles in this size range as a campaign average. The primary focus of FIN-03 was the measurement of INP concentrations via immersion freezing at temperatures > −33 °C. Additionally, some measurements were made in the deposition nucleation regime at these same temperatures, representing one of the first efforts to include both mechanisms within a field campaign. INP concentrations via immersion freezing agreed within factors ranging from nearly 1 to 5 times on average between matched (time and temperature) measurements, and disagreements only rarely exceeded 1 order of magnitude for sampling times coordinated to within 3 h. Comparisons were restricted to temperatures lower than −15 °C due to the limits of detection related to sample volumes and very low INP concentrations. Outliers of up to 2 orders of magnitude occurred between −25 and −18 °C; a better agreement was seen at higher and lower temperatures. Although the 5–10 factor agreement of INP measurements found in FIN-03 aligned with the results of the FIN-02 laboratory comparison phase, giving confidence in progress of this measurement field, this level of agreement still equates to temperature uncertainties of 3.5 to 5 °C that may not be sufficient for numerical cloud modeling applications that utilize INP information. INP activity in the immersion-freezing mode was generally found to be an order of magnitude or more, making it more efficient than in the deposition regime at 95 %–99 % water relative humidity, although this limited data set should be augmented in future efforts. To contextualize the study results, an assessment was made of the composition of INPs during the late-summer to early-fall period of this study inferred through comparison to existing ice nucleation parameterizations and through measurement of the influence of thermal and organic carbon digestion treatments on immersion-freezing ice nucleation activity. Consistent with other studies in continental regions, biological INPs dominated at temperatures of > −20 °C and sometimes colder, while arable dust-like or other organic-influenced INPs were inferred to dominate below −20 °C.

54 ENVIRONMENTAL SCIENCES↗

Electron energy distributions in the extended gas nebulae associated with high-z AGN: Maxwell–Boltzmann versus κ distributions

ABSTRACT Emission line observations together with photoionization models provide important information about the ionization mechanisms, densities, temperatures, and metallicities in active galactic nuclei (AGN)-ionized gas. Photoionization models usually assume Maxwell–Boltzmann (M–B) electron energy distributions (EED), but it has been suggested that using κ distributions may be more appropriate and could potentially solve the discrepancies in temperatures and abundances found in HII regions and planetary nebulae. We consider the impact of the presence of κ distributions in photoionized nebulae associated with AGN and study how this might affect spectral modelling and abundance analyses for such regions. Using the photoionization code MAPPINGS 1e, we compute models adopting M–B and κ distributions of electron energies and compare the behaviour of emission line ratios for different values of κ, gas metallicity, density, ionization parameter, and SED slope. We find that the choice of EED can have a large impact on some UV and optical emission lines emitted by photoionized nebulae associated with AGN, and that the impact of adopting a κ distribution is strongly dependent on gas metallicity and ionization parameter. We compile a sample of line ratios for 143 type 2 AGN and compare our models against the observed line ratios. We find that for 98 objects, κ distributions provide a better fit to the observed line ratios than M–B distributions. In addition, we find that adopting κ-distributed electron energies results in significant changes in the inferred gas metallicity and ionization parameter in a significant fraction of objects.

Morais, S. G. (ORCID:0000000200927813)↗

Autotrophy of green non-sulphur bacteria in hot spring microbial mats: biological explanations for isotopically heavy organic carbon in the geological record

Inferences about the evidence of life recorded in organic compounds within the Earth's ancient rocks have depended on 13C contents low enough to be characteristic of biological debris produced by the well-known CO2 fixation pathway, the Calvin cycle. 'Atypically' high values have been attributed to isotopic alteration of sedimentary organic carbon by thermal metamorphism. We examined the possibility that organic carbon characterized by a relatively high 13C content could have arisen biologically from recently discovered autotrophic pathways. We focused on the green non-sulphur bacterium Chloroflexus aurantiacus that uses the 3-hydroxypropionate pathway for inorganic carbon fixation and is geologically significant as it forms modern mat communities analogous to stromatolites. Organic matter in mats constructed by Chloroflexus spp. alone had relatively high 13C contents (-14.9%) and lipids diagnostic of Chloroflexus that were also isotopically heavy (-8.9% to -18.5%). Organic matter in mats constructed by Chloroflexus in conjunction with cyanobacteria had a more typical Calvin cycle signature (-23.5%). However, lipids diagnostic of Chloroflexus were isotopically enriched (-15.1% to -24.1%) relative to lipids typical of cyanobacteria (-33.9% to -36.3%). This suggests that, in mats formed by both cyanobacteria and Chloroflexus, autotrophy must have a greater effect on Chloroflexus carbon metabolism than the photoheterotrophic consumption of cyanobacterial photosynthate. Chloroflexus cell components were also selectively preserved. Hence, Chloroflexus autotrophy and selective preservation of its products constitute one purely biological mechanism by which isotopically heavy organic carbon could have been introduced into important Precambrian geological features.

NASA Discipline Exobiology↗

Analysis of ISS Plasma Interaction

To date, the International Space Station (ISS) has been one of the largest objects flown in lower earth orbit (LEO). The ISS utilizes high voltage solar arrays (160V) that are negatively grounded leading to pressurized elements that can float negatively with respect to the plasma. Because laboratory measurements indicate a dielectric breakdown potential difference of 80V, arcing could occur on the ISS structure. To overcome the possibility of arcing and clamp the potential of the structure, two Plasma Contactor Units (PCUs) were designed, built, and flown. Also a limited amount of measurements of the floating potential for the present ISS configuration were made by a Floating Potential Probe (FPP), indicating a minimum potential of 24 Volts at the measurement location. A predictive tool, the ISS Plasma Interaction Model (PIM) has been developed accounting for the solar array electron collection, solar array mast wire and effective conductive area on the structure. The model has been used for predictions of the present ISS configuration. The conductive area has been inferred based on available floating potential measurements. Analysis of FPP and PCU data indicated distribution of the conductive area along the Russian segment of the ISS structure. A significant input to PIM is the plasma environment. The International Reference Ionosphere (IRI 2001) was initially used to obtain plasma temperature and density values. However, IRI provides mean parameters, leading to difficulties in interpretation of on-orbit data, especially at eclipse exit where maximum charging can occur. This limits our predicative capability. Satellite and Incoherent Scatter Radar (ISR) data of plasma parameters have also been collected. Approximately 130,000 electron temperature (Te) and density (Ne) pairs for typical ISS eclipse exit conditions have been extracted from the reduced Langmuir probe data flown aboard the NASA DE-2 satellite. Additionally, another 18,000 Te and Ne pairs of ISR data from several radar locations around the globe were used to assure consistency of the satellite data. PIM predictions for ISS charging made with this data correlated very well with FPP data, indicating that the general physics of spacecraft charging with high voltage solar arrays have been captured. The predictions also provided the probabilities of occurrences for ISS charging. These probabilities give a numerical measure of the number of times when the ISS will approach or exceed the vehicle plasma hazard conditions for each configuration. In this paper we shall present the interaction mechanisms between the ISS and the surrounding plasma and give an overview of the PIM components. PIM predictions are compared with available data followed by a discussion of the variability of plasma parameters and the conductive area on the ISS. The ISS PIM will be further tested and verified as data from the Floating Potential Measurement Unit become available, and construction of the ISS continues.

Reddell, Brandon↗

Advances in Entry Modeling for Impact Risk Assessment

A summary of recent advancements in the detailed modeling of asteroid atmospheric entry processes made through NASA’s Asteroid Threat Assessment Project (ATAP) is presented. Understanding, and accurately modeling these processes and their associated uncertainties is critical to predicting all downstream impact effects, such as blast wave and thermal damage footprints. Furthermore, there is (perhaps thankfully) a dearth of empirical data for large impactors of the kind that would pose a threat to human populations, on which to anchor and/or validate models used in risk assessments. Therefore, we must rely heavily on detailed theoretical and numerical modeling to develop robust assessments for decision makers. To that end, ATAP has made some significant progress in advancing the capabilities in this area. Two areas in particular are highlighted in the present work: meteoroid ablation mechanisms, and bolide luminosity. The first research area – meteoroid ablation mechanisms – has focused on performing novel high-enthalpy wind tunnel experiments on meteorites and meteorite analogs, and utilizing the resulting data to develop high-fidelity models for impactor mass loss at scale. These efforts have resulted in several insights. Of note, these data suggest a differential vaporization process where volatiles are liberated preferentially when the asteroidal material is subject to high heat, while refractory components remain in the molten layer on the surface. A numerical model has been developed which considers this phenomena, and its effect on the bulk impactors effective heat of ablation is examined. The second research thrust that is discussed is focused on accurate modeling of impactor radiation phenomena. While another submission to the conference will discuss the application of our approach to thermal ground damage modeling, here, we present an overview of our extensive efforts to utilize available ground- and space-based observations of large bolides (~1m diameter, and above) to inform and validate our detailed modeling approaches. These methods have been shown to accurately reconstruct the detailed spectra for the Benesov bolide, as well as approximate the burn footprint for the Tunguska event. Recent effort has focused on reconciling light curve data from multiple sources (all-sky camera networks, GLM, US government sensors), using our validated model, and providing a calibrated model for luminous efficiency which can then be utilized to infer impactor properties such as shape, mass, and composition. This work will be demonstrated through an exemplar case study focusing on a large bolide event with observational data (e.g. Chelyabinsk, Flensburg, etc.). Finally, our team’s assessment on the current maturity of atmospheric entry models, and priorities for future research will be provided.

ATAP↗

Magnetoanisotropic thermal transport in a topological Josephson junction

Highlights: • Thermal conductance is sensitive to orientation and magnitude of the exchange field. • The NDTC is exhibited and tuned by the magnitude of the exchange field. • The interplay of zero-energy ABSs and thermal transport has been well inferred. • Our findings can apply to the energy control in hybridized mesoscopic systems. The thermal transport in Joesphson junctions, particularly, based on topological insulators (TIs), has not been extensively studied so far. Herein, we theoretically investigate the magnetoanisotropic thermal transport in a TI-based Josephson hybrid structure with a ferromagnetic insulator sandwiched between two superconductors (SCs). It is shown that the thermal conductance and current both are considerably sensitive to the orientation and magnitude of the exchange field. In particular, the negative differential thermal conductance is exhibited, which is modulated by the exchange field. More interestingly, zero-energy Andreev bound states are confirmed with the help of the thermal conductance, which in turn demonstrates itself to be of significance in the thermal transport. It is expected that these findings could have potential applications for the energy control in various hybridized mesoscopic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multi-fidelity information fusion with concatenated neural networks

Recently, computational modeling has shifted towards the use of statistical inference, deep learning, and other data-driven modeling frameworks. Although this shift in modeling holds promise in many applications like design optimization and real-time control by lowering the computational burden, training deep learning models needs a huge amount of data. This big data is not always available for scientific problems and leads to poorly generalizable data-driven models. This gap can be furnished by leveraging information from physics-based models. Exploiting prior knowledge about the problem at hand, this study puts forth a physics-guided machine learning (PGML) approach to build more tailored, effective, and efficient surrogate models. For our analysis, without losing its generalizability and modularity, we focus on the development of predictive models for laminar and turbulent boundary layer flows. In particular, we combine the self-similarity solution and power-law velocity profile (low-fidelity models) with the noisy data obtained either from experiments or computational fluid dynamics simulations (high-fidelity models) through a concatenated neural network. We illustrate how the knowledge from these simplified models results in reducing uncertainties associated with deep learning models applied to boundary layer flow prediction problems. The proposed multi-fidelity information fusion framework produces physically consistent models that attempt to achieve better generalization than data-driven models obtained purely based on data. While we demonstrate our framework for a problem relevant to fluid mechanics, its workflow and principles can be adopted for many scientific problems where empirical, analytical, or simplified models are prevalent. In line with grand demands in novel PGML principles, this work builds a bridge between extensive physics-based theories and data-driven modeling paradigms and paves the way for using hybrid physics and machine learning modeling approaches for next-generation digital twin technologies.

42 ENGINEERING↗

Enhanced electro-osmosis in propylene carbonate salt solutions

Properties of solid–liquid interfaces and surface charge characteristics mediate ionic and molecular transport through porous systems, affecting many processes such as separations. Herein, we report experiments designed to probe the electrochemical properties of solid–liquid interfaces using a model system of a single polyethylene terephthalate (PET) pore in contact with aqueous and propylene carbonate solutions of LiClO 4 . First, the existence and polarity of surface charges were inferred from current–voltage curves recorded when a pore was placed in contact with a LiClO4 concentration gradient. Second, the electro-osmotic transport of uncharged polystyrene particles through the PET pore provided information on the polarity and the magnitude of the pore walls’ zeta potential. Our experiments show that the PET pores become effectively positively charged when in contact with LiClO 4 solutions in propylene carbonate, even though in aqueous LiClO 4 , the same pores are negatively charged. Additionally, the electro-osmotic velocity of the particles revealed a significantly higher magnitude of the positive zeta potential of the pores in propylene carbonate compared to the magnitude of the negative zeta potential in water. The presented methods of probing the properties of solid–liquid interfaces are expected to be applicable to a wide variety of solid and liquid systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantifying the effect of late bombardment on terrestrial zircons

The first 500 million years of Earth history is thought to be a period of intense planetary bombardment, but the timing and flux of this meteorite bombardment is poorly understood. In particular, on the basis of an inferred lunar impact history, some workers have hypothesized a ~3.9 Ga terminal cataclysm (TC) in which there was marked increase in the impact flux affecting the Moon, the Earth and possibly other terrestrial planets. Minerals that survived this enigmatic period offer a way to test early planetary bombardment models as they may contain telltale micro- to nanoscale shock features. Here, we present results from a numerical modeling calculation that assesses the probability that a zircon residing in the crust would escape shock melting or shock deformation during a TC bombardment event. Even with conservative pressure estimates for zircon shock deformation and intermediate bombardment intensities, we find that only ~6% of ≥4.0 Ga crust would be expected to survive a 3.9 Ga cataclysm without experiencing either complete melting or zircon shock metamorphism. In this study, we couple this modeling with a search for shock effects in the oldest zircons from the Acasta Gneiss Complex, which would have been present in the Earth's crust during a putative 3.9 Ga TC. Spatially correlated electron and NanoSIMS ion microscopy of 4.02 Ga igneous zircons from Acasta reveals no evidence of ancient shock. These data, together with similar results from other Hadean zircon suites, confirm that a post-Hadean TC is unlikely to have occurred. We suggest that the dearth of pre-3.9 Ga terrestrial crust and zircons is instead best explained by endogenic processes related to the mechanisms of early crust formation. Our modeling allows us to evaluate bombardment scenarios from the terrestrial zircon record by applying probabilistic interpretations to zircon shock deformation data. This approach will be valuable for other planetary bodies, allowing broader conclusions to be drawn from geographically limited datasets.

58 GEOSCIENCES↗

Atomic-resolution imaging as a mechanistic tool for studying single-site heterogeneous catalysis

Heterogeneous catalysts dominate the chemical industry but typically feature diverse, incompletely defined active sites. Thus, describing structure-activity relationships, unlike homogeneous catalysts, remains challenging. In contrast, molecularly defined single-site heterogeneous catalysts (SSHCs), using appropriate tools, are poised to address these challenges and provide new avenues for catalysis research and development. The present study explores eco-friendly H 2 production mediated by discrete MoO 2 sites supported on carbon nanohorns (CNHs) and active for alcohol dehydrogenation. While informative, detailed ensemble EXAFS/XANES/XPS, kinetic measurements, and DFT analysis alone cannot provide a full molecular picture of the reaction pathway. Here, using single-molecule atomic-resolution time-resolved electron microscopy (SMART-EM), we propose the identification of four key catalytic intermediates anchored to CNHs and uncover a new reaction pathway involving alkoxide/hemiacetal equilibration and acetal oligomerization. Furthermore, these intermediates are inferred through a combination of theory and SMART-EM, showcasing the potential of SMART-EM as a complementary tool for exploring mechanistic hypotheses in catalysis.

36 MATERIALS SCIENCE↗

Identifying Bayesian optimal experiments for uncertain biochemical pathway models

Abstract Pharmacodynamic (PD) models are mathematical models of cellular reaction networks that include drug mechanisms of action. These models are useful for studying predictive therapeutic outcomes of novel drug therapies in silico. However, PD models are known to possess significant uncertainty with respect to constituent parameter data, leading to uncertainty in the model predictions. Furthermore, experimental data to calibrate these models is often limited or unavailable for novel pathways. In this study, we present a Bayesian optimal experimental design approach for improving PD model prediction accuracy. We then apply our method using simulated experimental data to account for uncertainty in hypothetical laboratory measurements. This leads to a probabilistic prediction of drug performance and a quantitative measure of which prospective laboratory experiment will optimally reduce prediction uncertainty in the PD model. The methods proposed here provide a way forward for uncertainty quantification and guided experimental design for models of novel biological pathways.

97 MATHEMATICS AND COMPUTING↗

The ghost in the radiation: robust encodings of the black hole interior

We reconsider the black hole firewall puzzle, emphasizing that quantum error-correction, computational complexity, and pseudorandomness are crucial concepts for understanding the black hole interior. We assume that the Hawking radiation emitted by an old black hole is pseudorandom, meaning that it cannot be distinguished from a perfectly thermal state by any efficient quantum computation acting on the radiation alone. We then infer the existence of a subspace of the radiation system which we interpret as an encoding of the black hole interior. This encoded interior is entangled with the late outgoing Hawking quanta emitted by the old black hole, and is inaccessible to computationally bounded observers who are outside the black hole. Specifically, efficient operations acting on the radiation, those with quantum computational complexity polynomial in the entropy of the remaining black hole, commute with a complete set of logical operators acting on the encoded interior, up to corrections which are exponentially small in the entropy. Thus, under our pseudorandomness assumption, the black hole interior is well protected from exterior observers as long as the remaining black hole is macroscopic. On the other hand, if the radiation is not pseudorandom, an exterior observer may be able to create a firewall by applying a polynomial-time quantum computation to the radiation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A physics-constrained neural network for multiphase flows

The present study develops a physics-constrained neural network (PCNN) to predict sequential patterns and motions of multiphase flows (MPFs), which includes strong interactions among various fluid phases. To predict the order parameters, which locate individual phases in the future time, a neural network (NN) is applied to quickly infer the dynamics of the phases by encoding observations. The multiphase consistent and conservative boundedness mapping algorithm (MCBOM) is next implemented to correct the predicted order parameters. This enforces the predicted order parameters to strictly satisfy the mass conservation, the summation of the volume fractions of the phases to be unity, the consistency of reduction, and the boundedness of the order parameters. Then, the density of the fluid mixture is updated from the corrected order parameters. Finally, the velocity in the future time is predicted by another NN with the same network structure, but the conservation of momentum is included in the loss function to shrink the parameter space. The proposed PCNN for MPFs sequentially performs (NN)-(MCBOM)-(NN), which avoids nonphysical behaviors of the order parameters, accelerates the convergence, and requires fewer data to make predictions. Numerical experiments demonstrate that the proposed PCNN is capable of predicting MPFs effectively.

Mechanics↗