Bridging molecular-scale interfacial science with continuum-scale models.

Solid-water interfaces are crucial for clean water, conventional and renewable energy, and effective nuclear waste management. However, reflecting the complexity of reactive interfaces in continuum-scale models is a challenge, leading to oversimplified representations that often fail to predict real-world behavior. This is because these models use fixed parameters derived by averaging across a wide physicochemical range observed at the molecular scale. Recent studies have revealed the stochastic nature of molecular-level surface sites that define a variety of reaction mechanisms, rates, and products even across a single surface. To bridge the molecular knowledge and predictive continuum-scale models, we propose to represent surface properties with probability distributions rather than with discrete constant values derived by averaging across a heterogeneous surface. This conceptual shift in continuum-scale modeling requires exponentially rising computational power. By incorporating our molecular-scale understanding of solid-water interfaces into continuum-scale models we can pave the way for next generation critical technologies and novel environmental solutions.

Steric Effects on Single-Molecule Conductance in Flat-Lying Phenanthrene.

A previous combined experimental and theoretical study found that the position of anchoring groups on a phenanthrene (PHE) backbone played a large role in determining the single-molecule conductance of the PHE derivative. However, a consistent 0.1 G0 feature was found across all PHE derivatives. To understand this, the previously investigated PHE derivatives were placed flat on a simulated Au substrate with a scanning tunneling microscope (STM) tip over PHE and conductance was calculated using the non-equilibrium Green's function technique in conjunction with density functional theory (NEGF-DFT). The location of the tip was varied to find the most conductive and most energetically favorable arrangements, which did not coincide. Furthermore, the variation in conductance found in erect junctions was not present when PHE derivatives were lying flat, with all derivatives calculated to have conductance values around 0.1 G0.

Reversible solvent interactions with UiO-67 metal-organic frameworks.

The utility of UiO-67 Metal-Organic Frameworks (MOFs) for practical applications requires a comprehensive understanding of intermolecular host-guest MOF-analyte interactions. To investigate intermolecular interactions between UiO-67 MOFs and complex molecules, it is useful to evaluate the interactions with simple polar and non-polar analytes. This problem is approached by investigating the interactions of polar (acetone and isopropanol) and non-polar (n-heptane) molecules with functionalized UiO-67 MOFs via temperature programmed desorption mass spectrometry and temperature programmed Fourier transform infrared spectroscopy. We find that isopropanol, acetone, and n-heptane bind reversibly and non-destructively to UiO-67 MOFs, where MOF and analyte functionality influence relative binding strengths (n-heptane ≈ isopropanol > acetone). During heating, all three analytes diffuse into the internal pore environment and directly interact with the µ3-OH groups located within the tetrahedral pores, evidenced by the IR response of ν(µ3-OH). We observe nonlinear changes in the infrared cross sections of the ν(CH) modes of acetone, isopropanol, and n-heptane following diffusion into UiO-67. Similarly, acetone's ν(C=O) infrared cross section increases dramatically when diffused into UiO-67. Ultimately, this in situ investigation provides insights into how individual molecular functional groups interact with UiO MOFs and enables a foundation where MOF interactions with complex molecular systems can be evaluated.

Cation Modifies Interfacial Water Structures on Platinum during Alkaline Hydrogen Electrocatalysis.

The molecular details of an electrocatalytic interface play an essential role in the production of sustainable fuels and value-added chemicals. Many electrochemical reactions exhibit strong cation-dependent activities, but how cations affect reaction kinetics is still elusive. We report the effect of cations (K+, Li+, and Ba2+) on the interfacial water structure using second-harmonic generation (SHG) and classical molecular dynamics (MD) simulation. The second- (χH2O(2)) and third-order (χH2O(3)) optical susceptibilities of water on Pt are smaller in the presence of Ba2+ compared to those of K+, suggesting that cations can affect the interfacial water orientation. MD simulation reproduces experimental SHG observations and further shows that the competition between cation hydration and interfacial water alignment governs the net water orientation. The impact of cations on interfacial water supports a cation hydration-mediated mechanism for hydrogen electrocatalysis; i.e., the reaction occurs via water dissociation followed by cation-assisted hydroxide/water exchange on Pt. Our study highlights the role of interfacial water in electrocatalysis and how innocent additives (such as cations) can affect the local electrochemical environment.

