Structure, Polarization, and Sum Frequency Generation Spectrum of Interfacial Water on Anatase TiO2

Water Structure in the First Layers on TiO2: A Key Factor for Boosting Solar-Driven Water-Splitting Performances

The water hydrogen-bonded network is strongly perturbed in the first layers in contact with the semiconductor surface. Even though this aspect influences the outer-sphere electron transfer, it was not recognized that it is a crucial factor impacting the solar-driven water-splitting performances. To fill this gap, we have selected two TiO2 anatase samples (with and without B-doping), and by extensive experimental and computational investigations, we have demonstrated that the remarkable 5-fold increase in water-splitting photoactivity of the B-doped sample cannot be ascribed to effects typically associated to enhanced photocatalytic properties, such as band gap, heterojunctions, crystal facets, and other aspects. Studying these samples by combining FTIR measurements under controlled humidity with first-principles simulations sheds light on the role and nature of the first-layer water structure in contact with the photocatalyst surfaces. It turns out that the doping hampers the percolation of tetrahedrally coordinated water molecules while enhancing the population of topological H-bond defects forming approximately linear H-bonded chains. This work unveils how doping the semiconductor surface affects the local electric field, determining the water splitting rate by influencing the H-bond topologies in the first water layers. This evidence opens new prospects for designing efficient photocatalysts for water splitting.

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.

Adsorption of Glycine on TiO2 in Water from On-the-fly Free-Energy Calculations and In Situ Electrochemical Impedance Spectroscopy

We report here an experimental-computational study of hydrated TiO2 anatase nanoparticles interacting with glycine, where we obtain quantitative agreement of the measured adsorption free energies. Ab initio simulations are performed within the tight binding and density functional theory in combination with enhanced free-energy sampling techniques, which exploit the thermodynamic integration of the unbiased mean forces collected on-the-fly along the molecular dynamics trajectories. The experiments adopt a new and efficient setup for electrochemical impedance spectroscopy measurements based on portable screen-printed gold electrodes, which allows fast and in situ signal assessment. The measured adsorption free energy is -30 kJ/mol (both from experiment and calculation), with preferential interaction of the charged NH3+ group which strongly adsorbs on the TiO2 bridging oxygens. This highlights the importance of the terminal amino groups in the adsorption mechanism of amino acids on hydrated metal oxides. The excellent agreement between computation and experiment for this amino acid opens the doors to the exploration of the interaction free energies for other moderately complex bionano systems.

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.

Long-range proton and hydroxide ion transfer dynamics at the water/CeO2 interface in the nanosecond regime: reactive molecular dynamics simulations and kinetic analysis

The structural properties, dynamical behaviors, and ion transport phenomena at the interface between water and cerium oxide are investigated by reactive molecular dynamics (MD) simulations employing neural network potentials (NNPs). The NNPs are trained to reproduce density functional theory (DFT) results, and DFT-based MD (DFT-MD) simulations with enhanced sampling techniques and refinement schemes are employed to efficiently and systematically acquire training data that include diverse hydrogen-bonding configurations caused by proton hopping events. The water interfaces with two low-index surfaces of (111) and (110) are explored with these NNPs, and the structure and long-range proton and hydroxide ion transfer dynamics are examined with unprecedented system sizes and long simulation times. Various types of proton hopping events at the interface are categorized and analyzed in detail. Furthermore, in order to decipher the proton and hydroxide ion transport phenomena along the surface, a counting analysis based on the semi-Markov process is formulated and applied to the MD trajectories to obtain reaction rates by considering the transport as stochastic jump processes. Through this model, the coupling between hopping events, vibrational motions, and hydrogen bond networks at the interface are quantitatively examined, and the high activity and ion transport phenomena at the water/CeO2 interface are unequivocally revealed in the nanosecond regime.

Mesoscopic Model for Reversible Adsorption Stage of Albumin and Fibrinogen on TiO2 Surface

The competitive behavior of proteins in the reversible adsorption stage plays a crucial role in determining the composition of the protein layer and the subsequent biological responses to the biomaterial. However, such competitive adsorption is a mesoscopic process at physiological protein concentration, and neither a macroscopic experiment nor microscopic MD (molecular dynamics) simulation is suitable to clarify it. Here, we proposed a mesoscopic DPD (dissipative particle dynamics) model to illustrate the competitive process of albumin and fibrinogen on TiO2 surface with its parameters deduced from our previous MD simulation, and proved the model well retained the diffusion and adsorption properties of proteins in the competitive adsorption on the plane surface. We then applied the model to the competitive adsorption on the surfaces with different nanostructures and observed that when the nanostructure size is much larger than that of protein, the increase in surface area is the main influencing factor; when the nanostructure size is close to that of protein, the coordination between the nanostructure and the size and shape of protein significantly affects the competitive adsorption process. The model has revealed many mechanical phenomena observed in previous experimental studies and has the potential to contribute to the development of high-performance biomaterials.

