Results for “proxy application” — Page 10

DOE OSTI2015

Thread-level parallelization and optimization of NWChem for the Intel MIC architecture

In the multicore era it was possible to exploit the increase in on-chip parallelism by simply running multiple MPI processes per chip. Unfortunately, manycore processors' greatly increased thread- and data-level parallelism coupled with a reduced memory capacity demand an altogether different approach. In this paper we explore augmenting two NWChem modules, triples correction of the CCSD(T) and Fock matrix construction, with OpenMP in order that they might run efficiently on future manycore architectures. As the next NERSC machine will be a self-hosted Intel MIC (Xeon Phi) based supercomputer, we leverage an existing MIC testbed at NERSC to evaluate our experiments. In order to proxy the fact that future MIC machines will not have a host processor, we run all of our experiments in native mode. We found that while straightforward application of OpenMP to the deep loop nests associated with the tensor contractions of CCSD(T) was sufficient in attaining high performance, significant e ort was required to safely and efeciently thread the TEXAS integral package when constructing the Fock matrix. Ultimately, our new MPI+OpenMP hybrid implementations attain up to 65× better performance for the triples part of the CCSD(T) due in large part to the fact that the limited on-card memory limits the existing MPI implementation to a single process per card. Additionally, we obtain up to 1.6× better performance on Fock matrix constructions when compared with the best MPI implementations running multiple processes per card.

CCSD(T)↗

DOE OSTI2021

Time to Use Dendrohydrological Data in Water Resources Management?

Despite a wide consensus on the potential importance of streamflow reconstructions to water resources planning and management activities, applications of dendrohydrology have been piecemeal. This review article describes barriers to widespread use of dendrohydrological data and suggests research needs for operationalizing streamflow records reconstructed from these proxy data.

42 ENGINEERING↗

DOE OSTI2023

L1 0 Ordering in MnAl and FeNi Influenced by Magnetic Field and Strain

Due to various materials supply chain challenges, magnets free of constrained elements are attracting increasing interest. Magnetic materials such as rare-earth free FeNi and MnAl have been receiving considerable attention due to the high magnetocrystalline anisotropy and other associated magnetic properties derived from their unique chemically ordered tetragonal crystal structure, denoted as the L1 0 structure. However, synthesis of L1 0 FeNi has had limited success due to the extremely low atomic mobilities of Fe and Ni. In this work, isostructural MnAl was first studied as proxy to understand the L1 0 ordering process. Here, evidence of L1 0 ordering in FeNi derived from TEM studies is presented, where ordering was facilitated by the application of strain and magnetic field provided during thermal treatment of a severely plastically deformed FeNi alloy.

36 MATERIALS SCIENCE↗

DOE OSTI2014

Thread-Level Parallelization and Optimization of NWChem for the Intel MIC Architecture

In the multicore era it was possible to exploit the increase in on-chip parallelism by simply running multiple MPI processes per chip. Unfortunately, manycore processors' greatly increased thread- and data-level parallelism coupled with a reduced memory capacity demand an altogether different approach. In this paper we explore augmenting two NWChem modules, triples correction of the CCSD(T) and Fock matrix construction, with OpenMP in order that they might run efficiently on future manycore architectures. As the next NERSC machine will be a self-hosted Intel MIC (Xeon Phi) based supercomputer, we leverage an existing MIC testbed at NERSC to evaluate our experiments. In order to proxy the fact that future MIC machines will not have a host processor, we run all of our experiments in tt native mode. We found that while straightforward application of OpenMP to the deep loop nests associated with the tensor contractions of CCSD(T) was sufficient in attaining high performance, significant effort was required to safely and efficiently thread the TEXAS integral package when constructing the Fock matrix. Ultimately, our new MPI OpenMP hybrid implementations attain up to 65x better performance for the triples part of the CCSD(T) due in large part to the fact that the limited on-card memory limits the existing MPI implementation to a single process per card. Additionally, we obtain up to 1.6x better performance on Fock matrix constructions when compared with the best MPI implementations running multiple processes per card.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Challenges in quantifying unparameterized spatial uncertainties in deep geologic repositories for nuclear waste