COF-supported zirconium oxyhydroxide as a versatile heterogeneous catalyst for Knoevenagel condensation and nerve agent hydrolysis.

A composite of catalytic Lewis acidic zirconium oxyhydroxides (8 wt %) and a covalent organic framework (COF) was synthesized. X-ray diffraction and infrared (IR) spectroscopy reveal that COF's structure is preserved after loading with zirconium oxyhydroxides. Electron microscopy confirms a homogeneous distribution of nano- to sub-micron-sized zirconium clusters in the COF. 3D X-ray tomography captures the micron-sized channels connecting the well-dispersed zirconium clusters on the COF. The crystalline ZrOx(OH)y@COF's nanostructure was model-optimized via simulated annealing methods. Using 0.8 mol % of the catalyst yielded a turnover number of 100-120 and a turnover frequency of 160-360 h-1 for Knoevenagel condensation in aqueous medium. Additionally, 2.2 mol % of catalyst catalyzes the hydrolysis of dimethyl nitrophenyl phosphate, a simulant of nerve agent Soman, with a conversion rate of 37% in 180 min. The hydrolytic detoxification of the live agent Soman is also achieved. Our study unveils COF-stabilized ZrOx(OH)y as a new class of zirconium-based Lewis + Bronsted-acid catalysts.

Modulation of Charge Transport through Single Molecules Induced by Solvent-Stabilized Intramolecular Charge Transfer.

The modulation of charge transport through single molecules can be established by using the intrinsic characteristics of molecules and the physical properties of their environment. Therefore, the impact of the solvent on the electronic properties of molecules in the junction and their charge transport behavior are of great interest. Here, for the first time, we focused on charge transport through dimethylaminobenzonitrile (DMABN). This molecule shows unique behavior, specifically noticeable electronic structure modulations in bulk solvents, e.g., dual fluorescence in a polar environment. Using the scanning tunneling microscopy break junction (STM-BJ) technique, we find an order of magnitude increase in conductance along with a second conductance value in polar solvents over nonpolar solvents. Inspired by the twisted intramolecular charge transfer (TICT) explanation of the famous dual fluorescence of DMABN in polar solvents, we hypothesize stabilization of twisted DMABN molecules in the junction in more polar solvents. Ab initio molecular dynamics (AIMD) simulations using density functional theory (DFT) show that DMABN can twist in the junction and have a larger dipole moment compared to planar DMABN junction geometries, supporting the hypothesis. The nonequilibrium Green's function with the DFT approach (NEGF-DFT) is used to calculate the conductance throughout the AIMD trajectory, finding a significant change in the frontier orbitals and transmission function at large internal twisting angles, which can explain the dual conductance in polar solvents in STM-BJ experiments.

Oxide- and Silicate-Water Interfaces and Their Roles in Technology and the Environment.

Interfacial reactions drive all elemental cycling on Earth and play pivotal roles in human activities such as agriculture, water purification, energy production and storage, environmental contaminant remediation, and nuclear waste repository management. The onset of the 21st century marked the beginning of a more detailed understanding of mineral aqueous interfaces enabled by advances in techniques that use tunable high-flux focused ultrafast laser and X-ray sources to provide near-atomic measurement resolution, as well as by nanofabrication approaches that enable transmission electron microscopy in a liquid cell. This leap into atomic- and nanometer-scale measurements has uncovered scale-dependent phenomena whose reaction thermodynamics, kinetics, and pathways deviate from previous observations made on larger systems. A second key advance is new experimental evidence for what scientists hypothesized but could not test previously, namely, interfacial chemical reactions are frequently driven by "anomalies" or "non-idealities" such as defects, nanoconfinement, and other nontypical chemical structures. Third, progress in computational chemistry has yielded new insights that allow a move beyond simple schematics, leading to a molecular model of these complex interfaces. In combination with surface-sensitive measurements, we have gained knowledge of the interfacial structure and dynamics, including the underlying solid surface and the immediately adjacent water and aqueous ions, enabling a better definition of what constitutes the oxide- and silicate-water interfaces. This critical review discusses how science progresses from understanding ideal solid-water interfaces to more realistic systems, focusing on accomplishments in the last 20 years and identifying challenges and future opportunities for the community to address. We anticipate that the next 20 years will focus on understanding and predicting dynamic transient and reactive structures over greater spatial and temporal ranges as well as systems of greater structural and chemical complexity. Closer collaborations of theoretical and experimental experts across disciplines will continue to be critical to achieving this great aspiration.

