Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ground Processing”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Supersonically expanded sodium metal-dilute halogen gas interactions. The importance of reaction populated and energy storing reservoir states and population inversion created amplification in Na2

The reactions of Cl, Br, and I with Nan=2,3 produced in a supersonic expansion form Na2* and Na* excited states extending across the visible and ultraviolet regions. Emission in the region extending from 410 to 600 nm indicates selectively formed excited state Na2 emission features. Experimental evidence suggests that this emission is associated with Na3 + X reactions. Broadband (0.5 cm−1) laser measurements demonstrate gain (population inversion) for select features at∼524–528(1%), ∼492(0.3%), and ∼458.7–461(0.8%) nm. Single mode (0.007 cm−1) measurements extending from 528.03 to 527.63 nm demonstrate amplification involving five to six individual rovibronic levels with a maximum gain close to 3% recorded at 527.9 nm. The observed gain is associated with select transitions from levels of the Na2 11Πu state populated, via identified curve crossings, through collision induced transfer from long-lived Na2 21Σg+ and 11Πg reservoir states. Collision induced population buildup in the lowest vibrational levels of these reservoir states and collision induced transfer to the Na2 11Πu state create a population inversion in transitions to the X 1Σg+ state of Na2. The observed amplification is aided by rapid vibrational and rotational relaxation in both the Na2 ground and excited reservoir states producing amplifiers in the visible region like the HF amplifier in the infrared. This study suggests the importance of reaction populated and energy storing long-lived reservoir states in small sodium molecule combustion processes and indicates the potential for providing new short wavelength visible and ultraviolet amplifiers for future laser-based chemical propulsion concepts.

Chemistry↗
DOE OSTI2023

A robust estimator of mutual information for deep learning interpretability

Abstract We develop the use of mutual information (MI), a well-established metric in information theory, to interpret the inner workings of deep learning (DL) models. To accurately estimate MI from a finite number of samples, we present GMM-MI (pronounced ‘Jimmie’), an algorithm based on Gaussian mixture models that can be applied to both discrete and continuous settings. GMM-MI is computationally efficient, robust to the choice of hyperparameters and provides the uncertainty on the MI estimate due to the finite sample size. We extensively validate GMM-MI on toy data for which the ground truth MI is known, comparing its performance against established MI estimators. We then demonstrate the use of our MI estimator in the context of representation learning, working with synthetic data and physical datasets describing highly non-linear processes. We train DL models to encode high-dimensional data within a meaningful compressed (latent) representation, and use GMM-MI to quantify both the level of disentanglement between the latent variables, and their association with relevant physical quantities, thus unlocking the interpretability of the latent representation. We make GMM-MI publicly available in this GitHub repository.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Insights of warm-cloud biases in Community Atmospheric Model 5 and 6 from the single-column modeling framework and Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) observations

There has been a growing concern that most climate models predict precipitation that is too frequent, likely due to lack of reliable subgrid variability and vertical variations in microphysical processes in low-level warm clouds. In this study, the warm-cloud physics parameterizations in the singe-column configurations of NCAR Community Atmospheric Model version 6 and 5 (SCAM6 and SCAM5, respectively) are evaluated using ground-based and airborne observations from the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) field campaign near the Azores islands during 2017–2018. The 8-month single-column model (SCM) simulations show that both SCAM6 and SCAM5 can generally reproduce marine boundary layer cloud structure, major macrophysical properties, and their transition. The improvement in warm-cloud properties from the Community Atmospheric Model 5 and 6 (CAM5 to CAM6) physics can be found through comparison with the observations. Meanwhile, both physical schemes underestimate cloud liquid water content, cloud droplet size, and rain liquid water content but overestimate surface rainfall. Modeled cloud condensation nuclei (CCN) concentrations are comparable with aircraft-observed ones in the summer but are overestimated by a factor of 2 in winter, largely due to the biases in the long-range transport of anthropogenic aerosols like sulfate. We also test the newly recalibrated autoconversion and accretion parameterizations that account for vertical variations in droplet size. Compared to the observations, more significant improvement is found in SCAM5 than in SCAM6. This result is likely explained by the introduction of subgrid variations in cloud properties in CAM6 cloud microphysics, which further suppresses the scheme's sensitivity to individual warm-rain microphysical parameters. The predicted cloud susceptibilities to CCN perturbations in CAM6 are within a reasonable range, indicating significant progress since CAM5 which produces an aerosol indirect effect that is too strong. The present study emphasizes the importance of understanding biases in cloud physics parameterizations by combining SCM with in situ observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Quantifying uncertainty in Pareto estimates of global lake area

