Search NASA⌕ Search

SEARCH · Search NASA

Results for “PHS”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Chemical Functional Groups Regulate Ion Concentrations and pHs in Nanopores

Understanding ion behaviors in functionalized nanopores is essential to deciphering reactions in both natural and engineered systems, such as sediments, biological ion channels, and membranes. While many efforts have shown the modified ion behaviors in the functionalized nanopores, a direct measurement and analysis to show how chemical functional groups affect ion concentrations in nanopores are critically needed. In this work, we present a plasmonic nanosensor that can measure the local concentrations of protons, anions (phosphate, nitrate, sulfate, and arsenate), and cations (mercury, lead, and copper) in functionalized nanopores, and we compare their concentrations in nanopores with the corresponding bulk concentrations. Notably, chemical functional groups induced ion concentrations differently in nanopores. In pristine nanopores and methyl- and phenyl-functionalized nanopores, we discovered an unexpected concurrence of an enhanced anion concentration and a suppressed cation concentration. In addition, the nanopore pH is dependent on bulk solution compositions and can be lower by 2.5 units, even when the bulk solution is well-buffered. In contrast, for hydrophilic (amine, thiol, and carboxyl) nanopores, pH depended on the p K a of the functional groups, and the heavy metal concentrations depended on chemical interactions with the functional groups. Our findings provide a better understanding of water chemistry in nanopores and can help precisely control ions in nanopores to benefit the design of membrane-based desalination techniques, CO 2 storage, and porous catalysts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Refining water and carbon fluxes modeling in terrestrial ecosystems via plant hydraulics integration

Plant hydraulics substantially affects terrestrial water and carbon cycles by modulating water transport and carbon assimilation. Despite improved drought simulations in certain ecosystems through their integration into land surface models (LSMs), the broader application of plant hydraulics in diverse ecosystems and hydroclimates is still underexplored. Here, in this study, we implemented the recently developed Noah-Multiparameterization Land Surface Model (Noah-MP LSM) equipped with a plant hydraulics scheme (Noah-MP-PHS) across 40 FLUXNET sites globally. Employing the Shuffled Complex Evolution-University of Arizona (SCE-UA) auto-calibration algorithm, we optimized key plant hydraulics parameters for these sites spanning eight vegetation types in both arid and humid climates. Noah-MP-PHS significantly improves the simulation of evapotranspiration (ET) and gross primary production (GPP) by better representing atmospheric and soil water stress compared to traditional soil hydraulic schemes (SHSs, such as Noah and CLM). The augmented Noah-MP-PHS models reduce surface flux overestimation and underestimation, exhibiting an average increase of 0.14 and 0.15 in Kling-Gupta Efficiency (KGE) compared to Noah and CLM, respectively. The explicit consideration of plant capacitance in PHS reveals substantial deep-layer and nocturnal root water uptake especially under dry conditions. We employed eXplainable Machine learning (XML) to quantify the model’s relative sensitivity to newly introduced leaf-, stem and root-related parameters in PHS. The sensitivity analysis reveals a rise in root parameter importance and a decline in leaf and stem parameters as conditions shift from humid to arid. These findings indicate that as aridity states vary, the most influential parameters affecting surface fluxes variation may change in parameter calibration for PHS applications. Our findings underscore the importance of incorporating plant hydraulics into LSMs to enhance simulations of terrestrial water and carbon dynamics. These findings are crucial for understanding ecosystem responses to global climate changes and guide the broader application of PHS at larger scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Enhancing SWAT with mechanistic plant hydraulics: development and application in the Hanjiang River Basin

