Computation of Overhauser dynamic nuclear polarization processes reveals fundamental correlation between water dynamics, structure, and solvent restructuring entropy.

Hydration water dynamics, structure, and thermodynamics are crucially important to understand and predict water-mediated properties at molecular interfaces. Yet experimentally and directly quantifying water behavior locally near interfaces at the sub-nanometer scale is challenging, especially at interfaces submerged in biological solutions. Overhauser dynamic nuclear polarization (ODNP) experiments measure equilibrium hydration water dynamics within 8-15 angstroms of a nitroxide spin probe on instantaneous timescales (10 picoseconds to nanoseconds), making ODNP a powerful tool for probing local water dynamics in the vicinity of the spin probe. As with other spectroscopic techniques, concurrent computational analysis is necessary to gain access to detailed molecular level information about the dynamic, structural, and thermodynamic properties of water from experimental ODNP data. We chose a model system that can systematically tune the dynamics of water, a water-glycerol mixture with compositions ranging from 0 to 0.3 mole fraction glycerol. We demonstrate the ability of molecular dynamics (MD) simulations to compute ODNP spectroscopic quantities, and show that translational, rotational, and hydrogen bonding dynamics of hydration water align strongly with spectroscopic ODNP parameters. Moreover, MD simulations show tight correlations between the dynamic properties of water that ODNP captures and the structural and thermodynamic behavior of water. Hence, experimental ODNP readouts of varying water dynamics suggest changes in local structural and thermodynamic hydration water properties.

Bridging the Gap in Cryopreservation Mechanism: Unraveling the Interplay between Structure, Dynamics, and Thermodynamics in Cryoprotectant Aqueous Solutions.

Cryoprotectants play a crucial role in preserving biological material, ensuring their viability during storage and facilitating crucial applications such as the conservation of medical compounds, tissues, and organs for transplantation. However, the precise mechanism by which cryoprotectants modulate the thermodynamic properties of water to impede the formation and growth of ice crystals, thus preventing long-term damage, remains elusive. This is evident in the use of empirically optimized recipes for mixtures that typically contain DMSO, glycerol, and various sugar constituents. Here, we use terahertz calorimetry, Overhauser nuclear polarization, and molecular dynamics simulations to show that DMSO exhibits a robust structuring effect on water around its methyl groups, reaching a maximum at a DMSO mole fraction of XDMSO = 0.33. In contrast, glycerol exerts a smaller water-structuring effect, even at higher concentrations (Scheme 1). These results potentially suggest that the wrapped water around DMSO's methyl group, which can be evicted upon ligand binding, may render DMSO a more surface-active cryoprotectant than glycerol, while glycerol may participate more as a viscogen that acts on the entire sample. These findings shed light on the molecular intricacies of cryoprotectant solvation behavior and have potentially significant implications for optimizing cryopreservation protocols.

Nanoscale water-polymer interactions tune macroscopic diffusivity of water in aqueous poly(ethylene oxide) solutions.

The separation and anti-fouling performance of water purification membranes is governed by both macroscopic and molecular-scale water properties near polymer surfaces. However, even for poly(ethylene oxide) (PEO) - ubiquitously used in membrane materials - there is little understanding of whether or how the molecular structure of water near PEO surfaces affects macroscopic water diffusion. Here, we probe both time-averaged bulk and local water dynamics in dilute and concentrated PEO solutions using a unique combination of experimental and simulation tools. Pulsed-Field Gradient NMR and Overhauser Dynamic Nuclear Polarization (ODNP) capture water dynamics across micrometer length scales in sub-seconds to sub-nanometers in tens of picoseconds, respectively. We find that classical models, such as the Stokes-Einstein and Mackie-Meares relations, cannot capture water diffusion across a wide range of PEO concentrations, but that free volume theory can. Our study shows that PEO concentration affects macroscopic water diffusion by enhancing the water structure and altering free volume. ODNP experiments reveal that water diffusivity near PEO is slower than in the bulk in dilute solutions, previously not recognized by macroscopic transport measurements, but the two populations converge above the polymer overlap concentration. Molecular dynamics simulations reveal that the reduction in water diffusivity occurs with enhanced tetrahedral structuring near PEO. Broadly, we find that PEO does not simply behave like a physical obstruction but directly modifies water's structural and dynamic properties. Thus, even in simple PEO solutions, molecular scale structuring and the impact of polymer interfaces is essential to capturing water diffusion, an observation with important implications for water transport through structurally complex membrane materials.