Abstract Size is a critical factor determining the rate and occurrence of specific lake processes such as carbon sequestration and greenhouse gas emissions and emerging evidence suggests that small lakes in particular have particularly large CO 2 flux rates. Because we do not have a complete census of all lakes, upscaling estimates of such processes to small lakes at broad spatial scales requires the use of lake size‐abundance distributions rather than empirical measurements of area. Existing lake census efforts are incomplete such that as lakes become smaller, they are more likely to be omitted either because they are too small to be resolved from remote sensing products or because of limited ground surveying effort (i.e., “censoring” of small lakes relative to large lakes). The present study explores one potential shortcoming of prior approaches estimating global lake area using lake size‐abundance distributions. Namely, that these prior approaches rely on frequentist curve fitting techniques combined with an ad‐hoc cutoff determination strategy (visual inspection to determine a likely censoring point). This yields an over‐exact lake area estimate that is typically reported with no uncertainty bounds. I show how these shortcomings can be addressed with a Bayesian model that produces larger estimates of lake area uncertainty relative to the typical approach. When used as part of a sensitivity analysis, such an approach has the potential to enable more robust intercomparisons among studies of aquatic processes upscaling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Frost Table Elevation across a Low-Centered and a High-Centered Polygon, Mapped using Ground Penetrating Radar, Utqiagvik (Barrow), Alaska, 2015

Ground Penetrating Radar, active layer thickness point measurements, and ground elevation surveys were conducted on a low-centered and a high-centered polygon in Utqiagvik (Barrow), AK (on October 2nd, 2015) to map frost table topography and evaluate its influence on hydrological flow. Data were collected to improve our understanding of the topography of the permafrost table across polygons which may have a major influence on hydrology. This dataset includes 2 zip files containing the raw data, 2 csv files of processed data, and a pdf user guide.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion quantum Monte Carlo

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing the accuracy. Advances in real space methods include techniques for accurate computation of bandgaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods, including GW and density functional based techniques. To provide an improved foundation for these calculations, we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Proton-Coupled Electron Transfer in a Ruthenium(II) Bipyrimidine Complex in Its Ground and Excited Electronic States

In this work, proton-coupled electron transfer (PCET) was studied for the ground and excited electronic states of a [Ru(terpy)(bpm)(OH 2 )(PF 6 ) 2 ] complex, Ru-bpm. Cyclic voltammetry measurements show that the Ru(II)-aqua moiety undergoes PCET to form a Ru(IV)-oxo moiety in the anodic region, while the bpm ligand undergoes PCET to form bpmH 2 in the cathodic region. The photophysical behavior of Ru-bpm was studied using steady-state and femtosecond transient UV–vis absorption spectroscopy, coupled with density functional theory (DFT) calculations. The lowest-lying excited state of Ru-bpm is described as a (Ru → bpm) metal-to-ligand charge-transfer (MLCT) state, while the metal-centered (MC) excited state was found computationally to be close in energy to the lowest-energy bright MLCT state (MC state was 0.16 eV above the MLCT state). The excited-state kinetics of Ru-bpm were found via transient absorption spectroscopy to be short-lived and were fit well to a biexponential function with lifetimes $τ_1$ = 4 ps and $τ_2$ = 65 ps in aqueous solution. Kinetic isotope effects of 1.75 ($τ_1$) and 1.61 ($τ_2$) were observed for both decay components, indicating that the solvent plays an important role in the excited-state dynamics of Ru-bpm. Based on the pH-dependent studies and the results from prior studies of similar Ru-complexes, we hypothesize that the 3 MLCT state forms an excited-state hydrogen-bond adduct with the solvent molecules and that this process occurs with a 4 ps lifetime. The formation of such a hydrogen-bond complex is consistent with the electronic density accumulation at the peripheral N atoms of the bpm moiety in the 3 MLCT state. The hydrogen-bonded state 3 MLCT decays to the ground state with a 65 ps lifetime. Such a short lifetime is likely associated with the efficient vibrational energy transfer from the 3 MLCT state to the solvent.

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

Layer-wise Imaging Dataset from Powder Bed Additive Manufacturing Processes for Machine Learning Applications (Peregrine v2022-10.1)