Plant transpiration plays a critical role in global water and energy cycles, requiring better process understanding as climate change intensifies drought stress and alters plant responses. Most hydrological models such as the widely-used SWAT lack representation of plant hydraulics, the mechanistic processes controlling plant water regulation and transpiration. Here, this study developed SWAT-PHS by integrating a plant hydraulics scheme (PHS) into SWAT hydrological model, enabling explicit simulation of root water uptake, sap flow, storage and transpiration at 30-minute timescales for watershed-scale modeling. In the Hanjiang River Basin, SWAT-PHS mitigated overestimation of runoff during the rainy season and underestimation during the dry season, reducing the overall simulation error by 29% across the entire simulation period. The model can simulate reasonable plant water dynamics, including diurnal transpiration patterns and drought responses showing declining transpiration flux, hydraulic buffering through stem water storage, and depth-dependent root water uptake strategies. Sensitivity analysis shows that SWAT-PHS captured mechanistic relationships between plant hydraulic traits and transpiration, with root distribution and stem capacitance positively affecting annual transpiration while vulnerability parameters showed negative effects. This work provides a pathway for improving hydrologic modeling and water resource management by better representing plant water regulation under climate change and expected intensifying water stress conditions.

China↗
DOE OSTI2025

Extreme Ultraviolet and Beyond Extreme Ultraviolet Lithography Using Amorphous Zeolitic Imidazolate Resists Deposited by Atomic/Molecular Layer Deposition

Amorphous zinc-imidazolate (aZnMIm) resists show potential to meet the demands for next-generation high-numerical aperture (high-NA) metal-containing extreme ultraviolet (EUV) resist materials, given their ease of deposition by atomic/molecular layer deposition (ALD/MLD) at thicknesses of 20 nm and below. Here, this study demonstrates that aZnMIm thin films, previously identified as high-resolution electron beam resists, can also function as negative-tone EUV photoresists. Water development achieves high sensitivity (5 mJ/cm 2 ) but leaves significant residue, while acetic acid development results in poor contrast. A hybrid approach─water followed by acetic acid─enables residue-free development with a sensitivity of 181 mJ/cm 2 . Dry development using 1,1,1,5,5,5-hexafluoroacetylacetone (hfacH) is also possible but shows lower sensitivity (375 mJ/cm 2 ) compared to wet development methods. EUV photoelectron spectroscopy (PES), reflectometry/EUV absorption, total electron yield (TEY), residual gas analysis (RGA), and time-of-flight secondary ion mass spectrometry (TOF-SIMS) were used to investigate the effects of EUV irradiation on aZnMIm resists. Reflectometry experiments reveal an aZnMIm EUV absorption coefficient of 6.2 μm –1 , while PES and TEY analyses show that, compared to poly(4-hydroxystyrene) (PHS), a polymer-based reference resist, aZnMIm emits more primary and secondary electrons but generates fewer slow electrons relative to its primary electron emission; its total electron yield is similar to that of poly(methyl methacrylate) (PMMA) resists. When exposed to EUV, aZnMIm predominantly outgasses H 2 , as determined by RGA. TOF-SIMS measurements demonstrate that high-dose EUV exposure only partially fragments the 2-methylimidazole (2MIm) organic linkers, unlike high-dose electron beam exposure, which is known to completely degrade them. Additionally, aZnMIm resists show promise for potential beyond EUV lithography (BEUVL) due to the presence of Zn, which provides higher sensitivity at a wavelength of 6.7 nm compared to other metal ions, such as Sn, that are currently used in the best-performing EUV metal–organic resists. TEY measurements demonstrate that aZnMIm emits nearly twice as many electrons as PHS at 6.7 nm. The BEUV TEY of aZnMIm also surpasses that of PMMA, poly(pentafluorostyrene), and poly(4-iodostyrene), with the latter two being known for their high EUV TEYs. This work provides insight into zeolitic imidazolate framework (ZIF)-based EUV and BEUV resists and highlights their potential for both wet and dry development.

lithography↗
DOE OSTI2026

Plutonium(III) versus uranium(III) and samarium(III) in small molecule activation chemistry