Spatially heterogeneous uncertainties are prevalent in geophysical modeling applications, such as probabilistic post-closure performance assessment (PA) of deep geologic repositories for nuclear waste. Such uncertainties are often highly influential to model outputs, so it is desirable to identify the most important mechanisms by which they influence model predictions. However, these uncertainties are often unparameterized in the sense that there is no set of parameters that can be specified to yield a particular realization of the uncertainty. Additionally, the uncertainty is not intrinsically endowed with a parameterization that captures a realization’s mechanistic influence on model outputs. Therefore, in this work we present a novel methodology to develop and assess a set of proxy variables that aim to represent this influence. We show how they can be computed, downselected, and incorporated into the construction of statistical surrogate models mapping model inputs to outputs. We present our methodology in the context of a motivating application problem in deep geologic repository PA and discuss the challenges in capturing the effects of these spatial heterogeneities in uncertainty analyses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2026

First Demonstration of Resonant Pitch-Angle Scattering of Relativistic Electrons by Externally Launched Helicon Waves

Helicon waves (a.k.a whistler waves) satisfying the normal wave-particle cyclotron resonance are observed to limit the growth and maximum energy of relativistic electrons (REs) in low-density Ohmic DIII-D tokamak plasmas. Following the application of helicon waves, pitch-angle scattering of high-energy REs causes an increase in both synchrotron and electron-cyclotron emissions. The hard x-ray emission, a proxy for the RE population, ceases to grow. Energy-resolved hard x-ray measurements also show a striking decrease in the number of high-energy REs (above the resonance at approximately 8MeV) to below the noise floor and an increase in low-energy (∼ 4 MeV) REs. This occurs despite the toroidal electric field remaining high enough to drive exponential RE growth in the absence of helicon waves. Furthermore, these results open new directions for limiting the maximum energy of RE populations in laboratory and fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Validation of (not-historical) large-event near-fault ground-motion simulations for use in civil engineering applications

Ground-motion simulations generated from physics-based wave propagation models are gaining increasing interest in the engineering community for their potential to inform the performance-based design and assessment of infrastructure residing in active seismic areas. A key prerequisite before the ground-motion simulations can be used with confidence for application in engineering domains is their comprehensive and rigorous investigation and validation. This article provides a four-step methodology and acceptance criteria to assess the reliability of simulated ground motions of not historical events, which includes (1) the selection of a population of real records consistent with the simulated scenarios, (2) the comparison of the distribution of Intensity Measures (IMs) from the simulated records, real records, and Ground-Motion Prediction Equations (GMPEs), (3) the comparison of the distribution of simple proxies for building response, and (4) the comparison of the distribution of Engineering Demand Parameters (EDPs) for a realistic model of a structure. Specific focus is laid on near-field ground motions (<10km) from large earthquakes (M w 7), for which the database of real records for potential use in engineering applications is severely limited. The methodology is demonstrated through comparison of (2490) near-field synthetic records with 5 Hz resolution generated from the Pitarka et al (2019) kinematic rupture model with a population of (38) pulse-like near-field real records from multiple events and, when applicable, with NGA-W2 GMPEs. Lastly, the proposed procedure provides an effective method for informing and advancing the science needed to generate realistic ground-motion simulations, and for building confidence in their use in engineering domains.

42 ENGINEERING↗

DOE OSTI2020

μ⋆ masses: weak-lensing calibration of the Dark Energy Survey Year 1 redMaPPer clusters using stellar masses