This release consists of six datasets which together include multi-modal layer-wise powder bed images from two different powder bed printing technologies. These datasets are designed primarily to facilitate the development and testing of new computer vision and machine learning based anomaly and defect detection algorithms. The authors provide both training data with corresponding ground truth pixel masks and evaluation data with corresponding baseline prediction pixel masks made by a trained neural network. The laser powder bed fusion (L-PBF) datasets are sourced from EOS M290 and AddUp FormUp 350 printers and the binder jet (BJ) dataset is sourced from an ExOne M-Flex printer. The materials represented in these datasets include 17-4 PH Stainless Steel, GammaPrint-700, Inconel 718, Maraging Steel, and H13 Steel. The sensor imaging modalities represented include visible-light (VL), temporally-integrated (i.e., long duration exposure) near-infrared (TI-NIR), and wide-band infrared (IR). To download the dataset: (1) Create a Globus account. (2) Create a Globus Endpoint on your computer. (3) Transfer the dataset from the OLCF DOI-DOWNLOADS Collection to your Collection. Common troubleshooting steps: (a) Confirm that the transfer is going from OLCF DOI-DOWNLOADS to your Collection. (b) Create an exception for Globus in your antivirus software so that it can create an Endpoint. (c) Manually create a Globus access directory (where the data will be downloaded) by going to the Preferences > Access tab.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Zn speciation and fate in soils and sediments along the ground transportation route of Zn ore to a smelter

We discuss assessment of Zn toxicity/mobility based on its speciation and transformations in soils as critical for maintaining human and ecosystem health. Zn-concentrate (56% Zn as ZnS, sphalerite) has been imported through a seaport and transported to a Zn-smelter for several decades, and smelting processes resulted in aerial deposition of Zn and sulfuric acids in two geochemically distinct territories around the smelter (mountain-slope and riverside). XAFS analysis showed that the mountain-slope soils contained franklinite (ZnFe 2 O 4 ) and amorphous (e.g., sorbed) species of Zn(II), whereas the riverside sediments contained predominantly hydrozincite [Zn 5 (OH) 6 (CO 3 ) 2 ], sphalerite, and franklinite. The mountain-slope soils had low pH and moderate levels of total Zn (~1,514 ppm), whereas the riverside sediments had neutral pH and higher total Zn (12,363 ppm). The absence of sphalerite and the predominance of franklinite in the mountain-slope soils are attributed to the susceptibility of sphalerite and the resistance of franklinite to dissolution at acidic pH. These results are compared to previous Zn analyses along the transportation routes, which showed that Zn-concentrate spilled along the roadside in dust and soils underwent transformation to various O-coordinated Zn species. Overall, Zn-concentrate dispersed in soils and sediments during transportation and smelting transforms into Zn phases of diverse stability and bioavailability during long-term weathering.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Application of unsupervised deep learning to image segmentation and in-situ contact angle measurements in a CO 2 -water-rock system

Rock surface wettability is a critical property that regulates multiphase flows in porous media, which can be quantified using the surface contact angle (CA). X-ray micro-computed tomography (μCT) provides an effective approach to in-situ measurements of surface CAs. However, the CA measurement accuracy depends significantly on the quality of CT image segmentation, which is the clustering of CT pixels into separate phases. Inspired by this, we developed a deep learning (DL)-based CA measurement workflow. Motivated by the recent tremendous progress in unsupervised learning techniques and aiming to avoid expensive manual data annotations, an unsupervised DL pipeline for CT image segmentation was proposed and implemented, which includes unsupervised model training and post-processing. The unsupervised model training was driven by a novel loss function constrained with feature similarity and spatial continuity and implemented by iterative forward and backward paths; the former clustered the pixel-wise feature vectors extracted by convolution neural networks, whereas the latter updated the parameters using gradient descent. An over-segmentation strategy was adopted for model training. The post-processing steps based on agglomerative hierarchical clustering (AHC) were implemented to further merge the over-segmented model output to the desired cluster number, which is intended to improve the efficiency of image segmentation. The developed unsupervised DL pipeline was compared with other commonly-used image segmentation methods using pixel-wise and physics-based evaluation metrics on a synthetic raw-image dataset, which had a known ground truth. The unsupervised DL pipeline showed the best performance. Next, the segmented images were input to an automatic CA measurement tool, and the results were validated by comparisons with manual measurements. The CA values from the manual and automatic measurements showed similar distributions and statistical properties. The automatic measurement demonstrated a wider spectrum because of the much larger number of measurement data points. The primary novelty of the unsupervised DL pipeline developed in this study lies in the novel loss function and the over-segmentation strategy associated with AHC post-processing. Finally, the workflow has been proven an efficient tool for pore-scale wettability characterization, which has a wide range of applications in fundamental studies of multiphase flows in natural porous media, which have critical implications to geological carbon sequestration, hydrocarbon energy recovery, and contaminant transport in groundwater.

