Insights into the structure and dynamics of K+ ions at the muscovite-water interface from machine learning potential simulations.

The surfaces of many minerals are covered by naturally occurring cations that become partially hydrated and can be replaced by hydronium or other cations when the surface is exposed to water or an aqueous solution. These ion exchange processes are relevant to various chemical and transport phenomena, yet elucidating their microscopic details is challenging for both experiments and simulations. In this work, we make a first step in this direction by investigating the behavior of the native K+ ions at the interface between neat water and the muscovite mica (001) surface with ab-initio-based machine learning molecular dynamics and enhanced sampling simulations. Our results show that the desorption of the surface K+ ions in pure ion-free water has a significant free energy barrier irrespective of their local surface arrangement. In contrast, facile K+ diffusion between mica's ditrigonal cavities characterized by different Al/Si orderings is observed. This behavior suggests that the K+ ions may favor a dynamic disordered surface arrangement rather than complete desorption when exposed to deionized water.

Correlated electron-nuclear dynamics of photoinduced water dissociation on rutile TiO2.

Elucidating the mechanism of photoinduced water splitting on TiO2 is important for advancing the understanding of photocatalysis and the ability to control photocatalytic surface reactions. However, incomplete experimental information and complex coupled electron-nuclear motion make the microscopic understanding challenging. Here we analyse the atomic-scale pathways of photogenerated charge carrier transport and photoinduced water dissociation at the prototypical water-rutile TiO2(110) interface using first-principles dynamics simulations. Two distinct mechanisms are observed. Field-initiated electron migration leads to adsorbed water dissociation via proton transfer to a surface bridging oxygen. In the other pathway, adsorbed water dissociation occurs via proton donation to a second-layer water molecule coupled to photoexcited-hole transfer promoted by in-plane surface lattice distortions. Two stages of non-adiabatic in-plane lattice motion-expansion and recovery-are observed, which are closely associated with population changes in Ti3d orbitals. Controlling such highly correlated electron-nuclear dynamics may provide opportunities for boosting the performance of photocatalytic materials.

Aqueous Titania Interfaces.

Water-metal oxide interfaces are central to many phenomena and applications, ranging from material corrosion and dissolution to photoelectrochemistry and bioengineering. In particular, the discovery of photocatalytic water splitting on TiO2 has motivated intensive studies of water-TiO2 interfaces for decades. So far, a broad understanding of the interaction of water vapor with several TiO2 surfaces has been obtained. However, much less is known about liquid water-TiO2 interfaces, which are more relevant to many practical applications. Probing these complex systems at the molecular level is experimentally challenging and is sometimes possible only through computational studies. This review summarizes recent advances in the atomistic understanding, mostly through computational simulations, of the structure and dynamics of interfacial water on TiO2 surfaces. The main focus is on the nature, molecular or dissociated, of water in direct contact with low-index defect-free crystalline surfaces. The hydroxyls resulting from water dissociation are essential in the photooxidation of water and critically affect the surface chemistry of TiO2.

A first-principles machine-learning force field for heterogeneous ice nucleation on microcline feldspar.

The formation of ice in the atmosphere affects precipitation and cloud properties, and plays a key role in the climate of our planet. Although ice can form directly from liquid water under deeply supercooled conditions, the presence of foreign particles can aid ice formation at much warmer temperatures. Over the past decade, experiments have highlighted the remarkable efficiency of feldspar minerals as ice nuclei compared to other particles present in the atmosphere. However, the exact mechanism of ice formation on feldspar surfaces has yet to be fully understood. Here, we develop a first-principles machine-learning model for the potential energy surface aimed at studying ice nucleation at microcline feldspar surfaces. The model is able to reproduce with high-fidelity the energies and forces derived from density-functional theory (DFT) based on the SCAN exchange and correlation functional. Our training set includes configurations of bulk supercooled water, hexagonal and cubic ice, microcline, and fully-hydroxylated feldspar surfaces exposed to a vacuum, liquid water, and ice. We apply the machine-learning force field to study different fully-hydroxylated terminations of the (100), (010), and (001) surfaces of microcline exposed to a vacuum. Our calculations suggest that terminations that do not minimize the number of broken bonds are preferred in a vacuum. We also study the structure of supercooled liquid water in contact with microcline surfaces, and find that water density correlations extend up to around 10 Å from the surfaces. Finally, we show that the force field maintains a high accuracy during the simulation of ice formation at microcline surfaces, even for large systems of around 30 000 atoms. Future work will be directed towards the calculation of nucleation free-energy barriers and rates using the force field developed herein, and understanding the role of different microcline surfaces in ice nucleation.

