SEARCH · Search NASA

Results for “ANALOG SIMULATION”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Shape Anisotropy-Dependent Leaking in Magnetic Neurons for Bio-Mimetic Neuromorphic Computing

Spiking neural networks seek to emulate biological computation through interconnected artificial neuron and synapse devices. Spintronic neurons can leverage magnetization physics to mimic biological neuron functions, such as integration tied to magnetic domain wall (DW) propagation in a patterned nanotrack and firing tied to the resistance change of a magnetic tunnel junction (MTJ), captured in the domain wall-magnetic tunnel junction (DW-MTJ) device. Leaking, relaxation of a neuron when it is not under stimulation, is also predicted to be implemented based on DW drift as a DW relaxes to a low energy position, but it has not been well explored or demonstrated in device prototypes. Here, in this work, we study DW-MTJ artificial neurons capable of leaky integrate-and-fire (LIF) behavior and demonstrate geometry-dependent leaking dynamics that results in repeatable, tunable LIF operation. Studying the behavior of five different device designs, we show tuning the geometry, stimulating fields and currents, and location of electrical contacts results in a wide range of neuron behavior. Additionally, implementation of an asymmetric notch allows for nonlinear pinning which increased expressivity without sacrificing leaking. The measured behavior is implemented in a simulated spiking neural network that outperforms a 1D model of continuous DW motion and approaches the performance of an ideal LIF activation function. The results show that the analog LIF capability of DW-MTJ neurons combines many desirable neuron functions into a single device, which can result in varied forms of multifunctional neuromorphic computing.

42 ENGINEERING↗

DOE OSTI2025

Analcime-wairakite formation during experimental cement-bentonite alteration at 200–300 °C

Engineered barrier system materials in nuclear waste repositories may undergo hydrothermal alteration in response to groundwater saturation and heating events over their long operational timescales. Hydrothermal interactions between engineered materials (e.g., bentonite buffers, cements, waste canister materials) and the host rock environment may drive alteration processes that affect the advantageous properties of some barrier materials. However, such alteration will also promote the formation of zeolites, hydrothermally formed minerals that may themselves act to isolate radionuclides. Understanding the environmental conditions that control radionuclide-sorbing properties of zeolites, such as Si/Al ratio or Na content, will be valuable for assessing the changing properties of engineered barrier materials in the case of an in situ heating event. Here, we present experimental work characterizing the formation of zeolites during hydrothermal interactions between generic clay and cement barrier materials, with a focus on the analcime-wairakite zeolite series, which has known radionuclide sorption and exchange properties. We present the results of hydrothermal experiments combining uncured ordinary Portland cement powder with steel and Wyoming bentonite in the presence of Opalinus clay and a synthesized Opalinus clay groundwater to simulate water-saturated conditions in an argillaceous rock repository. The experiments were conducted isothermally at 200 or 300 °C for 8 to 24 weeks. In contrast to analogous studies that did not include cement reactants, we observed the formation of analcime-group minerals in all experiments. The addition of cement resulted in lower Si/Al ratios in the zeolites compared to similar studies that did not include cement. The presence of cement as a reactant was interpreted to promote analcime formation at 200 °C. At 300 °C, we observed higher calcium and silica concentrations in the aqueous solutions as well as increased wairakite formation and decreased analcime formation compared to the experiments at 200 °C. These results show a fully realized analcime-wairakite solid solution that falls between Si/Al = 2 (ideal) and a trend of analcime minerals that have increasing Si/Al ratios with increasing Na/(Na+Ca). Our results show predictable relationships between the Si/Al ratio and the analcime-wairakite content under hydrothermal conditions and illustrate that repository material interactions may promote the formation of zeolites in the analcime-wairakite solid-solution series during heating events in the subsurface.

58 GEOSCIENCES↗

DOE OSTI2025

Electrochemically Mediated Au–C(sp 2 ) Anchors for Molecular Electronics