Thermal transport across TiO2-H2O interface involving water dissociation: Ab initio-assisted deep potential molecular dynamics

Water dissociation on TiO2 surfaces has been known for decades and holds great potential in various applications, many of which require a proper understanding of thermal transport across the TiO2-H2O interface. Molecular dynamics (MD) simulations play an important role in characterizing complex systems' interfacial thermal transport properties. Nevertheless, due to the imprecision of empirical force field potentials, the interfacial thermal transport mechanism involving water dissociation remains to be determined. To cope with this, a deep potential (DP) model is formulated through the utilization of ab initio datasets. This model successfully simulates interfacial thermal transport accompanied by water dissociation on the TiO2 surfaces. The trained DP achieves a total energy accuracy of ∼238.8 meV and a force accuracy of ∼197.05 meV/Å. The DPMD simulations show that water dissociation induces the formation of hydrogen bonding networks and molecular bridges. Structural modifications further affect interfacial thermal transport. The interfacial thermal conductance estimated by DP is ∼8.54 × 109 W/m2 K, smaller than ∼13.17 × 109 W/m2 K by empirical potentials. The vibrational density of states (VDOS) quantifies the differences between the DP model and empirical potentials. Notably, the VDOS disparity between the adsorbed hydrogen atoms and normal hydrogen atoms demonstrates the influence of water dissociation on heat transfer processes. This work aims to understand the effect of water dissociation on thermal transport at the TiO2-H2O interface. The findings will provide valuable guidance for the thermal management of photocatalytic devices.

High-Throughput Condensed-Phase Hybrid Density Functional Theory for Large-Scale Finite-Gap Systems: The SeA Approach

High-throughput electronic structure calculations (often performed using density functional theory (DFT)) play a central role in screening existing and novel materials, sampling potential energy surfaces, and generating data for machine learning applications. By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal DFT and furnish a more accurate description of the underlying electronic structure, albeit at a computational cost that often prohibits such high-throughput applications. To address this challenge, we have constructed a robust, accurate, and computationally efficient framework for high-throughput condensed-phase hybrid DFT and implemented this approach in the PWSCF module of Quantum ESPRESSO (QE). The resulting SeA approach (SeA = SCDM + exx + ACE) combines and seamlessly integrates: (i) the selected columns of the density matrix method (SCDM, a robust noniterative orbital localization scheme that sidesteps system-dependent optimization protocols), (ii) a recently extended version of exx (a black-box linear-scaling EXX algorithm that exploits sparsity between localized orbitals in real space when evaluating the action of the standard/full-rank V^xx operator), and (iii) adaptively compressed exchange (ACE, a low-rank V^xx approximation). In doing so, SeA harnesses three levels of computational savings: pair selection and domain truncation from SCDM + exx (which only considers spatially overlapping orbitals on orbital-pair-specific and system-size-independent domains) and low-rank V^xx approximation from ACE (which reduces the number of calls to SCDM + exx during the self-consistent field (SCF) procedure). Across a diverse set of 200 nonequilibrium (H₂O)₆₄ configurations (with densities spanning 0.4-1.7 g/cm³), SeA provides a 1-2 order-of-magnitude speedup in the overall time-to-solution, i.e., ≈8-26× compared to the convolution-based PWSCF(ACE) implementation in QE and ≈78-247× compared to the conventional PWSCF(Full) approach, and yields energies, ionic forces, and other properties with high fidelity. As a proof-of-principle high-throughput application, we trained a deep neural network (DNN) potential for ambient liquid water at the hybrid DFT level using SeA via an actively learned data set with ≈8,700 (H₂O)₆₄ configurations. Using an out-of-sample set of (H₂O)₅₁₂ configurations (at nonambient conditions), we confirmed the accuracy of this SeA-trained potential and showcased the capabilities of SeA by computing the ground-truth ionic forces in this challenging system containing >1,500 atoms.

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.

pH- and Facet-Dependent Surface Chemistry of TiO2 in Aqueous Environment from First Principles

TiO₂ is a relevant catalytic material, and its chemistry in aqueous environment is a challenging aspect to address. Also, the morphology of TiO₂ particles at the nanoscale is often complex, spanning from faceted to spherical. In this work, we study the pH- and facet-dependent surface chemistry of TiO₂/water interfaces by performing ab initio molecular dynamics simulations with the grand canonical formulation of species in solution. We first determined the acid-base equilibrium constants at the interface, which allows us to estimate the pH at the point of zero charge, an important experimental observable. Then, based on simulated equilibrium constants, we predict the amount of H⁺, OH^-, and adsorbed H₂O species present on the surfaces as a function of the pH, a relevant aspect for water splitting semi-reactions. We approximated the complex morphology of TiO₂ particles by considering the rutile (110) and (011), and anatase (101), (001), and (100) surfaces.