42 ENGINEERING↗
DOE OSTI2024

Non-destructive, whole-plant phenotyping reveals dynamic changes in water use efficiency, photosynthesis, and rhizosphere acidification of sorghum accessions under osmotic stress

Noninvasive phenotyping can quantify dynamic plant growth processes at higher temporal resolution than destructive phenotyping and can reveal phenomena that would be missed by end-point analysis alone. Additionally, whole-plant phenotyping can identify growth conditions that are optimal for both above- and below-ground tissues. However, noninvasive, whole-plant phenotyping approaches available today are generally expensive, complex, and non-modular. We developed a low-cost and versatile approach to noninvasively measure whole-plant physiology over time by growing plants in isolated hydroponic chambers. We demonstrate the versatility of our approach by measuring whole-plant biomass accumulation, water use, and water use efficiency every two days on unstressed and osmotically stressed sorghum accessions. We identified relationships between root zone acidification and photosynthesis on whole-plant water use efficiency over time. Our system can be implemented using cheap, basic components, requires no specific technical expertise, and should be suitable for any non-aquatic vascular plant species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Signatures of the exciton gas phase and its condensation in monolayer 1T-ZrTe 2

The excitonic insulator (EI) is a Bose-Einstein condensation (BEC) of excitons bound by electron-hole interaction in a solid, which could support high-temperature BEC transition. The material realization of EI has been challenged by the difficulty of distinguishing it from a conventional charge density wave (CDW) state. In the BEC limit, the preformed exciton gas phase is a hallmark to distinguish EI from conventional CDW, yet direct experimental evidence has been lacking. Here we report a distinct correlated phase beyond the 2×2 CDW ground state emerging in monolayer 1T-ZrTe 2 and its investigation by angle-resolved photoemission spectroscopy (ARPES) and scanning tunneling microscopy (STM). The results show novel band- and energy-dependent folding behavior in a two-step process, which is the signatures of an exciton gas phase prior to its condensation into the final CDW state. Our findings provide a versatile two-dimensional platform that allows tuning of the excitonic effect.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Absolute timing calibration of the Resolve microcalorimeter spectrometer on the X-ray Imaging and Spectroscopy Mission

The Resolve microcalorimeter spectrometer on the X-ray Imaging and Spectroscopy Mission (XRISM) is designed to have a good timing capability with the mission-level requirement of 1 ms as the absolute time tagging accuracy to suffice the needs for observatory science. In the ground calibration campaign, the absolute and relative timing offsets were measured using pulsed X-rays from the modulated X-ray sources. These offsets were used to determine calibration parameters, the Resolve timing coefficients, which are used in the offline correction of event times. In the orbit, reevaluation of the timing coefficients was carried out using the Crab pulsar, which was observed in two periods, one in the performance verification phase and the other in the guest observation phase. We report the absolute timing calibration of XRISM/Resolve using the ground and in-orbit data. Although the requirement is likely satisfied with the ground parameters, the timing coefficients have been refined using the in-orbit calibration to improve the timing accuracy. In addition, for the first time, we present the interpretation of the absolute timing offset originating from the analog and digital processing of X-ray events unique to a microcalorimeter spectrometer and present a complete view of the timing error of the Resolve instrument.

Astronomy and AstroPhysics↗
DOE OSTI2020

Ground-based_transient_electromagnetic_data_and_resistivity_models_for_SubTER_BelleCreek_2017-2018

Ground-based transient electromagnetic (TEM) data were acquired in selected locations around Belle Creek, Montana to define the resistivity structure of the near surface. Data were acquired using the ABEM WalkTEM system (Guideline Geo Ab, Sundbyberg, Sweden). TEM data were processed and numerically inverted to derive one-dimensional resistivity structure using SPIA (Aarhus Geosoftware, Aarhus, Denmark). This release includes raw and processed TEM data and resistivity models for each data location using the following structure: Ground-based_transient_electromagnetic_data_and_resistivity_models_for_SubTER_BelleCreek_2017-2018.xml SBC_TEM_2019_SoundingLocationsModels.gdb.zip - file geodatabase with all models and plots attached. This is a mirror of the contents of the model directory. Ground-based_transient_electromagnetic_data_and_resistivity_models_for_SubTER_BelleCreek_2017-2018.zip Data SBC_USF - raw data USF_format_description.pdf - manual describing the USF file format SBC_SoundingID-USF_lookup.csv - table describing the USF file naming convention and channel mapping *.usf files - raw, unprocessed data files SBC_TEM - processed data AarhusInv_manual.pdf - manual describing the TEM file format *.tem files - processed data files Models SBC_Model_Plots - plots of the data and the models for each sounding location TEM_model_data_dictionary.csv - data dictionary describing tabulated models SBC_TEM_models.csv - tabulated models