Terminal anchor groups play a key role in the stability and electronic properties of molecular junctions. Single molecule junctions typically consist of two preinstalled terminal anchors linking organic molecules to metal electrodes. Here, in this work, we show that p -terphenyl derivatives containing only a single terminal anchor show conductance features similar to junctions with two preinstalled terminal anchors. A set of p -terphenyl derivatives with one terminal anchor was prepared using automated chemical synthesis and characterized using single molecule electronics experiments, molecular dynamics (MD) simulations, bulk electrochemistry and spectroscopy, and nonequilibrium Green’s function-density functional theory (NEGF-DFT) calculations. Our results show that 4-amino- p -terphenyl (PPP) and related analogs exhibit a well-defined high conductance state that is diminished or absent in other p-terphenyl derivatives lacking a preinstalled amine terminal anchor or fluorine or methyl substitutions at the terminal para position. However, a low conductance state is observed in all amino- p -terphenyl derivatives with one preinstalled anchor due to molecular junctions formed by noncovalent dimeric π–π stacking interactions. The observed high conductance state diminishes upon the addition of reducing agents and is restored upon the addition of an oxidizing agent. Our results suggest that the high conductance state arises due to Au–C(sp 2 ) bond formation facilitated by a single electron oxidation event at the electrode surface. A series of control experiments with different anchor groups shows that primary amines play a key role in forming Au–C bonds for molecular junctions. Overall, these results suggest that Au–C bond formation gives rise to high conductance pathways in organic molecules containing only one preinstalled terminal anchor. Insights from this work can be leveraged in the design of molecular electronic devices, particularly in understanding the mechanisms of molecular binding and junction formation.

charge transport↗

DOE OSTI2025

Collisionless conduction in a high-beta plasma: a collision operator for whistler turbulence

The regulation of electron heat transport in high-β, weakly collisional, magnetized plasma is investigated. A temperature gradient oriented along a mean magnetic field can induce a kinetic heat-flux-driven whistler instability (HWI), which back-reacts on the transport by scattering electrons and impeding their flow. Previous analytical and numerical studies have shown that the heat flux for the saturated HWI scales as β$^{-1}_{e}$. These numerical studies, however, had limited scale separation and consequently large fluctuation amplitudes, which calls into question their relevance at astrophysical scales. To this end, we perform a series of particle-in-cell simulations of the HWI across a range of β$_e$ and temperature-gradient length scales under two different physical set-ups. The saturated heat flux in all of our simulations follows the expected β$^{-1}_{e}$ scaling, supporting the robustness of the result. We also use our simulation results to develop and implement several methods to construct an effective collision operator for whistler turbulence. The results point to an issue with the standard quasi-linear explanation of HWI saturation, which is analogous to the well-known 90° scattering problem in the cosmic-ray community. Despite this limitation, the methods developed here can serve as a blueprint for future work seeking to characterize the effective collisionality caused by kinetic instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Porous and Amorphous Mn x Mo 3 S 13 Chalcogel Electrode for High-Capacity Conversion-Based Lithium-Ion Batteries

While Li-ion batteries (LIBs) are a leading energy storage technology, their energy densities are limited by the low capacity of conventional intercalation cathodes, driving interest in high energy-density Li–S batteries that make use of conversion chemistry. Achieving high capacity, reversibility, and cycle stability, and controlling volume changes in conversion batteries during the charge–discharge process, however, remains challenging. Here, in this study, we present a porous, amorphous, sulfide-based Mn x Mo 3 S 13 chalcogel, which concurrently offers high capacity and cycle stability. The solution-processable room temperature synthesized Mn x Mo 3 S 13 (x = 0.25) chalcogel exhibits a local structure that resembles the Mo 3 S 13 cluster with Mn 2+ distributed across the Mo 3 S 13 matrix, as determined by synchrotron X-ray pair distribution function (PDF) and extended X-ray absorption fine structure (EXAFS). Ab initio molecular dynamics (AIMD) simulations reveal that Mn 2+ incorporation shortens the polysulfide chain in the gel matrix compared to the Mo 3 S 13 chalcogel, while forming a coordination environment with disulfide groups, analogous to the experimental findings. A Li/Mn 0.25 Mo 3 S 13 half-cell delivers 897 mAh g –1 capacity during the first discharge and retains 571 mAh g –1 capacity after 100 cycles at a C/3 rate. Distribution of relaxation time (DRT) unveils a stable solid–electrolyte interphase (SEI) formation upon cycling that enables charge–discharge reversibility. Here, the enhanced capacity retention and cycle stability compared to those of the Li/Mo 3 S 13 cell are attributed to the reduced dissolution of active mass into the electrolyte, facilitated by the formation of shorter polysulfide chains within the Mn 0.25 Mo 3 S 13 structure and the strong affinity of Lewis-acidic Mn 2+ for polysulfide anions generated during the charge–discharge process of the Li/Mn 0.25 Mo 3 S 13 cell. Thus, this work illustrates a design principle of material for high-capacity and cycle-stable Li-metal sulfide batteries.

