Interfacial Effects on the Band Edges of Functionalized Si Surfaces in Liquid Water

Water effect on the band edges of anatase TiO2 surfaces: A theoretical study on charge migration across surface heterojunctions and facet-dependent photoactivity

The charge migration mechanism across the surface heterojunction constructed on an anatase TiO2 nanocrystal is still under debate. To solve this longstanding question, we present a systematic study of the band edges (vs. standard hydrogen electrode, SHE) of aqueous TiO2 interfaces with anatase (101), (100) and (001) surfaces, using a combination of density functional theory-based molecular dynamics (DFTMD) and efficient computational SHE (cSHE) methods. Our calculations show that the conduction band minimum (CBM) of the (101) surface is lower than that of (001) and (100) surfaces, which is thermodynamically favorable for electrons migrating to the (101) surface through the surface heterojunction, while the hole preferentially accumulates on the (100) surface due to its highest valence band minimum (VBM). In addition, we qualitatively explore the facet-dependent photocatalytic activity of anatase TiO2. Due to the possession of both the beneficial atomic structure (with 100% undercoordinated Ti5c atoms at the surface) and electronic structure (more strongly oxidizing holes in the VBM and efficient electron-hole spatial separation separation), the (001) surface exhibits the most efficient photocatalytic performance for water oxidation. Furthermore, it is confirmed that the use of simplified theoretical models neglecting the detailed atomic structures of water at the aqueous interface is inadequate to predict the band alignment of semiconductors relative to water redox potentials, so that it may result in substantial errors in evaluating the photocatalytic performance of materials to be used for water splitting.

In situ electrochemical Raman spectroscopy and ab initio molecular dynamics study of interfacial water on a single-crystal surface

The dynamics and chemistry of interfacial water are essential components of electrocatalysis because the decomposition and formation of water molecules could dictate the protonation and deprotonation processes on the catalyst surface. However, it is notoriously difficult to probe interfacial water owing to its location between two condensed phases, as well as the presence of external bias potentials and electrochemically induced reaction intermediates. An atomically flat single-crystal surface could offer an attractive platform to resolve the internal structure of interfacial water if advanced characterization tools are developed. To this end, here we report a protocol based on the combination of in situ Raman spectroscopy and ab initio molecular dynamics (AIMD) simulations to unravel the directional molecular features of interfacial water. We present the procedures to prepare single-crystal electrodes, construct a Raman enhancement mode with shell-isolated nanoparticle, remove impurities, eliminate the perturbation from bulk water and dislodge the hydrogen bubbles during in situ electrochemical Raman experiments. The combination of the spectroscopic measurements with AIMD simulation results provides a roadmap to decipher the potential-dependent molecular orientation of water at the interface. We have prepared a detailed guideline for the application of combined in situ Raman and AIMD techniques; this procedure may take a few minutes to several days to generate results and is applicable to a variety of disciplines ranging from surface science to energy storage to biology.

GPU Acceleration of Large-Scale Full-Frequency GW Calculations

Many-body perturbation theory is a powerful method to simulate electronic excitations in molecules and materials starting from the output of density functional theory calculations. By implementing the theory efficiently so as to run at scale on the latest leadership high-performance computing systems it is possible to extend the scope of GW calculations. We present a GPU acceleration study of the full-frequency GW method as implemented in the WEST code. Excellent performance is achieved through the use of (i) optimized GPU libraries, e.g., cuFFT and cuBLAS, (ii) a hierarchical parallelization strategy that minimizes CPU-CPU, CPU-GPU, and GPU-GPU data transfer operations, (iii) nonblocking MPI communications that overlap with GPU computations, and (iv) mixed precision in selected portions of the code. A series of performance benchmarks has been carried out on leadership high-performance computing systems, showing a substantial speedup of the GPU-accelerated version of WEST with respect to its CPU version. Good strong and weak scaling is demonstrated using up to 25 920 GPUs. Finally, we showcase the capability of the GPU version of WEST for large-scale, full-frequency GW calculations of realistic systems, e.g., a nanostructure, an interface, and a defect, comprising up to 10 368 valence electrons.

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.

Crystallographic Effects of GaN Nanostructures in Photoelectrochemical Reaction