Controlling Water Delivery to an Electrochemical Interface with Surfactants

Most electrochemical reactions require delivery of protons, often from water, to surface-adsorbed species. However, water also acts as a competitor to many such processes by directly reacting with the electrode, which necessitates using water in small amounts. Controlling the water content and structure near the surface is an important frontier in directing the reactivity and selectivity of electrochemical reactions. Surfactants accumulate near surfaces, and therefore, they can be used as agents to control interfacial water. Using mid-IR spectro-electrochemistry, we show that a modest concentration (1 mM) of the cationic surfactant CTAB in mixtures of 10 M water in an organic solvent (dDMSO) has a large effect on the interfacial water concentration, changing it by up to ∼35% in the presence of an applied potential. The major cause of water content change is displacement due to the accumulation or depletion of surfactants driven by potential. Two forces drive the surfactants to the electrode: the applied potential and the hydrophobic interactions with the water in the bulk. We have quantified their competition by varying the water content in the bulk. To our knowledge, for the first time, we have identified the electrochemical equivalent of the hydrophobic drive. For our system, a change in applied potential of 1 V has the same effect as adding a 0.55 mole fraction of water to the bulk. This work illustrates the significance of surfactants in the partitioning of water between the bulk and the surface and paves the way toward engineering interfacial water structures for controlling electrochemical reactions.

Water structure, dynamics and reactivity on a TiO2-nanoparticle surface: new insights from ab initio molecular dynamics

Water structure, dynamics and reactivity at the surface of a small TiO₂-nanoparticle fully immersed in water was investigated by an ab initio molecular dynamics simulation. Several modes of water binding were identified by assigning each atom to an atom type, representing a distinct chemical environment in the ab initio ensemble, and then computing radial distribution functions between the atom types. Surface reactivity was investigated by monitoring how populations of atom types change during the simulation. In order to acquire further insight, electron densities for a set of representative system snapshots were analyzed using an atoms-in-molecules approach. Our results reveal that water dissociation, where a water molecule splits at a bridging oxygen site to form a hydroxyl group and a protonated oxygen bridge, can occur by a mechanism involving transfer of a proton over several water molecules. The hydroxyl group and protonated oxygen bridge formed in the process persist (on a 10 ps time scale) and the hydroxyl group undergoes exchange using a mechanism similar to the one responsible for water dissociation. Rotational and translational dynamics of water molecules around the nanoparticle were analyzed in terms of reorientational time correlation functions and mean square displacement. While reorientation of water O-H vectors decreases quickly in the proximity of the nanoparticle surface, translational diffusion slows down more gradually. Our results give new insight into water structure, dynamics and reactivity on TiO₂-nanoparticle surfaces and suggest that water dissociation on curved TiO₂-nanoparticle surfaces can occur via more complex mechanisms than those previously identified for flat defect-free surfaces.

Improving the Accuracy of Atomistic Simulations of the Electrochemical Interface

Atomistic simulation of the electrochemical double layer is an ambitious undertaking, requiring quantum mechanical description of electrons, phase space sampling of liquid electrolytes, and equilibration of electrolytes over nanosecond time scales. All models of electrochemistry make different trade-offs in the approximation of electrons and atomic configurations, from the extremes of classical molecular dynamics of a complete interface with point-charge atoms to correlated electronic structure methods of a single electrode configuration with no dynamics or electrolyte. Here, we review the spectrum of simulation techniques suitable for electrochemistry, focusing on the key approximations and accuracy considerations for each technique. We discuss promising approaches, such as enhanced sampling techniques for atomic configurations and computationally efficient beyond density functional theory (DFT) electronic methods, that will push electrochemical simulations beyond the present frontier.

Theoretical Modeling of Electrochemical Proton-Coupled Electron Transfer

Proton-coupled electron transfer (PCET) plays an essential role in a wide range of electrocatalytic processes. A vast array of theoretical and computational methods have been developed to study electrochemical PCET. These methods can be used to calculate redox potentials and pK_a values for molecular electrocatalysts, proton-coupled redox potentials and bond dissociation free energies for PCET at metal and semiconductor interfaces, and reorganization energies associated with electrochemical PCET. Periodic density functional theory can also be used to compute PCET activation energies and perform molecular dynamics simulations of electrochemical interfaces. Various approaches for maintaining a constant electrode potential in electronic structure calculations and modeling complex interactions in the electric double layer (EDL) have been developed. Theoretical formulations for both homogeneous and heterogeneous electrochemical PCET spanning the adiabatic, nonadiabatic, and solvent-controlled regimes have been developed and provide analytical expressions for the rate constants and current densities as functions of applied potential. The quantum mechanical treatment of the proton and inclusion of excited vibronic states have been shown to be critical for describing experimental data, such as Tafel slopes and potential-dependent kinetic isotope effects. The calculated rate constants can be used as input to microkinetic models and voltammogram simulations to elucidate complex electrocatalytic processes.