Relationships between Molecular Structural Order Parameters and Equilibrium Water Dynamics in Aqueous Mixtures.

Water's unique thermophysical properties and how it mediates aqueous interactions between solutes have long been interpreted in terms of its collective molecular structure. The seminal work of Errington and Debenedetti [Nature 2001, 409, 318-321] revealed a striking hierarchy of relationships among the thermodynamic, dynamic, and structural properties of water, motivating many efforts to understand (1) what measures of water structure are connected to different experimentally accessible macroscopic responses and (2) how many such structural metrics are adequate to describe the collective structural behavior of water. Diffusivity constitutes a particularly interesting experimentally accessible equilibrium property to investigate such relationships because advanced NMR techniques allow the measurement of bulk and local water dynamics in nanometer proximity to molecules and interfaces, suggesting the enticing possibility of measuring local diffusivities that report on water structure. Here, we apply statistical learning methods to discover persistent structure-dynamic correlations across a variety of simulated aqueous mixtures, from alcohol-water to polypeptoid-water systems. We investigate a variety of molecular water structure metrics and find that an unsupervised statistical learning algorithm (namely, sequential feature selection) identifies only two or three independent structural metrics that are sufficient to predict water self-diffusivity accurately. Surprisingly, the translational diffusivity of water across all mixed systems studied here is strongly correlated with a measure of tetrahedral order given by water's triplet angle distribution. We also identify a separate small number of structural metrics that well predict an important thermodynamic property, the excess chemical potential of an idealized methane-sized hydrophobe in water. Ultimately, we offer a Bayesian method of inferring water structure by using only structure-dynamics linear regression models with experimental Overhauser dynamic nuclear polarization (ODNP) measurements of water self-diffusivity. This study thus quantifies the relationships among several distinct structural order parameters in water and, through statistical learning, reveals the potential to leverage molecular structure to predict fundamental thermophysical properties. In turn, these findings suggest a framework for solving the inverse problem of inferring water's molecular structure using experimental measurements such as ODNP studies that probe local water properties.

End-to-End Distance Probability Distributions of Dilute Poly(ethylene oxide) in Aqueous Solution.

We introduce a powerful, widely applicable approach to characterizing polymer conformational distributions, specifically the end-to-end distance distributions, P(Ree), accessed through double electron-electron resonance (DEER) spectroscopy in conjunction with molecular dynamics (MD) simulations. The technique is demonstrated on one of the most widely used synthetic, disordered, water-soluble polymers: poly(ethylene oxide) (PEO). Despite its widespread importance, no systematic experimental characterization of PEO's Ree conformational landscape exists. The evaluation of P(Ree) is particularly important for short polymers or (bio)polymers with sequence complexities that deviate from simple polymer physics scaling laws valid for long chains. In this study, we characterize the Ree landscape by measuring P(Ree) for low molecular weight (MW: 0.22-2.6 kDa) dilute PEO chains. We use DEER with end-conjugated spin probes to resolve Ree populations from ∼2-9 nm and compare them with full distributions from MD. The P( Ree)'s from DEER and MD show remarkably good agreement, particularly at longer chain lengths where populations in the DEER-unresolvable range (<1.5 nm) are low. Both the P(Ree) and the root-mean-square R̃ee indicate that aqueous PEO is a semiflexible polymer in a good solvent, with the latter scaling linearly with molecular weight up to its persistence length (lp ∼ 0.48 nm), and rapidly transitioning to excluded volume scaling above lp. The R̃ee scaling is quantitatively consistent with that from experimental scattering data on high MW (>10 kDa) PEO and the P(Ree)'s crossover to the theoretical distribution for an excluded volume chain.

Water Structure and Properties at Hydrophilic and Hydrophobic Surfaces.

The properties of water on both molecular and macroscopic surfaces critically influence a wide range of physical behaviors, with applications spanning from membrane science to catalysis to protein engineering. Yet, our current understanding of water interfacing molecular and material surfaces is incomplete, in part because measurement of water structure and molecular-scale properties challenges even the most advanced experimental characterization techniques and computational approaches. This review highlights progress in the ongoing development of tools working to answer fundamental questions on the principles that govern the interactions between water and surfaces. One outstanding and critical question is what universal molecular signatures capture the hydrophobicity of different surfaces in an operationally meaningful way, since traditional macroscopic hydrophobicity measures like contact angles fail to capture even basic properties of molecular or extended surfaces with any heterogeneity at the nanometer length scale. Resolving this grand challenge will require close interactions between state-of-the-art experiments, simulations, and theory, spanning research groups and using agreed-upon model systems, to synthesize an integrated knowledge of solvation water structure, dynamics, and thermodynamics.