Tuning the surface structure of the photoelectrode provides one of the most effective ways to address the critical challenges in artificial photosynthesis, such as efficiency, stability, and product selectivity, for which gallium nitride (GaN) nanowires have shown great promise. In the GaN wurtzite crystal structure, polar, semipolar, and nonpolar planes coexist and exhibit very different structural, electronic, and chemical properties. Here, through a comprehensive study of the photoelectrochemical performance of GaN photocathodes in the form of films and nanowires with controlled surface polarities we show that significant photoelectrochemical activity can be observed when the nonpolar surfaces are exposed in the electrolyte, whereas little or no activity is measured from the GaN polar c-plane surfaces. The atomic origin of this fundamental difference is further revealed through density functional theory calculations. This study provides guideline on crystal facet engineering of metal-nitride photo(electro)catalysts for a broad range of artificial photosynthesis chemical reactions.

Boosting the photocatalytic H2 evolution activity of type-II g-GaN/Sc2CO2 van der Waals heterostructure using applied biaxial strain and external electric field

Two-dimensional (2D) van der Waals (vdW) heterostructures are a new class of materials with highly tunable bandgap transition type, bandgap energy and band alignment. Herein, we have designed a novel 2D g-GaN/Sc2CO2 heterostructure as a potential solar-driven photocatalyst for the water splitting process and investigate its catalytic stability, interfacial interactions, and optical and electronic properties, as well as the effects of applying an electric field and biaxial strain using first-principles calculation. The calculated lattice mismatch and binding energy showed that g-GaN and Sc2CO2 are in contact and may form a stable vdW heterostructure. Ab initio molecular dynamics and phonon dispersion simulations show thermal and dynamic stability. g-GaN/Sc2CO2 has an indirect bandgap energy with appropriate type-II band alignment relative to the water redox potentials. Meanwhile, the interfacial charge transfer from g-GaN to Sc2CO2 can effectively separate electron-hole pairs. Moreover, a potential drop of 3.78 eV is observed across the interface, inducing a built-in electric field pointing from g-GaN to Sc2CO2. The heterostructure shows improved visible-light optical absorption compared to the isolated g-GaN and Sc2CO2 monolayers. Our study demonstrates that tunable electronic and structural properties can be realised in the g-GaN/Sc2CO2 heterostructure by varying the electric field and biaxial strain. In particular, the compressive strain and negative electric field are more effective for promoting hydrogen production performance. Since it is challenging to tune the electric field and biaxial strain experimentally, our research provides strategies to boost the performance of MXene-based heterojunction photocatalysts in solar harvesting and optoelectronic devices.

Electronic structure of strain-tunable Janus WSSe–ZnO heterostructures from first-principles

The electronic structure of semiconducting 2D materials such as monolayer transition metal dichalcogenides (TMDs) are known to be tunable via environment and external fields, and van der Waals (vdW) heterostructures consisting of stacks of distinct types of 2D materials offer the possibility to further tune and optimize the electronic properties of 2D materials. In this work, we use density functional theory (DFT) calculations to calculate the structure and electronic properties of a vdW heterostructure of Janus monolayer WSSe with monolayer ZnO, both of which possess out of plane dipole moments. The effects of alignment, biaxial and uniaxial strain, orientation, and electric field on dipole moments and band edge energies of this heterostructure are calculated and examined. We find that the out of plane dipole moment of the ZnO monolayer is highly sensitive to strain, leading to the broad tunability of the heterostructure band edge energies over a range of experimentally-relevant strains. The use of strain-tunable 2D materials to control band offsets and alignment is a general strategy applicable to other vdW heterostructures, one that may be advantageous in the context of clean energy applications, including photocatalytic applications, and beyond.

Using strain engineering to optimize novel heterostructure materials to produce hydrogen from water.

Probing the electronic properties of the electrified silicon/water interface by combining simulations and experiments

Silicon (Si) is broadly used in electrochemical and photoelectrochemical devices, where the capacitive and Faradaic reactions at the Si/water interfaces are critical for signal transduction or noise generation. However, probing the electrified Si/water interface at the microscopic level remains a challenging task. Here we focus on hydrogenated Si surfaces in contact with water, relevant to transient electronics and photoelectrochemical modulation of biological cells and tissues. We show that by carrying out first-principles molecular dynamics simulations of the Si(100)/water interface in the presence of an electric field we can realistically correlate the computed flat-band potential and tunneling current images at the interface with experimentally measured capacitive and Faradaic currents. Specifically, we validate our simulations in the presence of bias by performing pulsed chronoamperometry measurements on Si wafers in solution. Consistent with prior experiments, our measurements and simulations indicate the presence of voltage-dependent capacitive currents at the interface. We also find that Faradaic currents are weakly dependent on the applied bias, which we relate to surface defects present in newly prepared samples.

Electrification of water interface