Electrocatalysis in Alkaline Media and Alkaline Membrane-Based Energy Technologies

Hydrogen energy-based electrochemical energy conversion technologies offer the promise of enabling a transition of the global energy landscape from fossil fuels to renewable energy. Here, we present a comprehensive review of the fundamentals of electrocatalysis in alkaline media and applications in alkaline-based energy technologies, particularly alkaline fuel cells and water electrolyzers. Anion exchange (alkaline) membrane fuel cells (AEMFCs) enable the use of nonprecious electrocatalysts for the sluggish oxygen reduction reaction (ORR), relative to proton exchange membrane fuel cells (PEMFCs), which require Pt-based electrocatalysts. However, the hydrogen oxidation reaction (HOR) kinetics is significantly slower in alkaline media than in acidic media. Understanding these phenomena requires applying theoretical and experimental methods to unravel molecular-level thermodynamics and kinetics of hydrogen and oxygen electrocatalysis and, particularly, the proton-coupled electron transfer (PCET) process that takes place in a proton-deficient alkaline media. Extensive electrochemical and spectroscopic studies, on single-crystal Pt and metal oxides, have contributed to the development of activity descriptors, as well as the identification of the nature of active sites, and the rate-determining steps of the HOR and ORR. Among these, the structure and reactivity of interfacial water serve as key potential and pH-dependent kinetic factors that are helping elucidate the origins of the HOR and ORR activity differences in acids and bases. Additionally, deliberately modulating and controlling catalyst-support interactions have provided valuable insights for enhancing catalyst accessibility and durability during operation. The design and synthesis of highly conductive and durable alkaline membranes/ionomers have enabled AEMFCs to reach initial performance metrics equal to or higher than those of PEMFCs. We emphasize the importance of using membrane electrode assemblies (MEAs) to integrate the often separately pursued/optimized electrocatalyst/support and membranes/ionomer components. Operando/in situ methods, at multiscales, and ab initio simulations provide a mechanistic understanding of electron, ion, and mass transport at catalyst/ionomer/membrane interfaces and the necessary guidance to achieve fuel cell operation in air over thousands of hours. We hope that this Review will serve as a roadmap for advancing the scientific understanding of the fundamental factors governing electrochemical energy conversion in alkaline media with the ultimate goal of achieving ultralow Pt or precious-metal-free high-performance and durable alkaline fuel cells and related technologies.

Room temperature bilayer water structures on a rutile TiO2(110) surface: hydrophobic or hydrophilic?

The lack of understanding of the molecular-scale water adsorbed on TiO₂ surfaces under ambient conditions has become a major obstacle for solving the long-time scientific and applications issues, such as the photo-induced wetting phenomenon and designing novel advanced TiO₂-based materials. Here, with the molecular dynamics simulation, we identified an ordered water bilayer structure with a two-dimensional hydrogen bonding network on a rutile TiO₂(110) surface at ambient temperature, corroborated by vibrational sum-frequency generation spectroscopy. The reduced number of hydrogen bonds between the water bilayer and water droplet results in a notable water contact angle (25 ± 5°) of the pristine TiO₂ surface. This surface hydrophobicity can be enhanced by the adsorption of the formate/acetate molecules, and diminishes with dissociated H₂O molecules. Our new physical framework well explained the long-time controversy on the origin of the hydrophobicity/hydrophilicity of the TiO₂ surface, thus help understanding the efficiency of TiO₂ devices in producing electrical energy of solar cells and the photo-oxidation of organic pollutants.

An ordered water bilayer structure was identified on a rutile TiO₂(110) surface at ambient temperature by combining VSFG experiments and MD simulations, which well explained the long-time controversy on the wetting behaviors of the TiO₂ surface.

Effects of Stearyl Alcohol Monolayer on the Structure, Dynamics and Vibrational Sum Frequency Generation Spectroscopy of Interfacial Water

The structure, dynamics and vibrational spectroscopy of water surface covered by a monolayer of stearyl alcohol (STA) are investigated by means of molecular dynamics simulations and vibrational sum frequency generation...

Interface Sensitivity in Electron/Ion Yield X-ray Absorption Spectroscopy: The TiO2-H2O Interface

To understand corrosion, energy storage, (electro)catalysis, etc., obtaining chemical information on the solid-liquid interface is crucial but remains extremely challenging. Here, X-ray absorption spectroscopy (XAS) is used to study the solid-liquid interface between TiO₂ and H₂O. A thin film (6.7 nm) of TiO₂ is deposited on an X-ray-transparent SiN_x window, acting as the working electrode in a three-electrode flow cell. The spectra are collected based on the electron emission resulting from the decay of the X-ray-induced core-hole-excited atoms, which we show is sensitive to the solid-liquid interface within a few nm. The drain currents measured at the working and counter electrodes are identical but of opposite sign. With this method, we found that the water layer next to anatase is spectroscopically similar to ice. This result highlights the potential of electron-yield XAS to obtain chemical and structural information with a high sensitivity for the species at the electrode-electrolyte interface.