We report the PuIII complex, [PuIII(CpMe4)3] (1-Pu), and demonstrate its differences in small molecule reactivity compared to the UIII and SmIII analogs, [UIII(CpMe4)3] (1-U) and [SmIII(CpMe4)3] (1-Sm), respectively. 1-Pu reductively cleaves the small molecule (PhS)2, affording a PuIII complex, [{PuIII(CpMe4)2}2(μ-SPh)2] (2-Pu), while retaining the PuIII center and eliminating (CpMe4)2 as a by-product, a fingerprint of a sterically induced reduction (SIR) reaction. Sm is often used as a surrogate for Pu, but the analogous [SmIII(CpMe4)3], (1-Sm), is unreactive. The (PhS)2 cleavage by 1-U proceeds solely via a metal-based oxidation (i.e., UIII → UIV), to form [UIV(CpMe4)3(SPh)] (3-U). Only 1-U reacts with (PhHN)2, affording the reductive cleavage product, [UIV(CpMe4)3(NHPh)] (4-U). The difference in reactivity of 1-Pu compared to complexes 1-Sm and 1-U was unexpected, and since SIR chemistry can enable complexes to participate in otherwise impossible reductive transformation of substrates, this reinforces the importance of studying small molecule reactivity with the transuranic elements.

Keener, Megan↗
DOE OSTI2025

Acidities of MgO surface sites: implications for the formation mechanism of Mg(OH) 2

The hydroxylation of periclase (MgO) to brucite (Mg(OH) 2 ) is thought to be an important intermediate step when using MgO to capture CO 2 from the atmosphere. However, the mechanism of hydroxylation of MgO to form Mg(OH) 2 is poorly understood. In this work, we used atomic-scale density functional tight binding simulations coupled with the metadynamics rare event method to analyze the surface chemistry of MgO and the acid dissociation equilibrium constants (pK a ) of its surface sites. The method and parameters were validated by calculating the pK a for hydroxylation of the first shell water bound to aqueous Mg 2+ ion. The pK a value derived using a probabilistic method was 12.3, which is in fair agreement with the accepted value of 11.4, with the difference between them equal to a ∼5 kJ mol −1 error in the calculations. We then extended these pK a calculations to probe the hydroxylation reactions of the surface sites of the MgO(100)–water interface, arriving at pK a s of 5.4 to deprotonate terminal water molecules bound to the surface magnesium sites (η-OH 2 or 〉MgOH 2 ), and 13.9 to deprotonate hydroxylated bridging oxygen sites (μ 5 -oxo or 〉O). Hydroxide (OH − ) adsorption on the surface was also probed and found to be less thermodynamically favorable than deprotonation of the terminal water molecule. The plausibility of the computed pK a s was verified using an activity-based speciation model and compared to pH measurements of water equilibrated with MgO nanoparticles and single crystals. The model predicted a solution pH of 7.1 when surface sites buffered and the pH of 12.0 when MgO dissolution dominated. These are close to the experimental initial solution pHs of 7–7.5 and the long term pHs of ∼10.5. The similarity suggests that the calculated pK a values from the DFTB+/metadynamics simulations are plausible and that these methods can be a useful tool to probe reaction mechanisms involving covalent bonds.

