Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface level”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Observations of surface energy fluxes and meteorology in the seasonally snow-covered high-elevation East River watershed during SPLASH, 2021–2023

Abstract. From autumn 2021 through summer 2023, scientists from the National Oceanic and Atmospheric Administration (NOAA) and partners conducted the Study of Precipitation, the Lower Atmosphere, and Surface for Hydrometeorology (SPLASH) campaign in the East River watershed of Colorado. One objective of SPLASH was to observe the transfer of energy between the atmosphere and the surface, which was done at several locations. Two remote sites were chosen that did not have access to power utilities. These were along the valley floor near the East River in the vicinity of the unincorporated town of Gothic, Colorado. Energy balance measurements were made at these locations using autonomous, single-level flux towers referred to as atmospheric surface flux stations (ASFSs). The ASFSs were deployed on 28 September 2021 at the Kettle Ponds Annex site and on 12 October 2021 at the Avery Picnic site and operated until 19 July and 21 June 2023, respectively. Measurements included basic meteorology; upward and downward longwave and shortwave radiative fluxes and subsurface conductive flux, each at 1 min resolution; 3-D winds from a sonic anemometer and H2O/CO2 from an open-path gas analyzer, both at 20 Hz from which sensible, latent heat, and CO2 fluxes were derived; and profiles of soil properties in the upper 0.5 m (both sites) and temperature profiles through the snow (at Avery Picnic), each reported between 10 min and 6 h. The system uptime was 97 % (Kettle Ponds) and 90 % (Avery Picnic), and collectively 1184 d of data was obtained between the stations. The purpose of this article is to document the ASFS deployment at SPLASH, to document the data acquisition and post-processing of measurements, and to serve as a guide for interested users of the data sets, which are archived at Zenodo (https://doi.org/10.5281/zenodo.10313363, Cox et al., 2023b; https://doi.org/10.5281/zenodo.10327409, Cox et al., 2023c; https://doi.org/10.5281/zenodo.10313894, Cox et al., 2023d; https://doi.org/10.5281/zenodo.10307825, Cox et al., 2023e; https://doi.org/10.5281/zenodo.10310520, Cox et al., 2023f) with the Creative Commons Attribution 4.0 International license.

Cox, Christopher J. (ORCID:0000000322037173)↗
DOE OSTI2022

Can fractal dimensions objectivize gastropod shell morphometrics? A case study from Lake Lugu (SW China)

1. Morphometrics are fundamental for the analysis of size and shape in fossils, particularly because soft parts or DNA are rarely preserved and hard parts such as shells are commonly the only source of information. Geometric morphometrics, that is, landmark analysis, is well established for the description of shape but it exhibits a couple of shortcomings resulting from subjective choices during landmarking (number and position of landmarks) and from difficulties in resolving shape at the level of micro-sculpture. 2. With the aid of high-resolution 3D scanning technology and analyses of fractal dimensions, we test whether such shortcomings of linear and landmark morphometrics can be overcome. As a model group, we selected a clade of modern viviparid gastropods from Lake Lugu, with shells that show a high degree of sculptural variation. Linear and landmark analyses were applied to the same shells in order to establish the fractal dimensions. 3. The genetic diversity of the gastropod clade was assessed. The genetic results suggest that the gastropod clade represents a single species. The results of all morphometric methods applied are in line with the genetic results, which is that no specific morphotype could be delimited. Apart from this overall agreement, landmark and fractal dimension analyses do not correspond to each other but represent data sets with different information. Generally, the fractal dimension values quantify the roughness of the shell surface, the resolution of the 3D scans determining the level. In our approach, we captured the micro-sculpture but not the first-order sculptural elements, which explains that fractal dimension and landmark data are not in phase. 4. We can show that analyzing fractal dimensions of gastropod shells opens a window to more detailed information that can be considered in evolutionary and ecological contexts. We propose that using low-resolution 3D scans may successfully substitute landmark analyses because it overcomes the subjective landmarking. Analyses of 3D scans with higher resolution than used in this study will provide surface roughness information at the mineralogical level. We suggest that fractal dimension analyses of a combination of differently resolved 3D models will significantly improve the quality of shell morphometrics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Quantitative Effects of Disorder on Chemically Modified Amorphous Carbon Electrodes

Real materials are disordered. This disorder influences the properties of these materials and the chemical processes that occur at their interfaces. Gaining a molecular-level understanding of the underlying physical manifestations caused by disordered materials is crucial to unraveling and ultimately controlling the efficiency and performance of these materials in a range of energy-related devices. This understanding necessitates measurement techniques through which disorder can be detected, quantified, and monitored. However, such quantitative measurements are notoriously difficult, as effects often average out in ensemble measurements. In this work, we describe how a combination of electrochemical and spatially resolved surface spectroscopy measurements illuminate a molecular-level picture of disorder in materials. Using amorphous carbon as an intrinsically disordered material, we covalently attached a monolayer of ferrocene. Interfacial electron transfer across the amorphous carbon–ferrocene interface is highly sensitive to disruptions of order. Here, by systematically varying linker properties and surface loadings, the influence of lateral interactions between nonuniformly distributed ferrocene headgroups on ensemble electrochemical measurements is demonstrated. Electrochemical and imaging data collectively indicate that conformational flexibility of the ferrocene moieties provides a mechanism to elude repulsive and unbalanced lateral interactions, while rigid linkages provide direct information about the underlying disorder of the material. This study is the first of its kind to quantify and visualize molecular disorder and heterogeneity with an experimental model accessed through ensemble measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Temporal Variability in Reservoir Surface Area Is an Important Source of Uncertainty in GHG Emission Estimates

Ebullitive methane (CH 4 ) emissions in lentic ecosystems tend to concentrate at river-lake interfaces and within shallow littoral zones. However, inconsistent definitions of the littoral zone and static representations of the lake or reservoir surface area contribute to major uncertainties in greenhouse gas (GHG) emissions estimates, particularly in reservoirs with large water-level fluctuations. This study examines temporal variation in littoral and total surface areas of US reservoirs and demonstrates how different methods and data sources lead to discrepencies in reservoir GHG emissions at large scales and over time. We also explore variability in remotely sensed water occurrence according to maximum surface area, reservoir purposes, and hydrologic regions. Notably, the largest relative variability in surface area is exhibited by small reservoirs with a maximum surface area <1 km 2 and non-hydroelectric reservoirs. Additionally, we use a case study of measured CH 4 emissions from the southeastern United States (Douglas Reservoir) to illustrate the effects of varying surface area on reservoir-wide GHG estimates. Upscaled CH 4 emissions in Douglas Reservoir differed by nearly two-fold depending on the source of total surface area data and whether estimates accounted for seasonal fluctuations in surface area. During seasonal drawdown in Douglas Reservoir, relative littoral area varies non-linearly; periods of lower pool elevation (and thus larger relative littoral area) likely contribute disproportionately high CH 4 emission rates compared to the commonly sampled summer season when water levels are at full-pool elevation. Improved GHG monitoring and upscaling techniques require accounting for temporal variability in reservoir surface extent and littoral area.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Neural Network Sampling of the Free Energy Landscape for Nitrogen Dissociation on Ruthenium

In heterogeneous catalysis, free energy profiles of reactions govern the mechanisms, rates, and equilibria. Energetics are conventionally computed using the harmonic approximation (HA), which requires determination of critical states a priori. Here, we use neural networks to efficiently sample and directly calculate the free energy surface (FES) of a prototypical heterogeneous catalysis reaction—the dissociation of molecular nitrogen on ruthenium—at density-functional-theory-level accuracy. We find that the vibrational entropy of surface atoms, often neglected in HA for transition metal catalysts, contributes significantly to the reaction barrier. The minimum free energy path for dissociation reveals an “on-top” adsorbed molecular state prior to the transition state. While a previously reported flat-lying molecular metastable state can be identified in the potential energy surface, it is absent in the FES at relevant reaction temperatures. Finally, these findings demonstrate the importance of identifying critical points self-consistently on the FES for reactions that involve considerable entropic effects.

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

Remote Influence of Andean Convection on Amazonian Rainfall and Its Mechanisms

Models from Coupled Model Intercomparison Project Phase 6 produce too much precipitation over the Andes but too little over the Amazon or the Wet Andes-Dry Amazon (WADA) bias pattern. Unlike the conventional view that convection parameterization and land model deficiencies can contribute to Amazonian rainfall biases, we approach this long-standing biased model behavior through the lens of Andean convection. Using Community Earth System Model v1.1 and focusing on the wet season, our mechanism-denial experiments demonstrate that Andean convection notably reduces precipitation over the Amazon during austral summer. The Andean forced Amazonian response operates on weather timescale. Furthermore, the reduction of Amazonian rainfall is detectable within a few hours after initial Andean forcing. The precipitation response is primarily driven by variations in the moisture budget and is moderated by changes in convective available potential energy over the Amazon. Changes in the total advection of moisture over the Amazon are dominated by the vertical advection term and can be attributed to discrepancies in the dynamic omega field. In the experiments, the Andean east flank region is scrutinized where the vertical velocity and moisture fields play an intermediary role for the Andean driven WADA connection. The Andean forcing induces descending anomalies on the Andean east flank. The disturbances of wind and geopotential fields over the Andean east flank propagate eastward via Kelvin waves. Over the Amazon, descending anomalies and advective drying lead to reduction of mid-to-high level cloud, increase of shortwave cloud forcing and surface net radiation, and enhancement of themodynamic stability and rainfall reduction.

58 GEOSCIENCES↗
DOE OSTI2021

The Influence of Ocean Topography on the Upwelling of Carbon in the Southern Ocean

The physical circulation of the Southern Ocean sets the surface concentration and thus air-sea exchange of CO 2 . However, we have a limited understanding of the three-dimensional circulation that brings deep carbon-rich waters to the surface. Here, we introduce and analyze a novel high-resolution ocean model simulation with active biogeochemistry and online Lagrangian particle tracking. We focus our attention on a subset of particles with high dissolved inorganic carbon (DIC) that originate below 1,000 m and eventually upwell into the near-surface layer (upper 200 m). We find that 71% of the DIC-enriched water upwelling across 1,000 m is concentrated near topographic features, which occupy just 33% of the Antarctic Circumpolar Current. Once particles upwell to the near-surface layer, they exhibit relatively uniform pCO 2 levels and DIC decorrelation timescales, regardless of their origin. Furthermore, our results show that Southern Ocean bathymetry plays a key role in delivering carbon-rich waters to the surface.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Variations in boundary layer stability across Antarctica: a comparison between coastal and interior sites

The range of boundary layer stability profiles, from the surface to 500 m a.g.l. (above ground level), present in radiosonde observations from two continental-interior (South Pole Station and Dome Concordia Station) and three coastal (McMurdo Station, Georg von Neumayer Station III, and Syowa Station) Antarctic sites, is examined using the self-organizing maps (SOMs) neural network algorithm. A wide range of potential temperature profiles is revealed, from shallow boundary layers with strong near-surface stability to deeper boundary layers with weaker or near-neutral stability, as well as profiles with weaker near-surface stability and enhanced stability aloft, above the boundary layer. Boundary layer regimes were defined based on the range of profiles revealed by the SOM analysis; 20 boundary layer regimes were identified to account for differences in stability near the surface as well as above the boundary layer. Strong, very strong, or extremely strong stability, with vertical potential temperature gradients of 5 to in excess of 30 K per 100 m, occurred more than 80 % of the time at South Pole and Dome Concordia in the winter. Weaker stability was found in the winter at the coastal sites, with moderate and strong stability (vertical potential temperature gradients of 1.75 to 15 K per 100 m) occurring 70 % to 85 % of the time. Even in the summer, moderate and strong stability is found across all five sites, either immediately near the surface or aloft, just above the boundary layer. While the mean boundary layer height at the continental-interior sites was found to be approximately 50 m, the mean boundary layer height at the coastal sites was deeper, around 110 m. Further, a commonly described two-stability-regime system in the Arctic associated with clear or cloudy conditions was applied to the 20 boundary layer regimes identified in this study to understand if the two-regime behavior is also observed in the Antarctic. It was found that moderate and strong stability occur more often with clear- than cloudy-sky conditions, but weaker stability regimes occur almost equally for clear and cloudy conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Deep level defects and cation sublattice disorder in ZnGeN 2

The III-nitrides have revolutionized lighting technology and power electronics. Expanding the nitride semiconductor family to include the heterovalent ternary nitrides opens up new and exciting opportunities for device design that may help overcome some of the limitations of the binary nitrides as well as new applications. However, the more complex cation sublattice also gives rise to new complexities with both native and extrinsic defects, such as disorder on the cation sublattice. Here, depth-resolved cathodoluminescence spectroscopy and surface photovoltage spectroscopy measurements of defect energy levels in ZnGeN 2 combined with transmission microscopy and crystal lattice diffraction reveal optical signatures of mid-gap states that can be associated with cation sublattice disorder. The energies of these characteristic optical signatures in ZnGeN 2 thin films grown by metal-organic chemical vapor deposition are in good agreement with energy levels predicted by density functional theory for specific thermodynamically stable native point defects. To test the dependence of these defect signatures on crystal atomic composition, we correlated spatially-resolved optical and atomic composition measurements using spatially resolved x-ray photoelectron spectroscopy with systematically varied growth conditions on the same ZnGeN 2 films. The resultant elemental maps versus defect spectral energies and intensities suggest that cation antisite complexes (Zn Ge -Ge Zn ) versus isolated native point defects form preferentially and introduce mid-gap states that dominate electron-hole pair recombination. Complexing of Zn Ge and Ge Zn antisites manifests as disorder in the cation sub-lattice and leads to the formation of wurtzitic ZnGeN 2 as indicated by transmission electron microscopy diffraction patterns and x-ray diffraction (XRD) reciprocal space maps. Furthermore, these findings emphasize the importance of growth conditions to control cation place exchange.

42 ENGINEERING↗
DOE OSTI2025

Dual roles of stepped surfaces: Catalytic initiators and stabilizers of oxygen-induced reconstructions

Understanding surface restructuring under reactive conditions is crucial for designing next-generation catalysts with enhanced activity and selectivity. Here, we employ in situ transmission electron microscopy to directly observe the dynamic behavior of Cu(100) and Cu(410) surfaces under both oxidizing and vacuum annealing conditions, revealing a complex interplay among surface crystallography, local oxygen coverage, Cu atom mobility, and step-edge reactivity. The stepped Cu(410) surface acts as an active site for O 2 dissociation, triggering the oscillatory transformation of the c(2 × 2)–O phase into the more stable (2$\sqrt2$ ×$\sqrt2$)R45°–O missing-row (MR) structure on the adjacent flat Cu(100) terrace. Under subsequent vacuum annealing, this same Cu(410) facet exhibits remarkable structural resilience, preserving the MR reconstruction and chemisorbed oxygen. In contrast, the Cu(100) surface undergoes reversible transitions from the MR structure back to the c(2 × 2)–O phase. These results highlight the critical role of surface morphology in directing both the formation and stability of oxygen-induced reconstructions, demonstrating that stepped surfaces serve dual roles as both catalytic initiators and structural stabilizers. Furthermore, this work offers atomic-level insights into the environment-responsive behavior of copper surfaces, establishing a mechanistic basis for designing Cu-based catalysts through facet-specific control of surface reactivity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Copper migration and surface oxidation of CuxBi2Se3 in ambient pressure environments

Abstract Chemical modifications such as intercalation can be used to modify surface properties or to further functionalize the surface states of topological insulators (TIs). Using ambient pressure x-ray photoelectron spectroscopy, we report copper migration in C u x B i 2 S e 3 , which occurs on a timescale of hours to days after initial surface cleaving. The increase in near-surface copper proceeds along with the oxidation of the sample surface and large changes in the selenium content. These complex changes are further modeled with core-level spectroscopy simulations, which suggest a composition gradient near the surface which develops with oxygen exposure. Our results shed light on a new phenomenon that must be considered for intercalated TIs—and intercalated materials in general—that surface chemical composition can change when specimens are exposed to ambient conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Surfactants Used in Colloidal Synthesis Modulate Ni Nanoparticle Surface Evolution for Selective CO 2 Hydrogenation

Colloidal chemistry holds promise to prepare uniform and size-controllable pre-catalysts; however, it remains a challenge to unveil the atomic-level transition from pre-catalysts to active catalytic surfaces under the reaction conditions to enable the mechanistic design of catalysts. In this work, we report an ambient-pressure X-ray photoelectron spectroscopy study, coupled with in situ environmental transmission electron microscopy, infrared spectroscopy, and theoretical calculations, to elucidate the surface catalytic sites of colloidal Ni nanoparticles for CO 2 hydrogenation. We show that Ni nanoparticles with phosphine ligands exhibit a distinct surface evolution compared with amine-capped ones, owing to the diffusion of P under oxidative (air) or reductive (CO 2 + H 2 ) gaseous environments at elevated temperatures. The resulting NiP x surface leads to a substantially improved selectivity for CO production, in contrast to the metallic Ni, which favors CH 4 . The further elimination of surface metallic Ni sites by designing multi-step P incorporation achieves unit selectivity of CO in high-rate CO 2 hydrogenation.

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

Antiferromagnetic topological insulating state in $\mathrm{Tb_{0.02}Bi_{1.08}Sb_{0.9}Te_2S}$ single crystals

Topological insulators are emerging materials with insulating bulk and symmetry protected nontrivial surface states. One of the most fascinating transport behaviors in a topological insulator is the quantized anomalous Hall insulator, which has been observed in magnetic-topological-insulator-based devices. Here in this work, we report a successful doping of rare earth element Tb into Bi 1.08 Sb 0.9 Te 2 S topological insulator single crystals, in which the Tb moments are antiferromagnetically ordered below ~10 K. Benefiting from the in-bulk-gap Fermi level, transport behavior dominant by the topological surface states is observed below ~150 K. At low temperatures, strong Shubnikov–de Haas oscillations are observed, which exhibit 2D-like behavior. The topological insulator with long range magnetic ordering in rare earth doped Bi 1.08 Sb 0.9 Te 2 S single crystal provides an ideal platform for quantum transport studies and potential applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Ab-initio simulation of spin-vibronic spectra of methoxy radical

Despite the fact that experimental and theoretical work on the spectrum of methoxy has stretched from the microwave to the ultraviolet and proceeded for nearly 50 years, parts of the spectrum have remained a challenge to simulate theoretically and make reliable line-by-line assignments. The spectral complexity arises because the radical has a non-zero electron spin and significant vibronic coupling between the two elec- tronic components of the ground state due to the presence of a conical intersection. This work describes a completely ab initio effort to understand and assign the spin- vibronic levels of the X 2E state from 0 to above 3000 cm−1, a region that includes the fundamental transitions of the C-H symmetric and asymmetric stretches that have not previously been identified uniquely. A potential energy surface for methoxy was calculated at the EOM-CCSDT/ANO1 level of theory. Subsequently this potential energy surface was fit to a quartic power series expansion of all nine vibrational nor- mal coordinates (as determined at the minimum of the conical intersection) by the use of a machine-learning-based algorithm. After the addition of spin-orbit coupling, the spin-vibronic problem was solved using both the Krylov-Schur and Lanczos algorithms with the SOCJT3 software to converge eigenvalues up to 3500 cm−1 and their eigen- vectors. The latter were used, in conjunction with the calculated dipole moment and its derivatives (calculated using finite differences at EOM-CCSDT/ANO1 level), to determine spectral intensities for the spin-vibronic spectra. The calculated transition frequencies and intensities were used to simulate and assign the observed transitions of the spin-vibronic spectra of the radical. The credibility of the assignments and their significance is discussed in detail.

Sharma, Ketan [University of Florida, Gainesville,↗
DOE OSTI2023

Deep and Shallow Gap States in Reduced and n-Type Doped m -ZrO 2

Monoclinic zirconium dioxide (m-ZrO 2 ) is a wide-band-gap functional oxide with many applications. Understanding the role of intrinsic defects and dopants is important for improving the properties of m-ZrO 2 relevant to its applications. In this work, we characterize the electronic states and energy levels of oxygen vacancies (VO’s) and n-type Nb, Ta, Rb, H, and F dopants in the bulk and at the majority (1̅11) surface of m-ZrO 2 , using accurate dielectric-dependent hybrid functional calculations. Our results show that VO’s generally behave as deep donors with the excess electrons localized at the vacant O sites in the form of F centers, in agreement with previous studies. In contrast, surface VO’s at two-fold coordinated oxygen sites show a rather shallow (1+/0) transition level. Shallow levels are also obtained for adsorbed hydrogen (H ads ) and fluorine substituting an oxygen atom (FO) either in the bulk or at the surface, whereas Nb and Ta form deep levels. Altogether, our results suggest that FO is the best candidate for realizing n-type conductivity in bulk m-ZrO 2 , while surface VO’s and Hads provide an efficient way to control the surface chemistry of m-ZrO 2 .

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

A novel hardmask-to-substrate pattern transfer method for creating 3D, multi-level, hierarchical, high aspect-ratio structures for applications in microfluidics and cooling technologies

Abstract This letter solves a major hurdle that mars photolithography-based fabrication of micro-mesoscale structures in silicon. Conventional photolithography is usually performed on smooth, flat wafer surfaces to lay a 2D design and subsequently etch it to create single-level features. It is, however, unable to process non-flat surfaces or already etched wafers and create more than one level in the structure. In this study, we have described a novel cleanroom-based process flow that allows for easy creation of such multi-level, hierarchical 3D structures in a substrate. This is achieved by introducing an ultra-thin sacrificial silicon dioxide hardmask layer on the substrate which is first 3D patterned via multiple rounds of lithography. This 3D pattern is then scaled vertically by a factor of 200–300 and transferred to the substrate underneath via a single shot deep etching step. The proposed method is also easily characterizable—using features of different topographies and dimensions, the etch rates and selectivities were quantified; this characterization information was later used while fabricating specific target structures. Furthermore, this study comprehensively compares the novel pattern transfer technique to already existing methods of creating multi-level structures, like grayscale lithography and chip stacking. The proposed process was found to be cheaper, faster, and easier to standardize compared to other methods—this made the overall process more reliable and repeatable. We hope it will encourage more research into hybrid structures that hold the key to dramatic performance improvements in several micro-mesoscale devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

The Role of Snowmelt and Subsurface Heterogeneity in Headwater Hydrology of a Mountainous Catchment in Colorado: A Model‐Data Integration Approach

Mountainous headwater streams are sustained by both snowmelt‐driven streamflow and groundwater discharge in the Upper Colorado River Basin. However, predicting headwater stream discharge magnitude and peak flow timing is challenging in mountainous terrains, where snowmelt rates vary with vegetation type and elevation, and heterogeneous subsurface physical properties influence groundwater storage and its release. We used a model‐data integration approach to investigate the roles of snowmelt and subsurface structure in stream discharge and groundwater level. We ran an ensemble of 100 integrated surface‐subsurface hydrologic models for a mountainous headwater catchment near Crested Butte, Colorado, USA. We also evaluated and calibrated these models against observed data sets, including snow depth measurements using distributed temperature probes, stream discharge, and groundwater levels. Calibration with multiple data sources using neural density estimators has further constrained uncertainty in subsurface properties and snowmelt rates. Results indicated that observed slower snowmelt rates in evergreen forests delayed the peak flow and baseflow onset. In upstream areas with lower subsurface permeability, water was stored within the subsurface but was not released as interflow or shallow groundwater flow, and thereby not contributing to downstream streamflow during recession limb periods. Double peaks in groundwater occurred in areas with spatial subsurface heterogeneity, in our case due to the contrast between granodiorite and Mancos shale. These process‐based insights into groundwater and snowmelt dynamics in mountainous headwaters will help improve predictions of headwater hydrology.

Wang, Lijing [University of Connecticut, Storrs, C↗
DOE OSTI2024

Prediction of O and OH Adsorption on Transition Metal Oxide Surfaces from Bulk Descriptors

In the search for stable and active catalysts, density functional theory and machine learning (ML) based models can accelerate the screening of materials. While stability is conveniently addressed on the bulk level of computation, the modelling of catalytic activity requires expensive surface simulations. Here, in this work, we develop models for the surface adsorption energy of O and OH intermediates across a consistent and extensive dataset of pure transition metal oxide surfaces. We show that adsorption energies across metal oxidation states of +2 to +6 are well captured from the metal-oxygen bond strength extracted from the bulk level calculation. Specifically, we calculate the integrated crystal orbital Hamiltonian population (ICOHP) of the metal-oxygen bond in the bulk oxide and employ a simple normalization scheme to obtain a strong correlation with adsorption energetics. By combining our ICOHP descriptor with non DFT features in a Gaussian Process regression (GPR) model, we achieve high model accuracy with mean absolute errors of 0.166 and 0.219 eV for OH and O adsorption, respectively. By targeting the O-OH adsorption energy difference with our GPR model, we predict the the oxygen evolution reaction (OER) activity from bulk descriptors only. Furthermore, we utilize the strong correlation between the COHP and metal oxygen bond lengths to rapidly predict adsorption energetics and catalytic activity from the optimized bulk geometry. Our approach can enable an efficient search for active catalysts by eliminating the need for surface calculations in the initial screening phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