Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface waters”

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

Advanced Anti-Fouling Coatings to Improve the Efficiency of Coal Power Plants

As the total cost of carbon to generate energy has become a global concern, operators are increasingly looking at all parts of the generation cycle to find areas where efficiency gains may be found. It has been long identified that fouling of heat exchangers is a persistent cause of up to 2.5% of global CO 2 emissions. Unfortunately, practice has also demonstrated that unless a powerful economic driver exists to encourage preemptive mitigation of fouling, there will always be a strong tendency for operators to minimize any form of intervention due to high costs and challenges in scheduling downtime. The objective of this proposed research effort was to demonstrate how existing power plants could lower their carbon emissions and significantly improve heat transfer efficiency using new surface treatment materials to control fouling in a variety of heat exchange equipment. The surface treatment material which was optimized and deployed in this effort is now known commercially as HeatX. It is a low-surface energy, water- and oil-repellent, abrasion resistant material which can be applied in-situ to a wide variety of previously worn/used/in-service substrates. Once applied, it provides a barrier against corrosion, scale deposit formation, and biofilm adhesion on the circulating water-containing tube-side. Alternatively, if applied to the tube exterior, the non-wetting nature of the surface was demonstrated to promote dropwise condensation, subsequently lowering condenser backpressure and increasing overall plant efficiency. As part of this cooperative effort, the Department of Energy’s support was crucial to de-risk and demonstrate the concept of HeatX, while validating both the performance and economic benefit in multiple pilot field studies. The HeatX material properties were optimized in this effort for full field applicability to heat exchangers and condensers, and Oceanit developed the necessary procedures and protocols to provide enough material to support extended length, multi-year demonstrations in the power generation, desalination, and refining industries, making this technology broadly applicable and ready for commercial transition. Field deployment case studies at thermal power plants have shown that the HeatX treatment can provide economic savings of up to $15,000 per day for an operator based on avoiding maintenance costs and lowering fuel usage. The complete mitigation of fouling effects can increase the efficiency of equipment by up to 7%, in a field where gains of 0.5% are generally seen as operationally significant. The HeatX treatment has also demonstrated exceptional lifetime and compatibility with a wide variety of seawater and hydrocarbon environments, further increasing both the return on investment and the effective emission reduction. Such efficiency improvements correlate to massive carbon savings. For every 1 GW of capacity, operators can see carbon emissions reductions of 300,000 tons of CO 2 per year. When looking at the bigger picture, improved condenser function across the U.S. has the potential to prevent 221.3 million tons of CO 2 emissions, equivalent to the sequestration capacity of 129 million acres of forest. If applied on a global scale, 1.26 billion metric tons of CO 2 could be averted from the atmosphere, the same amount of carbon sequestrated by 1.5 billion acres of forest annually or 262,000 wind turbines operating annually. As businesses across multiple industries take a more active role in focusing on environmental, social, and governance (ESG) solutions as part of their core business operations, HeatX will be an attractive technology for commercial investment.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Multivariate environmental and trait-based controls of transpiration in the Central Amazon Rainforest