Superhydrophilicity of α-alumina surfaces results from tight binding of interfacial waters to specific aluminols.

HYPOTHESIS: Understanding the microscopic driving force of water wetting is challenging and important for design of materials. The relations between structure, dynamics and hydrogen bonds of interfacial water can be investigated using molecular dynamics simulations. EXPERIMENTS AND SIMULATIONS: Contact angles at the alumina (0001) and (112‾0) surfaces are studied using both classical molecular dynamics simulations and experiments. To test the superhydrophilicity, the free energy cost of removing waters near the interfaces are calculated using the density fluctuations method. The strength of hydrogen bonds is determined by their lifetime and geometry. FINDINGS: Both surfaces are superhydrophilic and the (0001) surface is more hydrophilic. Interactions between surfaces and interfacial waters promote a templating effect whereby the latter are aligned in a pattern that follows the underlying lattice of the surfaces. Translational and rotational dynamics of interfacial water molecules are slower than in bulk water. Hydrogen bonds between water and both surfaces are asymmetric, water-to-aluminol ones are stronger than aluminol-to-water ones. Molecular dynamics simulations eliminate the impacts of surface contamination when measuring contact angles and the results reveal the microscopic origin of the macroscopic superhydrophilicity of alumina surfaces: strong water-to-aluminol hydrogen bonds.

Identifying UiO-67 Metal-Organic Framework Defects and Binding Sites through Ammonia Adsorption.

Ammonia is a widely used toxic industrial chemical that can cause severe respiratory ailments. Therefore, understanding and developing materials for its efficient capture and controlled release is necessary. One such class of materials is 3D porous metal-organic frameworks (MOFs) with exceptional surface areas and robust structures, ideal for gas storage/transport applications. Herein, interactions between ammonia and UiO-67-X (X: H, NH2 , CH3 ) zirconium MOFs were studied under cryogenic, ultrahigh vacuum (UHV) conditions using temperature-programmed desorption mass spectrometry (TPD-MS) and in-situ temperature-programmed infrared (TP-IR) spectroscopy. Ammonia was observed to interact with µ3 -OH groups present on the secondary building unit of UiO-67-X MOFs via hydrogen bonding. TP-IR studies revealed that under cryogenic UHV conditions, UiO-67-X MOFs are stable towards ammonia sorption. Interestingly, an increase in the intensity of the C-H stretching mode of the MOF linkers was detected upon ammonia exposure, attributed to NH-π interactions with linkers. These same binding interactions were observed in grand canonical Monte Carlo simulations. Based on TPD-MS, binding strength of ammonia to three MOFs was determined to be approximately 60 kJ mol-1 , suggesting physisorption of ammonia to UiO-67-X. In addition, missing linker defect sites, consisting of H2 O coordinated to Zr4+ sites, were detected through the formation of nNH3 ⋅H2 O clusters, characterized through in-situ IR spectroscopy. Structures consistent with these assignments were identified through density functional theory calculations. Tracking these bands through adsorption on thermally activated MOFs gave insight into the dehydroxylation process of UiO-67 MOFs. This highlights an advantage of using NH3 for the structural analysis of MOFs and developing an understanding of interactions between ammonia and UiO-67-X zirconium MOFs, while also providing directions for the development of stable materials for efficient toxic gas sorption.