The surface charge of a water interface determines many fundamental processes in physical chemistry and interface science, and it has been intensively studied for over a hundred years. We summarize experimental methods to characterize the surface charge densities developed so far: electrokinetics, double-layer force measurements, potentiometric titration, surface-sensitive nonlinear spectroscopy, and surface-sensitive mass spectrometry. Then, we elucidate physical ion adsorption and chemical electrification as examples of electrification mechanisms. In the end, novel effects on surface electrification are discussed in detail. We believe that this clear overview of state of the art in a charged water interface will surely help the fundamental progress of physics and chemistry at interfaces in the future.

Review on Interfacial Bonding Mechanism of Functional Polymer Coating on Glass in Atomistic Modeling Perspective

Atomistic modeling methods are successfully applied to understand interfacial interaction in nanoscale size and analyze adhesion mechanism in the organic-inorganic interface. In this paper, we review recent representative atomistic simulation works, focusing on the interfacial bonding, adhesion strength, and failure behavior between polymer film and silicate glass. The simulation works are described under two categories, namely non-bonded and bonded interaction. In the works for non-bonded interaction, three main interactions, namely van der Waals interaction, polar interaction, and hydrogen bonds, are investigated, and the contributions to interfacial adhesion energy are analyzed. It is revealed that the most dominant interaction for adhesion is hydrogen bonding, but flexibility of the polymer film and modes of adhesion measurement test do affect adhesion and failure behavior. In the case of bonded interactions, the mechanism of covalent silane bond formation through condensation and hydrolysis process is reviewed, and surface reactivity, molecular density, and adhesion properties are calculated with an example of silane functionalized polymer. Besides interfacial interactions, effects of external conditions, such as surface morphology of the glass substrate and relative humidity on the adhesion and failure behavior, are presented, and modeling techniques developed for building interfacial system and calculating adhesion strengths are briefly introduced.

A Review of Recent Developments in Molecular Dynamics Simulations of the Photoelectrochemical Water Splitting Process

In this review, we provide a short overview of the Molecular Dynamics (MD) method and how it can be used to model the water splitting process in photoelectrochemical hydrogen production. We cover classical non-reactive and reactive MD techniques as well as multiscale extensions combining classical MD with quantum chemical and continuum methods. Selected examples of MD investigations of various aqueous semiconductor interfaces with a special focus on TiO2 are discussed. Finally, we identify gaps in the current state-of-the-art where further developments will be needed for better utilization of MD techniques in the field of water splitting.

Origin of Asymmetric Electric Double Layers at Electrified Oxide/Electrolyte Interfaces

The structure of electric double layers (EDLs) dictates the chemistry and the physics of electrified interfaces, and the differential capacitance is the key property for characterizing EDLs. Here we develop a theoretical model for computing the differential Helmholtz capacitance C_H of oxide-electrolyte interfaces using density functional theory-based finite-field molecular dynamics simulations. It is found that the dipole of interfacial adsorbed groups (i.e., water molecule, hydroxyl ion, and proton) at the electrified SnO₂(110)/NaCl interfaces significantly modulates the double layer potential which leads to the asymmetric distribution of C_H. We also find that the dissociative water adsorption prefers the inner sphere binding of counterions, which in turn leads to a higher Helmholtz capacitance, compared with that of the nondissociative case at the interface. This work provides a molecular interpretation of asymmetric EDLs seen experimentally in a range of metal oxides/hydroxides.

Electronic, optical, and thermoelectric properties of Janus In-based monochalcogenides

Inspired by the successfully experimental synthesis of Janus structures recently, we systematically study the electronic, optical, and electronic transport properties of Janus monolayers In₂XY(X/Y= S, Se, Te withX≠Y) in the presence of a biaxial strain and electric field using density functional theory. Monolayers In₂XYare dynamically and thermally stable at room temperature. At equilibrium, both In₂STe and In₂SeTe are direct semiconductors while In₂SSe exhibits an indirect semiconducting behavior. The strain significantly alters the electronic structure of In₂XYand their photocatalytic activity. Besides, the indirect-direct gap transitions can be found due to applied strain. The effect of the electric field on optical properties of In₂XYis negligible. Meanwhile, the optical absorbance intensity of the Janus In₂XYmonolayers is remarkably increased by compressive strain. Also, In₂XYmonolayers exhibit very low lattice thermal conductivities resulting in a high figure of meritZT, which makes them potential candidates for room-temperature thermoelectric materials.

Gold-Decorated Silicon Nanowire Photocatalysts for Intracellular Production of Hydrogen Peroxide

