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.

Transition from Sequential to Concerted Proton-Coupled Electron Transfer of Water Oxidation on Semiconductor Photoanodes

Accelerating proton transfer has been demonstrated as key to boosting water oxidation on semiconductor photoanodes. Herein, we study proton-coupled electron transfer (PCET) of water oxidation on five typical photoanodes [i.e., α-Fe2O3, BiVO4, TiO2, plasmonic Au/TiO2, and nickel-iron oxyhydroxide (Ni1-xFexOOH)-modified silicon (Si)] by combining the rate law analysis of H2O molecules with the H/D kinetic isotope effect (KIE) and operando spectroscopic studies. An unexpected and universal half-order kinetics is observed for the rate law analysis of H2O, referring to a sequential proton-electron transfer pathway, which is the rate-limiting factor that causes the sluggish water oxidation performance. Surface modification of the Ni1-xFexOOH electrocatalyst is observed to break this limitation and exhibits a normal first-order kinetics accompanied by much enhanced H/D KIE values, facilitating the turnover frequency of water oxidation by 1 order of magnitude. It is the first time that Ni1-xFexOOH is found to be a PCET modulator. The rate law analysis illustrates an effective strategy for modulating PCET kinetics of water oxidation on semiconductor surfaces.

Precisely Constructed Metal Sulfides with Localized Single-Atom Rhodium for Photocatalytic C-H Activation and Direct Methanol Coupling to Ethylene Glycol

Although there are many studies on photocatalytic environmental remediation, hydrogen evolution, and chemical transformations, less success has been achieved for the synthesis of industrially important and largely demanded bulk chemicals using semiconductor photocatalysis, which holds great potential to drive unique chemical reactions that are difficult to implement by the conventional heterogeneous catalysis. The performance of semiconductors used for photochemical synthesis is, however, usually unsatisfactory due to limited efficiencies in light harvesting, charge-carrier separation, and surface reactions. The precise construction of heterogeneous photocatalysts to facilitate these processes is an attractive but challenging goal. Here, single-atom rhodium-doped metal sulfide nanorods composed of alternately stacked wurtzite/zinc-blende segments are successfully designed and fabricated, which demonstrate record-breaking efficiencies for visible light-driven preferential activation of C-H bond in methanol to form ethylene glycol (EG), a key bulk chemical used for the production of polyethylene terephthalate (PET) polymer. The wurtzite/zinc-blende heterojunctions lined regularly in one dimension accelerate the charge-carrier separation and migration. Single-atom rhodium selectively deposited onto the wurtzite segment with photogenerated holes accumulated facilitates methanol adsorption and C-H activation. The present work paves the way to harnessing photocatalysis for bulk chemical synthesis with structure-defined semiconductors.

Accelerated interfacial proton transfer for promoting the electrocatalytic activity

Interfacial pH is critical to electrocatalytic reactions involving proton-coupled electron transfer (PCET) processes, and maintaining an optimal interfacial pH at the electrochemical interface is required to achieve high activity. However, the interfacial pH varies inevitably during the electrochemical reaction owing to slow proton transfer at the interfacial layer, even in buffer solutions. It is therefore necessary to find an effective and general way to promote proton transfer for regulating the interfacial pH. In this study, we propose that promoting proton transfer at the interfacial layer can be used to regulate the interfacial pH in order to enhance electrocatalytic activity. By adsorbing a bifunctional 4-mercaptopyridine (4MPy) molecule onto the catalyst surface via its thiol group, the pyridyl group can be tethered on the electrochemical interface. The pyridyl group acts as both a good proton acceptor and donor for promoting proton transfer at the interfacial layer. Furthermore, the pK_a of 4MPy can be modulated with the applied potentials to accommodate the large variation of interfacial pH under different current densities. By in situ electrochemical surface-enhanced Raman spectroscopy (in situ EC-SERS), we quantitatively demonstrate that proton transfer at the interfacial layer of the Pt catalyst coated with 4MPy (Pt@4MPy) remains ideally thermoneutral during the H⁺ releasing electrocatalytic oxidation reaction of formic acid (FAOR) at high current densities. Thus, the interfacial pH is controlled effectively. In this way, the FAOR apparent current measured from Pt@4MPy is twice that measured from a pristine Pt catalyst. This work establishes a general strategy for regulating interfacial pH to enhance the electrocatalytic activities.