Tropical forest tree mortality is increasing due to more severe droughts, yet our understanding of how tree traits and life strategies are linked to drought stress has been limited by measurement scarcity. The BIONTE (BIOmass and NuTrient Experiment) near Manaus, Brazil hosts one of the world’s largest sap flow installations, with sensors in 90 canopy trees across a wood density gradient monitored since June 2022. The 2023 El Niño drought provided a unique opportunity to evaluate how water availability impacts tree transpiration. An interpretable machine learning framework was used to study the complex interactions between transpiration and multiple environmental variables such as soil water availability and vapor pressure deficit (VPD), and how these interactions vary with wood density and individual trees. We found varying responses of transpiration from different trees during the El Niño drought. Transpiration generally increased with temperature, with stronger effects in wetter areas and in trees with low to medium wood density. However, this response was modulated by stomatal sensitivity to VPD, which constrained transpiration under high atmospheric demand, particularly in intermediate-moisture area. The inflection in transpiration rate at high temperatures (>32°C) underscores the role of stomatal and hydraulic regulation in limiting water loss and protecting trees from excessive evaporative demand. Analysis of soil water contribution to transpiration revealed unimodal patterns in wetter area, with peak contributions near 0.45 cm 3 cm -3 of surface soil water and declining or flat responses beyond that threshold, suggesting a shift from water- to energy-limited transpiration. In contrast, drier areas exhibited limited transpiration sensitivity to soil water conditions and minimal trait-based variation in VPD responses, indicating supply-limited conditions. Despite higher wood density trees being generally more resilient, this study shows diverse tree drought resilience, prompting further investigation into the specific traits and dynamics between environmental variables in regulating transpiration and other physiological processes in trees.

Drought↗
DOE OSTI2023

Insights on Simulating Summer Warming of the Great Lakes: Understanding the Behavior of a Newly Developed Coupled Lake-Atmosphere Modeling System

The Laurentian Great Lakes are the world's largest freshwater system and regulate the climate of the Great Lakes region, which has been increasingly experiencing climatic, hydrological, and ecological changes. An accurate mechanistic representation of the Great Lakes thermal structure in Regional Climate Models (RCMs) is paramount to studying the climate of this region. Currently, RCMs have primarily represented the Great Lakes through coupled one-dimensional (1D) column lake models; this approach works well for small inland lakes but is unable to resolve the realistic hydrodynamics of the Great Lakes and leads to inaccurate representations of lake surface temperature (LST) that influence regional climate and weather patterns. This work overcomes this limitation by developing a fully two-way coupled modeling system using the Weather Research and Forecasting model and a three-dimensional (3D) hydrodynamic model. The coupled model system resolves the interactive physical processes between the atmosphere, lake, and surrounding watersheds; and validated against a range of observational data. The model is then used to investigate the potential impacts of lake-atmosphere coupling on the simulated summer LST of Lake Superior. By evaluating the difference between our two-way coupled modeling system and our observation-driven modeling system, we find that coupled-lake atmosphere dynamics can lead to a higher LST during June-September through higher net surface heat flux entering the lake in June and July and a lower net surface heat flux entering the lake in August and September. The unstratified water in June distributes the entering surface heat flux throughout the water column leading to a minor LST increase, while the stratified waters of July create a conducive thermal structure for the water surface to warm rapidly under the higher incoming surface heat flux. This research provides insight into the coupled modeling system behavior, which is critical for enhancing our predictive understanding of the Great Lakes climate system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Relationships between Water’s Structure and Solute Affinity at Polypeptoid Brush Surfaces

Excellent antifouling surfaces are generally thought to create a tightly bound layer of water that resists solute adsorption, and highly hydrophilic surfaces such as those with zwitterionic functionalities are of significant current interest as antifoulant strategies. Furthermore, despite significant proofs-of-concept, we still lack a fundamental understanding of how the nanoscopic structure of this hydration layer translates to reduced fouling, how surface chemistry can be tuned to achieve antifouling through hydration water, and why, in particular, zwitterionic surfaces seem so promising. Here, we use molecular dynamics simulations and free energy calculations to investigate the molecular relationships among surface chemistry, hydration water structure, and surface–solute affinity across a variety of surface-decorated chemistries. Specifically, we consider polypeptoid-decorated surfaces that display well-known experimental antifouling capabilities and that can be synthesized sequence specifically, with precise backbone positioning of, e.g., charged groups. Through simulations, we calculate the affinities of a range of small solutes to polypeptoid brush surfaces of varied side-chain chemistries. We then demonstrate that measures of the structure of surface hydration water in response to a particular surface chemistry signal solute–surface affinity; specifically, we find that zwitterionic chemistries produce solute–surface repulsion through highly coordinated hydration water while suppressing tetrahedral structuring around the solute, in contrast to uncharged surfaces that show solute–surface affinity. Based on the relationship of this structural perturbation to the affinity of small-molecule solutes, we propose a molecular mechanism by which zwitterionic surface chemistries enhance solute repulsion, with broader implications for the design of antifouling surfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Decontamination of urban surfaces contaminated with radioactive materials and consequent onsite recycling of the waste water