Hydrogen peroxide (H₂O₂) plays diverse biological roles, and its effects in part depend on its spatiotemporal presence, in both intra- and extracellular contexts. A full understanding of the physiological effects of H₂O₂ in both healthy and disease states is hampered by a lack of tools to controllably produce H₂O₂. Here, we address this issue by showing visible-light-induced production of exogenous H₂O₂ by free-standing, gold-decorated silicon nanowires internalized in human umbilical vein endothelial cells. We further show that the photocatalytic production of H₂O₂ is a general phenomenon of gold-silicon hybrid materials and is enhanced upon annealing.

Machine learning dielectric screening for the simulation of excited state properties of molecules and materials

Accurate and efficient calculations of absorption spectra of molecules and materials are essential for the understanding and rational design of broad classes of systems. Solving the Bethe-Salpeter equation (BSE) for electron-hole pairs usually yields accurate predictions of absorption spectra, but it is computationally expensive, especially if thermal averages of spectra computed for multiple configurations are required. We present a method based on machine learning to evaluate a key quantity entering the definition of absorption spectra: the dielectric screening. We show that our approach yields a model for the screening that is transferable between multiple configurations sampled during first principles molecular dynamics simulations; hence it leads to a substantial improvement in the efficiency of calculations of finite temperature spectra. We obtained computational gains of one to two orders of magnitude for systems with 50 to 500 atoms, including liquids, solids, nanostructures, and solid/liquid interfaces. Importantly, the models of dielectric screening derived here may be used not only in the solution of the BSE but also in developing functionals for time-dependent density functional theory (TDDFT) calculations of homogeneous and heterogeneous systems. Overall, our work provides a strategy to combine machine learning with electronic structure calculations to accelerate first principles simulations of excited-state properties.

van der Waals heterostructures based on MSSe (M = Mo, W) and graphene-like GaN: enhanced optoelectronic and photocatalytic properties for water splitting

The geometric structure, electronic, optical and photocatalytic properties of MSSe-g-GaN (M = Mo, W) van der Waals (vdW) heterostructures are investigated by performing first-principles calculations. We find that the MoSSe-g-GaN heterostructure exhibits type-II band alignment for all stacking patterns. While the WSSe-g-GaN heterostructure forms the type-II or type-I band alignment for the stacking model-I or model II, respectively. The average electrostatic potential shows that the potential of g-GaN is deeper than the MSSe monolayer, leading to the formation of an electrostatic field across the interface, causing the transfer of photogenerated electrons and holes. Efficient interfacial formation of interface and charge transfer reduce the work function of MSSe-g-GaN vdW heterostructures as compared to the constituent monolayer. The difference in the carrier mobility for electrons and holes suggests that these heterostructures could be utilized for hole/electron separation. Absorption spectra demonstrate that strong absorption from infrared to visible light in these vdW heterostructures can be achieved. Appropriate valence and conduction band edge positions with standard redox potentials provide enough force to drive the photogenerated electrons and holes to dissociate water into H+/H2 and O2/H2O at pH = 0.

Electronic properties and enhanced photocatalytic performance of van der Waals heterostructures of ZnO and Janus transition metal dichalcogenides

Vertical stacking of two-dimensional materials into layered van der Waals heterostructures has recently been considered as a promising candidate for photocatalytic and optoelectronic devices because it can combine the advantages of the individual 2D materials. Janus transition metal dichalcogenides (JTMDCs) have emerged as an appealing photocatalytic material due to the desirable electronic properties. Hence, in this work, we systematically investigate the geometric features, electronic properties, charge density difference, work function, band alignment and photocatalytic properties of ZnO-JTMDC heterostructures using first-principles calculations. Due to the different kinds of chalcogen atoms on both sides of JTMDC monolayers, two different possible stacking patterns of ZnO-JTMDC heterostructures have been constructed and considered. We find that all these stacking patterns of ZnO-JTMDC heterostructures are dynamically and energetically feasible. Moreover, both ZnO-MoSSe and ZnO-WSSe heterostructures are indirect band gap semiconductors and present type-I and type-II band alignments for model-I and model-II, respectively. The Rashba spin polarization of the ZnO-WSSe heterostructure for model-I is greater than that in the others. Furthermore, valence (conduction) band edge potentials are calculated to understand the photocatalytic behavior of these systems. Energetically favorable band edge positions in ZnO-Janus heterostructures make them suitable for water splitting at zero pH. We found that the ZnO-Janus heterostructures are promising candidates for water splitting with conduction and valence band edges positioned just outside of the redox interval.

The electrochemical interface in first-principles calculations