Probing the self-ionization of liquid water with ab initio deep potential molecular dynamics.

The chemical equilibrium between self-ionized and molecular water dictates the acid-base chemistry in aqueous solutions, yet understanding the microscopic mechanisms of water self-ionization remains experimentally and computationally challenging. Herein, Density Functional Theory (DFT)-based deep neural network (DNN) potentials are combined with enhanced sampling techniques and a global acid-base collective variable to perform extensive atomistic simulations of water self-ionization for model systems of increasing size. The explicit inclusion of long-range electrostatic interactions in the DNN potential is found to be crucial to accurately reproduce the DFT free energy profile of solvated water ion pairs in small (64 and 128 H2O) cells. The reversible work to separate the hydroxide and hydronium to a distance [Formula: see text] is found to converge for simulation cells containing more than 500 H2O, and a distance of [Formula: see text] 8 Å is the threshold beyond which the work to further separate the two ions becomes approximately zero. The slow convergence of the potential of mean force with system size is related to a restructuring of water and an increase of the local order around the water ions. Calculation of the dissociation equilibrium constant illustrates the key role of long-range electrostatics and entropic effects in the water autoionization process.

Elucidating the water-anatase TiO2(101) interface structure using infrared signatures and molecular dynamics.

The structure and dynamics of water on solid surfaces critically affect the chemistry of materials in ambient and aqueous environments. Here, we investigate the hydrogen bonding network of water adsorbed on the majority (101) surface of anatase TiO2, a widely used photocatalyst, using polarization- and azimuth-resolved infrared spectroscopy combined with neural network potential molecular dynamics simulations. Our results show that one monolayer of water saturates the undercoordinated titanium (Ti5c) sites, forming one-dimensional chains of molecule hydrogen bonded to surface undercoordinated bridging oxygen (O2c) atoms. As the coverage increases, water adsorption on O2c sites leads to significant restructuring of the water monolayer and the formation of a two-dimensional hydrogen bond network characterized by tightly bound pairs of water molecules on adjacent Ti5c and O2c sites. This structural motif likely persists at ambient conditions, influencing the reactions occurring there. The results reported here provide critical details of the structure of the water-anatase (101) interface that were previously hypothesized but unconfirmed experimentally.

Acid-Base Chemistry of a Model IrO2 Catalytic Interface.

Iridium oxide (IrO2) is one of the most efficient catalytic materials for the oxygen evolution reaction (OER), yet the atomic scale structure of its aqueous interface is largely unknown. Herein, the hydration structure, proton transfer mechanisms, and acid-base properties of the rutile IrO2(110)-water interface are investigated using ab initio based deep neural-network potentials and enhanced sampling simulations. The proton affinities of the different surface sites are characterized by calculating their acid dissociation constants, which yield a point of zero charge in agreement with experiments. A large fraction (≈80%) of adsorbed water dissociation is observed, together with a short lifetime (≈0.5 ns) of the resulting terminal hydroxy groups, due to rapid proton exchanges between adsorbed H2O and adjacent OH species. This rapid surface proton transfer supports the suggestion that the rate-determining step in the OER may not involve proton transfer across the double layer into solution, as indicated by recent experiments.

Modeling the aqueous interface of amorphous TiO2 using deep potential molecular dynamics.

Amorphous titanium dioxide (a-TiO2) is widely used as a coating material in applications such as electrochemistry and self-cleaning surfaces where its interface with water has a central role. However, little is known about the structures of the a-TiO2 surface and aqueous interface, particularly at the microscopic level. In this work, we construct a model of the a-TiO2 surface via a cut-melt-and-quench procedure based on molecular dynamics simulations with deep neural network potentials (DPs) trained on density functional theory data. After interfacing the a-TiO2 surface with water, we investigate the structure and dynamics of the resulting system using a combination of DP-based molecular dynamics (DPMD) and ab initio molecular dynamics (AIMD) simulations. Both AIMD and DPMD simulations reveal that the distribution of water on the a-TiO2 surface lacks distinct layers normally found at the aqueous interface of crystalline TiO2, leading to an ∼10 times faster diffusion of water at the interface. Bridging hydroxyls (Ti2-ObH) resulting from water dissociation decay several times more slowly than terminal hydroxyls (Ti-OwH) due to fast Ti-OwH2 → Ti-OwH proton exchange events. These results provide a basis for a detailed understanding of the properties of a-TiO2 in electrochemical environments. Moreover, the procedure of generating the a-TiO2-interface employed here is generally applicable to studying the aqueous interfaces of amorphous metal oxides.