ABSTRACT We present the weak-lensing mass calibration of the stellar-mass-based μ⋆ mass proxy for redMaPPer galaxy clusters in the Dark Energy Survey Year 1. For the first time, we are able to perform a calibration of μ⋆ at high redshifts, z > 0.33. In a blinded analysis, we use ∼6000 clusters split into 12 subsets spanning the ranges 0.1 ≤ z < 0.65 and μ⋆ up to ${\sim} 5.5 \times 10^{13} \, \mathrm{M}_{\odot }$, and infer the average masses of these subsets through modelling of their stacked weak-lensing signal. In our model, we account for the following sources of systematic uncertainty: shear measurement and photometric redshift errors, miscentring, cluster-member contamination of the source sample, deviations from the Navarro–Frenk–White halo profile, halo triaxiality, and projection effects. We use the inferred masses to estimate the joint mass–μ⋆–z scaling relation given by $\langle M_{200c} | \mu _{\star },z \rangle = M_0 (\mu _{\star }/5.16\times 10^{12} \, \mathrm{M_{\odot }})^{F_{\mu _{\star }}} ((1+z)/1.35)^{G_z}$. We find $M_0= (1.14 \pm 0.07) \times 10^{14} \, \mathrm{M_{\odot }}$ with $F_{\mu _{\star }}= 0.76 \pm 0.06$ and Gz = −1.14 ± 0.37. We discuss the use of μ⋆ as a complementary mass proxy to the well-studied richness λ for: (i) exploring the regimes of low z, λ < 20 and high λ, z ∼ 1; and (ii) testing systematics such as projection effects for applications in cluster cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Extraction and biomolecular analysis of dermal interstitial fluid collected with hollow microneedles

Dermal interstitial fluid (ISF) is an underutilized information-rich biofluid potentially useful in health status monitoring applications whose contents remain challenging to characterize. Here, we present a facile microneedle approach for dermal ISF extraction with minimal pain and no blistering for human subjects and rats. Extracted ISF volumes were sufficient for determining transcriptome, and proteome signatures. We noted similar profiles in ISF, serum, and plasma samples, suggesting that ISF can be a proxy for direct blood sampling. Dynamic changes in RNA-seq were recorded in ISF from induced hypoxia conditions. Finally, we report the first isolation and characterization, to our knowledge, of exosomes from dermal ISF. The ISF exosome concentration is 12–13 times more enriched when compared to plasma and serum and represents a previously unexplored biofluid for exosome isolation. This minimally invasive extraction approach can enable mechanistic studies of ISF and demonstrates the potential of ISF for real-time health monitoring applications.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2017

Decomposing Fuel Economy and Greenhouse Gas Regulatory Standards in the Energy Conversion Efficiency and Tractive Energy Domain

The three foundational elements that determine mobile source energy use and tailpipe carbon dioxide (CO2) emissions are the tractive energy requirements of the vehicle, the on-cycle energy conversion efficiency of the propulsion system, and the energy source. The tractive energy requirements are determined by the vehicle's mass, aerodynamic drag, tire rolling resistance, and parasitic drag. Oncycle energy conversion of the propulsion system is dictated by the tractive efficiency, non-tractive energy use, kinetic energy recovery, and parasitic losses. The energy source determines the mobile source CO2 emissions. For current vehicles, tractive energy requirements and overall energy conversion efficiency are readily available from the decomposition of test data. For future applications, plausible levels of mass reduction, aerodynamic drag improvements, and tire rolling resistance can be transposed into the tractive energy domain. Similarly, by combining thermodynamic, mechanical efficiency, and kinetic energy recovery fundamentals with logical proxies, achievable levels of energy conversion efficiency can be established to allow for the evaluation of future powertrain requirements. Combining the plausible levels of tractive energy and on-cycle efficiency provides a means to compute sustainable vehicle and propulsion system scenarios that can achieve future regulations. Using these principles, the regulations established in the United States (U.S.) for fuel consumption and CO2 emissions are evaluated. Fleet-level scenarios are generated and compared to the technology deployment assumptions made during rule-making. When compared to the rule-making assumptions, the results indicate that a greater level of advanced vehicle and propulsion system technology deployment will be required to achieve the model year 2025 U.S. standards for fuel economy and CO2 emissions.

CAFE↗

DOE OSTI2022

Tailoring polyvinyl alcohol-sodium alginate (PVA-SA) hydrogel beads by controlling crosslinking pH and time