First-principles predictions play an important role in understanding chemistry at the electrochemical interface. Electronic structure calculations are straightforward for vacuum interfaces, but do not easily account for the interfacial fields and solvation that fundamentally change the nature of electrochemical reactions. Prevalent techniques for first-principles prediction of electrochemical processes range from expensive explicit solvation using ab initio molecular dynamics, through a hierarchy of continuum solvation techniques, to neglecting solvation and interfacial field effects entirely. Currently, no single approach reliably captures all relevant effects of the electrochemical double layer in first-principles calculations. This review systematically lays out the relation between all major approaches to first-principles electrochemistry, including the key approximations and their consequences for accuracy and computational cost. Focusing on ab initio methods for thermodynamic properties of aqueous interfaces, we first outline general considerations for modeling electrochemical interfaces, including solvent and electrolyte dynamics and electrification. We then present the specifics of various explicit and implicit models of the solvent and electrolyte. Finally, we discuss the compromise between computational efficiency and accuracy, and identify key outstanding challenges and future opportunities in the wide range of techniques for first-principles electrochemistry.

Computational Approaches to Photoelectrode Design through Molecular Functionalization for Enhanced Photoelectrochemical Water Splitting

Photoelectrochemical water splitting is a promising carbon-free approach to produce hydrogen from water. A photoelectrochemical cell consists of a semiconductor photoelectrode in contact with an aqueous electrolyte. Its performance is sensitive to properties of the photoelectrode/electrolyte interface, which may be tuned through functionalization of the photoelectrode surface with organic molecules. This can lead to improvements in the photoelectrode's properties. This Minireview summarizes key computational investigations on using molecular functionalization to modify photoelectrode stability, barrier height, and catalytic activity. It is discussed how first-principles density functional theory, first-principles molecular dynamics, and device modeling simulations can provide predictive insights and complement experimental investigations of functionalized photoelectrodes. Challenges and future directions in the computational modeling of functionalized photoelectrode/electrolyte interfaces within the context of experimental studies are also highlighted.

Interfacial Effects on the Band Edges of Ta3N5 Photoanodes in an Aqueous Environment: A Theoretical View

Ta₃N₅, as a fascinating photoanode for solar hydrogen production, is expected to split water without any bias, because its band edge potentials straddle H₂O redox potentials. Unfortunately, Ta₃N₅ photoanodes can split water only when a bias of at least 0.6–0.9 V is applied. It means that they exhibit an onset potential as high as 0.6–0.9 V_RHE (reversible hydrogen electrode). In this study, density functional theory calculations show that the band edge potentials of Ta₃N₅ have a shift of approximately −0.42 eV relative to vacuum level when exposed to water. The increased ratio of dissociated water at Ta₃N₅-water interface will further make the band edge potentials shift −0.85 eV relative to vacuum level, implying the anodic shifts of the onset potential for water oxidation. The findings may reveal the mystery of the unexpectedly high onset potential of Ta₃N₅, as high as 0.6–0.9 V_RHE.

•

We have studied interfacial effects on the band edges of Ta₃N₅ in an aqueous environment

•

Both water and the hydroxylated surface promote the formation of the interface dipole

•

High onset potentials of Ta₃N₅ may be ascribed to negative shift of band edge potentials

Optoelectronic and solar cell applications of Janus monolayers and their van der Waals heterostructures

Janus monolayers and their van der Waals heterostuctures are investigated by hybrid density functional theory calculations.

The role of defects and excess surface charges at finite temperature for optimizing oxide photoabsorbers

Computational screening of materials for solar to fuel conversion technologies has mostly focused on bulk properties, thus neglecting the structure and chemistry of surfaces and interfaces with water. We report a finite temperature study of WO₃, a promising anode for photoelectrochemical cells, carried out using first-principles molecular dynamics simulations coupled with many-body perturbation theory. We identified three major factors determining the chemical reactivity of the material interfaced with water: the presence of surface defects, the dynamics of excess charge at the surface, and finite temperature fluctuations of the surface electronic orbitals. These general descriptors are essential for the understanding and prediction of optimal oxide photoabsorbers for water oxidation.

Probing the Ion-Specific Effects at the Water/Air Interface and Water-Mediated Ion Pairing in Sodium Halide Solution with Ab Initio Molecular Dynamics