Observation of electron orbital signatures of single atoms within metal-phthalocyanines using atomic force microscopy.

Resolving the electronic structure of a single atom within a molecule is of fundamental importance for understanding and predicting chemical and physical properties of functional molecules such as molecular catalysts. However, the observation of the orbital signature of an individual atom is challenging. We report here the direct identification of two adjacent transition-metal atoms, Fe and Co, within phthalocyanine molecules using high-resolution noncontact atomic force microscopy (HR-AFM). HR-AFM imaging reveals that the Co atom is brighter and presents four distinct lobes on the horizontal plane whereas the Fe atom displays a "square" morphology. Pico-force spectroscopy measurements show a larger repulsion force of about 5 pN on the tip exerted by Co in comparison to Fe. Our combined experimental and theoretical results demonstrate that both the distinguishable features in AFM images and the variation in the measured forces arise from Co's higher electron orbital occupation above the molecular plane. The ability to directly observe orbital signatures using HR-AFM should provide a promising approach to characterizing the electronic structure of an individual atom in a molecular species and to understand mechanisms of certain chemical reactions.

Activating Bi p-orbitals in Dispersed Clusters of Amorphous BiOx for Electrocatalytic Nitrogen Reduction.

Catalytic strategies based on main group metals are significantly less advanced than those of transition metal catalysis, leaving untapped areas of potentially fruitful research. We here demonstrate an effective approach for the modulation of Bi 6p energy levels during the construction of atomically dispersed clusters of amorphous BiOx . Bi oxidation state is proposed to strongly affects the nitrogen fixation activity, with the half-occupied pz orbitals of the Bi2+ ions being highly efficient toward electron injection into the inert N2 molecule. With sufficient catalytic sites to adsorb and activate N2 , the bonding between N2 and catalyst is able to be in situ identified. The catalyst shows an outstanding Faraday efficiency (≈30 %) and high yield (≈113 µg h-1 mg-1 cat ) in NH3 production, outperforming most of the existing catalysts in aqueous solution. These results lay the basis for developing the potential of p-block elements for catalysis of multi-electron reactions.

Water dissociation at the water-rutile TiO2(110) interface from ab initio-based deep neural network simulations.

The interaction of water with TiO2 surfaces is of crucial importance in various scientific fields and applications, from photocatalysis for hydrogen production and the photooxidation of organic pollutants to self-cleaning surfaces and bio-medical devices. In particular, the equilibrium fraction of water dissociation at the TiO2-water interface has a critical role in the surface chemistry of TiO2, but is difficult to determine both experimentally and computationally. Among TiO2 surfaces, rutile TiO2(110) is of special interest as the most abundant surface of TiO2's stable rutile phase. While surface-science studies have provided detailed information on the interaction of rutile TiO2(110) with gas-phase water, much less is known about the TiO2(110)-water interface, which is more relevant to many applications. In this work, we characterize the structure of the aqueous TiO2(110) interface using nanosecond timescale molecular dynamics simulations with ab initio-based deep neural network potentials that accurately describe water/TiO2(110) interactions over a wide range of water coverages. Simulations on TiO2(110) slab models of increasing thickness provide insight into the dynamic equilibrium between molecular and dissociated adsorbed water at the interface and allow us to obtain a reliable estimate of the equilibrium fraction of water dissociation. We find a dissociation fraction of 22 ± 6% with an associated average hydroxyl lifetime of 7.6 ± 1.8 ns. These quantities are both much larger than corresponding estimates for the aqueous anatase TiO2(101) interface, consistent with the higher water photooxidation activity that is observed for rutile relative to anatase.

Modeling the Solvation and Acidity of Carboxylic Acids Using an Ab Initio Deep Neural Network Potential.

Formic and acetic acid constitute the simplest of carboxylic acids, yet they exhibit fascinating chemistry in the condensed phase such as proton transfer and dimerization. The go-to method of choice for modeling these rare events have been accurate but expensive ab initio molecular dynamics simulations. In this study, we present a deep neural network potential trained using accurate ab initio data that can be used in tandem with enhanced-sampling methods to perform an efficient exploration of the free-energy surface of aqueous solutions of weak carboxylic acids. In particular, we show that our model captures proton dissociation and provides a good estimate of the pKa, as well as the dimerization of formic and acetic acid. This provides a suitable starting point for applications in different research areas where computational efficiency coupled with the accuracy of ab initio methods is required.

Realizing Two-Electron Transfer in Ni(OH)2 Nanosheets for Energy Storage.