25 ENERGY STORAGE↗

DOE OSTI2026

Algorithm to extract direction in 2D discrete distributions and a continuous Frobenius norm

In this study, we present a novel algorithm for determining directionality in 2D distributions of discrete data. We compare a reference dataset with a known direction to a measured dataset with an unknown direction by the Frobenius norm of the difference (FND) to find the unknown direction. To generalize this concept, we develop a continuous Frobenius norm of the difference (CFND) as a continuous analog of the FND and derive its analytical expression. By relating fitted and normalized 2D Gaussian distributions, we show that the CFND approximates the FND, and we validate this relationship with computer simulations. We find that a first-order approximation of the CFND between two similar Gaussian distributions takes the form of an absolute sine function, offering a simple analytical form with potential applications in specialized areas such as segmented inverse beta decay neutrino detectors, astronomy, machine learning, and more. Our methodology consists of modeling a 2D Gaussian distribution, binning the data into a histogram, and encoding it as a square matrix. Rotating this matrix around its geometric center and comparing it to a measured dataset using the FND gives us rotational data that we fit with an absolute sine function. The location of the minimum of this fit is the angle closest to the true angle of the direction in the measured dataset. We present the derivation and discuss initial applications of the CFND in our novel algorithm, demonstrating its success in approximating directionality in 2D distributions.

Physics↗

DOE OSTI2025

Algorithm to extract direction in 2D discrete distributions and a continuous Frobenius norm

In this study, we present a novel algorithm for determining directionality in 2D distributions of discrete data. We compare a reference dataset with a known direction to a measured dataset with an unknown direction by the Frobenius norm of the difference (FND) to find the unknown direction. To generalize this concept, we develop a continuous Frobenius norm of the difference (CFND) as a continuous analog of the FND and derive its analytical expression. By relating fitted and normalized 2D Gaussian distributions, we show that the CFND approximates the FND, and we validate this relationship with computer simulations. We find that a first-order approximation of the CFND between two similar Gaussian distributions takes the form of an absolute sine function, offering a simple analytical form with potential for specialized applications in segmented inverse beta decay (IBD) neutrino detectors, astronomy, machine learning, and more. Although this method may easily extend to 3D scalar fields, our focus here is on 2D real-valued fields as it directly applies to directionality. Our methodology consists of modeling a 2D Gaussian distribution, binning the data into a histogram, and encoding it as a square matrix. Rotating this matrix around its geometric center and comparing it to a measured dataset using the FND gives us rotational data that we fit with an absolute sine function. The location of the minimum of this fit is the angle closest to the true angle of the direction in the measured dataset. We present the derivation and discuss initial applications of the CFND in our novel algorithm, demonstrating its success in approximating directionality in 2D distributions.