Abstract Hydrogel-encapsulated catalysts are an attractive tool for low-cost intensification of (bio)-processes. Polyvinyl alcohol-sodium alginate hydrogels crosslinked with boric acid and post-cured with sulfate (PVA-SA-BS) have been applied in bioproduction and water treatment processes, but the low pH required for crosslinking may negatively affect biocatalyst functionality. Here, we investigate how crosslinking pH (3, 4, and 5) and time (1, 2, and 8 h) affect the physicochemical, elastic, and process properties of PVA-SA-BS beads. Overall, bead properties were most affected by crosslinking pH. Beads produced at pH 3 and 4 were smaller and contained larger internal cavities, while optical coherence tomography suggested polymer cross-linking density was higher. Optical coherence elastography revealed PVA-SA-BS beads produced at pH 3 and 4 were stiffer than pH 5 beads. Dextran Blue release showed that pH 3-produced beads enabled higher diffusion rates and were more porous. Last, over a 28-day incubation, pH 3 and 4 beads lost more microspheres (as cell proxies) than beads produced at pH 5, while the latter released more polymer material. Overall, this study provides a path forward to tailor PVA-SA-BS hydrogel bead properties towards a broad range of applications, such as chemical, enzymatic, and microbially catalyzed (bio)-processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Transforming microseismic clouds into near real-time visualization of the growing hydraulic fracture

SUMMARY Microseismic observations during unconventional reservoir stimulation are typically seen as a proxy for clusters of hydraulic fractures and the extent of the stimulated reservoir. Such straightforward interpretation is often misleading and fails to provide a physically reasonable image of the fracturing process. This paper demonstrates the application of a physics-based machine learning algorithm which enables a rapid and accurate fracture mapping from the microseismic data. Our training and validation data set relies on a history-matched geomechanical modelling workflow implemented in GEOS software for the Hydraulic Fracturing Test Site 1 (HFTS-1) project. For this study we augmented the simulated fracture growth through geostatistical modelling of induced seismicity, so that the synthetic microseismic catalogue matches the main statistical properties of the field observations. We formulated the problem of mapping the actual fracture in the clutter of events to parallel common video segmentation workflows: several past video frames (microseismic density snapshots) are passed through a deep convolutional network to classify whether a given voxel is associated with a fracture or intact rock. We found that for accurate fracture mapping, the network’s input and architecture must be augmented to incorporate the fluid injection parameters (pressure, rate, concentration of proppant, and location of the perforation within the cluster). The error rate for the network reached as little as 10 per cent of the fracture area, while a conventional microseismic interpretation approach yielded ∼300 per cent. Our approach also yields must faster predictions than conventional methods (minutes instead of weeks), and could enable engineers to make rapid decisions regarding engineering parameters (pumping rate, viscosity) in real time during stimulation.

58 GEOSCIENCES↗

DOE OSTI2026

Machine learning enables reconstruction of past fire regimes from charcoal-derived fire intensity and fuel composition