The theoretical capacity of a given electrode material is ultimately determined by the number of electrons transferred in each redox center. The design of multi-electron transfer processes could break through the limitation of one-electron transfer and multiply the total capacity but is difficult to achieve because multiple electron transfer processes are generally thermodynamically and kinetically more complex. Here, we report the discovery of two-electron transfer in monolayer Ni(OH)2 nanosheets, which contrasts with the traditional one-electron transfer found in multilayer materials. First-principles calculations predict that the first oxidation process Ni2+ → Ni3+ occurs easily, whereas the second electron transfer in Ni3+ → Ni4+ is strongly hindered in multilayer materials by both the interlayer hydrogen bonds and the domain H structure induced by the Jahn-Teller distortion of the Ni3+ (t2g6eg1)-centered octahedra. In contrast, the second electron transfer can easily occur in monolayers because all H atoms are fully exposed. Experimentally, the as-prepared monolayer is found to deliver an exceptional redox capacity of ∼576 mA h/g, nearly 2 times the theoretical capacity of one-electron processes. In situ experiments demonstrate that monolayer Ni(OH)2 can transfer two electrons and most Ni ions transform into Ni4+ during the charging process, whereas bulk Ni(OH)2 can only be transformed partially. Our work reveals a new redox reaction mechanism in atomically thin Ni(OH)2 nanosheets and suggests a promising path toward tuning the electron transfer numbers to multiply the capacity of the relevant energy storage materials.

Pathways for Electron Transfer at MgO-Water Interfaces from Ab Initio Molecular Dynamics.

The nature of electron transfer across metal oxide-water interfaces depends significantly on the band gap of the oxide and its band edge energies relative to the potentials of relevant aqueous redox couples. Here we focus on the water interface with MgO, a prototypical wide band gap oxide whose conduction band edge is close in energy to that of water. We investigate the behavior of an excess electron at and out of equilibrium near the interface using ab initio molecular dynamics based on hybrid density functional theory. Our simulations show that under equilibrium conditions the excess electron (donated by an Al impurity in MgO) localizes to a midgap defect state comparable in energy and shape to a hydrated electron in bulk water. To characterize the electron transfer from the conduction band of MgO to interfacial product states, we dope near-equilibrium configurations of the pristine MgO-water system with Al and run short trajectories of these instantaneously out-of-equilibrium systems. We observe two distinct products associated with the excess electron: a surface-localized electron (esurf-) and an aqueous hydrogen radical (H•). The H• pathway exhibits a much higher activation barrier despite being more exoergic, making esurf- the kinetic product. Our characterization of the pathways on the basis of Marcus theory is consistent with the poor observed utility of MgO for water radiolysis. Moreover, we anticipate that the computational framework employed here will be broadly applicable to assessing electron transfer mechanisms at aqueous, photocatalytic interfaces.

Breaking a dative bond with mechanical forces.

Bond breaking and forming are essential components of chemical reactions. Recently, the structure and formation of covalent bonds in single molecules have been studied by non-contact atomic force microscopy (AFM). Here, we report the details of a single dative bond breaking process using non-contact AFM. The dative bond between carbon monoxide and ferrous phthalocyanine was ruptured via mechanical forces applied by atomic force microscope tips; the process was quantitatively measured and characterized both experimentally and via quantum-based simulations. Our results show that the bond can be ruptured either by applying an attractive force of ~150 pN or by a repulsive force of ~220 pN with a significant contribution of shear forces, accompanied by changes of the spin state of the system. Our combined experimental and computational studies provide a deeper understanding of the chemical bond breaking process.

Increasing Iridium Oxide Activity for the Oxygen Evolution Reaction with Hafnium Modification.

Synthesis and implementation of highly active, stable, and affordable electrocatalysts for the oxygen evolution reaction (OER) is a major challenge in developing energy efficient and economically viable energy conversion devices such as electrolyzers, rechargeable metal-air batteries, and regenerative fuel cells. The current benchmark electrocatalyst for OER is based on iridium oxide (IrOx) due to its superior performance and excellent stability. However, large scale applications using IrOx are impractical due to its low abundance and high cost. Herein, we report a highly active hafnium-modified iridium oxide (IrHfxOy) electrocatalyst for OER. The IrHfxOy electrocatalyst demonstrated ten times higher activity in alkaline conditions (pH = 11) and four times higher activity in acid conditions (pH = 1) than a IrOx electrocatalyst. The highest intrinsic mass activity of the IrHfxOy catalyst in acid conditions was calculated as 6950 A gIrOx-1 at an overpotential (η) of 0.3 V. Combined studies utilizing operando surface enhanced Raman spectroscopy (SERS) and DFT calculations revealed that the active sites for OER are the Ir-O species for both IrOx and IrHfxOy catalysts. The presence of Hf sites leads to more negative charge states on nearby O sites, shortening of the bond lengths of Ir-O, and lowers free energies for OER intermediates that accelerate the OER process.

