Search NASA⌕ Search

SEARCH · Search NASA

Results for “uncertainty characterization”

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 397 records · Page 22

Exploration of multifidelity UQ sampling strategies for computer network applications.

Network modeling is a powerful tool to enable rapid analysis of complex systems that can be challenging to study directly using physical testing. Here, two approaches are considered: emulation and simulation. The former runs real software on virtualized hardware, while the latter mimics the behavior of network components and their interactions in software. Although emulation provides an accurate representation of physical networks, this approach alone cannot guarantee the characterization of the system under realistic operative conditions. Operative conditions for physical networks are often characterized by intrinsic variability (payload size, packet latency, etc.) or a lack of precise knowledge regarding the network configuration (bandwidth, delays, etc.); therefore uncertainty quantification (UQ) strategies should be also employed. UQ strategies require multiple evaluations of the system with a number of evaluation instances that roughly increases with the problem dimensionality, i.e., the number of uncertain parameters. It follows that a typical UQ workflow for network modeling based on emulation can easily become unattainable due to its prohibitive computational cost. In this paper, a multifidelity sampling approach is discussed and applied to network modeling problems. The main idea is to optimally fuse information coming from simulations, which are a low-fidelity version of the emulation problem of interest, in order to decrease the estimator variance. By reducing the estimator variance in a sampling approach it is usually possible to obtain more reliable statistics and therefore a more reliable system characterization. Several network problems of increasing difficulty are presented. For each of them, the performance of the multifidelity estimator is compared with respect to the single fidelity counterpart, namely, Monte Carlo sampling. For all the test problems studied in this work, the multifidelity estimator demonstrated an increased efficiency with respect to MC.

97 MATHEMATICS AND COMPUTING↗

Thermal property characterization of phase change materials in building applications: A systematic review of fundamentals, recent progress, and future directions

Phase change materials (PCMs) can reduce building peak loads and enable demand-responsive thermal energy storage (TES), but their deployment depends on reliable measurement and interpretation of thermal properties across laboratory, intermediate, and application scales. Here, this review systematically examines characterization methods, testing protocols, and recent advances for neat PCMs and PCM composites, emphasizing thermal conductivity, enthalpy-related properties (phase change temperature, latent heat, specific heat), and cycling stability. For thermal conductivity, we compare steady-state and transient techniques and note limitations when phase transition and contact resistance affect measurements. For enthalpy–temperature characterization, we discuss differential scanning calorimetry together with intermediate- and bulk-scale methods, including T-history, heat flow meter testing, and three-layer calorimetry (3LC), to generate application-relevant enthalpy–temperature profiles. Cycling stability is organized into four experimental families: thermoelectric–air, fully thermoelectric, water-bath, and in situ chamber approaches, with attention to separating reversible supercooling from true degradation such as phase segregation. We highlight emerging noncontact diagnostics, including infrared thermography and embedded sensing, for spatially resolved validation and multiscale interpretation. Finally, we review the growing use of AI and machine learning for property prediction, inverse characterization from experimental signals, and real-time state estimation in building-integrated TES. Key needs include harmonized protocols, interlaboratory benchmarking, uncertainty reporting, and metadata-rich datasets to accelerate reproducible PCM qualification for grid-flexible buildings.

AI↗

Predicting Thermal Performance of an Enhanced Geothermal System From Tracer Tests in a Data Assimilation Framework

Abstract Predicting the thermal performance of an enhanced geothermal system (EGS) requires a comprehensive characterization of the underlying fracture flow patterns from practically available data such as tracer data. However, due to the inherent complexities of subsurface fractures and the generally insufficient geological/geophysical data, interpreting tracer data for fracture flow characterization and thermal prediction remains a challenging task. The present study aims to tackle the challenge by leveraging a data assimilation method to maximize the utilization of information inherently contained in tracer data, and meanwhile maintain the flexibility to handle various uncertainties. A tracer data interpretation framework was proposed with the following three components integrated: (a) We use principal component analysis (PCA) to reduce the dimensionality of model parameter space. (b) We use ES‐MDA (ensemble smoother with multiple data assimilation) to invert for fracture aperture/flow fields and obtain posterior model ensembles for uncertainty quantification. Various data types are assimilated jointly to improve the predictive ability of the posterior ensemble. (c) The inverted fracture aperture fields are then incorporated into reservoir models to predict thermal performance. We developed a field‐scale EGS model to verify the ability of the framework to characterize highly heterogeneous fracture aperture/flow fields and predicting thermal performance. We also applied the framework to a mesoscale field experiment to demonstrate its potential application in real‐world geothermal reservoirs. The results indicate that the proposed framework can effectively retrieve fracture flow information from tracer data for thermal prediction and uncertainty quantification, and thus provide informative guidance for EGS optimization and risk management.