Enhancing rapid remediation strategies is paramount for recovery after a large-scale nuclear contamination event in an urban environment. Some current strategies recommend use of readily available equipment, materials, and facilities to expedite recovery. For example, applying pressurized water to contaminated surfaces may effectively remove radioactive contamination. In this study, a commercial power washer removes soluble forms of 152 Eu 3 + , 85 Sr 2 + , and 137 Cs + contamination from common porous building materials, and computer simulations characterize the recycling of the resultant contaminated wash water. Pressure washing the porous building materials under spray conditions typical with do-it-yourself units improved decontamination factors (DFs) for 152 Eu compared to low-pressure application of tap water (majority of two-tailed t-test p-values < 0.1), but pressure did not improve DFs for 137 Cs or 85 Sr. For both pressurized and low-pressure applications, adding potassium ions (K + ) to promote ion exchange reactions produced significantly higher DFs for tested radionuclides on asphalt, brick, and concrete. The resultant contaminated wash water can be processed through self-prepared chemical filtration beds of clay and sand. Modeled in a prior study, the beds yielded linear trends (R 2 > 0.98) in sensitivity analyses between most bed configuration variables and bed performance variables, permitting flexible ad-hoc bed design. The experimental and simulation results led to estimates of the remediation rate and waste generated after cleaning 250 m 2 of cesium-contaminated concrete from the combined deployment of a power washer and two different mobile treatment beds. Furthermore, the first treatment bed was designed to reduce treatment time and processed 1900 L of wash solution in 70 min using 880 kg of clay/sand infill material. Designed to reduce the solid waste generated, the second bed processed the same solution volume in 1040 min (17 h) using 170 kg of clay/sand infill material. The results of this analysis warrant further investigation of power washing with recycled salt solution as an effective rapid decontamination method with manageable waste.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

DECOVALEX-2023: Task F1 Final Report

DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75 th percentile for all benchmarks. At the 90 th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Daily water stable isotopes, transpiration, and matrix potential data for an aspen and engelmann stand in the East River Watershed (version 2)

We provide daily stable isotope (2H & 18O) ratios in soil water and xylem (plant stem) water, as well as the sap flow (transpiration) and the soil's matric potential at a forested site near Gothic, Colorado, in the East River catchment. We measured the stable isotopic composition of the transpiration and the daily transpiration flux sum of three aspen and three engelmann spruce. In both forest stands, we installed a soil profile and measured the soil matric potential at 15, 30, and 60 cm depth as well as the stable isotopes of soil pore water at 5, 10, 30, 60, and 90 cm depths. All isotope measurements were done in situ via vapor probes connected to a cavity ring down spectrometer (Picarro L1240i).We further report the daily meteorological data observed at billy barr near our study site. We also provide for each tree the relative share of root water uptake derived from the isotope measurements via a Bayesian mixing model (MixSIAR).The daily data is provided as a time series in "Iso_MP_Sap_DataDaily_ESSDiveUpload.csv" and the units are provided in "dd.csv"; the location of the instrumented trees and soil profiles are given as latitude and longitude coordinates saved as CSV and KMZ files; and a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata.The data was gathered to investigate the short-term changes of the water sources (i.e., variation of root water uptake from different soil depths) of the studied subalpine trees.Update 07/16/2025: The relative and absolute plant water uptake depths were grouped to ensure that the MixSIAR model was applied with endmembers that differed in their d2H value by at least 3 permill and at least 1 permill for d18O. Whenever the difference between observed d2H values for two or more probes at neighboring depths was less than 3 permill, we used the average value for the source water endmember. For days at which probes that were not next to each other measurements did not differ at least 3 permill, the average of all probes between these two depths was used as the water source endmember representing the depth range between these two probes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Fine Aerosol Acidity and Water during Summer in the Eastern North Atlantic