Hydrogen Bonds and H3O+ Formation at the Water Interface with Formic Acid Covered Anatase TiO2.

Carboxylic acid-modified TiO2 surfaces in aqueous environment are of widespread interest, yet atomic-scale understanding of their structure is limited. We here investigate formic acid (FA) on anatase TiO2 (101) (A-101) in contact with water using density functional theory (DFT) and ab initio molecular dynamics (AIMD). Isolated FA molecules adsorbed in a deprotonated bridging bidentate (BD) form on A-101 are found to remain stable at the interface with water, with the acid proton transferred to a surface oxygen to form a surface bridging hydroxyl (ObrH). With increasing FA coverage, adsorbed monolayers of only BD and successively of alternating monodentate (MD) and BD species give rise to a higher concentration of surface ObrH's. Simulations of these adsorbed monolayers in water environment show that some protons are released from the surface ObrH's to water resulting in a negatively charged surface with nearby solvated H3O+ ions. These results provide insight into the complex acid-base equilibrium between an oxide surface, adsorbates and water and can also help obtain a better understanding of the wetting properties of chemically modified TiO2 surfaces.

Effects of applied voltage on water at a gold electrode interface from ab initio molecular dynamics.

Electrode-water interfaces under voltage bias demonstrate anomalous electrostatic and structural properties that are influential in their catalytic and technological applications. Mean-field and empirical models of the electrical double layer (EDL) that forms in response to an applied potential do not capture the heterogeneity that polarizable, liquid-phase water molecules engender. To illustrate the inhomogeneous nature of the electrochemical interface, Born-Oppenheimer ab initio molecular dynamics calculations of electrified Au(111) slabs interfaced with liquid water were performed using a combined explicit-implicit solvent approach. The excess charges localized on the model electrode were held constant and the electrode potentials were computed at frequent simulation times. The electrode potential in each trajectory fluctuated with changes in the atomic structure, and the trajectory-averaged potentials converged and yielded a physically reasonable differential capacitance for the system. The effects of the average applied voltages, both positive and negative, on the structural, hydrogen bonding, dynamical, and vibrational properties of water were characterized and compared to literature where applicable. Controlled-potential simulations of the interfacial solvent dynamics provide a framework for further investigation of more complex or reactive species in the EDL and broadly for understanding electrochemical interfaces in situ.

Polaron-Adsorbate Coupling at the TiO2(110)-Carboxylate Interface.

Understanding how adsorbates influence polaron behavior is of fundamental importance in describing the catalytic properties of TiO2. Carboxylic acids adsorb readily at TiO2 surfaces, yet their influence on polaronic states is unknown. Using UV photoemission spectroscopy (UPS), two-photon photoemission spectroscopy (2PPE), and density functional theory (DFT) we show that dissociative adsorption of formic and acetic acids has profound, yet different, effects on the surface density, crystal field, and photoexcitation of polarons in rutile TiO2(110). We also show that these variations are governed by the contrasting electrostatic properties of the acids, which impacts the extent of polaron-adsorbate coupling. The density of polarons in the surface region increases more in formate-terminated TiO2(110) relative to acetate. Consequently, increased coupling gives rise to new photoexcitation channels via states 3.83 eV above the Fermi level. The onset of this process is 3.45 eV, likely adding to the catalytic photoyield.

Hydration structure of flat and stepped MgO surfaces.

We investigate the solvation structure of flat and stepped MgO(001) in neutral liquid water using ab initio molecular dynamics based on a hybrid density functional with dispersion corrections. Our simulations show that the MgO surface is covered by a densely packed layer of mixed intact and dissociated adsorbed water molecules in a planar arrangement with strong intermolecular H-bonds. The water dissociation fractions in this layer are >20% and >30% on the flat and stepped surfaces, respectively. Slightly above the first water layer, we observe metastable OH groups perpendicular to the interface, similar to those reported in low temperature studies of water monolayers on MgO. These species receive hydrogen bonds from four nearby water molecules in the first layer and have their hydrophobic H end directed toward bulk water, while their associated protons are bound to surface oxygens. The formation of these OH species is attributed to the strong basicity of the MgO surface and can be relevant for understanding various phenomena from morphology evolution and growth of (nano)crystalline MgO particles to heterogeneous catalysis.