The ion-specific effects at the water/air interface represent a fundamentally essential topic of research, and high-level ab initio simulations are still demanding to reveal the microscopic picture of the interactions between ions and water at the solvation interface. In this work, we present a fragment-based ab initio molecular dynamics (AIMD) simulation of sodium halide solution droplet (in a neutral mixture of Na⁺, F^-, Cl^-, and Br^- ions) at the MP2/aug-cc-pVDZ level. We show that the studied halide ions exhibit surface preference in the order (F^- < Cl^- < Br^-) which is in accordance with the experimental observation. The resulting potential of mean force (PMF) for Br^- produces a distinct minimum at the water/air interface, while the minimum of the PMF for F^- appears in the bulk region. The ion-pairing interactions between halide anions and Na⁺ cations are characterized, and it reveals that the specific solvent-separated ion pairs (SIPs) are more preferred than the direct contact ion pairs (CIPs). The transition between different types of SIPs is observed. Other structural and dynamical properties of ions and ion-hydration shells are investigated. These results provide broader and new physical insights for understanding the ion-specific behavior in interfacial solvation at the atomistic level.

Where do photogenerated holes at the g-C3N4/water interface go for water splitting: H2O or OH-?

Graphitic carbon nitride (g-C3N4), a metal-free two-dimensional photocatalyst, has drawn increasing attention due to its application in photocatalytic water splitting. However, its quantum efficiency is limited by the poor performance of the oxygen evolution reaction (OER). Therefore, it is important to clarify the behavior of photogenerated holes in the OER. In this work, we investigate the energy level alignment using the GW method and the exciton properties using the Bethe-Salpeter equation within the ab initio many-body Green's function theory at the g-C3N4/water interface. We found that the g-C3N4 substrate can elevate energy levels of OH- and H2O molecules at the interface by up to 0.6 eV. This effect can make the electronic levels of OH- surpass the valence band maximum (VBM) of g-C3N4. However, orbital energies of H2O molecules remain far below the VBM of g-C3N4. This indicates that a photogenerated hole after exciting g-C3N4 can relax to OH- instead of neutral H2O. Moreover, OH- could be directly oxidized through electron transfer from OH- to g-C3N4 by light near the optical absorption edge of g-C3N4, which is beneficial for efficient carrier separation at the interface.

Interface Engineering of Monolayer MoS2/GaN Hybrid Heterostructure: Modified Band Alignment for Photocatalytic Water Splitting Application by Nitridation Treatment

Interface engineering is a key strategy to deal with the two-dimensional (2D)/three-dimensional (3D) hybrid heterostructure, since the properties of this atomic-layer-thick 2D material can easily be impacted by the substrate environment. In this work, the structural, electronic, and optical properties of the 2D/3D heterostructure of monolayer MoS₂ on wurtzite GaN surface without and with nitridation interfacial layer are systematically investigated by first-principles calculation and experimental analysis. The nitridation interfacial layer can be introduced into the 2D/3D heterostructure by remote N₂ plasma treatment to GaN sample surface prior to stacking monolayer MoS₂ on top. The calculation results reveal that the 2D/3D integrated heterostructure is energetically favorable with a negative formation energy. Both interfaces demonstrate indirect band gap, which is a benefit for longer lifetime of the photoexcited carriers. Meanwhile, the conduction band edge and valence band edge of the MoS₂ side increases after nitridation treatment. The modification to band alignment is then verified by X-ray photoelectron spectroscopy measurement on MoS₂/GaN heterostructures constructed by a modified wet-transfer technique, which indicates that the MoS₂/GaN heterostructure without nitridation shows a type-II alignment with a conduction band offset (CBO) of only 0.07 eV. However, by the deployment of interface nitridation, the band edges of MoS₂ move upward for ∼0.5 eV as a result of the nitridized substrate property. The significantly increased CBO could lead to better electron accumulation capability at the GaN side. The nitridized 2D/3D heterostructure with effective interface treatment exhibits a clean band gap and substantial optical absorption ability and could be potentially used as practical photocatalyst for hydrogen generation by water splitting using solar energy.

Design Strategy for the Molecular Functionalization of Semiconductor Photoelectrodes: A Case Study of p-Si(111) Photocathodes for H2 Generation

Functionalization of semiconductor photoelectrodes is actively pursued as an approach to improve the efficiency of photoelectrochemical reactions by modulating the semiconductor's barrier height, but the selection of molecules for functionalization remains largely empirical. We propose a simple but effective design strategy for the organic functionalization of photocathodes for high-efficiency hydrogen generation based on first-principles density functional theory (DFT) calculations. The surface dipole of the functionalized photocathode determines its barrier height, which can be optimized to enhance charge separation at the semiconductor-electrolyte interface. Focusing on p-Si(111) photocathodes functionalized with different mixed aryl/methyl monolayers, we use DFT to systematically investigate the effect of - the presence and distribution of pi bonds, binding atom (the atom in the functional group that bonds with the Si surface), functional group length, and electrophilic substituent groups - on the surface dipole and charge rearrangement at the functionalized surface. We find that the most important factors affecting the surface dipole are the intrinsic molecular dipole moment of the organic moiety, the presence of electrophilic substituent groups, and the binding atom. Using these findings, a three-step design strategy is proposed for the experimental realization of high-performing functionalized p-Si(111) photocathodes by maximizing the surface dipole.