Background Fire is a foundational ecological process that shapes ecosystem structure, diversity, and resilience. Quantifying paleofire regime attributes such as frequency, severity, and intensity is essential for understanding the historical range of variability in fire behavior and its ecological effects. While frequency and severity are often reconstructed in paleofire studies, quantitative reconstructions of fire intensity remain limited. Recent work has shown that maximum pyrolysis temperature—a proxy for fire intensity—and plant species type can be inferred from charcoal using transmission Fourier-transform infrared (FTIR) spectroscopy. However, the sample preparation for transmission FTIR is destructive and time-consuming, limiting application and reuse of materials for other analyses. We evaluated reflectance FTIR spectroscopy as a non-destructive alternative for reconstructing combustion temperature and plant species from laboratory-generated charcoal. We also examined the influence of contrasting airflow environments (ambient air versus nitrogen-rich) on pyrolysis temperature and plant species reconstruction prediction accuracies and compared predictive performance between a novel, neural network–based deep learning model with the traditional modern analogue technique (MAT) using k-nearest neighbor functions. As proof of concept, we apply our enhanced methodology to ancient charcoal to demonstrate applicability at improving long-term fire regime reconstructions and the ability to link paleofire records with contemporary fire ecology. Results Our analysis shows that transmission and reflectance FTIR spectra yield comparable spectral profiles. However, sample preparation for reflectance FTIR is minimal and non-destructive, unlike transmission FTIR which is destructive. We demonstrate that oxygen environments improved reconstruction accuracy relative to nitrogen-rich conditions. Finally, our deep learning neural network (DL) achieved testing accuracies of 98.7% for temperature and 96.2% for species identification, outperforming MAT’s k-NN approach (89.8% and 65.9%, respectively). A Shapley importance analysis identified 5 key spectral regions that greatly influenced the model’s temperature or species categorization. When applied to ancient charcoal, our results show historic fires from the most recent past primarily burned at low intensities (400–500 °C), reflective of natural fire regimes in ponderosa pine forests. Our results corroborate charcoal morphology data that suggests all ancient charcoal originated from burned woody plant types. Conclusions By combining reflectance FTIR spectroscopy with a deep learning approach, we provide the first accuracies high enough to confidently identify both species and temperature from laboratory-produced charcoal, improving quantitative reconstructions of fire intensity and fuel composition from paleofire records. This opens a wide range of research into the link between fire and larger drivers (i.e., climate or human) and greater ecological understanding of fire regimes beyond that of burn scars or recent observations. These methodological improvements have direct relevance for fire management by improving interpretation of historical fire behavior, informing fuel–fire relationships, and providing a scalable analytical framework applicable to both long-term ecological studies and contemporary fire science.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2018

Retrievals of Riming and Snow Density From Vertically Pointing Doppler Radars

Retrievals of ice and snow are made from Ka- and W-band zenith-pointing Doppler radars at Hyytiälä, Finland, during the snow experiment component of the Biogenic Aerosols: Effects on Clouds and Climate (2014) field campaign. In a novel optimal estimation retrieval, mean Doppler velocity is exploited to retrieve a density factor parameter, which modulates the mass, shape, terminal velocity, and backscatter cross sections of ice particles. In a case study including aggregate snow and graupel we find that snow rate and ensemble mean ice density can be retrieved to within 50% of in situ measurements at the surface using dual-frequency Doppler radar retrievals. While Doppler measurements are essential to the retrieval of particle density, the dual-frequency ratio provides a strong constraint on particle size. The retrieved density factor is strongly correlated with liquid water path, indicating that riming is the primary process by which the density factor is modulated. Using liquid water path as a proxy for riming, profiles classified as unrimed snow, rimed snow, and graupel exhibit distinct features characteristic of aggregation and riming processes, suggesting the potential to make estimates of process rates from these retrievals. We discuss the potential application of the technique to future satellite missions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Heat transport in liquid water from first-principles and deep neural network simulations

In this work, we compute the thermal conductivity of water within linear response theory from equilibrium molecular dynamics simulations, by adopting two different approaches. In one, the potential energy surface (PES) is derived on the fly from the electronic ground state of density functional theory (DFT) and the corresponding analytical expression is used for the energy flux. In the other, the PES is represented by a deep neural network (DNN) trained on DFT data, whereby the PES has an explicit local decomposition and the energy flux takes a particularly simple expression. By virtue of a gauge invariance principle, established by Marcolongo, Umari, and Baroni, the two approaches should be equivalent if the PES were reproduced accurately by the DNN model. We test this hypothesis by calculating the thermal conductivity, at the GGA (PBE) level of theory, using the direct formulation and its DNN proxy, finding that both approaches yield the same conductivity, in excess of the experimental value by approximately 60%. Besides being numerically much more efficient than its direct DFT counterpart, the DNN scheme has the advantage of being easily applicable to more sophisticated DFT approximations, such as meta-GGA and hybrid functionals, for which it would be hard to derive analytically the expression of the energy flux. We find in this way that a DNN model, trained on meta-GGA (SCAN) data, reduces the deviation from experiment of the predicted thermal conductivity by about 50%, leaving the question open as to whether the residual error is due to deficiencies of the functional, to a neglect of nuclear quantum effects in the atomic dynamics, or, likely, to a combination of the two.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Recent Developments in Hydrogeologic Applications for Strain Tensor Analyses