The Synthesis Science of Targeted Vapor-Phase Metal-Organic Framework Postmodification.

The postmodification of metal organic frameworks (MOFs) affords exceedingly high surface area materials with precisely installed chemical features, which provide new opportunities for detailed structure-function correlation in the field of catalysis. Here, we significantly expand upon the number of vapor-phase postmodification processes reported to date through screening a library of atomic layer deposition (ALD) precursors, which span metals across the periodic table and which include ligands from four distinct precursor classes. With a large library of precursors and synthesis conditions, we discern trends in the compatibility of precursor classes for well-behaved ALD in MOFs (AIM) and identify challenges and solutions to more precise postsynthetic modification.

Selective Methane Oxidation to Methanol on Cu-Oxo Dimers Stabilized by Zirconia Nodes of an NU-1000 Metal-Organic Framework.

Mononuclear and dinuclear copper species were synthesized at the nodes of an NU-1000 metal-organic framework (MOF) via cation exchange and subsequent oxidation at 200 °C in oxygen. Copper-exchanged MOFs are active for selectively converting methane to methanol at 150-200 °C. At 150 °C and 1 bar methane, approximately a third of the copper centers are involved in converting methane to methanol. Methanol productivity increased by 3-4-fold and selectivity increased from 70% to 90% by increasing the methane pressure from 1 to 40 bar. Density functional theory showed that reaction pathways on various copper sites are able to convert methane to methanol, the copper oxyl sites with much lower free energies of activation. Combining studies of the stoichiometric activity with characterization by in situ X-ray absorption spectroscopy and density functional theory, we conclude that dehydrated dinuclear copper oxyl sites formed after activation at 200 °C are responsible for the activity.

Quantum confinement in few layer SnS nanosheets.

Orthorhombic tin monosulfide (SnS) consists of layers of covalently bound Sn and S atoms held together by weak van der Waals forces and is a stable two-dimensional material with potentially useful properties in emerging applications such as valleytronics. Large-scale sustainable synthesis of few-layer (e.g., 1-10 layers) SnS is a challenge, which also slows progress in understanding their properties as a function of number of layers. Herein we describe solvothermal synthesis of SnS in water or ethylene glycol. The latter yields a flower-like morphology where the petals are SnS nanoplates and sonication and separation of these flowers via differential centrifugation yields 1-10 layer SnS nanoplates. The direct optical absorption edges of these SnS nanoplates blue-shift due to quantum confinement from 1.33 to 1.88 eV as the thickness (number of layers) is decreased from ∼5 nm (10 layers) to ∼2 nm (4 layers).

Well-Defined Rhodium-Gallium Catalytic Sites in a Metal-Organic Framework: Promoter-Controlled Selectivity in Alkyne Semihydrogenation to E-Alkenes.

Promoters are ubiquitous in industrial heterogeneous catalysts. The wider roles of promoters in accelerating catalysis and/or controlling selectivity are, however, not well understood. A model system has been developed where a heterobimetallic active site comprising an active metal (Rh) and a promoter ion (Ga) is preassembled and delivered onto a metal-organic framework (MOF) support, NU-1000. The Rh-Ga sites in NU-1000 selectively catalyze the hydrogenation of acyclic alkynes to E-alkenes. The overall stereoselectivity is complementary to the well-known Lindlar's catalyst, which generates Z-alkenes. The role of the Ga in promoting this unusual selectivity is evidenced by the lack of semihydrogenation selectivity when Ga is absent and only Rh is present in the active site.

Application and Limitations of Nanocasting in Metal-Organic Frameworks.

