Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface interactions”

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

Synoptic and Mesoscale Modulation of Dynamic and Thermodynamic Impacts on Central Arctic Sea Ice During MOSAiC

The primary objectives of this proposed project are to determine 1) to what extent dynamic and thermodynamic air‐ice interactions in the Central Arctic are driven by atmospheric synoptic and mesoscale phenomena, and 2) the key processes involved in these interactions. These phenomena include large‐scale features such as fronts, high‐ and low‐pressure systems, low‐level jets, and cloud bands, but may also include smaller‐scale phenomena such as mesoscale instabilities (e.g., symmetric instability or slantwise convection), transient meso‐lows and gravity waves. Dynamic air‐ice interactions include direct atmospheric stress and atmospheric stress divergence/deformation; thermodynamic air‐ice interactions involve the net surface energy budget, especially in association with non‐stratocumulus clouds generally occurring with stronger dynamical forcing.

58 GEOSCIENCES↗
DOE OSTI2025

Operando synchrotron X-ray analysis of melt pool dynamics in an Al-Sn immiscible alloy

The melt flow in an Al-50vol.% Sn immiscible alloy, produced by single-track laser melting of Al and Sn elemental powders, was studied in real time. High-speed synchrotron X-ray imaging was used to track Al and Sn liquids’ movements, and also to examine elemental distributions in the laser tracks, complimented by electron microscopy after solidification. Key aspects, including melt pool geometry, keyhole instability, and flow dynamics (flow pattern and velocity), were examined using digital image analysis. Relatively deeper melt pools formed at 400 W and 300 mm/s exhibited greater stability, with smooth surfaces, consistent outward flow, and minor vortices near the keyhole. In contrast, shallower pools produced at higher scanning speeds (>500 mm/s) demonstrated greater instability with increased surface waviness, and stronger velocity fluctuations, leading to numerous micro-vortices and increased Al-Sn heterogeneity. Velocity scale estimations, supported by experimental observations, examined the roles of vapour pressure, Marangoni effect, buoyancy, inertial, and surface tension forces in the flow. The results revealed that vapour pressure and mechanical waves dominated at high scanning speeds (shallow pools), while Marangoni forces were equally significant in deep pools at lower speeds (300 mm/s). Buoyancy was found to have minimal impact in both cases. Furthermore, the interaction between inertial and surface tension forces played a critical role in determining the degree of waviness of the pools’ surfaces. These findings offer valuable insights into melt pool dynamics during laser processing of immiscible alloys and other metallic systems using elemental powders, and provide guidance for developing high-fidelity computational fluid dynamics models.

immiscible alloys↗
DOE OSTI2025

Structure and Dynamics of CO 2 at the Air–Water Interface from Classical and Neural Network Potentials

The accurate description of the structure and dynamics of CO 2 at the instantaneous air–water interface, along with the effects of surface fluctuations on the CO 2 -transport processes, is essential for the development of negative emission technologies aimed at minimizing climate change. In this study, we performed molecular dynamics simulations of CO 2 at the air–water interface using neural network potentials (NNPs) trained on ab initio data generated through density-functional-theory-based molecular dynamics simulations. We compared these results with classical force fields to assess their performance in modeling interfacial CO 2 behavior. Our findings revealed that the asymmetric interactions, coupled with thermal surface fluctuations at the air–water interface, significantly influence CO 2 transport into the aqueous phase. The simulations demonstrate that classical force fields underestimate both the free energy of CO 2 transport and the strength of its interactions at the interface compared with the neural network potentials. In conclusion, the free energy and the interfacial dynamics of CO 2 are primarily influenced by the distribution of water within the instantaneous interfacial water layer, responsible for creating an asymmetric intermolecular interaction environment within the interfacial region.

Ab initio molecular dynamics↗
DOE OSTI2025

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗
DOE OSTI2025

Relativistic approach to manipulating angular distribution of charged particles via kinetic equations

Deflection angles of charged particles interacting with materials play a critical role in various plasma applications. The development of a mathematically well-posed kinetic collision operator that accounts for deflection angles of strong Coulomb interactions remains a fundamental open problem. This paper presents a relativistic method for modifying the electromagnetic field in an anisotropic and adjustable manner to manipulate a system of charged particles, specifically by the transfer of angular momentum from a superluminal wave source to particles at specific times and locations. The method provides a mechanism to influence the scattering outcomes of strong interactions by manipulating the angular distribution of particles, and thus the deflection angles of their interactions with a solid surface, without requiring detailed knowledge of the kinetic collision operator. To this end, we demonstrate how a specific type of singularity, generated by Maxwell's equations for a superluminal wave source at the boundary of the plasma, can modify the electromagnetic field in a highly directional manner. The proposed method can lead to the development of novel approaches for controlling interactions of charged particles with a material in plasma systems. Published by the American Physical Society 2025