Changes in fluid pressure deform porous media and this effect occurs in a variety of hydrogeologic processes, from the change in storage during pumping or injection to fluctuations in water levels caused by barometric pressure. We have developed instruments for measuring small strains in porous media, and we have used the resulting strain data to evaluate well testing, hydraulic fracturing, manual loading at the ground surface, and ambient hydrologic processes, like rainfall and evaporation. A particularly important application is the use of strain tensor data measured at shallow depths to analyze well tests or hydraulic fractures conducted at much greater depths. An early demonstration of this technique was conducted at the North Avant Field north of Tulsa, Oklahoma, where Pennsylvanian sandstone creates a confined aquifer and oil reservoir at a depth of 530m. We have showed that the strains caused by injecting into the aquifer could be measured at a depth of 30m and used to evaluate the properties of the aquifer. We recently expanded the array of strainmeters at the North Avant Field by deploying three more instruments at shallow depth (30m) along with a deep instrument at 520m depth in the winter, 2021. The deep instrument is deployed in shale caprock slightly above the aquifer. To our knowledge, the deep strainmeter at the North Avant Field is the deepest strainmeter ever deployed and this required refining methods originally developed for shallow deployments. The instrument was lowered to depth on oil field tubing and cemented in place using techniques and materials developed for use in oil wells. Optical fiber used to communicate with the instrument was cemented in the annulus on the outside of tubing. This is significant because the techniques we used could readily be extended to greater depth, suggesting that strainmeters can be deployed over a wide range of depths for monitoring critical subsurface processes. For example, it suggests that strainmeters could be deployed through the caprock to monitor for leaks from underlying CO2 storage reservoirs. The strainmeter array at the North Avant Field has been used to characterize deformation during a series of injection tests in the spring and summer, 2021. All the new strainmeters respond to pumping, and the strainmeters we installed earlier also responded. To our knowledge, this is the first application of well testing in a deep aquifer that was monitored by an array of strainmeters—our earlier work used strainmeters at a single location. We are currently analyzing the strain data using an analytical solution, a proxy-based Bayesian inversion algorithm, and other methods. Strainmeter data has also been used to characterize periodic pumping tests by us and Riley Blais. A periodically varying pumping rate causes both hydraulic head and strain signals that vary with the same period as the pumping. The peaks and troughs of the head in monitoring wells lag behind the peaks and troughs of the head in the pumping well, and this lag time increases with distance from the pumping well. The lag time of the pressure and the distance to the monitoring well can be used in a simple analysis to estimate the hydraulic diffusivity of the aquifer. The lag time determined from strain data can be used to estimate aquifer properties using the same analysis that works for the pressure only for strain data measured at particular locations. That is because the strain field in a confining unit advances upward, laterally and then downward even though the pressure in the underlying aquifer only advances laterally, according to our recent simulations. We have field data showing that a small periodic signal superimposed on an injection rate at the North Avant Field will create a periodic strain signal at shallow strainmeters. The field data and the recent simulations suggest that including a periodic component to injection or pumping and then monitoring the resulting strain signal could be a way to monitor the subsurface.

Murdoch, Larry↗

DOE OSTI2023

Recent Developments in Hydrogeologic Applications for Strain Tensor Analyses