Adapa, Sai Krishna Reddy [Oak Ridge National Labor↗
DOE OSTI2026

Contrasting effects of glutamate and branched-chain amino acid metabolism on acid tolerance in a Castellaniella isolate from acidic groundwater

Groundwater acidification co-occurring with nitrate pollution is a common, global environmental health hazard. Denitrifying bacteria have been leveraged for the in situ removal of nitrate in groundwater. However, co-existing stressors—such as low pH—reduce the efficacy of biological removal processes. Castellaniella sp. str. MT123 is a complete denitrifier that was isolated from acidic, nitrate-contaminated groundwater. The strain grows robustly by nitrate respiration at pH < 6.0, completely reducing nitrate to dinitrogen gas. Genomic analyses of MT123 revealed few previously characterized acid tolerance genes. Thus, we utilized a combination of proteomics, metabolomics, and competitive mutant fitness to characterize the genetic mechanisms of MT123 acclimation to growth under mildly acidic conditions. We found that glutamate accumulation is critical in the acid acclimation of MT123, possibly through consumption of intracellular protons via glutamate decarboxylation to GABA. This is despite the fact that MT123 lacks the canonical glutamate decarboxylase-glutamate/GABA antiporter system implicated in acid tolerance in other bacteria. In contrast, branched-chain amino acid (BCAA) accumulation was detrimental to cell growth at lower pHs, possibly through indirect mechanisms impacting the cellular glutamate pool. Genetic analysis previously linked MT123 to a population of Castellaniella that bloomed—concurrent to nitrate removal—during a biostimulation effort to reduce groundwater nitrate concentrations at MT123’s location of origin. Thus, our analyses provide novel insight into mechanisms of acclimation to acidic conditions in a strain with significant potential for nitrate bioremediation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Cation exchange separation of radium and actinium using lactic acid DTPA buffer

The separation of actinium ( 228 Ac) from radium ( 223,228 Ra) on cation exchange resin with a diethylenetriaminpentaacetic acid-lactate buffer solution is demonstrated with series of columns. High yield, high radiopurity (~ 100%) separations of Ac from Ra are feasible with small columns in biologically compatible conditions and pHs. As Ac is eluted before Ra on these columns, further studies were performed to determine if this separation system could be applied to Ra/Ac isotope generators, but these were not successful. In conclusion, the separations presented in this work may be relevant for radiopharmaceutical purifications of 225 Ac which is typically obtained from its 225 Ra parent isotope.

and nuclear chemistry↗
DOE OSTI2026

Demonstration of a 228 Ra/ 228 Ac isotope generator

An isotope generator to produce 228 Ac from 228 Ra was developed using a cation exchange resin column eluted with an acetate-diethylenetriaminpentaacetic acid buffer. Here, the elution behavior of 228 Ac on the column was studied with solutions of various pHs to select the optimal conditions. Two identical isotope generators were eluted for 47 days with over 100 mL of eluant with no detectable 228 Ra breakthrough and high 228 Ac yields (∼95%). The separation requires only biocompatible reagents and is performed at pH 4.6, conditions suitable for radiopharmaceutical studies of Ac.

and nuclear chemistry↗
DOE OSTI2023

Environmental impact of solution pH on the formation and migration of iron colloids in deep subsurface energy systems

Deep subsurface stimulation processes often promote fluid-rock interactions that can lead to the formation of small colloidal particles that are suspected to migrate through the rock matrix, partially or fully clog pores and microfractures, and promote the mobilization of contaminants. Thus, the goal of this work is to understand the geochemical changes of the host rock in response to reservoir stimulation that promote the formation and migration of colloids. Two different carbonate-rich shales were exposed to different solution pHs (pH = 2 and 7). In this study, iron and other mineral transformations at the shale-fluid interface were first characterized by synchrotron-based XRF mapping. Then, colloids that were able to migrate from the shale into the bulk fluid were characterized by synchrotron-based extended X-ray absorption structure (EXAFS), scanning electron microscopy (SEM), and single-particle inductively coupled plasma time-of-flight mass spectrometry (sp-icpTOF-MS). When exposed to the pH = 2 solution, extensive mineral dissolution and secondary precipitation was observed; iron-(oxyhydr)oxide colloids colocated with silicates were observed by SEM at the fluid-shale interfaces, and the mobilization of chromium and nickel with these iron colloids into the bulk fluid was detected by sp-icpTOF-MS. Iron EXAFS spectra of the solution at the shale-fluid interface suggests the rapid (within minutes) formation of ferrihydrite-like nanoparticles. Thus, we demonstrate that the pH neutralization promotes the mobilization of existing silicate minerals and the rapid formation of new iron colloids. These Fe colloids have the potential to migrate through the shale matrix and mobilize other heavy metals (such as Cr and Ni, in this study) and impacting groundwater quality, as well produced waters from these hydraulic fracturing operations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Exploring viral particle, soil, and extraction buffer physicochemical characteristics and their impacts on extractable viral communities

Viruses are expected to be pivotal members of soil ecosystems, and recent advances in viral size fraction metagenomic (viromic) approaches have improved our ability to interrogate soil viral ecology. However, viromics relies on extraction buffers to effectively desorb viral particles from the soil matrix for downstream analysis, and viral extraction efficiency could be affected by the interplay between viral particles, soils, and extraction buffer chemistry. Here, we investigated whether extraction buffer chemistry affected extractable viral community composition measured by viromics from different soil types, for both biological (samples collected 1 m apart) and technical (subsamples from the same soil homogenate) replicates. We investigated pH manipulation of protein-supplemented phosphate-buffered saline buffer (PPBS, pHs 4.5, 5.5, 6.5, and 7.5) on forest, grassland and wetland soils that exhibited different soil edaphic properties, and we tested different buffer chemistries (PPBS, Carbonated Buffer, Glycine, and Saline Magnesium) on just the wetland soil. Spatial distance (i.e., biological replicate) was the primary driver of extractable viral community composition across all buffers and soils tested. Differences in viral community composition according to extraction buffer properties were only observed in the grassland technical replicates at PPBS buffer pH 4.5, and in both the wetland technical and biological replicates treated with different buffer chemistries, but the effects of buffer chemistry were secondary to spatial distance in the biological replicates. The lack of buffering capacity in the grassland soil technical replicates likely increased sorption of viral particles at pH 4.5, but neither protein composition nor isoelectric point explained this phenomenon. In conclusion, given that most soil viral ecological studies to date include sample collection over distances much farther apart than the 1-m distances considered here, results suggest that extraction buffer chemistry is likely of much lower importance than ecological considerations, such as spatial distance, in the design of future soil viral ecological investigations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Speciation and Aqueous Dissolution of Macronutrients in Fire Ash: Variation across Ecosystems and the Effects on Nutrient Cycling

This study investigated the speciation and aqueous dissolution of macronutrients in fire ash from diverse ecosystems and speciation of ash and smoke from laboratory burning, exploring the variations and their causes. The speciation of phosphorus (P), calcium (Ca), and potassium (K) in fire ash from five globally distributed ecosystems was characterized by using X-ray absorption spectroscopy and sequential fractionation. Aqueous dissolution of the macronutrients was measured by batch experiments at acidic and alkaline pHs. The results showed that P existed mainly as Ca phosphates, Ca as double carbonates, calcite, and sulfates, and most K was associated with Ca carbonates. Mineralogy and the relative abundance of the species were primarily controlled by elemental stoichiometry and fire temperature. Differences in Ca and P speciation existed between ash and smoke from laboratory burning, possibly caused by the temperature difference and/or mass fractionation during burning. The rates, extents, and pH dependencies of macronutrient dissolution differed among macronutrients and depended on their speciation, with K being highly soluble and the P and Ca regulated by solution pH. The variability in ash macronutrient chemistry and ecosystemspecific fire ash loads resulted in varying loads and availability of individual macronutrient from fire among ecosystems. This study provides a mechanistic understanding of how fires transform the chemistry of macronutrients and affect macronutrient returns to soils across different ecosystems, which is essential for evaluating the disturbance to ecosystem nutrient cycling by fires.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

pH-Driven Restructuring of Hydration Layers and Cation Adsorption at the Alumina–Water Interface

Oxide−water interfaces underpin ion separations, catalysis, and electrochemical energy technologies, where the electrical double layer (EDL) controls adsorption, transport, and reactivity. However, the molecular-scale links between pHdependent surface protonation, hydration-layer structure, and counterion adsorption remain poorly defined. Here, we combine in situ crystal truncation rod (CTR) and resonant anomalous X-ray reflectivity (RAXR) with streaming potential measurements and ab initio molecular dynamics (AIMD) simulations to resolve the chemical and structural evolution of the EDL at the single-crystal alumina (012)−water interface in 10 mM Rb + over pH 3−12. CTR measurements reveal two distinct adsorbed water layers at ∼2.2 and ∼3.5 Å above the surface. Each water layer shifts toward the substrate at transition pHs near 6.5 and 10.6, respectively, reflecting changes in primary hydration layer structure in response to the deprotonation of bridging and terminal aluminol groups. RAXR shows a 10-fold increase in Rb + coverage and a decrease in mean adsorption height from ∼3.5 to ∼2.7 Å with increasing pH, indicating enhanced counterion binding accompanied by Stern layer contraction. Streaming potential measurements demonstrate that the zeta potential, i.e., the potential at the hydrodynamic shear plane, is positive at pH 3 and becomes negative at pH ≥ 3.5. The negative charge magnitude increases with increasing pH, consistent with progressive surface deprotonation at higher pH. AIMD identifies inner- and outer-sphere Rb+ complexes whose adsorption heights and coordination geometries depend sensitively on the protonation state of surface oxygens, providing atomistic support for the experimentally inferred trends. These measurements establish two discrete, site-specific pH transitions in hydration-layer structure that track aluminol (de)protonation and quantitatively link them to a pH-driven contraction of the Stern layer (increasing Rb + coverage and decreasing adsorption height). This provides a direct structural basis for connecting surface acid−base chemistry to ion binding distances at an oxide−water interface.

Adsorption↗
DOE OSTI2025

Structural determinants for pH-dependent activation of a plant metacaspase

Arabidopsis thaliana metacaspase 9 (AtMC9) plays roles in clearing dead cells, forming xylem vessels, and regulating immunity and programmed cell death in plants. The protease's activation is controlled by pH levels, but the exact structural mechanism behind this has not been elucidated. In this work, high-resolution crystal structures for AtMC9 at active (pH 5.5 and pH 4.2) and inactive (pH 7.5) conditions are reported. The three structures are similar except for local conformations where their hydrogen bonding interactions with solvents are mediated through the protonation of specific titratable amino acid residues' side chains. By combining structural analysis, molecular dynamics simulations under constant pHs, and biochemical assays coupled with site-directed mutagenesis, we show that the regulation of AtMC9 activation involves multiple titratable glutamate and histidine residues across the three domains of p20, linker, and p10. Specifically, deprotonated Glu112, His193, and His208 can suppress AtMC9 proteolytic activity, while protonation of Glu255 and His307 at acidic pH may promote it. This study provides valuable insights into the pH-dependent activation of AtMC9 and could potentially lead to improving crops with enhanced immunity and controlled cell death, ultimately increasing agricultural productivity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Facile electrocatalytic proton reduction by a [Fe–Fe]-hydrogenase bio-inspired synthetic model bearing a terminal CN — ligand

An azadithiolate bridged CN — bound pentacarbonyl bis-iron complex, mimicking the active site of [Fe–Fe] H 2 ase is synthesized. The geometric and electronic structure of this complex is elucidated using a combination of EXAFS analysis, infrared and Mössbauer spectroscopy and DFT calculations. The electrochemical investigations show that complex 1 effectively reduces H + to H 2 between pH 0–3 at diffusion-controlled rates (10 11 M —1 s —1 ) i.e. 10 8 s —1 at pH 3 with an overpotential of 140 mV. Electrochemical analysis and DFT calculations suggests that a CN— ligand increases the p K a of the cluster enabling hydrogen production from its Fe(I)–Fe(0) state at pHs much higher and overpotential much lower than its precursor bis-iron hexacarbonyl model which is active in its Fe(0)–Fe(0) state. The formation of a terminal Fe–H species, evidenced by spectroelectrochemistry in organic solvent, via a rate determining proton coupled electron transfer step and protonation of the adjacent azadithiolate, lowers the kinetic barrier leading to diffusion controlled rates of H 2 evolution. The stereo-electronic factors enhance its catalytic rate by 3 order of magnitude relative to a bis-iron hexacarbonyl precursor at the same pH and potential.

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

Soil Water Retention and Hydraulic Conductivity Data and Model at Pump House in East River Watershed, Colorado 2019-2024

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Pump House at Mount Crested Butte in the East River Watershed. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format ER-X-Y, where ER refers to East River, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, ER-PHS, ER-LMC, ER-LMF, and ER-SMN are associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and ER-RBTn (upslope n=1) are sampling transects during the 2019 Rootball Campaign. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2026

pH-Driven Restructuring of Hydration Layers and Cation Ad-sorption at the Alumina-Water Interface

Oxide-water interfaces underpin ion separation, catalysis, and electrochemical energy technologies, where the electrical double layer (EDL) controls adsorption, transport, and reactivity. Yet, the molecular-scale link between pH-dependent surface protonation, hydration-layer structure, and counter-ion adsorption remains poorly defined. Here, we combine in situ crystal truncation rod (CTR) and resonant anomalous X-ray reflectivity (RAXR) with streaming potential measurements and ab initio molecular dynamics (AIMD) simulations to resolve the chemical and structural evolution of the EDL at the single-crystal alumina (012)-water interface in 10 mM Rb+ over pH 3-12. CTR measurements reveal two distinct adsorbed water layers at ~2.2 and ~3.5 Å above the surface that each shift toward the substrate at transition pHs near 6.5 and 10.6, respectively, directly reflecting changes in primary hydration layer structure in response to the deprotonation of bridging and terminal aluminol groups. RAXR shows a 10-fold increase in Rb+ coverage and a decrease in mean adsorption height from ~3.5 to ~2.7 Å with increasing pH, indicating enhanced counter-ion binding accompanied by Stern layer contraction. Streaming potential measurements demonstrate that the zeta potential, i.e., potential at the hydrodynamic shear plane, is positive at pH 3 and becomes negative at pH ≥3.5. This negative charge magnitude increases with increasing pH, consistent with progressive surface deprotonation at higher pH. AIMD identifies inner- and outer-sphere Rb+ complexes whose adsorption heights and coordination geometries depend sensitively on the protonation state of surface oxygens, providing atomistic support for the experimentally inferred trends. These measurements establish two discrete, site-specific pH transitions in hydration-layer structure that track aluminol (de)protonation and quantitatively link them to a pH-driven contraction of the Stern layer (increasing Rb+ coverage and decreasing adsorption height). This provides a direct structural basis for connecting surface acid-base chemistry to ion binding distances at an oxide-water interface.

Electrical double layer (EDL), Surface protonation↗
DOE OSTI2024

Structural characterization of ligand binding and pH-specific enzymatic activity of mouse Acidic Mammalian Chitinase

Chitin is an abundant biopolymer and pathogen-associated molecular pattern that stimulates a host innate immune response. Mammals express chitin-binding and chitin-degrading proteins to remove chitin from the body. One of these proteins, Acidic Mammalian Chitinase (AMCase), is an enzyme known for its ability to function under acidic conditions in the stomach but is also active in tissues with more neutral pHs, such as the lung. Here, we used a combination of biochemical, structural, and computational modeling approaches to examine how the mouse homolog (mAMCase) can act in both acidic and neutral environments. We measured kinetic properties of mAMCase activity across a broad pH range, quantifying its unusual dual activity optima at pH 2 and 7. We also solved high-resolution crystal structures of mAMCase in complex with oligomeric GlcNAcn, the building block of chitin, where we identified extensive conformational ligand heterogeneity. Leveraging these data, we conducted molecular dynamics simulations that suggest how a key catalytic residue could be protonated via distinct mechanisms in each of the two environmental pH ranges. These results integrate structural, biochemical, and computational approaches to deliver a more complete understanding of the catalytic mechanism governing mAMCase activity at different pH. Engineering proteins with tunable pH optima may provide new opportunities to develop improved enzyme variants, including AMCase, for therapeutic purposes in chitin degradation.

59 BASIC BIOLOGICAL SCIENCES↗