geophysical surveying↗
DOE OSTI2021

Revolutionizing observations and predictability of Arctic system dynamics through next-generation dense, heterogeneous and intelligent wireless sensor networks with embedded AI

The Arctic region is a complex and dynamic system, where changes in temperature and partitioning between water in solid, liquid and vapor phase are fundamentally reshaping hydro-biogeochemical fluxes from the bedrock to the atmosphere. The increasing temperature and frequency of extreme heat, precipitation and fire events (Wang et al., 2017; Meredith et al., 2019) have critical consequences for the Arctic ecosystems and people, as well as for the global climate system. Although many observational and remote sensing activities have been on-going, detailed ecosystem processes – particularly related to controls on and impacts of permafrost thaw – are difficult to detect or measure. To improve the understanding and prediction of these extreme events and their consequences, we need a revolution in the way we collect ground-based observations from hillslope to pan-Arctic scales and couple them with satellite imagery and models to deliver rich, actionable, and scale-appropriate data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Direct cross section measurement of 102 Pd ⁢(𝛾,𝑝) and 102 Pd ⁢(𝛾,𝛼) for the astrophysical 𝑝 process

Background: A handful of neutron-deficient stable nuclei, known as the “p nuclei,” cannot be produced through astrophysical neutron capture processes. Instead, some of these nuclei are proposed to be produced by 𝛾-induced reactions on existing r- and s-process seeds. The specific astrophysical site or sites are not yet identified, however, with uncertainties in the cross sections of these 𝛾-induced reactions playing a role. Databases of reaction rates for astrophysical simulations often rely on theoretical statistical model calculations, such as Hauser-Feshbach, for rates where no experimental information is known. However, reasonable variations in the choice of parametrizations of various nuclear properties can create order-of-magnitude variations in the final predicted cross sections and reaction rates, which are then propagated through the models to the predicted final abundances. Purpose: To better constrain these statistical model calculations and ultimately reduce the uncertainties from the nuclear physics on our understanding of the p nuclei, a measurement of the cross sections of 𝛾-induced reactions on the p-nucleus 102 Pd was undertaken. This work represents the first measurement of its kind, using segmented silicon detectors to measure prompt charged particle emission from 𝛾-induced reactions. Methods: Quasimonoenergetic gamma beams from the High Intensity 𝛾 Source facility bombarded an enriched 102 Pd target. A segmented silicon array was arranged to detect the particles emitted from (𝛾,𝑝) and (𝛾,𝛼) reactions. Results: Reaction cross sections were deduced at multiple 𝛾-beam energies between 10 and 19 MeV, and compared to statistical model calculations using talys-1.96. The 102 Pd ⁢(𝛾,𝑝)⁢ 101 Rh reaction cross section was reasonably well reproduced by a subset of photon strength functions and level densities, though the strength to the ground state of 101 Rh was underestimated at higher incident gamma energies. The 102 Pd ⁢(𝛾,𝛼)⁢ 98 Ru was in general overpredicted by the various alpha-nucleus optical model potentials. Conclusions: While the theoretical cross sections used to model the (𝛾,𝑝) reactions for the p process may be reasonable, a more careful approach is needed in the case of (𝛾,𝛼). Further work to probe gamma-induced reaction cross sections at and near the p nuclei is warranted.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.

ATS↗
DOE OSTI2021

Quantum Solver of Contracted Eigenvalue Equations for Scalable Molecular Simulations on Quantum Computing Devices

The accurate computation of ground and excited states of many-fermion quantum systems is one of the most consequential, contemporary challenges in the physical and computational sciences whose solution stands to benefit significantly from the advent of quantum computing devices. Existing methodologies using phase estimation or variational algorithms have potential drawbacks such as deep circuits requiring substantial error correction or non-trivial high-dimensional classical optimization. In this work, we introduce a quantum solver of contracted eigenvalue equations, the quantum analogue of classical methods for the energies and reduced density matrices of ground and excited states. The solver does not require deep circuits or difficult classical optimization and achieves an exponential speed-up over its classical counterpart. We demonstrate the algorithm though computations on both a quantum simulator and two IBM quantum processing units.

97 MATHEMATICS AND COMPUTING↗