Photo-induced hydrophilicity at the ZnO(112̄0) surface: an evolutionary algorithm-aided density functional theory study

Evolutionary algorithm-aided density functional theory calculations were utilized to determine the stable adsorption structures of H₂O at ZnO(112̄0) extensively under different coverages. By decomposing the adsorption energetics, we illustrate that H₂O dissociation to a large extent is actually hampered by the barrier for induced distortion of the ZnO surface, and once the surface becomes less difficult to be distorted it will exhibit higher hydrophilicity or even superhydrophilicity. Specifically, photo-stimulation modelling suggests that the surface Zn-O bonds can be weakened by photo-excitation, and the layer of fully dissociated H₂O can be then facilitated to form. Accordingly, a novel mechanism for photo-induced superhydrophilicity is proposed.

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

Carboxylic acid-modified TiO₂ 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 TiO₂ (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 (O_brH). 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 O_brH's. Simulations of these adsorbed monolayers in water environment show that some protons are released from the surface O_brH's to water resulting in a negatively charged surface with nearby solvated H₃O⁺ 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 TiO₂ surfaces.

Adsorption properties of albumin and fibrinogen on hydrophilic/hydrophobic TiO2 surfaces: A molecular dynamics study

In serval experimental researches, UV-induced hydrophilicity enabled better hemocompatibility in the TiO₂ surface, which was considered to be caused by the removal of the carboxylic acid contamination from the surface. In this paper, we altered the surface wetting property by applying the formate contamination on the rutile (110) surface, and systematically investigated the adsorption properties of albumin and fibrinogen on hydrophilic/hydrophobic TiO₂ surface. Unique contacts were found between the charged residues and the hydrophilic surface, anchoring the protein on the surface. The small size and the heart shape of albumin make it easy to cross the stable water layers near the surface. Besides, albumin has a higher proportion of charged residues, so it can form more unique contacts on the hydrophilic surface. Therefore, the albumin tends to adsorb on the hydrophilic surface. For the hydrophobic surface, the water layers near the surface are weakened, which helps the fibrinogen diffusing to the surface and adjusting its orientation. Although the hydrophobic surface cannot form the unique contacts, the larger size of fibrinogen can provide more residues to form enough ordinary contacts after adjusting, and then achieves stable adsorption. Therefore, fibrinogen tends to adsorb on the hydrophobic surface.

Disentangling Sum-Frequency Generation Spectra of the Water Bending Mode at Charged Aqueous Interfaces

The origin of the sum-frequency generation (SFG) signal of the water bending mode has been controversially debated in the past decade. Unveiling the origin of the signal is essential, because different assignments lead to different views on the molecular structure of interfacial water. Here, we combine collinear heterodyne-detected SFG spectroscopy at the water-charged lipid interfaces with systematic variation of the salt concentration. The results show that the bending mode response is of a dipolar, rather than a quadrupolar, nature and allows us to disentangle the response of water in the Stern and the diffuse layers. While the diffuse layer response is identical for the oppositely charged surfaces, the Stern layer responses reflect interfacial hydrogen bonding. Our findings thus corroborate that the water bending mode signal is a suitable probe for the structure of interfacial water.

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.

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.

Ab initio molecular dynamics of an aqueous electrode interface reveal the electrostatic, structural, and dynamic effects of quantifiable voltage biases on water.

Supercooled liquid-like dynamics in water near a fully hydrated titania surface: Decoupling of rotational and translational diffusion

We report an ab initio molecular dynamics (MD) simulation investigating the effect of a fully hydrated surface of TiO₂ on the water dynamics. It is found that the universal relation between the rotational and translational diffusion characteristics of bulk water is broken in the water layers near the surface with the rotational diffusion demonstrating progressive retardation relative to the translational diffusion when approaching the surface. This kind of rotation-translation decoupling has so far only been observed in the supercooled liquids approaching glass transition, and its observation in water at a normal liquid temperature is of conceptual interest. This finding is also of interest for the application-significant studies of the water interaction with fully hydrated nanoparticles. We note that this is the first observation of rotation-translation decoupling in an ab initio MD simulation of water.

Reactive molecular dynamics simulations of hydration shells surrounding spherical TiO2 nanoparticles: implications for proton-transfer reactions

In many potential applications, nanoparticles are typically in an aqueous medium. This has strong influence on the stability, optical properties and reactivity, in particular for their functionalization. Therefore, the understanding of the chemistry at the interface between the solvent and the nanoparticle is of utmost importance. In this work, we present a comparative ReaxFF reactive molecular dynamics investigation on spherical TiO2 nanoparticles (NSs) of realistic size, with diameters from 2.2 to 4.4 nm, immersed in a large drop of bulk water. After force field validation for its use for a curved anatase TiO2 surface/water interface, we performed several simulations of the TiO2 nanoparticles of increasing size in a water drop. We found that water can be adsorbed jointly in a molecular and dissociative way on the surface. A Langmuir isotherm indicating an adsorption/desorption mechanism of water on the NS is observed. Regarding the dissociative adsorption, atomistic details reveal two different mechanisms, depending on the water concentration around the NS. At low coverage, the first mechanism involves direct dissociation of a single water molecule, whereas, at higher water coverage, the second mechanism is a proton transfer reaction involving two water molecules, also known as Grotthuss-like mechanism. Thermal annealing simulations show that several water molecules remain on the surface in agreement with the experimental reports. The capacity of adsorption is higher for the 2.2 and 3.0 nm NSs than for the 4.4 nm NS. Finally, a comparative investigation with flat surfaces indicates that NSs present a higher water adsorption capacity (undissociated and dissociated) than flat surfaces, which can be rationalized considering that NSs present many more low-coordinated Ti atoms available for water adsorption.

Morphology, Surface Structure and Water Adsorption Properties of TiO2 Nanoparticles: A Comparison of Different Commercial Samples

Water is a molecule always present in the reaction environment in photocatalytic and biomedical applications of TiO₂ and a better understanding of its interaction with the surface of TiO₂ nanoparticles is crucial to develop materials with improved performance. In this contribution, we first studied the nature and the surface structure of the exposed facets of three commercial TiO₂ samples (i.e., TiO₂ P25, SX001, and PC105) by electron microscopy and IR spectroscopy of adsorbed CO. The morphological information was then correlated with the water adsorption properties, investigated at the molecular level, moving from multilayers of adsorbed H₂O to the monolayer, combining medium- and near-IR spectroscopies. Finally, we assessed in a quantitative way the surface hydration state at different water equilibrium pressures by microgravimetric measurements.

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car-Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson's equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. This is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500-1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

Site dependent reactivity of Pt single atoms on anatase TiO2(101) in an aqueous environment

The TiO₂-Pt-water interface is of great relevance in photocatalysis where Pt is widely used as a co-catalyst for enhancing hydrogen evolution in aqueous TiO₂. Using ab initio molecular dynamics, we investigated this interface focusing on Pt single atoms supported on anatase TiO₂(101) in a water environment. Based on recent experiments showing a broad distribution of Pt coordination sites in TiO₂, we examined six distinct single-Pt supported species with different nominal Pt oxidation states, namely: Pt, PtOH, and PtO₂ species adsorbed on the stoichiometric surface; Pt adsorbed at a surface oxygen vacancy (O_v); and Pt substituting a surface Ti cation (Pt_Ti), both without and with an accompanying O_v (Pt_Ti + O_v). As found for the pristine anatase surface, interfacial water remained intact in the presence of a nearly neutral Pt adatom within the time duration of our simulations (∼15 ps). Similarly, no (or only temporary) water dissociation was observed at the Pt_Ti + O_v and PtO₂ interfaces, due to the formation of very stable planar Pt coordination structures that interact only weakly with water. In contrast, water dissociated with OH^- (H⁺) on the Pt atom when this substituted a surface Ti (oxygen) ion as well as on PtOH. The significant proton affinity of Pt atoms at surface oxygen vacancies suggests that negatively charged Pt species are particularly efficient at catalyzing hydrogen evolution in aqueous TiO₂.

Self-interaction error overbinds water clusters but cancels in structural energy differences

We gauge the importance of self-interaction errors in density functional approximations (DFAs) for the case of water clusters. To this end, we used the Fermi-Löwdin orbital self-interaction correction method (FLOSIC) to calculate the binding energy of clusters of up to eight water molecules. Three representative DFAs of the local, generalized gradient, and metageneralized gradient families [i.e., local density approximation (LDA), Perdew-Burke-Ernzerhof (PBE), and strongly constrained and appropriately normed (SCAN)] were used. We find that the overbinding of the water clusters in these approximations is not a density-driven error. We show that, while removing self-interaction error does not alter the energetic ordering of the different water isomers with respect to the uncorrected DFAs, the resulting binding energies are corrected toward accurate reference values from higher-level calculations. In particular, self-interaction-corrected SCAN not only retains the correct energetic ordering for water hexamers but also reduces the mean error in the hexamer binding energies to less than 14 meV/[Formula: see text] from about 42 meV/[Formula: see text] for SCAN. By decomposing the total binding energy into many-body components, we find that large errors in the two-body interaction in SCAN are significantly reduced by self-interaction corrections. Higher-order many-body errors are small in both SCAN and self-interaction-corrected SCAN. These results indicate that orbital-by-orbital removal of self-interaction combined with a proper DFA can lead to improved descriptions of water complexes.

Decoding the molecular water structure at complex interfaces through surface-specific spectroscopy of the water bending mode

The structure of interfacial water determines atmospheric chemistry, wetting properties of materials, and protein folding. The challenge of investigating the properties of specific interfacial water molecules has frequently been confronted using surface-specific sum-frequency generation (SFG) vibrational spectroscopy using the O-H stretch mode. While perfectly suited for the water-air interface, for complex interfaces, a potential complication arises from the contribution of hydroxyl or amine groups of non-water species present at the surface, such as surface hydroxyls on minerals, or O-H and N-H groups contained in proteins. Here, we present a protocol to extract the hydrogen bond strength selectively of interfacial water, through the water bending mode. The bending mode vibrational frequency distribution provides a new avenue for unveiling the hydrogen bonding structure of interfacial water at complex aqueous interfaces. We demonstrate this method for the water-CaF₂ and water-protein interfaces. For the former, we show that this method can indeed single out water O-H groups from surface hydroxyls, and that with increasing pH, the hydrogen-bonded network of interfacial water strengthens. Furthermore, we unveil enhanced hydrogen bonding of water, compared to bulk water, at the interface with human serum albumin proteins, a prototypical bio-interface.

Molecular Structure and Modeling of Water-Air and Ice-Air Interfaces Monitored by Sum-Frequency Generation

From a glass of water to glaciers in Antarctica, water-air and ice-air interfaces are abundant on Earth. Molecular-level structure and dynamics at these interfaces are key for understanding many chemical/physical/atmospheric processes including the slipperiness of ice surfaces, the surface tension of water, and evaporation/sublimation of water. Sum-frequency generation (SFG) spectroscopy is a powerful tool to probe the molecular-level structure of these interfaces because SFG can specifically probe the topmost interfacial water molecules separately from the bulk and is sensitive to molecular conformation. Nevertheless, experimental SFG has several limitations. For example, SFG cannot provide information on the depth of the interface and how the orientation of the molecules varies with distance from the surface. By combining the SFG spectroscopy with simulation techniques, one can directly compare the experimental data with the simulated SFG spectra, allowing us to unveil the molecular-level structure of water-air and ice-air interfaces. Here, we present an overview of the different simulation protocols available for SFG spectra calculations. We systematically compare the SFG spectra computed with different approaches, revealing the advantages and disadvantages of the different methods. Furthermore, we account for the findings through combined SFG experiments and simulations and provide future challenges for SFG experiments and simulations at different aqueous interfaces.

Anharmonic calculations of vibrational spectra for molecular adsorbates: A divide-and-conquer semiclassical molecular dynamics approach

The vibrational spectroscopy of adsorbates is becoming an important investigation tool for catalysis and material science. This paper presents a semiclassical molecular dynamics method able to reproduce the vibrational energy levels of systems composed by molecules adsorbed on solid surfaces. Specifically, we extend our divide-and-conquer semiclassical method for power spectra calculations to gas-surface systems and interface it with plane-wave electronic structure codes. The Born-Oppenheimer classical dynamics underlying the semiclassical calculation is full dimensional, and our method includes not only the motion of the adsorbate but also those of the surface and the bulk. The vibrational spectroscopic peaks related to the adsorbate are accounted together with the most coupled phonon modes to obtain spectra amenable to physical interpretations. We apply the method to the adsorption of CO, NO, and H₂O on the anatase-TiO₂ (101) surface. We compare our semiclassical results with the single-point harmonic estimates and the classical power spectra obtained from the same trajectory employed in the semiclassical calculation. We find that CO and NO anharmonic effects of fundamental vibrations are similarly reproduced by the classical and semiclassical dynamics and that H₂O adsorption is fully and properly described in its overtone and combination band relevant components only by the semiclassical approach.

Structure and sum-frequency generation spectra of water on uncharged Q4 silica surfaces: a molecular dynamics study

The structural characteristics and sum-frequency generation (SFG) spectra of water near neutral Q₄ silica surfaces were investigated by molecular dynamics simulations. The interactions of water molecules with atoms of the solid were described by different potential models, in particular by the CLAYFF [Cygan et al., J. Phys. Chem. B, 2004, 108, 1255] and INTERFACE [Heinz et al. Langmuir, 2013, 29, 1754] force fields. The calculations of the contact angle of water have shown that the silica surface modeled with CLAYFF behaves as macroscopically hydrophilic, in contrast to the surface described with the INTERFACE model. The hydrophilicity of CLAYFF stems from too attractive electrostatic surface-water interactions. Regardless of the surface's affinity for water, the aqueous phase has a layered structure in the direction perpendicular to the surface with density fluctuations decaying within a distance of 10 Å from the surface. The orientational ordering of H₂O molecules was found to be more short-range than the density fluctuations, especially for the hydrophobic surfaces. Modeling the SFG spectra has shown that the spectra of all studied hydrophobic silica-water interfaces are similar and have features in common with the spectrum of the water-vapor interface. The spectra fairly agree with experimental results obtained for the silica-water interface at low pH conditions [Myalitsin et al., J. Phys. Chem. C, 2016, 120, 9357]. The spectral response for the hydrophobic interface was computed to primarily arise from the topmost molecules of the first layer of interfacial water. In contrast, the SFG signal from the hydrophilic silica-water interface is accumulated over a greater distance extending for several water layers due to more long-range perturbation of the structure by the surface.

Free energy of proton transfer at the water-TiO2 interface from ab initio deep potential molecular dynamics

TiO₂ is a widely used photocatalyst in science and technology and its interface with water is important in fields ranging from geochemistry to biomedicine. Yet, it is still unclear whether water adsorbs in molecular or dissociated form on TiO₂ even for the case of well-defined crystalline surfaces. To address this issue, we simulated the TiO₂-water interface using molecular dynamics with an ab initio-based deep neural network potential. Our simulations show a dynamical equilibrium of molecular and dissociative adsorption of water on TiO₂. Water dissociates through a solvent-assisted concerted proton transfer to form a pair of short-lived hydroxyl groups on the TiO₂ surface. Molecular adsorption of water is ΔF = 8.0 ± 0.9 kJ mol^-1 lower in free energy than the dissociative adsorption, giving rise to a 5.6 ± 0.5% equilibrium water dissociation fraction at room temperature. Due to the relevance of surface hydroxyl groups to the surface chemistry of TiO₂, our model might be key to understanding phenomena ranging from surface functionalization to photocatalytic mechanisms.

Structure and reactivity of a water-covered anatase TiO2(001) surface

We systematically studied water adsorption and oxidation on the unreconstructed TiO₂(001) surface using first-principles calculations. Water first adsorbs on the surface in a dissociative state and then in a molecular state, as water coverage increases. The geometric properties of all adsorption structures suggest that the dissociative water molecules can induce stress release of the (001) surface at low coverage, reducing reactivity of the surface and thus leading to molecular adsorption of water on the surface at high coverage. The adsorption energy (or the surface energy) monotonously increases (or decreases) with the increase of the coverage, which further confirms that water, irrespective of its dissociative or molecular state, can improve the stability of the (001) surface and reduce its activity. We deeply investigated the mechanism of the oxygen evolution reaction (OER) on the water-covered (001) surface. A new water-assisted OER pathway is identified on the (001) surface, which includes the sequential transfer of protons from molecular water and surface hydroxyls, and O-O coupling processes. During the OER pathway, the O-O coupling step exhibits the largest thermodynamic energy and highest energy barrier, clarifying that it is the rate-determining step in the whole pathway. Our findings provide new insights into the strong dependence of water adsorption modes on coverage for the anatase TiO₂(001) surface and may explain the high oxidation activity of the TiO₂(001) surface in aqueous environments typical of TiO₂ photocatalysis.

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 pK_a) 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 pK_a 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 pK_a is only 1.0 pK_a 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.