Accuracy of dielectric-dependent hybrid functionals in the prediction of optoelectronic properties of metal oxide semiconductors: a comprehensive comparison with many-body GW and experiments

Understanding the electronic structure of metal oxide semiconductors is crucial to their numerous technological applications, such as photoelectrochemical water splitting and solar cells. The needed experimental and theoretical knowledge goes beyond that of pristine bulk crystals, and must include the effects of surfaces and interfaces, as well as those due to the presence of intrinsic defects (e.g. oxygen vacancies), or dopants for band engineering. In this review, we present an account of the recent efforts in predicting and understanding the optoelectronic properties of oxides using ab initio theoretical methods. In particular, we discuss the performance of recently developed dielectric-dependent hybrid functionals, providing a comparison against the results of many-body GW calculations, including G ₀ W ₀ as well as more refined approaches, such as quasiparticle self-consistent GW. We summarize results in the recent literature for the band gap, the band level alignment at surfaces, and optical transition energies in defective oxides, including wide gap oxide semiconductors and transition metal oxides. Correlated transition metal oxides are also discussed. For each method, we describe successes and drawbacks, emphasizing the challenges faced by the development of improved theoretical approaches. The theoretical section is preceded by a critical overview of the main experimental techniques needed to characterize the optoelectronic properties of semiconductors, including absorption and reflection spectroscopy, photoemission, and scanning tunneling spectroscopy (STS).

Integrating Ab Initio Simulations and X-ray Photoelectron Spectroscopy: Toward A Realistic Description of Oxidized Solid/Liquid Interfaces

Many energy storage and conversion devices rely on processes that take place at complex interfaces, where structural and chemical properties are often difficult to probe under operating conditions. A primary example is solar water splitting using high-performance photoelectrochemical cells, where surface chemistry, including native oxide formation, affects hydrogen generation. In this Perspective, we discuss some of the challenges associated with interrogating interface chemistry, and how they may be overcome by integrating high-level first-principles calculations of explicit interfaces with ambient pressure X-ray photoelectron spectroscopy and direct spectroscopic simulations. We illustrate the benefit of this combined approach toward insights into native oxide chemistry at prototypical InP/water and GaP/water interfaces. This example suggests a more general roadmap for obtaining a realistic and reliable description of the chemistry of complex interfaces by combining state-of-the-art computational and experimental techniques.

Electrochemical Capacitance of CO-Terminated Pt(111) Dominated by the CO-Solvent Gap

The distribution of electric fields within the electrochemical double layer depends on both the electrode and electrolyte in complex ways. These fields strongly influence chemical dynamics in the electrode-electrolyte interface but cannot be measured directly with submolecular resolution. We report experimental capacitance measurements for aqueous interfaces of CO-terminated Pt(111). By comparing these measurements with first-principles density functional theory (DFT) calculations, we infer microscopic field distributions and decompose contributions to the inverse capacitance from various spatial regions of the interface. We find that the CO is strongly electronically coupled to the Pt and that most of the interfacial potential difference appears across the gap between the terminating O and water and not across the CO molecule, as previously hypothesized. This "gap capacitance" resulting from hydrophobic termination lowers the overall capacitance of the aqueous Pt-CO interface and makes it less sensitive to electrolyte concentration compared to the bare metal.

Theoretical Determination of Band Edge Alignments at the Water-CuInS2(112) Semiconductor Interface

Knowledge of a semiconductor electrode's band edge alignment is essential for optimizing processes that occur at the semiconductor/electrolyte interface. Photocatalytic processes are particularly sensitive to such alignments, as they govern the transfer of photoexcited electrons or holes from the surface to reactants in the electrolyte solution. Reconstructions of a semiconductor surface during operation, as well as its interaction with the electrolyte solution, must be considered when determining band edge alignment. Here, we employ density functional theory + U theory to assess the stability of reconstructed CuInS₂ surfaces, a system which has shown promise for the active and selective photoelectrocatalytic reduction of CO₂ to CH₃OH. Using many-body Green's function theory combined with calculations of surface work functions, we determine band edge positions of explicitly solvated, reconstructed CuInS₂ surfaces. We find that there is a linear relationship between band edge position and net surface dipole, with the most stable solvent/surface structures tending to minimize this dipole because of generally weak interactions between the surface and solvating water molecules. We predict a conduction band minimum (CBM) of the solvated, reconstructed CuInS₂ surface of -2.44 eV vs vacuum at the zero-dipole intercept of the dipole/CBM trendline, in reasonable agreement with the experimentally reported CBM position at -2.64 eV vs vacuum. This methodology offers a simplified approach for approximating the band edge positions at complex semiconductor/electrolyte interfaces.

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.