Aerosol pH governs many important atmospheric processes that occur in the marine boundary layer, including regulating halogen and sulfur chemistries, and nutrient fertilization of surface ocean waters. In this study, we investigated the acidity of PM1 over the eastern North Atlantic during the Aerosol and Cloud Experiments in Eastern North Atlantic (ACE-ENA) aircraft campaign. The ISORROPIA-II thermodynamic model was used to predict PM1 pH and water. We first investigated the sensitivities of PM1 pH and water predictions to gas-phase NH3 and HNO3 concentrations. Our sensitivity analysis indicated that even though NH3 and HNO3 were present at very low concentrations in the eastern North Atlantic during the campaign, PM1 pH calculations can still be sensitive to NH3 concentrations. Specifically, NH3 was needed to constrain the pH of populations of PM1 that had low mass concentrations of NH4+ and non-volatile cations (NVCs). We next assumed that gas-phase NH3 and HNO3 concentrations during the campaign were 0.15 and 0.09 µg m−3, respectively, based on previous measurements conducted in the eastern North Atlantic. Using the assumption that PM1 were internally mixed (i.e., bulk PM1), we determined that PM1 pH ranged from 0.3–8.6, with a mean pH of 5.0 ± 2.3. The pH depended on both Hair+ and Wi. Hair+ was controlled primarily by the NVCs/SO42− molar ratio, while Wi was controlled by the SO42− mass concentration and RH. Changes in pH with altitude were driven primarily by changes in SO42−. Since aerosols in marine atmospheres are rarely internally mixed, the scenario where non-sea salt species and sea-salt species were present in two separate aerosol modes in the PM1 (i.e., completely externally mixed) was also considered. Smaller pH values were predicted for the aerosol mode comprised only of non-sea salt species compared to the bulk PM1 (difference of around 1 unit on average). This was due to the exclusion of sea-salt species (especially hygroscopic alkaline NVCs) in this aerosol mode, which led to increases in Hair+ values and decreases in Wi values. This result demonstrated that assumptions of aerosol mixing states can impact aerosol pH predictions substantially, which will have important implications for evaluating the nature and magnitude of pH-dependent atmospheric processes that occur in the marine boundary layer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Molecular Dynamics Simulation and Cryo-Electron Microscopy Investigation of AOT Surfactant Structure at the Hydrated Mica Surface

Structural properties of the anionic surfactant dioctyl sodium sulfosuccinate (AOT or Aerosol-OT) adsorbed on the mica surface were investigated by molecular dynamics simulation, including the effect of surface loading in the presence of monovalent and divalent cations. The simulations confirmed recent neutron reflectivity experiments that revealed the binding of anionic surfactant to the negatively charged surface via adsorbed cations. At low loading, cylindrical micelles formed on the surface, with sulfate head groups bound to the surface by water molecules or adsorbed cations. Cation bridging was observed in the presence of weakly hydrating monovalent cations, while sulfate groups interacted with strongly hydrating divalent cations through water bridges. The adsorbed micelle structure was confirmed experimentally with cryogenic electronic microscopy, which revealed micelles approximately 2 nm in diameter at the basal surface. At higher AOT loading, the simulations reveal adsorbed bilayers with similar surface binding mechanisms. Adsorbed micelles were slightly thicker (2.2–3.0 nm) than the corresponding bilayers (2.0–2.4 nm). Upon heating the low loading systems from 300 K to 350 K, the adsorbed micelles transformed to a more planar configuration resembling bilayers. The driving force for this transition is an increase in the number of sulfate head groups interacting directly with adsorbed cations.