15 GEOTHERMAL ENERGY↗

Unveiling Feedstock Variability: Insights into Corn Stover Conversion - Part I: Physicochemical Properties and Self-Degradation

Transforming agricultural waste into biofuels and bioproducts is crucial to advancing a low-carbon bioeconomy. However, the inherent variability in the composition and quality introduces uncertainties in the conversion efficiency and poses challenges in process development. Through integrating a high-throughput conversion system, material characterization techniques, and advanced data analysis tools, this study investigates the variability of corn stover and its subsequent impacts on carbohydrate conversion. The findings reveal that indoor storage substantially reduces the moisture and ash content and soil contamination, while other properties remain largely unchanged. Self-degradation due to microbial activity during storage decreases the carbohydrate content of corn stover but enhances glucose and xylose yields. A negative correlation is observed between sugar yields and lignin content across samples with varying ash and moisture content. The inhibitory effect of lignin diminishes in self-degraded samples likely due to the disrupted cell wall structure. Although self-degradation slightly increases cellulose crystallinity, no strong correlation was observed between the crystallinity and sugar yield. Hot water pretreatment under mild conditions effectively mitigates inherent variability, consistently improving the sugar yield from corn stover by up to 50%. By elucidating the feedstock variability and its impact on convertibility, these findings offer valuable insights into appropriate feedstock handling and management, highlighting potential strategies to address variability challenges.

09 BIOMASS FUELS↗

Femtosecond pulse parameter estimation from photoelectron momenta using machine learning

Abstract Deep learning models have provided huge interpretation power for image-like data. Specifically, convolutional neural networks (CNNs) have demonstrated incredible acuity for tasks such as feature extraction or parameter estimation. Here we test CNNs on strong-field ionization photoelectron spectra, training on theoretical data sets to ‘invert’ experimental data. Pulse characterization is used as a ‘testing ground’, specifically we retrieve the laser intensity, where ‘traditional’ measurements typically lead to 20% uncertainty. We report on crucial data augmentation techniques required to successfully train on theoretical data and return consistent results from experiments, including accounting for detector saturation. The same procedure can be repeated to apply CNNs in a range of scenarios for strong-field ionization. Using a predictive uncertainty estimation, reliable laser intensity uncertainties of a few percent can be extracted, which are consistently lower than those given by traditional techniques. Using interpretability methods can reveal parts of the distribution that are most sensitive to laser intensity, which can be directly associated with holographic interferences. The CNNs employed provide an accurate and convenient ways to extract parameters, and represent a novel interpretational tool for strong-field ionization spectra.

Szołdra, Tomasz (ORCID:0000000228970506)↗

Borehole Disposal of Radioactive Waste in Israel - 20388

Israel is assessing borehole disposal of radioactive waste. With limited geological options for disposal, intermediate-depth borehole disposal is being considered in the arid Yamin Plain region of the northeastern Negev desert at depths of several hundred meters below ground surface in the vadose zone. Unlike deep borehole disposal of several kilometers, which relies on emplacement below the depth of recirculating groundwater, the safety case for intermediate-depth borehole disposal relies more on the aridity of the vadose zone, the robustness of the waste package, and the other engineered barriers (e.g., seals, backfill materials) in the disposal borehole. As a tradeoff, the shallower depth may allow the use of a larger diameter borehole than would be possible for a deeper borehole. To study the suitability of the Yamin Plain region for borehole disposal, a small-diameter characterization borehole is being planned to retrieve core samples and to better understand the vadose zone geomechanical and hydrogeochemical properties and percolation flux. The information from the characterization borehole will inform the safety case, and together with performance assessment analyses will help to identify key areas of uncertainty and guide future research and development activities aimed at demonstrating the feasibility of the intermediate-depth borehole disposal concept in Israel. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