Moini, Nima (ORCID:0009000929568824)↗
DOE OSTI2024

CO–induced roughening of Cu(111): formation and detection of reactive nanoclusters on metal surfaces

The formation of nanoclusters on metal surfaces in the presence of reactive environments is a phenomenon with important implications for catalysis. These nanoclusters are composed of atoms ejected from undercoordinated sites such as step edges, and their presence alters the catalytic properties of solid materials. We perform density functional theory (DFT) and kinetic Monte Carlo (KMC) simulations to investigate the formation and reactivity of copper clusters on Cu(111). Our results indicate a considerably higher reactivity of small copper nanoclusters, with up to seven atoms in size on roughened copper surfaces than on pristine Cu(111) and Cu(211). Regarding the restructuring events that give rise to nanoclusters under CO atmospheres, we determine that the ejection of Cu atoms from step edges and their migration therefrom to adjacent Cu(111) terraces are, by and large, driven by CO coverage effects. By means of KMC simulations, which account for CO–CO lateral interactions and CO–induced surface restructuring, we show that temperature programmed desorption (TPD) holds promise for the detection of highly reactive nanoclusters. Furthermore, our approach showcases how surface restructuring and surface–adsorbate bond breaking can be combined when modeling surface reactions and contributes to the development of an advanced understanding of the nature of active site under reaction conditions.

catalyst dynamic restructuring↗
DOE OSTI2026

First-principles DFT modeling of nitrobenzene adsorption on the Ag(111) surface at varying monolayer coverages

Here, the adsorption behavior of nitrobenzene on the Ag(111) surface as a function of coverage is investigated using density functional theory. Adsorption energies and optimized geometries are analyzed together with isolated intermolecular interaction calculations and electronic structure analysis, including Bader charge partitioning and projected density of states, to disentangle the roles of adsorbate–surface bonding and through-space adsorbate–adsorbate interactions. At low and intermediate coverages (θ = 1/9 and θ = 2/9), bidentate adsorption at top sites is favored due to strong adsorbate–surface interactions, with additional stabilization at θ = 2/9 arising from favorable intermolecular separations. At higher coverage (θ = 1/3), bidentate adsorption is destabilized by strong intermolecular repulsion, and monodentate adsorption becomes energetically preferred as rotational freedom allows more favorable intermolecular spacing. Charge density differences, Bader charge, and density of states analyses show that charge transfer from the Ag surface is localized primarily on the nitro group and increases with coverage and adsorption denticity, although this increase does not directly correlate with adsorption strength at high coverage due to competing intermolecular interactions. These results demonstrate that surface coverage can induce a transition in preferred adsorption denticity driven by intermolecular interactions, highlighting the importance of adsorbate packing in organic molecule adsorption at metal interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Investigating Aerosol Hygroscopicity in the Subcloud Transition Zone and at the Surface in the Southern Great Plains

Aerosols beneath a cloud base, a subcloud transition zone (SCTZ), are key to understand both the aerosol-cloud interaction and aerosol-radiation interactions. Lidars have been the primary means of observing aerosols in the SCTZ by virtue of enhanced light scattered by aerosol particles. The enhanced light maybe caused by several factors: the aerosol swelling effect due to hygroscopicity under high relative humidity, cloud 3-dimensional (3D) effect, aerosol nucleation into cloud droplets, etc. While each factor and process has been known, their relative contributions are much poorly quantified. This study explores the hygroscopicity and optical properties of aerosols in the SCTZ and at ground level in the Southern Great Plains (SGP) region. Utilizing comprehensive observational data from the U.S. Department of Energy's Atmospheric Radiation Measurement at the Oklahoma SGP site, including ground-based aerosol measurements and Raman lidar profiles from April 2021 to April 2022, this study extensively analyzes the influence of aerosol hygroscopic growth and cloud fragments on aerosol optical properties. Distinct seasonal variations in aerosol hygroscopic characteristics are revealed. At the ground level, aerosols in autumn and winter exhibit stronger hygroscopicity due to a higher proportion of inorganic content than summer. In the SCTZ, aerosols during summer show enhanced backscatter due to strong cloud fragmentation effects, with numerous cloud fragments elevating hygroscopicity beyond that observed in autumn and winter. These insights are crucial for understanding the interactions between aerosols at the surface and cloud layers, evaluating cloud condensation nuclei beneath clouds, and their implications for atmospheric radiation and climate modeling.