58 GEOSCIENCES↗
DOE OSTI2022

Reorganization Energies for Interfacial Proton-Coupled Electron Transfer to a Water Oxidation Catalyst

The reorganization energy (λ) for interfacial electron transfer (ET) and proton-coupled electron transfer (PCET) from a conductive metal oxide (In 2 O 3 :Sn, ITO) to a surface-bound water oxidation catalyst were extracted from kinetic data measured as a function of the thermodynamic driving force. Visible light excitation resulted in rapid excited-state injection (k inj > 10 8 s -1 ) to the ITO, which photo-initiated the two interfacial reactions of interest. The rate constants for both reactions increased with the driving force, −ΔG°, to a saturating limit, k max , with rate constants consistently larger for ET than for PCET. Marcus-Gerischer analysis of the kinetic data provided the reorganization energy for interfacial PCET (0.90 ± 0.02 eV) and ET (0.40 ± 0.02 eV), respectively. The magnitude of k max for PCET was found to decrease with pH, behavior that was absent for ET. Both the decrease in k max and the larger reorganization energy for an unwanted competing PCET reaction from the ITO to the oxidized catalyst showcases a significant kinetic advantage for driving solar water oxidation at high pH. Computational analysis revealed a larger inner-sphere reorganization energy contribution for PCET than for ET arising from a more significant change in the Ru–O bond length for PCET. Extending Marcus-Gerischer theory to PCET by including the excited electron-proton vibronic states and the proton donor-acceptor motion provided an apparent reorganization energy of 1.01 eV. Furthermore, this study demonstrates that the Marcus-Gerischer theory initially developed for ET can be reliably extended to PCET for quantifying and interpreting reorganization energies observed experimentally.

Catalysts Kinetic parameters↗
DOE OSTI2022

Superhydrophobic to superhydrophilic wettability transition of functionalized SiO 2 nanoparticles

The superhydrophobic and superhydrophilic surfaces and their transitions are of great interest for the production of self-cleaning, anti-biofouling, or corrosion-resistant materials. This research reports the wettability transition from superhydrophobic to superhydrophilic SiO 2 nanoparticles functionalized with 1H, 1H, 2H, 2H-perfluorooctyltriethoxysilane (POTS) and induced by temperature. The functionalization of these nanoparticles was confirmed by Fourier transform infrared spectroscopy, X-ray photoelectron spectroscopy (XPS), and transmission electron microscopy. The functionalization of SiO 2 nanoparticles with POTS resulted in superhydrophobic surfaces with water contact angles up to 157°. A sudden transition to superhydrophilic behavior with water contact angles (WCA) below 5° was observed when the sample was heat-treated at 500 °C, despite the presence of fluorine on the surface of these nanoparticles, as confirmed by XPS and transmission electron microscopy. XPS suggested that the transition was caused by the change in orientation of the fluoroalkyl molecules and its partial decomposition due to the loss of the –CF 3 group, resulting in shorter chains with a tail-end group with C–O bonds, which promoted the superhydrophilicity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Facet-dependent structure and dissociation of water at pristine IrO 2 /water interfaces

Understanding the microscopic structure of water at metal oxide interfaces is crucial for advancing electrocatalysis. IrO 2 , specifically, has shown exceptional activity for electrochemical water oxidation, but we currently lack a fundamental understanding of how the surface structure of IrO 2 impacts water reactivity. In this work, we developed a machine learning potential trained to first-principles accuracy for modeling IrO 2 /water interfaces across different facets: (110), (100), (101), and (001). Using extensive machine learning molecular dynamics simulations, we investigated the spontaneous dissociation of water molecules at these interfaces. Our results reveal a distinct dissociation probability trend: (110) > (100) ≈ (101) > (001), which we attribute primarily to the reaction thermodynamics of surface water dissociation. A strong correlation is observed between the surface Ir–O bond distances and the dissociation probabilities, highlighting the role of surface geometry in modulating reactivity. As a consequence, the interfacial solvation structures and hydrogen bonding environments are dynamically tuned by the varying water dissociation capabilities across facets. This work elucidates how water dissociation energetics depend on surface orientation and interfacial structure, offering atomistic insights into manipulating reaction chemistry at electrocatalytic interfaces.