Adsorbing 4MPy on Pt surface promotes proton transfer at the interfacial layer, maintaining an optimal interfacial pH and promotes electrocatalytic reactions involving proton-coupled electron transfer (PCET) processes.

Proton-Coupled Defects Impact O-H Bond Dissociation Free Energies on Metal Oxide Surfaces

Proton-coupled electron transfer (PCET) reactions on metal oxides require coupling between proton transfer at the solid-liquid interface and electron transfer involving defects at or near the band edge. Herein, hybrid functional periodic density functional theory is used to elucidate the impact of proton-coupled defects on the bond dissociation free energies (BDFEs) of O-H bonds on anatase TiO₂ surfaces. These O-H BDFEs are directly related to interfacial PCET thermochemistry. Comparison between geometrically similar O-H bonds associated with different defect types, namely conduction d-band electrons or valence p-band holes, reveals that the BDFEs differ by ∼81 kcal/mol (3.50 eV), comparable to the wide TiO₂ band gap. These differences are shown to be determined primarily by differences in electron transfer driving forces, which are analyzed by using band energies and inner-sphere reorganization energies within a Marcus theory framework. These fundamental insights about the impact of proton-coupled defects on PCET thermochemistry at semiconductor surfaces have broad implications for electrocatalysis.

Interplay of physically different properties leading to challenges in separating lanthanide cations - an ab initio molecular dynamics and experimental study

Lanthanide elements have well-documented similarities in their chemical behavior, which make the valuable trivalent lanthanide cations (Ln3+) particularly difficult to separate from each other in water. In this work, we apply ab initio molecular dynamics simulations to compare the free energies (ΔGads) associated with the adsorption of lanthanide cations to silica surfaces at a pH condition where SiO- groups are present. The predicted ΔGads for lutetium (Lu3+) and europium (Eu3+) are similar within statistical uncertainties; this is in qualitative agreement with our batch adsorption measurements on silica. This finding is remarkable because the two cations exhibit hydration free energies (ΔGhyd) that differ by >2 eV, different hydration numbers, and different hydrolysis behavior far from silica surfaces. We observe that the similarity in Lu3+ and Eu3+ ΔGads is the result of a delicate cancellation between the difference in Eu3+ and Lu3+ hydration (ΔGhyd), and their difference in binding energies to silica. We propose that disrupting this cancellation at the two end points, either for adsorbed or completely desorbed lanthanides (e.g., via nanoconfinment or mixed solvents), will lead to effective Ln3+ separation.

Electrocatalysis for CO2 conversion: from fundamentals to value-added products

This timely and comprehensive review mainly summarizes advances in heterogeneous electroreduction of CO₂: from fundamentals to value-added products.

Mechanisms for Tuning Engineered Nanomaterials to Enhance Radiation Therapy of Cancer

Engineered nanomaterials that produce reactive oxygen species on exposure to X- and gamma-rays used in radiation therapy offer promise of novel cancer treatment strategies. Similar to photodynamic therapy but suitable for large and deep tumors, this new approach where nanomaterials acting as sensitizing agents are combined with clinical radiation can be effective at well-tolerated low radiation doses. Suitably engineered nanomaterials can enhance cancer radiotherapy by increasing the tumor selectivity and decreasing side effects. Additionally, the nanomaterial platform offers therapeutically valuable functionalities, including molecular targeting, drug/gene delivery, and adaptive responses to trigger drug release. The potential of such nanomaterials to be combined with radiotherapy is widely recognized. In order for further breakthroughs to be made, and to facilitate clinical translation, the applicable principles and fundamentals should be articulated. This review focuses on mechanisms underpinning rational nanomaterial design to enhance radiation therapy, the understanding of which will enable novel ways to optimize its therapeutic efficacy. A roadmap for designing nanomaterials with optimized anticancer performance is also shown and the potential clinical significance and future translation are discussed.

Computing Surface Acidity Constants of Proton Hopping Groups from Density Functional Theory-Based Molecular Dynamics: Application to the SnO2(110)/H2O Interface