Changes in fluid pressure deform porous media and this effect occurs in a variety of hydrogeologic processes, from the change in storage during pumping or injection to fluctuations in water levels caused by barometric pressure. We have developed instruments for measuring small strains in porous media, and we have used the resulting strain data to evaluate well testing, hydraulic fracturing, manual loading at the ground surface, and ambient hydrologic processes, like rainfall and evaporation. A particularly important application is the use of strain tensor data measured at shallow depths to analyze well tests or hydraulic fractures conducted at much greater depths. An early demonstration of this technique was conducted at the North Avant Field north of Tulsa, Oklahoma, where Pennsylvanian sandstone creates a confined aquifer and oil reservoir at a depth of 530m. We have showed that the strains caused by injecting into the aquifer could be measured at a depth of 30m and used to evaluate the properties of the aquifer. We recently expanded the array of strainmeters at the North Avant Field by deploying three more instruments at shallow depth (30m) along with a deep instrument at 520m depth in the winter, 2021. The deep instrument is deployed in shale caprock slightly above the aquifer. To our knowledge, the deep strainmeter at the North Avant Field is the deepest strainmeter ever deployed and this required refining methods originally developed for shallow deployments. The instrument was lowered to depth on oil field tubing and cemented in place using techniques and materials developed for use in oil wells. Optical fiber used to communicate with the instrument was cemented in the annulus on the outside of tubing. This is significant because the techniques we used could readily be extended to greater depth, suggesting that strainmeters can be deployed over a wide range of depths for monitoring critical subsurface processes. For example, it suggests that strainmeters could be deployed through the caprock to monitor for leaks from underlying CO2 storage reservoirs. The strainmeter array at the North Avant Field has been used to characterize deformation during a series of injection tests in the spring and summer, 2021. All the new strainmeters respond to pumping, and the strainmeters we installed earlier also responded. To our knowledge, this is the first application of well testing in a deep aquifer that was monitored by an array of strainmeters—our earlier work used strainmeters at a single location. We are currently analyzing the strain data using an analytical solution, a proxy-based Bayesian inversion algorithm, and other methods. Strainmeter data has also been used to characterize periodic pumping tests by us and Riley Blais. A periodically varying pumping rate causes both hydraulic head and strain signals that vary with the same period as the pumping. The peaks and troughs of the head in monitoring wells lag behind the peaks and troughs of the head in the pumping well, and this lag time increases with distance from the pumping well. The lag time of the pressure and the distance to the monitoring well can be used in a simple analysis to estimate the hydraulic diffusivity of the aquifer. The lag time determined from strain data can be used to estimate aquifer properties using the same analysis that works for the pressure only for strain data measured at particular locations. That is because the strain field in a confining unit advances upward, laterally and then downward even though the pressure in the underlying aquifer only advances laterally, according to our recent simulations. We have field data showing that a small periodic signal superimposed on an injection rate at the North Avant Field will create a periodic strain signal at shallow strainmeters. The field data and the recent simulations suggest that including a periodic component to injection or pumping and then monitoring the resulting strain signal could be a way to monitor the subsurface.

Murdoch, Larry↗

DOE OSTI2025

Enabling hydrogen chemisorption on charged graphene

Two-dimensional (2D) materials, including graphitic derivatives, have long been of interest for hydrogen storage applications, due to their high theoretical storage capacity, low weight, and other useful properties. However, poor kinetics for hydrogen adsorption and surface diffusion as part of the proposed spillover process for hydrogenation have limited their technological potential. Here, we use first-principles calculations to study electronic doping as a means to improve hydrogen chemisorption on graphene, which we use here as a proxy for graphitic derivatives more broadly. We find that positively charged graphene sheets have vastly improved kinetics for hydrogen diffusion and adsorption, while they limit unwanted hydrogen desorption. This combination of effects should favor hydrogen chemisorption via spillover. We connect these trends to the C–H bond, which introduces states near the Fermi level. These states are depopulated as electrons are removed, thereby lowering the bond energy and permitting more facile movement of hydrogen. Our results suggest that spillover mechanisms for hydrogen chemisorption should be revisited if strategies to apply a large charge to graphitic systems can be realized. Moreover, switchable application of the charge may lead to the reversible chemisorption of hydrogen. While the large magnitude of charging required suggests that graphene itself may not be suitable for reversible hydrogen chemisorption, the factors we identify and discuss could significantly boost the prospects of graphitic derivatives and other 2D or layered materials for hydrogen storage applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