organic↗
DOE OSTI2025

Unusual Electrochemical Activity of Thin SiO 2 Layers Leads to Instability of Molecular Attachment in Hybrid Photoelectrodes

Hybrid photoelectrodes, comprised of a light-absorbing semiconductor and a surface-integrated molecular catalyst, are attractive for applications in artificial photosynthesis, since they combine the advantages of broadband semiconductor light absorption with the selectivity of molecular catalysis. A widely used class of hybrid photoelectrodes is based on Si substrates passivated by a thin (<3 nm) layer of silicon oxide, which is commonly prepared by controlled chemical or thermal oxidation, resulting in chemical oxide (ChO) or thermal oxide (ThO) layers, respectively. However, the electrochemical stability of these oxide layers, and the chemical stability of the semiconductor-molecule assembly in hybrid photoelectrodes, are not well understood, with evidence that covalently-bound molecules detach from the oxide surface upon application of cathodic bias. We have examined the intrinsic electrochemical reactivity of silicon oxide layers and how it affects the attachment of molecular monolayers. We determined that the surface of Si|ThO is primarily terminated with hydrophobic siloxane moieties, whereas that of Si|ChO contains a higher concentration of hydrophilic silanol groups. Initial high current densities for Si|ChO under applied bias up to -2 V vs. Ag/AgCl, decrease during repeated cyclic voltammetry scans, due to the consumption of surface-bound water. This is manifested by a reversible wave around -0.5 V in CH 3 CN solution, and a similar pH-dependent wave in water, revealing the pK a of the silanol groups to be ~4. Here, our combined observations support the electrochemically-induced dehydration of the SiO 2 surface, which converts silanol groups to siloxanes and proceeds through an H-atom intermediate that is most likely stabilized by pentavalent Si. We propose that similar reactivity is responsible for the electrochemical loss of alkylsiloxane-attached molecules under cathodic bias, which has important implications for the choice of catalyst attachment strategy in hybrid photoelectrodes.

14 SOLAR ENERGY↗
DOE OSTI2021

Structure and Stability of Pristine and Carboxylate-Covered Anatase TiO 2 (001) in Aqueous Environment

The interactions of carboxylic acids and water with TiO 2 surfaces are important in applications ranging from solar cells to biomedical devices. Here we focus on the aqueous interface with the minority (001) surface of anatase TiO 2 (A-001) and the behavior of adsorbed formic and acetic acid monolayers at this interface. We investigated the structure and stability of the pristine reconstructed and formic/acetic acid covered A-001 surfaces in contact with water using density functional theory (DFT) calculations and ab initio molecular dynamics (AIMD) simulations. The (1 × 4) reconstruction of the pristine surface is found to be stable in aqueous environment, within the time scale of our simulation. Here, carboxylic acids adsorb in deprotonated bidentate (BD) form on A-001, with the dissociated proton transferred to a surface oxygen to form a bridging hydroxyl. Of the two possible configurations, BD bridging and BD straddling, of the adsorbed species, the latter is found to rapidly transform to a monodentate structure during our simulations. Further investigation of mixed acetate–formate monolayers on A-001 in water indicates that also BD bridging species can become unstable at the boundaries between formate and acetate-covered regions, transforming to a monodentate form that does not prevent water adsorption on the TiO 2 surface.

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

Onset dynamics of air-water menisci on rock fracture surfaces

Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.

58 GEOSCIENCES↗
DOE OSTI2022

On the additivity of climate responses to the volcanic and solar forcing in the early 19th century