Tuning the Lewis acidity of metal-organic frameworks for enhanced catalysis.

The kinetics of hydrolysis of dimethyl nitrophenyl phosphate (DMNP), a simulant of the nerve agent Soman, was studied and revealed transition metal salts as catalysts. The relative rates of DMNP hydrolysis by zirconium and hafnium chlorides are in accordance with their Lewis acidity. In situ conversion of zirconium chloride to zirconium oxy-hydroxide was identified as the key step. We propose a precursor-MOF activity relationship.

Probing Heterogeneous Charge Distributions at the α-Al2O3(0001)/H2O Interface.

Unlike metal or semiconductor electrodes, the surface charge resulting from the protonation or deprotonation of insulating mineral oxides is highly localized and heterogeneous in nature. In this work the Stark active C≡N stretch of potassium thiocyanate is used as a molecular probe of the heterogeneity of the interfacial electrostatic potential at the α-Al2O3(0001)/H2O interface. Vibrational sum frequency generation (vSFG) measurements performed in the OH stretching region suggest that thiocyanate species organize interfacial water similarly to halide ions. Changes in the electrostatic potential are then tracked via Stark shifts of the vibrational frequency of the thiocyanate stretch. Our vSFG measurements show that we can simultaneously measure the vSFG response of SCN- ions experiencing charged and neutral surface sites. We assign local potentials of +308 and -154 mV to positively and negatively charged aluminol groups that are present at pH = 4 and pH = 10, respectively. Thiocyanate anions at positively charged surface sites and negatively charged surface sites and those participating in contact ion pairing adopt similar orientations and are oppositely oriented relative to thiocyanate ions near neutral surface sites. All four species followed Langmuir adsorption isotherms. Density functional theory-molecular dynamics (DFT-MD) simulations of SCN- near the neutral α-Al2O3(0001)/H2O interface show that the vSFG response in the C≡N stretch region originates from a SCN-H-O-Al complex, suggesting the surface site specificity of these experiments. To our knowledge this is the first spectroscopic measurement of local potentials associated with a heterogeneously charged surface. The ability to probe the evolution of local charges in situ could provide vital insight into many industrial, electrochemical, and geochemically relevant interfaces.

Ions Tune Interfacial Water Structure and Modulate Hydrophobic Interactions at Silica Surfaces.

The structure and ultrafast dynamics of the electric double layer (EDL) are central to chemical reactivity and physical properties at solid/aqueous interfaces. While the Gouy-Chapman-Stern model is widely used to describe EDLs, it is solely based on the macroscopic electrostatic attraction of electrolytes for the charged surfaces. Structure and dynamics in the Stern layer are, however, more complex because of competing effects due to the localized surface charge distribution, surface-solvent-ion correlations, and the interfacial hydrogen bonding environment. Here, we report combined time-resolved vibrational sum frequency generation (TR-vSFG) spectroscopy with ab initio DFT-based molecular dynamics simulations (AIMD/DFT-MD) to get direct access to the molecular-level understanding of how ions change the structure and dynamics of the EDL. We show that innersphere adsorbed ions tune the hydrophobicity of the silica-aqueous interface by shifting the structural makeup in the Stern layer from dominant water-surface interactions to water-water interactions. This drives an initially inhomogeneous interfacial water coordination landscape observed at the neat interface toward a homogeneous, highly interconnected in-plane 2D hydrogen bonding (2D-HB) network at the ionic interface, reminiscent of the canonical, hydrophobic air-water interface. This ion-induced transformation results in a characteristic decrease of the vibrational lifetime (T1) of excited interfacial O-H stretching modes from T1 ∼ 600 fs to T1 ∼ 250 fs. Hence, we propose that the T1 determined by TR-vSFG in combination with DFT-MD simulations can be widely used for a quantitative spectroscopic probe of the ion kosmotropic/chaotropic effect at aqueous interfaces as well as of the ion-induced surface hydrophobicity.