Hu, Rong [Beijing Normal University, Beijing (Chin↗
DOE OSTI2024

CO 2 Uptake and Stability Enhancement in Vinyltrimethoxysilane‐Treated SBA‐15 Solid Amine‐Based Sorbents

Abstract Silica‐supported amine absorbents, including materials produced by tethering aminosilanes or infusion of poly(ethyleneimine), represent a promising class of materials for CO 2 capture applications, including direct air and point source capture. Various silica surface treatments and functionalization strategies are explored to enhance stability and CO 2 uptake in amine‐based solid sorbent systems. Here, the synthesis and characterization of novel vinyltrimethoxysilane‐treated Santa Barbara Amorphous‐15 (SBA‐15) supports and the corresponding enhancement in CO 2 uptake compared to various SBA‐15‐based control supports are presented. The relationship between CO 2 diffusion and amine efficiency in these systems is explored using a previously reported kinetic model. The synthesized materials are characterized with CO 2 and H 2 O isotherms, diffuse reflectance infrared Fourier transform spectroscopy, 1 H T 1 – T 2 relaxation correlation NMR, and rapid thermal cycling experiments. The novel support materials are shown to enable high amine efficiencies, approaching a fourfold improvement over standard SBA‐15‐supported amines, while simultaneously exhibiting excellent stability when cycled rapidly under humid conditions. As the poly(ethyleneimine) loadings are held constant across the various samples, enhancements in CO 2 uptake are attributed to differences in the way the poly(ethyleneimine) interacts with the support surface.

Vallace, Anthony↗
DOE OSTI2025

Ultrafast Dynamics of Plasmon-Coupled Excitons in Semiconducting Nanoplatelets

Exciton-plasmon coupling in nanomaterials produces many relevant phenomena for photonics applications including increased light-matter interactions, enhanced radiative rates of quantum emitters, and coherent energy exchange. In the case of exciton coupling to surface plasmon polaritons (SPPs), dispersive interactions controlled by the wavevector of optical excitation create the opportunity for tunable optical emission. Strong temporal impacts on exciton lifetimes can also occur in coupled systems, creating the opportunity for ultrafast control of exciton lifetime via changes in electronic coupling magnitude to a dispersive SPP. The coupling strength can be impacted by the morphology of the nanomaterials. Here, in this study, we utilize colloidal semiconductor nanoplatelets deposited onto thin silver plasmonic films, and compare the results to semiconductor quantum dots deposited on the silver films. We map the dispersion of the coupled systems and measure the ultrafast transient absorption response of the coupled systems. Due to the larger interaction areas of the nanoplatelets that lie flat on the silver films, a greater degree of coupling is found for the nanoplatelets, and much faster temporal responses are found as compared to quantum dots. Fresnel theory calculations that incorporate heavy and light hole features can reproduce the dispersion of the nanoplatelet-silver film, and a simple three-state model is developed to provide insights into the nature of the coupling at different photon energies along the dispersion curve.

colloidal semiconductor nanoplatelets↗
DOE OSTI2025

Designing protein–material interfaces

This article addresses recent advances in using de novo protein design to create coherent interfaces between proteins and inorganic materials, either through protein self-assembly on crystal lattices or through directed nucleation and growth of crystals by protein scaffolds. Inspired by natural protein-crystal interfaces, we focus on a class of designed helical repeat proteins that present a repeating pattern of charged amino acid residues. We describe the use of in situ imaging and spectroscopic methods to investigate both the assembly of these proteins and their ability to direct crystal nucleation and growth. Furthermore, the findings reveal the importance of surface charge, facet-specific binding, solvent organization, and, more generally, the balance of protein-substrate-solvent interactions in determining how organized protein-materials interfaces emerge. Moreover, the results demonstrate the vast potential of protein design in materials science and elucidate the mechanisms by which interactions between biomolecules and inorganic surfaces lead to unique materials and morphologies.

Biomaterials-Proteins↗
DOE OSTI2026

Catalytic Ammonia Synthesis over Pure, Defective, and Metal-Doped Rutile TiO 2 : A Periodic DFT Study

In this work, we aim to describe the energetics associated with the formation of ammonia from N 2 interacting with doped hydroxylated rutile TiO 2 (110) surfaces with the vacant O 2c site, following the reaction N 2 + 3H 2 O → 2NH 3 + 3/2O 2 . The water molecules interact with the surface, creating exposed Ti−OH groups that can transfer hydrogen to the adsorbed N 2 molecule. Two metal dopants are evaluated: Mo and Ta. For both metals, calculations show a dramatic decrease in the energy of most intermediates during the entire mechanism, leading to more favorable reaction mechanisms. Nonetheless, it is worth noting that when the Ti 6c site of the vacant site is doped with either Mo or Ta, there is a stronger effect on the energetics than doping on the exposed Ti 5c sites. The effect of increasing the concentration of metal dopants on the vacant site was also investigated. In this case, calculations indicate that a higher percentage of the dopant on the surface results in a more substantial decrease in the energy of most intermediates, suggesting that increasing the dopant content could be beneficial for the catalytic process.

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

Niobium hydride formation in superconducting qubit thin films

The formation of nonsuperconducting hydrides in 160–170-nm-thick films of niobium is examined. We identify six elastically distinct orientation relationships ɛ−Nb 4 ⁢H 3 takes within the matrix, solid-solution 𝛼−NbH 𝑥 phase. We employ a phase field model to assess the impact of elastic energy induced by the strain of phase transformation on the morphology and transformation dynamics of ɛ precipitates within a thin film. We consider the dimensions of the thin film, crystallographic growth direction, and diffusion rates to predict the timescale of hydride evolution. Leveraging the finite element method, we predict two-dimensional and three-dimensional equilibrium shapes of ɛ-hydrides within a bulk sample and in a thin film that has a traction-free surface. Our results suggest that niobium hydrides migrate to the free surface of the film. Precipitates which reach the free surface coarsen, while precipitates within the film dissolve. Precipitates in both two dimensions and three dimensions experience a repulsive interaction force at the free surface, that is attractive in the bulk, shown in experiment and theory of previous studies.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Drought-induced changes in groundwater-surface water exchange at Lake Mead area

This study focuses on the Lake Mead region in the southwestern United States, a key water reservoir serving over 25 million people and agricultural lands across several states. The area has experienced recurring anthropogenic droughts since the early 2000s. We investigate the hydrological response of the coupled surface-groundwater system to the 2020–2022 drought, one of the most severe on record. To this end, we use Sentinel-1 Interferometric Synthetic Aperture Radar (InSAR) data to quantify vertical land motion caused by the elastic response of the crust to water-mass loss in the lake vicinity. Next, we apply an inverse elastic load modeling framework to quantify water loss. We further assess possible hydraulic connectivity between Lake Mead and adjacent groundwater reservoirs. We detect ground uplift of up to 8 mm/yr near the lake center, likely due to crustal rebound from reduced water-mass loading. We estimated the total water storage loss at 3.03 ± 0.25 km 3 /yr across a 3150 km 2 area surrounding Lake Mead. Groundwater accounts for approximately a third of that, being 0.94 ± 0.32 km 3 /yr. In addition, observed time lags of 6–98 days between lake and groundwater level responses, corresponding to a lateral diffusivity of 3.2–86 m 2 /s, suggest spatially variable connectivity between the lake and aquifers. These findings highlight that drought impacts propagate through the subsurface within interconnected systems, resulting in reduced buffering capacity of groundwater resources following droughts, and emphasizing the need for more integrated surface and groundwater management strategies to enhance resilience under climate and anthropogenic stressors.

58 GEOSCIENCES↗
DOE OSTI2024

Identifying recharge sources and their impacts on a North Central New Mexico shallow aquifer using unsupervised machine learning

In this article, shallow aquifers are important but highly variable resources in arid to semi-arid regions. Limited shallow aquifer volume results in high sensitivity to recharge fluctuations, which can impact the local fauna and flora, and transport of contaminants in the aquifer or vadose zone. Aquifer response to external forcing (e.g., precipitation) is usually solved by estimating aquifer parameters and running physics-based models to match known fluctuations of hydraulic head. However, this technique is time and computationally expensive. Furthermore, high aquifer complexity decreases precision in physics-based models. Alternatively supervised machine learning is used to predict aquifer dynamics. However, these techniques rely on input data and struggle to interpret aquifer response for missing sources (i.e., snowpack data). To counter these problems, we propose an unsupervised machine learning technique (NMFk) to estimate the impact of different sources on aquifer recharge. NMFk is used to understand the influence of external forcing on shallow aquifer recharge in the Pajarito Plateau (Los Alamos, NM, USA). The results show how NMFk can be used to reduce the data dimension in a complex field dataset to three recharge signals that cause fluctuations within the field data. Here, the source signals are interpreted as rainfall, snowmelt, and a delayed aquifer response to the previous two signals. These results evidence how heterogeneous aquifers delimited by canyons incised into the Pajarito Plateau respond in similar ways to the source signals identified by NMFk. Furthermore, results show the importance of the local geology where faults act as sinks, and anthropogenic disturbances can facilitate infiltration amplifying the interpreted signal.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

APOA2 increases cholesterol efflux capacity to plasma HDL by displacing the C-terminus of resident APOA1

The ability of high-density lipoprotein (HDL) to promote cellular cholesterol efflux is a more robust predictor of cardiovascular disease protection than HDL-cholesterol levels in plasma. Previously, we found that lipidated HDL containing both apolipoprotein A-I (APOA1) and A-II (APOA2) promotes cholesterol efflux via the ATP-binding cassette transporter (ABCA1). In the current study, we directly added purified, lipid-free APOA2 to human plasma and found a dose-dependent increase in whole plasma cholesterol efflux capacity. APOA2 likewise increased the cholesterol efflux capacity of isolated HDL with the maximum effect occurring when equal masses of APOA1 and APOA2 coexisted on the particles. Follow-up experiments with reconstituted HDL corroborated that the presence of both APOA1 and APOA2 were necessary for the increased efflux. Using limited proteolysis and chemical cross-linking mass spectrometry, we found that APOA2 induced a conformational change in the N- and C-terminal helices of APOA1. Using reconstituted HDL with APOA1 deletion mutants, we further showed that APOA2 lost its ability to stimulate ABCA1 efflux to HDL if the C-terminal domain of APOA1 was absent, but retained this ability when the N-terminal domain was absent. Based on these findings, we propose a model in which APOA2 displaces the C-terminal helix of APOA1 from the HDL surface which can then interact with ABCA1—much like it does in lipid-poor APOA1. These findings suggest APOA2 may be a novel therapeutic target given this ability to open a large, high-capacity pool of HDL particles to enhance ABCA1-mediated cholesterol efflux.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Unveiling the porosity effect of superbase ionic liquid-modified carbon sorbents in CO 2 capture from air

Direct air capture (DAC) of CO 2 represents one of the most promising technologies to achieve negative carbon emissions. In this work, the superbase ionic liquids (ILs)-modified carbon substrates were developed for DAC of CO 2 by harnessing the strong CO 2 binding capability of IL and the ordered porous channels of the carbon supports. Detailed porosity analysis revealed that the IL with an aromatic cation and an oxygenate anion preferred to fill the micropores, and a thin layer was created on the surface of the mesopores. Strong π-π interaction between the IL layer and the carbon surface was disclosed by wide-angle X-ray scattering (WAXS) analysis, leading to enhanced thermal stability of the IL phase. For the same lL coating amount, the DAC of CO 2 evaluation revealed that a larger mesopore size and pore volume in the carbon/IL composite materials led to higher CO 2 uptake capacity by exposing more active sites to integrate CO 2 from diluted sources. Further, the thermodynamic analysis confirmed the critical role of IL coating in providing strong chemisorption sites and significantly improved selectivity to enrich the diluted CO 2 from the air atmosphere. This work provides guidance on leveraging the scaffolds' surface properties and porosities of the scaffolds to optimize DAC of CO 2 behavior.

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

Morphology Effects on Free Energies of Proton-Coupled Electron Transfer in Polyoxotungstates

Polyoxotungstates have previously been established to facilitate the hydrogenation of small molecule substrates via hydrogen atom transfer from reactive hydroxyl groups formed at the assembly surface. Understanding structure−function relationships that dictate the thermochemistry and kinetics of protoncoupled electron transfer is key to controlling this chemistry. In this work, we combine comprehensive electrochemical experiments and density functional theory calculations to address how different polyoxotungstate morphologies, specifically W 6 O 19 −2 , W 10 O 32 −4 , SiW 12 O 40 −4 , and P 2 W 18 O 62 −6 , affect the bond dissociation free energies of surface hydroxides (BDFE(O−H)) formed upon reduction of the assembly in acidic media. Our results reveal increasing hydroxide bond strengths with increasing cluster size, and that anisotropic cluster geometries result in substantial thermodynamic differentiation of H-binding sites. We demonstrate an excellent agreement between theory and experiments on the reported BDFE(O−H) values and, importantly, we elucidate how cluster size and shape affect electronic properties (local charges and frontier molecular orbitals), giving rise to sites with increased preference for hydrogen binding, demonstrated in higher BDFE(O−H). Overall, this work aids the understanding and design of polyoxometalates exhibiting surface sites with tailored interaction strengths.

anions↗