Data Analysis, Statistics and Probability (physics↗

DOE OSTI2026

Hydrophobic Metal–Organic Frameworks Enable Superior High-Pressure Ammonia Storage through Geometric Design

Hydrophobic metal–organic frameworks (MOFs) are typically overlooked for ammonia storage due to weak host–guest interactions. Here, we demonstrate that four structurally analogous aluminum-based MOFs exhibit a counterintuitive behavior whereby framework geometry, rather than ligand hydrophilicity, determines high-pressure NH 3 adsorption performance. The hydrophobic CAU-23 achieved an exceptional capacity matching hydrophilic analogs despite its poor low-pressure uptake. This pressure-dependent enhancement stems from the unique 4-cis-4-trans geometry of CAU-23 compared to the purely cis arrangement of MIL-160 and KMF-1 and the alternating cis-trans configuration of MOF-303. Critically, CAU-23 retained 95% capacity over three high-pressure cycles, whereas hydrophilic MOFs suffered 39–46% irreversible losses due to strong NH 3 -framework interactions that compromise structural integrity. Grand canonical Monte Carlo simulations reveal that high pressure enables NH 3 clustering through intermolecular hydrogen bonding, bypassing the need for strong host–guest interactions. High-pressure powder X-ray diffraction measurements confirm the exceptional mechanical resilience of CAU-23, showing complete structural recovery upon decompression despite exhibiting the highest pressure sensitivity among the studied MOFs. An extended analog, HE-CAU-23, validates this design principle with further enhanced capacity. Furthermore, these findings reveal a paradigm shift toward hydrophobic MOFs with optimized geometry for high-performance and regenerable gas storage applications.

MOFs↗

DOE OSTI2025

On a Critical Acceleration Scale of Dark Matter in ΛCDM and Dynamical Dark Energy

Abstract Universal acceleration a 0 emerges in various empirical laws, yet its fundamental nature remains unclear. Using Illustris and Virgo N -body simulations, we focus on the velocity and acceleration fluctuations in collisionless dark matter involving long-range gravity. For comparison, in the kinetic theory of gases, molecules undergo random elastic collisions involving short-range interactions, where only velocity fluctuations are relevant. Hierarchical structure formation proceeds through the merging of smaller halos to form larger halos, which facilitates a continuous energy cascade from small to large halos at a constant rate ε u ≈ −10 −7 m 2 s −3 . Velocity fluctuations involve a critical velocity u c ∝ (1 + z ) −3/4 . Acceleration fluctuations involve a critical acceleration a c ∝ (1 + z ) 3/4 . Two critical quantities are related by the rate of energy cascade ε u ≈ − a c u c /[2(3 π ) 2 ], where factor 3 π is from the angle of incidence during merging. With critical velocity u c on the order of 300 km s −1 at z = 0, the critical acceleration is determined to be a c 0 ≡ a c ( z = 0) ≈ 10 −10 m s −2 , suggesting a c might explain the universal acceleration a 0 ≈ 10 −10 m s −2 in the empirical Tully–Fisher relation or modified Newtonian dynamics. The redshift evolution a c ∝ (1 + z ) 3/4 is in good agreement with Magneticum and EAGLE simulations and in reasonable agreement with limited observations. This suggests a larger a 0 at a higher redshift such that galaxies of fixed mass rotate faster at a higher redshift. Note that for dark energy (DE) density ρ DE 0 ≈ a c 0 2 / G = 1 0 − 10 J m −3 , we postulate an entropic origin of the DE from acceleration fluctuations of dark matter, analogous to the gas pressure from velocity fluctuations. This leads to a dynamical DE coupled to the structure evolution involving a relatively constant DE density followed by a slow weakening phase, suggesting possible deviations from the standard ΛCDM paradigm.

N-body simulations↗

DOE OSTI2025

SIDM Concerto: Compilation and Data Release of Self-interacting Dark Matter Zoom-in Simulations

We present SIDM Concerto: 14 cosmological zoom-in simulations in cold dark matter (CDM) and self-interacting dark matter (SIDM) models based on the Symphony and Milky Way-est suites. SIDM Concerto includes one Large Magellanic Cloud– (LMC-) mass system (host mass ∼10 11 M ⊙ ), two Milky Way (MW) analogs (∼10 12 M ⊙ ), two group-mass hosts (∼10 13 M ⊙ ), and one low-mass cluster (∼10 14 M ⊙ ). Each host contains ≈2 × 10 7 particles and is run in CDM and one or more strong, velocity-dependent SIDM models. Our analysis of SIDM (sub)halo populations over seven subhalo mass decades reveals that (1) the fraction of core-collapsed isolated halos and subhalos peaks at a maximum circular velocity corresponding to the transition of the SIDM cross section from a v −4 to v 0 scaling; (2) SIDM subhalo mass functions are suppressed by ≈50% relative to CDM in LMC, MW, and group-mass hosts but are consistent with CDM in the low-mass cluster host; (3) subhalos’ inner density profile slopes, which are more diverse in SIDM than in CDM, are sensitive to both the amplitude and shape of the SIDM cross section. Our simulations provide a benchmark for testing SIDM predictions with astrophysical observations of field and satellite galaxies, strong lensing systems, and stellar streams. Data products are publicly available at doi:10.5281/zenodo.14933624.

dark matter↗

DOE OSTI2026

Lieb-Mattis States for Robust Entangled Differential Phase Sensing

We explore a two-node, entanglement-enhanced sensor network for differential phase sensing that exploits decoherence-free subspaces to suppress common-mode noise, a primary limitation of many state-of-the-art quantum sensors. We identify a class of entangled states that, while not strictly optimal, achieve the same asymptotic sensitivity scaling as optimal states and can be prepared efficiently from initially unentangled atomic ensembles. Importantly, the preparation time decreases with increasing system size, which makes the states compatible with realistic noise processes in present-day quantum sensors that operate with large particle numbers but lack full error correction. We illustrate these ideas using two cavity-mediated preparation protocols: (i) coherent, unitary entanglement generation analogous to bosonic two-mode squeezing, yielding Heisenberg scaling, and (ii) dissipative preparation through collective emission into a shared cavity mode, providing a square-root improvement beyond the standard quantum limit. Numerical simulations show that both approaches remain effective at experimentally realistic cavity cooperativities, establishing a practical path toward scalable, quantum-enhanced differential phase sensing.

Kaubruegger, Raphael [National Institute of Standa↗

DOE OSTI2024

Computing with a Chemical Reservoir

Contemporary computation is expensive, with large language models and artificial intelligence becoming more common in daily life. However, high-performance computing is reaching the limits in speed and energy expenditure, and domain science requires ever-increasing computational capacity, with simulations and data analysis pipelines ever-growing in complexity. As we progress towards post-exascale computation, with the associated high energy costs, new methods of energy-conscious computation are required. Novel analog and hybrid digital-analog systems can overcome these challenges, and chemical reactions offer a promising avenue. Computers based on chemistry can provide compact desktop devices with immense computational power. These devices are readily scalable by considering greater reaction systems or vessels, meeting the high-performance requirements for scientific workflows. In this article, we present ChemComp, a compilation pipeline for the conversion of ordinary differential equations into implementable chemical reactions. We then demonstrate the solving capabilities of ChemComp by emulating a potential chemical reservoir device. We leverage the multi-layer intermediate representation (MLIR) compiler framework to implement an expressive chemical reaction abstraction and propose a path for chemical reaction networks (CRNs) to represent mathematical problems effectively. Combined, we demonstrate a potential workflow that can harness chemistry’s computing power to create energy-efficient, high-performance computation systems for contemporary computing needs.

artificial intelligence↗

DOE OSTI2025

Studying Open Quantum Systems Relevant to Chemistry on a Trapped-Ion Quantum Simulator (Final Technical Report)

This project advances the trapped-ion quantum simulator as a versatile platform for studying open quantum system phenomena. We aim to contribute to the emerging quantum simulation toolkits and enable simulation of nanoscale energy processes. Trapped-ion platforms offer unique capabilities: their vibrational motion can be precisely manipulated, measured, and coherently coupled to auxiliary qubits. The vibrational mode can function both as a highly sensitive quantum sensor and a programmable environment bath. Using this platform, we achieved three major outcomes. First, we demonstrated using the vibrational mode as an ultrasensitive probe for testing fundamental physics, including possible nonlinear quantum mechanics effects. Second, we established that these modes can act as controllable baths in which tunable noise and loss can enhance or modify energy-transfer dynamics, providing the experimental preparation toward studying mechanisms relevant to chemical reactions and light-harvesting systems. Third, by introducing controllable nonlinear gain and loss, we showed theoretically how simulations using trapped ions can model vibrationally-assisted energy transport in a non‐Hermitian quantum system comprising a chromophore dimer weakly coupled to a vibrational mode. Exploring the non‐Hermitian dynamics of the whole system including vibrations, we found that energy transfer accompanied by absorption of phonons from a vibrational mode can be significantly enhanced near an exceptional point. This theoretical work on simulation of energy transfer processes in driven non‐Hermitian quantum systems revealed an interesting novel path to study open quantum systems dynamics under conditions of gain and loss. We then further explored the benefits of controllable gain and loss with an experimental realization of quantum analogs of nonlinear oscillators, namely, the van der Pol oscillator. Here we observed mutual synchronization mediated by collective dissipation between two oscillators. In parallel, we explored related quantum networking protocols using the same trapped-ion platform, developing fast, high-fidelity schemes for distributing entanglement. Together, these achievements show that trapped-ion vibrational modes provide a highly programmable and high-fidelity platform for investigating complex dissipative quantum behavior, while enabling new approaches to remote quantum sensing, energy science, and nonlinear quantum dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Jamming Detection for Low-Resolution SC-FDE Systems: A Machine Learning Approach

Jammers interfere with communication between base stations (BSs) and legitimate users, leading to degradation of wireless system performance. Our study focuses on jamming detection for wideband single-carrier frequency domain equalization (SC-FDE) systems with low-resolution analog-to digital converters (ADCs). In such systems, jamming detection is challenging because traditional analytical approaches cannot be directly applied due to the delay dispersion in wideband channels and the non-linearity induced by low-resolution ADCs. We propose a machine learning (ML)-based jamming detection method that directly uses the quantized receive signals. Significantly, our ML-based detector can be integrated into existing standard frameworks, such as unique word (UW)-based SC-FDE systems, as it uses existing pilots without requiring additional pilots for jamming detection. Through numerical simulations, we show that two or more bits provide satisfactory performance compared to unquantized scenarios. Additionally, we demonstrate that using more and well-separated pilot symbols improves performance.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2026

Towards Generalizable and Efficient Circuit Topology Design: A Graph-Transformer-based Surrogate Model with Curriculum Learning

Unlike circuit parameter and sizing optimizations, the automated design of analog circuit topologies poses significant challenges for learning-based approaches. One challenge arises from the combinatorial growth of the topology space with circuit size, which limits the topology optimization efficiency. Moreover, traditional circuit evaluation methods are time-consuming, while the presence of data discontinuity in the topology space makes the accurate prediction of circuit performance exceptionally difficult for unseen topologies. To tackle these challenges, we design a novel Graph-Transformer-based Network (GTN) as the surrogate model for circuit evaluation, offering a substantial acceleration in the speed of circuit topology optimization without sacrificing performance. Our GTN model architecture is designed to embed voltage changes in circuit loops and current flows in connected devices, enabling accurate performance predictions for circuits with unseen topologies. To address the cold start problem when scaling GTN to large-scale circuits, we further introduce a curriculum learning strategy that progressively trains GTN from small-scale to large-scale circuits. This approach enables the model to first learn fundamental physical principles from simpler topologies and gradually adapt to complex configurations, effectively bridging the circuit complexity gap and improving prediction accuracy. Taking the power converter circuit design as an experimental task, our GTN model significantly outperforms an analytical approach and baseline methods directly utilizing graph neural networks. Furthermore, GTN achieves less than 5% relative error and 196× speed-up compared with high-fidelity simulation. Notably, our GTN surrogate model empowers an automatic circuit design framework to discover circuits of comparable quality to those identified through high-fidelity simulation while reducing the time required by up to 98.2%. With curriculum learning, the enhanced GTN achieves a 51% improvement for performance prediction of large-scale circuits compared to the GTN model without this strategy. These advancements establish GTN as a scalable framework for automated analog circuit design across varying circuit complexity levels.

Lu, Haoshu [New Jersey Institute of Technology (NJ↗

DOE OSTI2024

Characterization of Pliocene and Miocene Formations in the Wilmington Graben, Offshore Los Angeles, for Large-Scale Geologic Storage of CO2

The project Characterization of Pliocene and Miocene Formations in the Wilmington Graben, Offshore Los Angeles, for Large-Scale Geologic Storage of CO2 is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The Los Angeles Basin presents an opportunity for large-scale geologic CO2 storage. Due to its large population and historical and geologic setting as one of the most prolific oil and gas producing basins in the United States, the region is home to more than 12 major power plants and oil refineries that produce more than 5 million metric tons of fossil fuel-related CO2 emissions each year. GeoMechanics Technologies worked to characterize the Pliocene and Miocene sediments in the Wilmington Graben, offshore of Los Angeles, California, for high-volume CO2 storage. The Graben is located offshore of the Los Angeles and Long Beach Harbor area, making it accessible yet geologically isolated from the nearby Wilmington oilfield and onshore areas. These sediments span more than 5,000 feet of vertical interval with an estimated storage resource of more than 100 million metric tons of CO2. The project team analyzed and interpreted existing geologic data within the region, including detailed exploration well log data and 2-D and 3-D seismic data. New seismic lines were acquired to fill in current data gap areas and two new characterization wells were drilled and logged. This information was integrated with existing geologic interpretations for adjacent onshore areas to help characterize optimal areas for CO2 storage and seals to safely store CO2. Integrated 3-D geologic and geomechanical models for the Wilmington Graben were developed to simulate the fate and transport of injected CO2 in the subsurface and to assess risks. This project contributed to the understanding of injectivity, containment mechanisms, rate of dissolution and mineralization, and storage capacity of the Wilmington Graben and associated analogous basins. This effort also provided greater insight into the potential for offshore geologic formations to safely and permanently store CO2.

.las↗

DOE OSTI2026

Geomechanical properties of the Meletta sandstone - the high-temperature heat storage reservoir rock of DeepStor

The DeepStor project aims at storing excess heat at temperatures up to 140 °C in the depleted Leopoldshafen oil field at a depth of about 1300 m. In order to gain knowledge on the target horizons, the different layers of the Meletta sandstone, samples cored in a block retrieved in a quarry near Nussloch where the Meletta sandstone outcrops were studied. Several petrophysical properties were investigated including mineralogy, porosity, permeability, thermal conductivity, P- and S-wave velocities and quality factor. A mechanical study focused on compressive strength under uniaxial and triaxial stress conditions, tensile strength, critical pressure and stress-dependence of physical properties. Our results show that the Meletta sandstone is heterogeneous, anisotropic, mechanically weak, stress-sensitive and prone to water weakening. Petrophysical measurements on few Meletta sandstone cores retrieved in boreholes at about 1250 m depth showed that the outcrop samples are significantly less cohesive and weaker. Therefore, one should consider the Nussloch outcrop as a good reservoir analog only with great caution. Several models were applied to account for our experimental results on anisotropy, mechanical behavior and stress-dependence of the outcrop samples. Our study provides a valuable data set which can feed numerical models to simulate the behavior of the host formation during heat exchange operations in the DeepStor project.

Geomechanics↗

DOE OSTI2024

On the connection between least squares, regularization, and classical shadows

Classical shadows (CS) offer a resource-efficient means to estimate quantum observables, circumventing the need for exhaustive state tomography. Here, we clarify and explore the connection between CS techniques and least squares (LS) and regularized least squares (RLS) methods commonly used in machine learning and data analysis. By formal identification of LS and RLS ``shadows'' completely analogous to those in CS---namely, point estimators calculated from the empirical frequencies of single measurements---we show that both RLS and CS can be viewed as regularizers for the underdetermined regime, replacing the pseudoinverse with invertible alternatives. Through numerical simulations, we evaluate RLS and CS from three distinct angles: the tradeoff in bias and variance, mismatch between the expected and actual measurement distributions, and the interplay between the number of measurements and number of shots per measurement. Compared to CS, RLS attains lower variance at the expense of bias, is robust to distribution mismatch, and is more sensitive to the number of shots for a fixed number of state copies---differences that can be understood from the distinct approaches taken to regularization. Conceptually, our integration of LS, RLS, and CS under a unifying ``shadow'' umbrella aids in advancing the overall picture of CS techniques, while practically our results highlight the tradeoffs intrinsic to these measurement approaches, illuminating the circumstances under which either RLS or CS would be preferred, such as unverified randomness for the former or unbiased estimation for the latter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