Nanocasting can be a useful strategy to transfer the catalytic metal clusters in metal-organic frameworks (MOFs) to an all-inorganic support such as silica. The incorporation of silica in the MOF pores as a secondary support has the potential to extend the application of the highly tunable metal-based active sites in MOFs to high temperature catalysis. Here, we demonstrate the applicability of the nanocasting method to a range of MOFs that incorporate catalytically attractive hexazirconium, hexacerium, or pentanickel oxide-based clusters (UiO-66, (Ce)UiO-66, (Ce)UiO-67, (Ce)MOF-808, DUT-9, and In- and Ni-postmetalated NU-1000). We describe, in tutorial form, the challenges associated with nanocasting of MOFs that are related to their small pore size and to considerations of chemical and mechanical stability, and we provide approaches to overcome some of these challenges. Some of these nanocast materials feature the site-isolated clusters in a porous, thermally stable silica matrix, suitable for catalysis at high temperatures; in others, structural rearrangement of clusters or partial cluster aggregation occurs, but extensive aggregation can be mitigated by the silica skeleton introduced during nanocasting.

Sinter-Resistant Platinum Catalyst Supported by Metal-Organic Framework.

Single atoms and few-atom clusters of platinum are uniformly installed on the zirconia nodes of a metal-organic framework (MOF) NU-1000 via targeted vapor-phase synthesis. The catalytic Pt clusters, site-isolated by organic linkers, are shown to exhibit high catalytic activity for ethylene hydrogenation while exhibiting resistance to sintering up to 200 °C. In situ IR spectroscopy reveals the presence of both single atoms and few-atom clusters that depend upon synthesis conditions. Operando X-ray absorption spectroscopy and X-ray pair distribution analyses reveal unique changes in chemical bonding environment and cluster size stability while on stream. Density functional theory calculations elucidate a favorable reaction pathway for ethylene hydrogenation with the novel catalyst. These results provide evidence that atomic layer deposition (ALD) in MOFs is a versatile approach to the rational synthesis of size-selected clusters, including noble metals, on a high surface area support.

Role of a Modulator in the Synthesis of Phase-Pure NU-1000.

NU-1000 is a robust, mesoporous metal-organic framework (MOF) with hexazirconium nodes ([Zr6O16H16]8+, referred to as oxo-Zr6 nodes) that can be synthesized by combining a solution of ZrOCl2·8H2O and a benzoic acid modulator in N,N-dimethylformamide with a solution of linker (1,3,6,8-tetrakis(p-benzoic acid)pyrene, referred to as H4TBAPy) and by aging at an elevated temperature. Typically, the resulting crystals are primarily composed of NU-1000 domains that crystallize with a more dense phase that shares structural similarity with NU-901, which is an MOF composed of the same linker molecules and nodes. Density differences between the two polymorphs arise from the differences in the node orientation: in NU-1000, the oxo-Zr6 nodes rotate 120° from node to node, whereas in NU-901, all nodes are aligned in parallel. Considering this structural difference leads to the hypothesis that changing the modulator from benzoic acid to a larger and more rigid biphenyl-4-carboxylic acid might lead to a stronger steric interaction between the modulator coordinating on the oxo-Zr6 node and misaligned nodes or linkers in the large pore and inhibit the growth of the more dense NU-901-like material, resulting in phase-pure NU-1000. Side-by-side reactions comparing the products of synthesis using benzoic acid or biphenyl-4-carboxylic acid as a modulator produce structurally heterogeneous crystals and phase-pure NU-1000 crystals. It can be concluded that the larger and more rigid biphenyl-4-carboxylate inhibits the incorporation of nodes with an alignment parallel to the neighboring nodes already residing in the crystal.

Assembly of dicobalt and cobalt-aluminum oxide clusters on metal-organic framework and nanocast silica supports.

NU-1000, a mesoporous metal-organic framework (MOF) featuring hexazirconium oxide nodes and 3 nm wide channels, was infiltrated with a reactive dicobalt complex to install dicobalt active sites onto the MOF nodes. The anchoring of the dicobalt complex onto NU-1000 occurred with a nearly ideal stoichiometry of one bimetallic complex per node and with the cobalt evenly distributed throughout the MOF particle. To access thermally robust multimetallic sites on an all-inorganic support, the modified NU-1000 materials containing either the dicobalt complex, or an analogous cobalt-aluminum species, were nanocast with silica. The resulting materials feature Co2 or Co-Al bimetallated hexazirconium oxide clusters within a silica matrix. The cobalt-containing materials are competent catalysts for the selective oxidation of benzyl alcohol to benzaldehyde. Catalytic activity depends on the number of cobalt ions per node, but does not vary significantly between the NU-1000 and silica supports. Hence, the multimetallic oxide clusters remain site-isolated and substrate-accessible within the nanocast materials.