On the stability of hydroxyl groups on substituted titania

The present study reports the stability of hydroxyl groups involving the surface coordinated oxygens of Pd,C,N-doped, and Pd/C and Pd/N-codoped anatase TiO₂, probed using DFT calculations.

Hydrogen bonding and molecular orientations across thin water films on sapphire

HYPOTHESIS

Water vapor binding to metal oxide surfaces produces thin water films with properties controlled by interactions with surface hydroxo sites. Hydrogen bonding populations vary across films and induce different molecular orientations than at the surface of liquid water. Identifying these differences can open possibilities for tailoring film-mediated catalytic reactions by choice of the supporting metal oxide substrate.

EXPERIMENTS

The (0001) face of a single sapphire (α-Al₂O₃) sample exposed to water vapor and the surface of liquid water were probed by polarization dependent Sum Frequency Generation-Vibration Spectroscopy (SFG-VS). Molecular dynamics (MD) provided insight into the hydrogen bond populations and molecular orientations across films and liquid water.

FINDINGS

SFG-VS revealed a submonolayer film on sapphire exposed to 43% relative humidity (R.H.), and a multilayer film at 78% R.H. Polarization dependent SFG-VS spectra showed that median tilt angles of free OH bonds on the top of films are at ∼43° from the normal of the (0001) face but at 38° on neat liquid water. These values align with MD simulations, which also show that up to 36% of all OH bonds on films are free. This offers new means for understanding how interfacial reactions on sapphire-supported water films could contrast with those involving liquid water.