Electronic structure of aqueous solutions: Bridging the gap between theory and experiments

Predicting the electronic properties of aqueous liquids has been a long-standing challenge for quantum mechanical methods. However, it is a crucial step in understanding and predicting the key role played by aqueous solutions and electrolytes in a wide variety of emerging energy and environmental technologies, including battery and photoelectrochemical cell design. We propose an efficient and accurate approach to predict the electronic properties of aqueous solutions, on the basis of the combination of first-principles methods and experimental validation using state-of-the-art spectroscopic measurements. We present results of the photoelectron spectra of a broad range of solvated ions, showing that first-principles molecular dynamics simulations and electronic structure calculations using dielectric hybrid functionals provide a quantitative description of the electronic properties of the solvent and solutes, including excitation energies. The proposed computational framework is general and applicable to other liquids, thereby offering great promise in understanding and engineering solutions and liquid electrolytes for a variety of important energy technologies.

Impact of Interfacial Defects on the Properties of Monolayer Transition Metal Dichalcogenide Lateral Heterojunctions

We explored the impact of interfacial defects on the stability and optoelectronic properties of monolayer transition metal dichalcogenide lateral heterojunctions using a density functional theory approach. As a prototype, we focused on the MoS₂-WSe₂ system and found that even a random alloy-like interface with a width of less than 1 nm has only a minimal impact on the band gap and alignment compared to the defect-less interface. The largest impact is on the evolution of the electrostatic potential across the monolayer. Similar to defect-less interfaces, a small number of defects results in an electrostatic potential profile with a sharp change at the interface, which facilitates exciton dissociation. Differently, a large number of defects results in an electrostatic potential profile switching smoothly across the interface, which is expected to reduce the capability of the heterojunction to promote exciton dissociation. These results are generalizable to other transition metal dichalcogenide lateral heterojunctions.

Modelling heterogeneous interfaces for solar water splitting

The generation of hydrogen from water and sunlight offers a promising approach for producing scalable and sustainable carbon-free energy. The key of a successful solar-to-fuel technology is the design of efficient, long-lasting and low-cost photoelectrochemical cells, which are responsible for absorbing sunlight and driving water splitting reactions. To this end, a detailed understanding and control of heterogeneous interfaces between photoabsorbers, electrolytes and catalysts present in photoelectrochemical cells is essential. Here we review recent progress and open challenges in predicting physicochemical properties of heterogeneous interfaces for solar water splitting applications using first-principles-based approaches, and highlights the key role of these calculations in interpreting increasingly complex experiments.

Energy level alignment at semiconductor–water interfaces from atomistic and continuum solvation models

Efficient electronic energy level alignment at solid–liquid interfaces with continuum solvation models.

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%.

Catalytic Reduction of CO2 by Renewable Organohydrides

Dihydropyridines are renewable organohydride reducing agents for the catalytic reduction of CO2 to MeOH. Here we discuss various aspects of this important reduction. A centerpiece, which illustrates various general principles, is our theoretical catalytic mechanism for CO2 reduction by successive hydride transfers (HTs) and proton transfers (PTs) from the dihydropyridine PyH2 obtained by 1H(+)/1e(-)/1H(+)/1e(-) reductions of pyridine. The Py/PyH2 redox couple is analogous to NADP(+)/NADPH in that both are driven to effect HTs by rearomatization. High-energy radical intermediates and their associated high barriers/overpotentials are avoided because HT involves a 2e(-) reduction. A HT-PT sequence dictates that the reduced intermediates be protonated prior to further reduction for ultimate MeOH formation; these protonations are aided by biased cathodes that significantly lower the local pH. In contrast, cathodes that efficiently reduce H(+) such as Pt and Pd produce H2 and create a high interfacial pH, both obstructing dihydropyridine production and formate protonation and thus ultimately CO2 reduction by HTPTs. The role of water molecule proton relays is discussed. Finally, we suggest future CO2 reduction strategies by organic (photo)catalysts.

Solvation effects on the band edge positions of photocatalysts from first principles

Continuum solvation theories predict large shifts in band positions of photocatalysts from vacuum to solution, in agreement with experiment.