Proton transfer at metal oxide/water interfaces plays an important role in electrochemistry, geochemistry, and environmental science. The key thermodynamic quantity to characterize this process is the surface acidity constant. An ab initio method that combines density functional theory-based molecular dynamics (DFTMD) and free energy perturbation theory has been established for computing surface acidity constants. However, it involves a reversible proton insertion procedure in which frequent proton hopping, e.g., for strong bases and some oxide surfaces (e.g., SnO₂), can cause instability issues in electronic structure calculation. In the original implementation, harmonic restraining potentials are imposed on all O-H bonds (denoted by the V_rH scheme) to prevent proton hopping and thus may not be applicable for systems involving spontaneous proton hopping. In this work, we introduce an improved restraining scheme with a repulsive potential V_rep to compute the surface acidities of systems in which proton hopping is spontaneous and fast. In this V_rep scheme, a Buckingham-type repulsive potential V_rep is applied between the deprotonation site and all other protons in DFTMD simulations. We first verify the V_rep scheme by calculating the pK_a values of H₂O and aqueous HS^- solution (i.e., strong conjugate bases) and then apply it to the SnO₂(110)/H₂O interface. It is found that the V_rep scheme leads to a prediction of the point of zero charge (PZC) of 4.6, which agrees well with experiment. The intrinsic individual pK_a values of the terminal five-coordinated Sn site (Sn_5cOH₂) and bridge oxygen site (Sn₂O_brH⁺) are 4.4 and 4.7, respectively, both being almost the same as the PZC. The similarity of the two pK_a values indicates that dissociation of terminal water has almost zero free energy at this proton hopping interface (i.e., partial water dissociation), as expected from the acid-base equilibrium on SnO₂.

Mass-zero constrained molecular dynamics for electrode charges in simulations of electrochemical systems

Classical molecular dynamics simulations have recently become a standard tool for the study of electrochemical systems. State-of-the-art approaches represent the electrodes as perfect conductors, modeling their responses to the charge distribution of electrolytes via the so-called fluctuating charge model. These fluctuating charges are additional degrees of freedom that, in a Born-Oppenheimer spirit, adapt instantaneously to changes in the environment to keep each electrode at a constant potential. Here, we show that this model can be treated in the framework of constrained molecular dynamics, leading to a symplectic and time-reversible algorithm for the evolution of all the degrees of freedom of the system. The computational cost and the accuracy of the new method are similar to current alternative implementations of the model. The advantage lies in the accuracy and long term stability guaranteed by the formal properties of the algorithm and in the possibility to systematically introduce additional kinematic conditions of arbitrary number and form. We illustrate the performance of the constrained dynamics approach by enforcing the electroneutrality of the electrodes in a simple capacitor consisting of two graphite electrodes separated by a slab of liquid water.

Frontiers in molecular simulation of solvated ions, molecules and interfaces

This themed collection is a collection of articles on frontiers in molecular simulation of solvated ions, molecules and interfaces.

Theoretical study of kinetics of proton coupled electron transfer in photocatalysis

Photocatalysis induced by sunlight is one of the most promising approaches to environmental protection, solar energy conversion, and sustainable production of fuels. The computational modeling of photocatalysis is a rapidly expanding field that requires to adapt and to further develop the available theoretical tools. The coupled transfer of protons and electrons is an important reaction during photocatalysis. In this work, we present the first step of our methodology development in which we apply the existing kinetic theory of such coupled transfer to a model system, namely, methanol photodissociation on the rutile TiO₂(110) surface, with the help of high-level first-principles calculations. Moreover, we adapt the Stuchebrukhov-Hammes-Schiffer kinetic theory, where we use the Georgievskii-Stuchebrukhova vibronic coupling to calculate the rate constant of the proton coupled electron transfer reaction for a particular pathway. In particular, we propose a modified expression to calculate the rate constant, which enforces the near-resonance condition for the vibrational wave function during proton tunneling.

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.

Computational Ag/AgCl Reference Electrode from Density Functional Theory-Based Molecular Dynamics

We have developed a scheme to compute the standard potential of the Ag/AgCl reference electrode using density functional theory-based molecular dynamics, similar to the computational standard hydrogen electrode (SHE) developed by Cheng, Sulpizi, and Sprik [J. Chem. Phys. 2009, 131, 154504], with which our new computational reference electrode was compared. We have obtained a similar value of the potential of the Ag/AgCl electrode versus SHE to the experiment. The newly developed computational reference electrode will be extended to nonaqueous solvents in the future, where it will be used to predict standard equilibrium potentials to be compared with experimental data.

Analytic theory of finite-size effects in supercell modeling of charged interfaces