The early 19th century was the coldest period over the past 500 years, when strong tropical volcanic events and a solar minimum coincided. The 1809 unidentified eruption and the 1815 Tambora eruption happened consecutively during the Dalton minimum of solar irradiance; however, the relative role of the two forcing (volcano and solar) agents is still unclear. In this study, we examine the responses from a set of early 19th century simulations with combined and separated volcanic and solar forcing agents, as suggested in the protocol for the past1000 experiment of the Paleoclimate Modelling Intercomparison Project – Phase 4 (PMIP4). From 20-member ensemble simulations with the Max Planck Institute Earth system model (MPI-ESM1.2-LR), we find that the volcano- and solar-induced surface cooling is additive in the global mean/large scale, regardless of combining or separating the forcing agents. The two solar reconstructions (SATIRE (Spectral and Total Irradiance REconstruction-Millennia model) and PMOD (Physikalisch-Meteorologisches Observatorium Davos)) contribute to a cooling before and after 1815 of ~0.05 and ~0.15 K monthly average near-surface air cooling, respectively, indicating a limited solar contribution to the early 19th century cold period. The volcanic events provide the main cooling contributions, inducing a surface cooling that peaks at ~0.82 K for the 1809 event and ~1.35 K for Tambora. After the Tambora eruption, the temperature in most regions increases toward climatology largely within 5 years, along with the reduction of volcanic forcing. In the northern extratropical oceans, the temperature increases slowly at a constant rate until 1830, which is related to the reduction of seasonality and the concurrent changes in Arctic sea-ice extent. The albedo feedback of Arctic sea ice is found to be the main contributor to the Arctic amplification of the cooling signal. Several non-additive responses to solar and volcanic forcing happen on regional scales. In the atmosphere, the stratospheric polar vortex tends to strengthen when combining both volcano and solar forcing, even though the two forcing agents separately induce opposite-sign changes in stratospheric temperatures and zonal winds. In the ocean, when combining the two forcings, additional surface cold water propagates to the northern extratropics from the additional solar cooling in the tropics, which results in regional cooling along the propagation. Overall, this study not only quantifies the surface responses from combinations of the volcano and solar forcing, but also highlights the components that cannot be simply added from the responses of the individual forcing agents, indicating that a relatively small forcing agent (such as solar in early 19th century) can impact the response from the large forcing (such as the 1815 Tambora eruption) when considering regional climates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Molecular Structure and Thermodynamics of CO 2 and Water Adsorption on Mica

The adsorption of CO 2 and water on clay surfaces plays a key role in applications, such as gas storage in saline aquifers and depleted hydrocarbon reservoirs, but is not yet fully understood. Here, we study the adsorption of CO 2 and water vapor using Grand Canonical Monte Carlo and molecular dynamics simulations. At a bulk pressure of 100 bar, pure CO 2 adsorbs strongly on mica and forms extensive layers next to it. CO 2 adsorption is lowered substantially if introducing water vapor above mica and is largely eliminated when the relative humidity (RH) approaches about 60%. When pure water vapor is introduced above a mica surface, a subnanometer thick liquid water film develops on it to form apparent liquid–solid and liquid–vapor interfaces simultaneously. Using the identification of truly interfacial molecules (ITIM) analysis, we delineate how individual water layers develop in this film as RH increases. We highlight that the water film is spatially heterogeneous and the true liquid–vapor interface emerges only at an RH of 60–80%. Introducing 100 bar of CO 2 into the water vapor above the mica surface modulates water adsorption nonlinearly: at RH = 0.01%, the water adsorption is reduced by ∼30%; as RH increases, the reduction is weakened, and eventually, enhancement of water adsorption by about 7% occurs at RH = 90%. These variations are attributed to the interplay of film thinning by high-pressure CO 2 , competition of mica surface sites by CO 2 molecules, and energetic and entropic stabilization of interfacial water by CO 2 molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