StarDICE III: characterization of the photometric instrument with a collimated beam projector

The measurement of Type Ia supernovae magnitudes provides cosmological distances, which constrain dark energy parameters. Current and upcoming large photometric surveys require improved photometric calibration precision to reduce systematic uncertainties in cosmological constraints. The StarDICE experiment aims to establish accurate broad-band flux references for these surveys, targeting sub-percent precision in magnitude measurements. Achieving this requires precise filter bandpass measurements for both StarDICE and survey instruments with sub-nanometre accuracy. To this end, we developed the Collimated Beam Projector (CBP), an optical device for calibrating the throughput of astronomical telescopes and their filters. The CBP uses a tunable laser source and a reversed telescope to emit a parallel monochromatic light beam, continuously monitored in flux and wavelength. The CBP output flux is measured with a large-area photodiode calibrated relative to a NIST photodiode. Using CBP measurements, we derive the StarDICE telescope throughput and filter transmissions, anchoring them to NIST’s absolute calibration. After analysing systematic uncertainties, we achieved sub-nanometre accuracy for filter central wavelengths, measured filter transmission with ~0.5 per cent precision per 1 nm bin, and detected out-of-band leakages at a relative level of 10 –4 ⁠. Furthermore, we synthesized equivalent transmission for full pupil illumination from four sampled positions in the StarDICE telescope mirror, with ~0.2 nm accuracy for central wavelengths and 7 mmag for broad-band fluxes. This demonstrates our ability to characterize telescope throughput down to the millimagnitude, paving the way for future developments, such as the Rubin-CBP for measuring the LSST at Vera Rubin Observatory, and a portable CBP version for in-situ transmission monitoring.

Calibration↗

Carbon Dioxide Removal Measurement, Reporting, and Verification Simulation toolkit (CDR MRVSim) v0.1

CDR MRVSim is a statistical software toolkit for techno-economic analysis of measurement, reporting, and verification of multiple carbon dioxide removal technologies. The current version applies Monte Carlo simulation to existing datasets estimate the cost and uncertainty of measuring the soil organic carbon (SOC) content of agricultural land, accounting for multiple sources of uncertainty, using data from field measurements. We are planning to incorporate enhanced rock weathering and other technologies into the software. The current model leverages existing data to characterize underlying variability in soil organic carbon and related parameters, including bulk density, in an agricultural field attempting to increase its SOC. We then simulate baseline and post-intervention "measurement campaigns" in which some number of SOC and bulk density measurements are conducted using a selected technology. This enables estimates of the cost and accuracy of measuring changes in SOC in the simulated field. The primary advantages over similar software are: 1) ability to quantify tradeoffs between cost and uncertainty in MRV across a range of possible MRV approaches 2) focus on guiding development of novel sensors by determining desirable sets of characteristics