The Ewald3D sum with the tinfoil boundary condition (e3dtf) evaluates the electrostatic energy of a finite unit cell inside an infinitely periodic supercell. Although it has been used as a de facto standard treatment of electrostatics for simulations of extended polar or charged interfaces, the finite-size effect on simulated properties has yet to be fully understood. There is, however, an intuitive way to quantify the average effect arising from the difference between the e3dtf and Coulomb potentials on the response of mobile charges to contact surfaces with fixed charges and/or to an applied external electric field. Although any charged interface formed by mobile countercharges that compensate the fixed charges fluctuates upon a change in the acting electric field, the distance between a pair of oppositely charged interfaces is found to be nearly stationary, which allows an analytic finite-size correction to the amount of countercharges. Application of the theory to solvated electric double layers (insulator/electrolyte interfaces) predicts that the state of complete charge compensation is invariant with respect to solvent permittivities, which is confirmed by a proper analysis of simulation data in the literature.

Diffusion and reaction pathways of water near fully hydrated TiO2 surfaces from ab initio molecular dynamics

Ab initio molecular dynamics simulations are reported for water-embedded TiO₂ surfaces to determine the diffusive and reactive behavior at full hydration. A three-domain model is developed for six surfaces [rutile (110), (100), and (001), and anatase (101), (100), and (001)] which describes waters as "hard" (irreversibly bound to the surface), "soft" (with reduced mobility but orientation freedom near the surface), or "bulk." The model explains previous experimental data and provides a detailed picture of water diffusion near TiO₂ surfaces. Water reactivity is analyzed with a graph-theoretic approach that reveals a number of reaction pathways on TiO₂ which occur at full hydration, in addition to direct water splitting. Hydronium (H₃O⁺) is identified to be a key intermediate state, which facilitates water dissociation by proton hopping between intact and dissociated waters near the surfaces. These discoveries significantly improve the understanding of nanoscale water dynamics and reactivity at TiO₂ interfaces under ambient conditions.

Visible light-driven C-H activation and C-C coupling of methanol into ethylene glycol

The development of new methods for the direct transformation of methanol into two or multi-carbon compounds via controlled carbon–carbon coupling is a highly attractive but challenging goal. Here, we report the first visible-light-driven dehydrogenative coupling of methanol into ethylene glycol, an important chemical currently produced from petroleum. Ethylene glycol is formed with 90% selectivity and high efficiency, together with hydrogen over a molybdenum disulfide nanofoam-modified cadmium sulfide nanorod catalyst. Mechanistic studies reveal a preferential activation of C−H bond instead of O−H bond in methanol by photoexcited holes on CdS via a concerted proton–electron transfer mechanism, forming a hydroxymethyl radical (⋅CH₂OH) that can readily desorb from catalyst surfaces for subsequent coupling. This work not only offers an alternative nonpetroleum route for the synthesis of EG but also presents a unique visible-light-driven catalytic C−H activation with the hydroxyl group in the same molecule keeping intact.

Direct transformation of methanol into two- or multi-carbon compounds is extremely attractive but remains a challenge. Here, the authors report an efficient photocatalytic route to the transformation of methanol into ethylene glycol and hydrogen over a molybdenum disulfide nanofoam-modified cadmium sulfide nanorod catalyst.

Determining Potentials of Zero Charge of Metal Electrodes versus the Standard Hydrogen Electrode from Density-Functional-Theory-Based Molecular Dynamics

We develop a computationally efficient scheme to determine the potentials of zero charge (PZC) of metal-water interfaces with respect to the standard hydrogen electrode. We calculate the PZC of Pt(111), Au(111), Pd(111) and Ag(111) at a good accuracy using this scheme. Moreover, we find that the interface dipole potentials are almost entirely caused by charge transfer from water to the surfaces, the magnitude of which depends on the bonding strength between water and the metals, while water orientation hardly contributes at the PZC conditions.

Structure of a model TiO2 photocatalytic interface