First-Principles Calculation of Water pKa Using the Newly Developed SCAN Functional.

Acid/base chemistry is an intriguing topic that still constitutes a challenge for computational chemistry. While estimating the acid dissociation constant (or pKa) could shed light on many chemistry processes, especially in the fields of biochemistry and geochemistry, evaluating the relative stability between protonated and nonprotonated species is often very difficult. Indeed, a prerequisite for calculating the pKa of any molecule is an accurate description of the energetics of water dissociation. Here, we applied constrained molecular dynamics simulations, a noncanonical sampling technique, to investigate the water deprotonation process by selecting the OH distance as the reaction coordinate. The calculation is based on density functional theory and the newly developed SCAN functional, which has shown excellent performance in describing water structure. This first benchmark of SCAN on a chemical reaction shows that this functional accurately models the energetics of proton transfer reactions in an aqueous environment. After taking Coulomb long-range corrections and nuclear quantum effects into account, the estimated water pKa is only 1.0 pKa unit different from the target experimental value. Our results show that the combination of SCAN and constrained MD successfully reproduces the chemistry of water and constitutes a good framework for calculating the free energy of chemical reactions of interest.

Anisotropic Conductivity at the Single-Molecule Scale.

In most junctions built by wiring a single molecule between two electrodes, the electrons flow along only one axis: between the two anchoring groups. However, molecules can be anisotropic, and an orientation-dependent conductance is expected. Here, we fabricated single-molecule junctions by using the electrode potential to control the molecular orientation and access individual elements of the conductivity tensor. We measured the conductance in two directions, along the molecular plane as the benzene ring bridges two electrodes using anchoring groups (upright) and orthogonal to the molecular plane with the molecule lying flat on the substrate (planar). The perpendicular (planar) conductance is about 400 times higher than that along the molecular plane (upright). This offers a new method for designing a reversible room-temperature single-molecule electromechanical switch that controllably employs the electrode potential to orient the molecule in the junction in either "ON" or "OFF" conductance states.

Potential-Induced High-Conductance Transport Pathways through Single-Molecule Junctions.

Employing single molecules as electronic circuit building blocks is one promising approach to electronic device miniaturization. We report single-molecule junction formation where the orientation of molecules can be controlled externally by the working electrode potential. The scanning tunneling microscopy break junction (STM-BJ) method is used to bridge tetrafluoroterephthalic acid (TFTPA) and terephthalic acid (TPA) molecules between the Au(111) electrode and the STM tip to measure the single-molecule conductance through the junction. When the Au(111) electrode is at negative potentials (with respect to the zero-charge potential), a highly ordered and flat-oriented superstructure forms, allowing for direct contact between the π system of the benzene ring of the molecules and the Au(111) electrode, leading to junction formation with no anchoring group involvement. Our first-principles nonequilibrium Green's function (NEGF) computation shows a flat configuration yields a conductance that is 3 orders of magnitude larger than for a molecule vertically connected to the electrodes via anchoring groups. Conductances of 0.24 ± 0.04 and 0.22 ± 0.02 G0 are experimentally measured with the flat configurations of TFTPA and TPA, respectively. These values are at least 2 orders of magnitude higher than the experimental values previously reported for the conductance of TPA bridged through carboxylic acid anchoring groups (3.8 × 10-4-3.2 × 10-3 G0). In contrast, a positively charged surface triggers an order-disorder transition eliminating the high-conductance states, most likely because the formation of the flat-oriented junction is prevented. The dependence of TFTPA conductance on the electrode potential (electrode Fermi level) suggests a LUMO mediated transport mechanism. Calculation confirms the lack of an effect of the addition of an electron-withdrawing group are investigated.

Bond-Dependent Thole Model for Polarizability and Spectroscopy.