Sherwin, Evan [Lawrence Berkeley National Laborato↗

Deep Gaussian process-based cost-aware batch Bayesian optimization for complex materials design campaigns

The accelerating pace and expanding scope of materials discovery demand optimization frameworks that efficiently navigate vast design spaces with complex response surfaces while judiciously allocating limited evaluation resources. We present a cost-aware, batch Bayesian optimization scheme powered by deep Gaussian process (DGP) surrogates and a heterotopic querying strategy. Our DGP surrogate, formed by stacking GP layers, models complex hierarchical relationships among high-dimensional compositional features and captures correlations across multiple target properties, propagating uncertainty through successive layers. We integrate evaluation cost into an upper-confidence-bound acquisition extension, which, together with heterotopic querying, proposes small batches of candidates in parallel, balancing exploration of under-characterized regions with exploitation of high-mean, low-variance predictions across correlated properties. Applied to refractory high-entropy alloys for high-temperature applications, our framework converges to optimal formulations in fewer iterations with cost-aware queries than conventional GP-based BO, highlighting the value of deep, uncertainty-aware, cost-sensitive strategies in materials campaigns.

36 MATERIALS SCIENCE↗

Characterization of carbonaceous aerosols during TRACER-CAT

Absorbing aerosols (AA) have an important impact on the global radiation budget and cloud properties. The composition and properties of AA can vary substantially throughout the atmosphere, depending on the particle source and the influence of chemical aging. Uncertainties associated with the radiative effects of AA remain substantial. A key contributor to this uncertainty is understanding the extent to which coatings in general, and water uptake especially, alters absorption by AA particles and how this depends on particle composition. We deployed new and existing experimental tools during the Tracking Aerosol Convection Interactions Experiment (TRACER) campaign in Houston, TX as part of the Carbonaceous Aerosols Thrust (CAT) to provide detailed characterization of aerosol optical, chemical, and physical properties. Our TRACER-CAT measurements complemented and expanded on the planned TRACER instrumentation, allowing for more detailed characterization of aerosol properties of relevance to cloud development (a core focus of TRACER), such as the composition of particles that can act as cloud condensation nuclei, than would otherwise be available. Our measurements have allowed for assessment of the relationship(s) between AA optical properties (with a focus on absorption) and the chemical and physical characteristics (including the mixing state of black carbon (BC) containing particles). These field observations occurred in collaboration with Los Alamos National Laboratory in summer 2022 during the TRACER intensive operating period. The instrumentation we co-deployed provided for measurement of (i) multi-wavelength dry aerosol absorption, scattering, and extinction, (ii) the size-dependent composition and abundance of sub-micron aerosol, differentiating between those particles that do and do not contain BC, (iii) BC-specific concentrations and size distributions, (iv) particle size, and (v) the first field measurements at an ARM site of the influence of RH on multi-wavelength absorption by ambient AA. We have leveraged the natural variability of the atmosphere and of aerosol sources in the Houston region to (i) specifically disentangle contributions to light absorption from BC, absorbing organic carbon (brown carbon), and coatings on BC, (ii) characterize the mixing state of BC and assess the factors that give rise to compositional differences between BC-containing and BC-free aerosol, and (iii) establish how water uptake influences absorption and how any such effect depends on particle composition and BC mixing state. Overall, our study contributed to the mission of the Atmospheric System Research program in multiple ways. Through the deployment of complementary, advanced instrumentation for characterization of a wide range of aerosol properties our work helped to maximize the scientific impact of the TRACER campaign. Our work also allowed for development of new insights into the relationship(s) between aerosol composition, hygroscopicity, and the mixing state of BC with aerosol optical properties. Through this, our work has provided knowledge that can improve understanding and model representation of aerosol processes as they affect the Earth’s radiation budget.

54 ENVIRONMENTAL SCIENCES↗

DEEPEN 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. Part of the DEEPEN project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology is based on the method proposed by Poux et al. (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., 2022. For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. More information on this process is described in the following sections. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty. There are two different versions of the Leapfrog model and associated favorability models: - v1.0: The first release in June 2023 - v2.1: The second release, with improvements made to the earthquake catalog (included additional identified events, removed duplicate events), to the temperature model (fixed a deep BHT), and to the index models (updated the seismicity-heat source index models for supercritical and EGS, and the resistivity-insulation index models for all three play types). Also uses the jet color map rather than the magma color map for improved interpretability. - v2.1.1: Updated to include v2.0 uncertainty results (see below for uncertainty model versions) There are two different versions of the associated uncertainty models: - v1.0: The first release in June 2023 - v2.0: The se...

15 GEOTHERMAL ENERGY↗

TEDS Test Plan for Performance Characterization of Packed-bed Thermal Energy Storage Unit

This report proposes a set of test plans for performance measurement and characterization of the packed-bed thermal energy storage (TES) unit for TEDS. The main goals of this test plan are three-fold: (i) support the high-fidelity data measurement for TEDS, (ii) characterize thermal and structural performance of the TEDS TES unit, and (iii) produce high quality validation-oriented data. To achieve the goals, it is first suggested to revisit the instrument calibration and data sampling protocols applied during the startup and commissioning tests of TEDS, based on the observation from the recent CFD validation study using the TEDS TES operational data. The uncertainty quantification of the thermal and hydraulic measurements is also suggested. Lastly, test matrices are provided for charge, discharge, and standby modes of the packed-bed TES unit, respectively, to characterize the thermal and structural performances during the envisioned TEDS operations.

25 ENERGY STORAGE↗

A Systems Approach to Estimating the Uncertainty Limits of X-Ray Radiographic Metrology

Micro- and nanomanufacturing capabilities have rapidly expanded over the past decade to include complex three-dimensional (3D) structure fabrication; however, the metrology required to accurately assess these processes via part inspection and characterization has struggled to keep pace. X-ray computed tomography (CT) is considered an ideal candidate for providing the critically needed metrology on the smallest scales, especially internal features, or inaccessible regions. X-ray CT supporting micro- and nanomanufacturing often push against the poorly understood resolution and variation limits inherent to the machines, which can distort or hide fine structures. In this study, we have developed and experimentally verify a comprehensive analytical uncertainty propagation signal variation flow graph (SVFG) model for X-ray radiography in this work to better understand resolution and image variability limits on the small scale. The SVFG approach captures, quantifies, and predicts variations occurring in the system that limit metrology capabilities, particularly in the micro/nanodomain. This work is the first step to achieving full uncertainty modeling of CT reconstructions and provides insight into improving X-ray attenuation imaging systems. The SVFG methodology framework is applied to generate a complete basis set of functions describing the major sources of variation in radiographs. Five models are identified, covering variation in energy, intensity, length, blur, and position. Radiographic system experiments are defined to measure the parameters required by the SVFGs. Best practices are identified for these measurements. The SVFG models are confirmed via direct measurement of variation to predict variation within 30% on average.

47 OTHER INSTRUMENTATION↗

Infrared spectroscopy data- and physics-driven machine learning for characterizing surface microstructure of complex materials

There is a need to characterize complex materials and their dynamics under reaction conditions to accelerate materials design. Adsorbate vibrational excitations are selective to adsorbate/surface interactions and infrared (IR) spectra associated with activating adsorbate vibrational modes are accurate, capture details of most modes, and can be obtained operando. Current interpretation depends on heuristic peak assignments for simple spectra, precluding the possibility of obtaining detailed structural information. Here, we combine data-based approaches with chemistry-dependent problem formulation to develop physics-driven surrogate models that generate synthetic IR spectra from first-principles calculations. Using synthetic IR spectra of carbon monoxide on platinum, we implement multinomial regression via neural network ensembles to learn probability distributions functions (pdfs) that describe adsorption sites and quantify uncertainty. We use these pdfs to infer detailed surface microstructure from experimental spectra and extend this methodology to other systems as a first step towards characterizing complex interfaces and closing the materials gap.

Science & Technology - Other Topics↗

A phased workflow to define permit‐ready locations for large volume CO 2 injection and storage

Abstract To‐date, only two UIC Class VI permits have been issued by the US Environmental Protection Agency. We illustrate a four‐phase workflow to first identify regional storage resources and then down‐select sites to yield permit‐ready locations that can accept and store large volumes of CO 2 . Specific permit requirements should guide objectives and define deliverables of respective workflow phases. In the first phase we used available regional data and screened structure and injection zones to locate resources that match CO 2 volumes planned to be captured. Available data were also used to assess presence and depth of usable groundwater, the key resource being protected via permitting. We then used advanced, closed‐form, analytical solutions (EASiTool) to estimate CO 2 injectivity into each hydrologically connected injection compartment. In the second phase we acquired and conditioned additional wireline logs and leased available seismic datasets. We interpreted the depositional systems from wireline well‐log character and mapped sandbody geometry to interpolate injection and confining‐zone distribution. Using available data, we mapped faults and locations of freshwater and overpressure (or other capacity‐limiting geologic parameters) in more detail. In the third phase, we used the augmented geologic data to develop a static model for the selected area, extracted the areas of highest interest, and generated and ran dynamic (flow) models. In a fourth phase, we reduced major uncertainties identified in earlier phases. Our case study indicates that to complete preparation of a permit application requires (1) improved lithologic characterization information (thicknesses and horizontal and vertical connectivity) and (2) better definition of poorly defined local faults. © 2023 The Authors. Greenhouse Gases: Science and Technology published by Society of Chemical Industry and John Wiley & Sons Ltd.

58 GEOSCIENCES↗

Land use and turbine technology influences on wind potential in the United States

As clean energy ambitions have expanded, critically evaluating renewable energy supply has become increasingly important to the energy research community and stakeholders. This study examines the onshore wind resource potential for the conterminous United States and its sensitivity to siting constraints and turbine technology innovation. We compile localized regulatory information and use high-resolution data to present multiple siting regimes covering relatively constrained to unconstrained potentials. Our efforts reveal high sensitivity to these variables and sizable uncertainty in the overall wind energy resource potential. Specifically, we find that siting constraints may shift the total capacity available to commercial wind energy by 2.3–15.1 TW. Furthermore, our results illustrate that technology advancement could require larger setbacks from buildings and infrastructure, reducing the total available capacity potential by 20% relative to estimates using current technology, but that this reduction is largely offset by increased generation such that the net effect on generation is 1%. The observed sensitivity to and uncertainty resulting from the variables we analyze suggest there is value in continued study and development of increasingly sophisticated approaches to characterizing wind resource potential.

14 SOLAR ENERGY↗

Developing a heavy-duty vehicle activity database to estimate start and idle emissions

Heavy-duty vehicle start and idling activities were characterized from two datasets to improve emission estimates in the MOtor Vehicle Emission Simulator (MOVES): 1. Fleet DNA from the National Renewable Energy Laboratory (NREL) and 2. A dataset collected by the University of California, Riverside for the California Air Resources Board. Furthermore, the combined dataset includes 564 commercial vehicles, over 23,000 vehicle days of operation and covers seven of the nine heavy-duty source types defined by MOVES. The start and idle activities are characterized and illustrated across MOVES source types, vocations, fleets, days, and hours. This study provides the most comprehensive analysis yet made publicly available to characterize start and idle activity for heavy-duty vehicles within the United States. The results also show there is significant uncertainty in the average heavy-duty idle and start activity due to the large variation in activity across fleets and vocations, and sparsity of nation-wide vehicle population data by vocation.

33 ADVANCED PROPULSION SYSTEMS↗

Persistent Homology Metrics Reveal Quantum Fluctuations and Reactive Atoms in Path Integral Dynamics

Nuclear quantum effects (NQEs) are known to impact a number of features associated with chemical reactivity and physicochemical properties, particularly for light atoms and at low temperatures. In the imaginary time path integral formalism, each atom is mapped onto a “ring polymer” whose spread is related to the quantum mechanical uncertainty in the particle’s position, i.e., its thermal wavelength. A number of metrics have previously been used to investigate and characterize this spread and explain effects arising from quantum delocalization, zero-point energy, and tunneling. Many of these shape metrics consider just the instantaneous structure of the ring polymers. However, given the significant interest in methods such as centroid molecular dynamics and ring polymer molecular dynamics that link the molecular dynamics of these ring polymers to real time properties, there exists significant opportunity to exploit metrics that also allow for the study of the fluctuations of the atom delocalization in time. Here we consider the ring polymer delocalization from the perspective of computational topology, specifically persistent homology, which describes the 3-dimensional arrangement of point cloud data, (i.e. atomic positions). We employ the Betti sequence probability distribution to define the ensemble of shapes adopted by the ring polymer. The Wasserstein distances of Betti sequences adjacent in time are used to characterize fluctuations in shape, where the Fourier transform and associated principal components provides added information differentiating atoms with different NQEs based on their dynamic properties. We demonstrate this methodology on two representative systems, a glassy system consisting of two atom types with dramatically different de Broglie thermal wavelengths, and ab initio molecular dynamics simulation of an aqueous 4 M HCl solution where the H-atoms are differentiated based on their participation in proton transfer reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