The interaction of water with TiO₂ is crucial to many of its practical applications, including photocatalytic water splitting. Following the first demonstration of this phenomenon 40 years ago there have been numerous studies of the rutile single-crystal TiO₂(110) interface with water. This has provided an atomic-level understanding of the water-TiO₂ interaction. However, nearly all of the previous studies of water/TiO₂ interfaces involve water in the vapour phase. Here, we explore the interfacial structure between liquid water and a rutile TiO₂(110) surface pre-characterized at the atomic level. Scanning tunnelling microscopy and surface X-ray diffraction are used to determine the structure, which is comprised of an ordered array of hydroxyl molecules with molecular water in the second layer. Static and dynamic density functional theory calculations suggest that a possible mechanism for formation of the hydroxyl overlayer involves the mixed adsorption of O₂ and H₂O on a partially defected surface. The quantitative structural properties derived here provide a basis with which to explore the atomistic properties and hence mechanisms involved in TiO₂ photocatalysis.

Hydrogen evolution based on the electrons/protons stored on amorphous TiO2

The hydrogen evolution reaction (HER) using recyclable mediator is being actively pursued as a route for solar energy conversion.

First-Principles Free-Energy Barriers for Photoelectrochemical Surface Reactions: Proton Abstraction at TiO_{2}(110)

We explicitly calculate the free-energy barrier for the initial proton abstraction in the water splitting reaction at rutile TiO_{2}(110) through ab initio molecular dynamics. Combining solid-state embedding, an energy based reaction coordinate and state-of-the-art free-energy reconstruction techniques renders the calculation tractable at the hybrid density-functional theory level. The obtained free-energy barrier of approximately 0.2 eV, depending slightly on the orientation of the first acceptor water molecule, suggests a hindered reaction on the pristine rutile surface.

Modeling and Simulations in Photoelectrochemical Water Oxidation: From Single Level to Multiscale Modeling

This review summarizes recent developments, challenges, and strategies in the field of modeling and simulations of photoelectrochemical (PEC) water oxidation. We focus on water splitting by metal-oxide semiconductors and discuss topics such as theoretical calculations of light absorption, band gap/band edge, charge transport, and electrochemical reactions at the electrode-electrolyte interface. In particular, we review the mechanisms of the oxygen evolution reaction, strategies to lower overpotential, and computational methods applied to PEC systems with particular focus on multiscale modeling. The current challenges in modeling PEC interfaces and their processes are summarized. At the end, we propose a new multiscale modeling approach to simulate the PEC interface under conditions most similar to those of experiments. This approach will contribute to identifying the limitations at PEC interfaces. Its generic nature allows its application to a number of electrochemical systems.

Density Functional Theory Calculation of the Band Alignment of (101̅0) In(x)Ga(1-x)N/Water Interfaces

Conduction band edge (CBE) and valence band edge (VBE) positions of InxGa1-xN photoelectrodes were computed using density functional theory methods. The band edges of fully solvated GaN and InN model systems were aligned with respect to the standard hydrogen electrode using a molecular dynamics hydrogen electrode scheme applied earlier to TiO2/water interfaces. Similar to the findings for TiO2, we found that the Purdew-Burke-Ernzerhof (PBE) functional gives a VBE potential which is too negative by 1 V. This cathodic bias is largely corrected by application of the Heyd-Scuseria-Ernzerhof (HSE06) hybrid functional containing a fraction of Hartree-Fock exchange. The effect of a change of composition was investigated using simplified model systems consisting of vacuum slabs covered on both sides by one monolayer of H2O. The CBE was found to vary linearly with In content. The VBE, in comparison, is much less sensitive to composition. The data show that the band edges straddle the hydrogen and oxygen evolution potentials for In fractions less than 47%. The band gap was found to exceed 2 eV for an In fraction less than 54%.

Calculation of Electrochemical Energy Levels in Water Using the Random Phase Approximation and a Double Hybrid Functional

Understanding charge transfer at electrochemical interfaces requires consistent treatment of electronic energy levels in solids and in water at the same level of the electronic structure theory. Using density-functional-theory-based molecular dynamics and thermodynamic integration, the free energy levels of six redox couples in water are calculated at the level of the random phase approximation and a double hybrid density functional. The redox levels, together with the water band positions, are aligned against a computational standard hydrogen electrode, allowing for critical analysis of errors compared to the experiment. It is encouraging that both methods offer a good description of the electronic structures of the solutes and water, showing promise for a full treatment of electrochemical interfaces.

A computational study of the competing reaction mechanisms of the photo-catalytic reduction of CO2 on anatase(101)

Three reaction pathways for the photocatalytic reduction of carbon dioxide to methane are investigated with density functional theory calculations.