We present a new model for the calculation of molecular polarizabilities from effective atomic polarizabilities. This model is based on the Thole modified dipole interaction model for molecular polarizabilities, where the total polarizability is computed as a sum of effective atomic polarizabilities modified by dipole-dipole interactions. We extend this model by making the atomic polarizabilities explicit functions of the interatomic distances, scaling them by the radius of the volume an atom occupies in a molecule. We use the SCAN functional to show that this model, denoted TholeL, yields accurate molecular polarizabilities with little dependence on the training set. We also demonstrate that the TholeL model yields accurate polarizabilities for configurations far from the ground state structure for a wide range of molecules. Finally, we show that the TholeL model can be used to generate accurate Raman spectra for water, crystalline urea, and urea in water from ab initio molecular dynamics simulations.

Effect of Functional and Electron Correlation on the Structure and Spectroscopy of the Al2O3(001)-H2O Interface.

Oxide-water interfaces are ubiquitous, with many applications in industry and the environment, yet there is a great deal of controversy over their properties and microscopic structure. This controversy stems, in part, from the unique H-bond networks formed at different surface terminations and mineral compositions. Density functional theory simulations of these interfaces require an accurate description of both the oxide mineral and water in diverse H-bond environments. Thus, herein we simulate the Al2O3(001)-H2O interface using the PBE, PBE-TS, RPBE, SCAN, and HSE06-TS functionals to determine how calculated interfacial properties depend on the choice of functional. We find that the structure of the first few layers of water at the surface is determined by electron correlation in a way that cannot be approximated using semiemipirical van der Waals corrections. Of the functionals investigated, we find that SCAN yields the most accurate interfacial structure, dynamics, and sum frequency generation spectrum. Furthermore, SCAN leads to a reduction in the order of the 2D H-bond network of water at the alumina surface predicted by GGA functionals, leading to a significant decrease in the anisotropy of the diffusion coefficient at the surface. We emphasize the importance of using a functional which accurately describes electron correlation for more complex oxides, such as transition-metal oxides, where electron correlation may play an even greater role in determining the structure and dynamics of the oxide-water interface.

Ultrabroadband mid-infrared noncollinear difference frequency generation in a silver thiogallate crystal.

We report the generation of ultrabroadband mid-infrared (mid-IR) pulses by noncollinear difference frequency mixing. The signal and the idler output beams of an optical parametric amplifier are combined in a silver thiogallate crystal (AgGaS2) to generate mid-infrared radiation. We show that a noncollinear geometry facilitates broadband phase matching. Spectral bandwidths up to 1750 cm-1 were obtained at an external noncollinear angle of 4.2 deg, which is more than three times broader than in a collinear geometry. The broadband spectrum is tunable in the range of 1500-4500 cm-1. Pulse energies up to 1 µJ were achieved. The broadband pulses were used in sum frequency generation in ZnSe and in vibrational absorption spectroscopy experiments of liquid samples.

Systematic Doping of Cobalt into Layered Manganese Oxide Sheets Substantially Enhances Water Oxidation Catalysis.

The effect on the electrocatalytic oxygen evolution reaction (OER) of cobalt incorporation into the metal oxide sheets of the layered manganese oxide birnessite was investigated. Birnessite and cobalt-doped birnessite were characterized by X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), Raman spectroscopy, and conductivity measurements. A cobalt:manganese ratio of 1:2 resulted in the most active catalyst for the OER. In particular, the overpotential (η) for the OER was 420 mV, significantly lower than the η = 780 mV associated with birnessite in the absence of Co. Furthermore, the Tafel slope for Co/birnessite was 81 mV/dec, in comparison to a Tafel slope of greater than 200 mV/dec for birnessite. For chemical water oxidation catalysis, an 8-fold turnover number (TON) was achieved (h = 70 mmol of O2/mol of metal). Density functional theory (DFT) calculations predict that cobalt modification of birnessite resulted in a raising of the valence band edge and occupation of that edge by holes with enhanced mobility during catalysis. Inclusion of extra cobalt beyond the ideal 1:2 ratio was detrimental to catalysis due to disruption of the layered structure of the birnessite phase.