Integrating Ab Initio Simulations and X-ray Photoelectron

Dec 14, 2017 - He is the LLNL Laboratory Lead for the Hydrogen Storage Materials—Advanced Research Consortium (HyMARC), funded by the Fuel Cell Tech...
1 downloads 6 Views 493KB Size
Subscriber access provided by Binghamton University | Libraries

Perspective

Integrating Ab Initio Simulations and X-ray Photoelectron Spectroscopy: Towards A Realistic Description of Oxidized Solid/Liquid Interfaces Tuan Anh Pham, Xueqiang Zhang, Brandon C. Wood, David Prendergast, Sylwia Ptasinska, and Tadashi Ogitsu J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.7b01382 • Publication Date (Web): 14 Dec 2017 Downloaded from http://pubs.acs.org on December 18, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry Letters is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 17

Integrating Ab Initio Simulations and X-ray Photoelectron Spectroscopy: Towards A Realistic Description of Oxidized Solid/Liquid Interfaces Tuan Anh Pham,∗,†,⊥ Xueqiang Zhang,∗,‡,¶,⊥ Brandon C. Wood,∗,† David Prendergast,∗,§ Sylwia Ptasinska,∗,‡,k and Tadashi Ogitsu∗,† Quantum Simulations Group, Lawrence Livermore National Laboratory, Livermore, CA 94551, United States, Radiation Laboratory, University of Notre Dame, Notre Dame, IN 46556, United States, Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, IN 46556, United States, Molecular Foundry, Lawrence Berkeley National Lab, Berkeley, CA 94720, United States, and Department of Physics, University of Notre Dame, Notre Dame, IN 46556, United States E-mail: [email protected]; [email protected]; [email protected]; [email protected]; [email protected]; [email protected]

Abstract

calculations of explicit interfaces with ambientpressure X-ray photoelectron spectroscopy and direct spectroscopic simulations. We illustrate the benefit of this combined approach towards 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.

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

Graphical TOC Entry



To whom correspondence should be addressed Quantum Simulations Group, Lawrence Livermore National Laboratory, Livermore, CA 94551, United States ‡ Radiation Laboratory, University of Notre Dame, Notre Dame, IN 46556, United States ¶ Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, IN 46556, United States § Molecular Foundry, Lawrence Berkeley National Lab, Berkeley, CA 94720, United States k Department of Physics, University of Notre Dame, Notre Dame, IN 46556, United States ⊥ Contributed equally to this work †

XPS

Intensity (Arb. units)

H2O/InP(001)

H2O*

534

M-O*H-M

533

532

M-O*H

531

M-O*-M

530

529

Binding energy (eV)

7

Band edge analyses InP

|VBM| (eV)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

θ=1.5 ML

6 5

δ(2x4) θ=0 ML

4 3

Vacuum

ACS Paragon Plus Environment

1

Water

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Understanding and controlling the chemistry of complex solid/liquid interfaces is essential for discovering and optimizing emerging energy materials, including photoelectrochemical cells (PEC). These devices, which produce hydrogen via solar water splitting in a process long considered the “Holy Grail” of electrochemistry, 1–3 require efficient and stable operation. This remains a challenging task, as many of the most efficient photoabsorbing materials, such as IIIV semiconductors, are unstable under harsh operating conditions at the solid/liquid interface in a PEC device. 4 On the other hand, metal oxides, although being more stable, possess transport properties and band gaps that are unfavorable for achieving high solar-to-hydrogen conversion efficiency. 5,6 In recent years, significant efforts have been directed toward improving the stablity and efficiency of photoabsorbing materials by interfacing the photoelectrode absorbers with catalysts and/or protective layers. 5,7 One important aspect of interface chemistry that has been discussed at length in the context of certain PEC semiconductors is the possible role of oxidation. Native oxidation occurs spontaneously on the surfaces of many non-oxide PEC photoelectrodes under operation, and can also be introduced deliberately during processing to alter electrochemical properties. 5,8–14 The various negative and positive impacts of these oxides continue to be the subject of much debate. The insulating nature of oxides can inhibit charge transfer and introduce additional kinetic barriers to redox processes; 8 for instance, in the case of GaP, it has shown that water-splitting reactions may be limited by the formation of a thin Ga2 O3 oxide layer due to insufficient band bending and unfavorable band alignment at the GaP/Ga2 O3 interface. 15 Oxide-related defect states in III-V semiconductors might also act as traps that can induce corrosion processes. 9 At the same time, it has been suggested that oxides can aid hydrogen evolution, possibly due to improved kinetics for hydrogen transport. 10–12,16 Oxide films have also been investigated as a potential strategy for stabilizing surfaces of photocorrosive electrode materials. 5,13,14

Page 2 of 17

These often contradictory roles highlight the fact that much is left to be understood regarding the nature and specific impacts of the surface oxidation at complex semiconductor/water interfaces. Many of the key knowledge gaps relate to the electronic properties of the interfaces, most notably the effects on the absolute energetic alignment of the semiconductor valence band maximum (VBM) and conduction band minimum (CBM) with respect to the various redox processes (e.g., hydrogen evolution versus corrosion) that might occur in the device. 1 While it is generally understood that these band positions depend on the details of the semiconductor interface, 17 the specific relationship between the electronic properties of oxidized PEC surfaces and their chemical constituents remains enigmatic. Realization of this correlation would open up opportunities not only for tuning interfacial chemistry to achieve highly efficient PECs but also for quantitative predictions of photoelectrode/water interfacial structure based on macroscopic observations such as experimental band edge measurements. In principle, obtaining the necessary understanding of chemical and electronic properties of interfaces requires experimental probes that can achieve atomic sensitivity under operating conditions, which currently presents significant technical challenges. 18 Band edge measurements yield information on the electronic structure of the photoelectrodes, but provide little understanding on the chemistry at the photoelectrode/liquid interface. In this regard, recent development of ambient pressure X-ray photoelectron spectroscopy (APXPS) offers a promising way for investigating chemical species and intermediates of solid/liquid interfaces, as this technique directly probes chemical processes near surfaces at elevated vapor pressures, which approximate in situ conditions that cannot easily be achieved with standard measurements. 19,20 In practice, however, definitive XPS identification of chemically similar species can be extremely challenging, particularly for complex amorphous oxides where intrinsic configurational heterogeneity significantly broadens

ACS Paragon Plus Environment

2

Page 3 of 17 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

and overlaps key spectral features. For wellcharacterized surfaces in vacuum or gas-phase adsorbents, experimental standards may aid in the interpretation of photoemission spectra; however, such reference standards are likely to be far less reliable for these more complex interfacial systems and intermediate species. In addition to difficulties in the preparation of standard samples, definitive fingerprinting of interfacial XPS spectral signatures using experimental references tends to be limited to cases in which discrete changes in formal oxidation state are encountered. On the other hand, subtler changes in local chemical environment may produce distinctions that are unique to the surface chemistry and cannot be unambiguously identified using existing bulk or molecular standards. In this regard, experiments can benefit from additional information derived from simulation and modelling, since recent advances in highperformance computing and theoretical techniques now make it possible to computationally probe the chemistry and electronic properties of the photoelectrode/water interface more explicitly. 16,21–38 In particular, first-principles molecular dynamics (FPMD) simulations allow for accessing complex interfacial chemical processes, 39 while high-level electronic structure theory, such as many-body perturbation theory within the GW approximation, 40 offers a promising way to compute accurate properties of interfaces, due to its ability to treat the electronic structure of semiconductors and aqueous solutions on the same footing. 41–49 In recent years, these techniques have advanced to the point where quantitative accuracy is realistically achievable on modern computing platforms. Nevertheless, the length and time scales accessible to first-principles simulations continue to constrain the configurational complexity. Perhaps more importantly, the exact results obtained from first-principles techniques depend strongly on the precise details of the model. Such details are particularly difficult to intuit for systems with a high degree of structural and/or chemical complexity, such as kinetically stabilized surface oxides. For these systems, the predictive power of simulations becomes critically dependent on experimental in-

sights. The rapid developments in advanced computational and experimental methods, as well as the mutual interdependence of the two approaches, afford a unique opportunity to directly integrate advanced first-principles methods with APXPS experiments for interrogating complex interfaces. One way to achieve this is to compute XPS core-binding energies directly from first-principles based on proposed interface models. This allows for direct comparison between chemical processes of theoretical models and those of realistic photoelectrode/water interfaces, while providing a more complete interpretation of the experimental photoemission spectra. By taking advantage of a close feedback loop between computational and experimental tools, one can hope to fully elucidate the photoelectrode chemistry (e.g., oxidation), as well as its relationship to the observed electronic properties (e.g., band edges) during operation. In this perspective, we provide a roadmap for how high-level electronic structure theory may be combined with APXPS and other in situ experimental data as a first step towards probing surface chemistry at a complex solid/liquid interface. We illustrate the specific procedure by applying it to a demonstration PEC system, focusing on the nature and role of native surface oxides of GaP and InP in water. Notably, these materials can be alloyed for highly efficient PEC photocathodes but have persistent stability issues. 4,7 In addition to a better understanding of the interfacial chemistry, we show how the approach can offer insight into the specific correlation between surface composition and band edges of photoabsorbers in liquid water, an important consideration for functionality and stability of a PEC device. The theoretical framework is based on FPMD simulations in conjunction with GW calculations, and leverages recently developed algorithms that now make it possible to obtain reliable quantitative predictions of band edge positions within full interface simulations. In our demonstration system, model InP and GaP surfaces representing different degrees of oxidation are employed to investigate the influence of sur-

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry Letters

Page 4 of 17

-2

-3

Energy (eV)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1: Models of InP/water interfaces considered in this work: Pristine reconstructed δ(2 × 4)(001) (left) and hydroxylated surfaces (right).

Eg=2.44 eV

-4

Eg=1.25 eV

SHE -5

-6 (110) (001)

GaP

(110) (001)

InP

Figure 2: Band edge positions of pristine In(Ga)P(110) and (001) surfaces in vacuum, as computed with the G0 W0 approximation. The G0 W0 results were obtained with the wavefunctions obtained with DFT and the PBE exchange-correlation functional. Experimental results are presented by blue lines. 52,53 The potential of the standard hydrogen electrode (SHE) is indicated by the dashed line.

face oxide species on the band edges, which is shown to be significant due to the change in the surface polarization. These theoretical results are combined with experimental band edge measurements to probe and compare the general chemical properties of model semiconductor surfaces in water. Finally, we outline a procedure for how APXPS experiments can be closely integrated with direct first-principles core-binding energy calculations of semiconductor/water interfaces to provide a more complete interpretation of the surface chemistry, including the contribution of specific oxygen-derived surface species to the photoelectron spectrum. In the case of III-V semiconductors, the results reveal clear and interesting differences in the structure and chemistry of native oxides on GaP and InP photoelectrodes at the interface with water, demonstrating the power of the integrated theory-experiment approach. We conclude by highlighting opportunities for improving and extending this approach towards a comprehensive and realistic description of complex interfaces under device operation conditions. Before discussing the results, we briefly summarize the theoretical models employed in this work. Specifically, to capture the structure and chemistry of GaP and InP photoelectrodes under ideal and working conditions, we considered pristine GaP(001) and InP (001) surfaces in vacuum, as well as different models of oxidized surfaces in direct contact with liquid water. For the former, we focus on the mixeddimer δ(2×4) In- and Ga-rich (001) reconstruction that is generally preferred under vacuum conditions, 50 while also making comparisons to non-polar (110) surfaces in order to address the

role of surface polarization. For the surfaces fully in contact with liquid water, we considered two model systems that represent two extreme limits of oxidation extent, which were generated using FPMD simulations. The mixeddimer δ(2×4) reconstruction with zero oxygen coverage was used to simulate the low oxidation limit. For the high oxidation limit, we used a fully hydroxylated Ga- and In-rich (001) surface derived from our prior investigations. 16,36 This model is based on an oxygen coverage of θ = 1.5 ML, consisting of surface M –OH–M bridges along h110i with M –OH dangling bond atop structures (M ≡ In or Ga). As we will show below, the consideration of model systems representing extreme limits of the oxidation degree allows for infering how the realistic surface chemistry envolves under operation. Representative snapshots of mixed-dimer δ(2×4) and hydroxylated surfaces in contact with liquid water are shown in Fig. 1 for InP. Further details of the model systems can be found in the Supporting Information and our previous publications. 9,16,36,51 Relationship Between Surface Oxidation and Electronic Properties. The relation between electronic properties and surface oxidation can be established based on a combination of highlevel electronic structure theory and estab-

ACS Paragon Plus Environment

4

Page 5 of 17 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

lished experimental band edges. Experimentally, the measured absolute band edges (VBM and CBM) of pristine semiconductor surfaces vac in vacuum (EVBM/CBM ) correspond to the ionization potential and the electron affinity, and are typically measured using photoemission experiments. On the other hand, band edges wat of semiconductors in liquid water (EVBM/CBM ) are usually obtained from capacitance-voltage measurements of the semiconductor/water interface (see Supporting Information), where the wat vac difference between EVBM/CBM and EVBM/CBM represents the collective effect of water on the band edges of the semiconductor surface. 54 In vac our calculations, EVBM/CBM are computed using many-body perturbation theory within the G0 W0 approximation 55 following the procedure wat described in Ref. 33,56. To compute EVBM/CBM , the effects of liquid water were taken into account by considering the variation in the electrostatic potential at the interface, as described more fully in Ref. 33. This method has the advantage of treating the electronic structure of the semiconductor and liquid water on equal footing, 43–46 and has previously demonstrated predictive capability for interfaces. 33 Further details of band edge calculations are described in the Supporting Information. To baseline the discussion of surface oxidation in water and to validate the predictive accuracy of our theoretical approach, we first summarize the first-principles results for pristine InP and GaP surfaces in vacuum, for which the structural and electronic properties are established. 52,53 We find that the G0 W0 approach accurately describes the electronic properties of both systems, yielding band gaps of 2.44 eV and 1.25 eV for GaP and InP, respectively. These values are in good agreement with experimental values of 2.26 eV and 1.35 eV. 52,53 Furthermore, the predicted VBM/CBM for InP (110) and δ(2×4) (001) surfaces are consistent with available experiments, yielding errors less than 0.15 eV as demonstrated in Fig. 2. 52,53 Good agreement between theory and experiment is also found for the GaP (110) surface 52,53 (experimental band positions for δ(2×4) GaP (001) surface are not available for comparison). This is a significant improvement over the conven-

tional density functional theory (DFT) calculations with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional, 57 which produce sizable errors in the computed band gap (about 50%) and band edge positions (about 0.6–0.7 eV). The results highlight the importance of surface structure and composition in determining the photoelectrode band edges, even in the absence of liquid water. For instance, δ(2×4) (001) surfaces exhibit band edges closer to vacuum than (110) surfaces by 0.7–1.0 eV for both GaP and InP. This shift in the electronic energy levels can be qualitatively explained based on the linear variation of the surface potential associated with the change in the surface dipole moment: 17 ∆Vsurf = −4πµsurf /S0 ,

(1)

where µsurf is the surface dipole moment component normal to the surface plane for a unit cell of area S0 ; a negative µsurf implies a positive ∆Vsurf and band positions further from the vacuum, and vice versa. According to the electronegativity relation χP = 2.19 > χIn/Ga = 1.78/1.81, electrons are extracted from In/Ga, resulting in a positive µsurf for the In/Ga rich δ(2×4)(001) surfaces. The band edges therefore move closer toward the vacuum level than for the non-polar (110) surfaces. Accordingly, effects of surface polarization are critical in the determination of the electronic properties of pristine GaP/InP surfaces. We now turn to the relation between electronic properties and surface oxidation for InP and GaP in the presence of liquid water. We note that our computational approach probes the band edges in the bulk region of photoelectrode materials where the bulk band gap is recovered; accordingly, the computed CBM and VBM follow the same trend in the presence of liquid water, reflecting the band edge evolution obtained in differential capacitance measurements. 59 We find the effects of water on the band edges are significant, as demonstrated in Fig. 3. Focusing on the VBM, the level for pristine δ(2×4) GaP(001) and InP(001) surfaces is shifted toward the vacuum level by 0.8–1.0 eV

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry Letters 7

7

InP 5

GaP

θ=1.5 ML

6

|VBM| (eV)

|VBM| (eV)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 17

δ(2x4) θ=0 ML

4

θ=1.5 ML

6 5

δ(2x4) θ=0 ML

4

3

3

Vacuum

Water

Vacuum

Water

Figure 3: Absolute values of the valence band maximum of GaP/InP in vacuum (mixed-dimer δ(2 × 4) surfaces), and in the presence of liquid water, including mixed-dimer δ(2 × 4) (θ = 0, magenta solid line) and hydroxylated (θ = 1.5 monolayer (ML), green solid line) surfaces. Results computed for the corresponding surfaces with a monolayer of water molecules and hydroxyls in vacuum are presented by magenta and green dashed lines, respectively. All theoretical results were obtained with the G0 W0 approximation using the PBE wavefunction. Experimental results obtained for pH=0–14 and neutral water are represented by filled boxes and the solid blue lines, which were deduced from Wrighton et al., 58 as discussed in more detail in the Supporting Information. InP(001) δ(2×4) surfaces with a monolayer of adsorbed water molecules, as well as hydroxylated surfaces with water removed. These surfaces are broadly representative of configurations with a monolayer of molecularly adsorbed (undissociated) water molecules and a monolayer of dissociated water molecules, respectively. Calculated band edges of these models in vacuum are represented by the dashed magenta and green lines in Fig. 3, respectively. Interestingly, the overall trend in the band edge shifts is similar to what was found for the full interface with liquid water (solid lines, Fig. 3). Accordingly, the photoelectrode band edges are largely determined by these first monolayers of adsorbed water molecules or hydroxyl groups. This behavior can be understood by considering that adsorbed water molecules on the mixeddimer surfaces introduce a net positive surface dipole, resulting in VBMs closer to the vacuum. On the other hand, the dense oxygen coverage on the hydroxylated surfaces leads to a negative surface dipole and shifts band edges in the opposite direction by a similar amount. The presence of the full interface with liquid water preserves this overall picture, and consequently leads to similar trends in the band edge positions. At the same time, the effect of the full interface is also quantitatively important, contributing on the order of 0.5–0.8 eV to the band

when the effect of liquid water is taken into account, as represented by the difference between the red and magenta solid lines in Fig. 3 for each materials. The pristine GaP(001) and InP(001) surfaces in water can be compared with the corresponding fully hydroxylated surfaces in water with an oxygen coverage of θ = 1.5 ML (solid green lines) to assess the effect of surface oxidation. In the latter case, the band edges are again shifted with respect to the pristine surface in vacuum, but in the opposite direction. Note that the hydroxylated surfaces in water can be interpreted as the limit of full surface water dissociation, with all available sites covered with hydroxyl; on the other hand, the pristine surface can be interpreted as the limit of exclusively molecularly adsorbed water that remains intact. 36 We conclude that three factors are important for determining the band edge positions of a photoelectrode in liquid water: surface polarity, the interface with water, and the extent of surface oxidation. However, these effects are closely interconnected, and together form the basis of the relationship between surface composition and electronic structure. While our simulations include full solid/liquid interfaces, to gain additional insights into the effects of interfaces with water on the calculated band edges, we also considered GaP(001) and

ACS Paragon Plus Environment

6

Page 7 of 17 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

0.0 ML and 1.5 ML represented by the theoretical models. Our conclusion is also consistent with other XPS experiments 61,62 that find a mixture of non-dissociated and dissociated adsorbed water on the GaP(001) and InP(001) surfaces in the presence of water vapor. The G0 W0 calculations also indicates significant differences between GaP and InP in the hydrogen-bond network at their interfaces with liquid water. As shown in Fig. 3, when compared to the corresponding results obtained for the full interfaces with bulk water (solid green and magenta lines), band edges of the δ(2 × 4) surface with a monolayer of adsorbed water molecules (dashed magenta lines), as well as those of the hydroxylated surface (dashed green lines) consistently show that the deviations are different for GaP and InP. Specifically, these deviations are notably larger for GaP (by ∼0.3 eV). Such an observation suggests that GaP band edges are affected more strongly by interfacial water, and conversely, that the perturbation introduced by the semiconductor surfaces on interfacial water molecules is stronger for GaP. The above conclusion reflects key aspects of our previous FPMD simulations on GaP/water and InP/water interfaces. 16 In particular, we reported that water at the interface with GaP(001) shows significant density perturbations from the bulk limit, whereas the effect of InP(001) on water is much weaker. This difference was attributed to higher covalent character in the Ga–O surface bonds, which results in a more rigid hydrogen-bond network near the interface. 16 The relatively larger impact of GaP on the properties of interfacial water implies a stronger direct interaction with the semiconductor surface, which is reflected in the band edges of Fig. 3. Interfacial Chemistry from X-ray Photoelectron Spectroscopy. Although the band edge analysis provides insight into the overall extent of surface oxidation, an additional probe is required for more specific identification of oxygen-bearing surface chemical species. In previous studies, we found that surface oxygen atoms can exist in several coordination environments, which were found to interchange dynam-

edge shift. This is likely due to the partial mitigation of the surface dipole due to hydrogen bonding with the liquid water. Overall, our analysis highlights the importance of dipolar effects on electronic band edges of the photoelectrodes in water that stem from complex chemistry involving adsorbed and dissociated water molecules at the solid/liquid interface. We note that the importance of surface dipoles on the electronic band edges of photoelectrodes and, by extension, performance has also been discussed in the literature. 60–63 For example, it has been shown that the surface dipole leads to a downward movement of the GaP band edges when in contact with water vapor (by about 0.17 eV in water pressure of 10−5 mbar 61 ). This has important implications for performance, since it reduces the offset between the CBM and the water reduction potential, thereby benefiting the use of GaP as a photocathode for hydrogen evolution. Interfacial Chemistry from Band Edge Analysis. Having demonstrated the surface oxidation-band edge correlation, the G0 W0 calculations can be combined with experimental band edge determinations to predict chemistry of GaP/water and InP/water interfaces. In particular, band edges of GaP and InP in water is known to vary with respect to the pH following the Nernst equation, 58,64,65 resulting in a wide range of experimental VBMs, i.e., 4.72– 5.52 eV and 4.88–5.68 eV below vacuum in the pH range of 0–14 for GaP and InP surfaces, respectively, as represented by the filled boxes in Fig. 3. When compared to the experimental band edges at neutral pH that corresponds to the simulation condition, our calculations show that theoretical VBMs of δ(2×4) and hydroxylated surfaces straddle the experimental data, indicating that the realistic surfaces are likely to be a mixed of the two model surfaces that represent two extreme limits of oxidation extent. As a consequence, this observation suggests that non-dissociated and dissociated adsorbed water coexist on the GaP(001) and InP(001) surfaces, and that water dissociation leads to the formation of surface oxygen and hydroxyl groups with an overall surface oxygen coverage between the two extreme limits of

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry Letters Binding energy (eV)

InP(001)/H2O

534

533

H2O*

GaP(001)/H2O

Intensity (Arb. units)

Intensity (Arb. units)

Binding energy (eV)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 17

532

531

M-O*H-M M-O*H

530

529

534

533

532

H2O* M-O*H-M

M-O*-M

M-O*H

531

530

529

M-O*-M

Figure 4: The O 1s photoemission spectra (red circles) recorded for InP(001) (left) and GaP(001) (right) at a water vapor pressure of 5 mbar. Peak assignments (vertical lines) for different surface species were derived from first-principles: surface oxygen bridge bonds (M –O*–M , magenta), atop hydroxyl (M –O*H, blue), bridge hydroxyl (M –O*H–M , green), and molecularly adsorbed water (H2 O*, maroon), where M ≡ In or Ga, and ∗ indicates the oxygen atom from which the photoelectrons are emitted. The experimental spectra (red circles) were fitted with four components using Gaussian/Lorentzian hybrid functions and calculated binding energies (see Suporting Information for details of the fitting procedure). sion that the GaP(001)/water interface shows more rigid hydrogen-bond network than its InP counterpart, which was initially deduced from FPMD simulations and the combined theoryexperiment band edge analysis. This apparent discrepancy turns out to be easily understood once the O 1s spectra are properly interpreted with aid from first-principles theory, highlighting the strength of an integrated experimenttheory approach to studying interfacial chemistry. To assist the interpretation of the experimental photoemission spectra, first-principles calculations of oxygen binding energies (BEs) were carried out for the oxidized GaP(001) and InP(001) surfaces in water, particularly for the hydroxylated interfaces as they provide a more complete set of surface oxide configurations. In addition, to elucidate oxygen speciation, the individual oxygen atoms contributing to the computed BEs were categorized according to the same coordination scheme described before, consisting of bridge oxygens, bridge and atop hydroxyls, and adsorbed water molecules. We emphasize that our computational procedure is applicable to the full semiconductor/water interface, thereby allowing for direct comparison between the experimental spectra and the calculated BEs (see Supporting Information). Ac-

ically in FPMD simulations: M –O–M , M –OH– M , M –OH, and H2 O. 9,16 This suggests a more complex surface speciation at the semiconductor/water interface, which can be probed more directly using APXPS experiments. Specifically, we recorded the O1s spectra of the semiconductor surfaces at room temperature, with the water vapor pressure adjusted from ultrahigh vacuum (UHV) to 5 mbar. The reproducibility of the XPS spectra were carefully checked; specifically, each experimental conditions, we repeated XPS measurements for at least three times. The O1s signal measured for InP and GaP at the highest pressure is presented in Fig. 4. Clear differences are visible in the shape of the photoemission spectra obtained for InP(001)/water and GaP(001)/water interfaces. Most notably, the spectrum of GaP is much broader and exhibits more complex features within its overall shape. Such an observation is consistent with previous XPS studies of GaP(001) and InP(001) surfaces at the interface with water, 61,62 and implies fundamentally different chemistry occurring at the two interfaces, as we have suggested. However, at first glance, the relative breadth of the O 1s spectrum for GaP compared to InP appears counterintuitive to our earlier conclu-

ACS Paragon Plus Environment

8

Page 9 of 17 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

ing task in the case of complex systems. The breadth of the overall photoemission spectrum of GaP(001) in Fig. 4 is therefore not attributable to thermal broadening alone and does not necessarily indicate a more fluid interchange between surface species. In addition, we verified that the shape is not merely an artifact of the band bending effects, as examination of the spectral shift of GaP and InP before and after the introduction of water vapor indicates that the variations in the band bending are relatively small (∼0.1 eV). Rather, it reflects the much more significant chemical differences between oxygen species on GaP, which lead to a larger collective spread in the corresponding BEs. In particular, a large separation between photoemission peaks obtained for oxygen atoms belonging to atop and bridge hydroxyls in the GaP spectrum implies that M – O bond strengths in M –OH and M –OH–M are significantly different for the GaP surface, in contrast to InP where these bond strengths are much more similar. This indicates that the interchange of atop and bridge hydroxyls is likely to be more kinetically limited for GaP compared to InP. Accordingly, the core-binding energy calculations also support the interpretation of a more rigid, well-defined network structure at the GaP/water interface predicted by the G0 W0 calculations. Finally, it is worth discussing possible effects of surface structure orientation on the binding modes of water molecules. Our current and previous studies 70 of the GaP(001) and GaP(111) surfaces indicate that the binding mode of water molecules are not sensitive to surface orientation. For example, we find that the shape of the O 1s spectra are similar, showing a consistent position of the main peak at 531.7±0.1 eV. However, this similarity may stem from the fact that these are both polar surfaces; accordingly, further investigation of other surface orientations with different characters, e.g., non-polar surface, is needed in order to provide a solid conclusion on the surface orientation dependence of water behavior at the interface with III-V semiconductors. Conclusions and Outlook. The combination of experimental measurements with high-level

cordingly, this method provides standards for interpreting photoemission spectra that are far more representative of the in situ environment under which the measurements are performed, while simultaneously accounting for the configurational complexity of the dynamical interfaces. The computed BEs of photoelectrons emitted from oxygen atoms belonging to the bridge bonds (M –O*–M ), atop (M –O*H) and bridge (M –O*H–M ) hydroxyls, and molecularly adsorbed water (H2 O*) are presented by blue vertical lines in Fig. 4 for both InP/water and GaP/water interfaces, together with fitted experimental spectra using the calculated BEs. Our analysis shows that the contribution to the spectra at lower binding energies (530-532 eV) are attributed to water dissociation products, whereas the broad peak around 533 eV is attributed to molecular water molecules on the semiconductor surfaces. The contribution to this peak stems mostly from water molecules in the first monolayer, but may also originate from those belonging to water clusters on the surface. We find that the photoelectron BEs follow the order H2 O* > O*H > M –O*–M for both materials; such an ordering is found to be in agreement with the assignment obtained in previous in-situ XPS studies of water dissociative adsorption onto metal surfaces 19,66,67 and Ga-based surface compounds. 68–70 Interestingly, the calculated BEs patterns obtained for oxygen configurations on GaP and InP are significantly different, particularly in the relative BEs rising from different hydroxyl groups. Specifically, it is shown that the photoelectron BEs obtained for oxygen belonging to the atop and bridge hydroxyls are rather close in energy for InP (less than 0.2 eV), together contributing to the main experimental peak. In contrast, these two peaks are well separated by almost 0.8 eV for GaP, leading to a much wider BEs pattern. This is consistent with the experimental observation, showing that the GaP surface yields a broader spectrum with more complex features. Accordingly, BEs patterns obtained from calculations can greatly assist in interpretation and assignment of photoemission spectra obtained from experiments, which is a challeng-

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

first-principles methods can provide valuable insight into the role of surface oxidation on InP and GaP photoabsorbers. For instance, our results clearly indicate that surface polarity, the interface with water, and the extent of surface oxidation are all important factors in determining band edges. Theoretical calculations of XPS reference spectra and band edges can be combined with analogous experiments to establish these factors in the context of the specific surface chemistry. When combined with APXPS and band edge measurements, the theoretical calculations can be used to deduce that dissociated and non-dissociated adsorbed water coexist on the GaP and InP surfaces, altering surface dipolar properties. Notable differences between GaP and InP are also observed, with the band edges of GaP exhibiting much higher sensitivity to the presence of water. Differences between the two surfaces are also reflected in the O 1s XPS spectra in the presence of water, which are broader for GaP than for InP. Using simulations of model interfaces, we trace the origin of this behavior to the higher rigidity of the water hydrogen-bond network at the GaP/water interface. The importance of the intrinsic surface polarity in the determination of the photoelectrode band edges demonstrated here indicates that surface treatment is critical for improving PEC efficiency. In this regard, careful engineering of the photoelectrode surface may result in a well-adapted band alignment with their native oxide layers and the water oxidation/reduction potentials, thereby benefiting water-splitting performance. Beyond these specific insights into GaP and InP oxidation, the recipe illustrated here— using experimental insights to refine models, and in turn using modeling data to inform experimental interpretation—can be viewed as a more general roadmap for investigations of wider classes of interfaces. High-level firstprinciples methods are employed to understand the specific relationship between surface chemistry and key electronic (e.g., band edge) and spectroscopic (e.g., XPS) features at the interface, based on model systems that represent features present in the real system. The reference electronic structures and spectra of these

Page 10 of 17

models are then compared with in situ experimental probes of these same quantities at the complex interface, with the similarities and differences used to infer how the surface chemistry evolves under operation. For example, APXPS experiments can be used to probe not only the chemistry but also the electronic properties of solid/water interfaces, as demonstrated in recent studies, 71,72 thereby providing a direct comparison with theoretical calculations. This combined theory-experiment approach offers a promising way for investigating surface morphology and chemical composition, which is important to guide future strategies for improving device performance through interface engineering. As first-principles methods advance to access larger and longer simulations, and as resolution and sensitivity of X-ray probes simultaneously continue to improve, these methods will become capable of probing interfaces with increasing complexity. In this event, the combination of computational and experimental approaches will certainly play an increasingly important role in interfacial science. Whereas historically, solid/liquid simulations have often focused on simplified models using singlemolecule or monolayer interface descriptions, computational developments are also likely to make explicit interface simulations more routine. 73 This can offer a qualitatively or quantitatively different view of surface chemistry and electronic properties, as we will illustrate in the example covered within this perspective. On the experimental side, equally important is the role of in situ probes, which are beginning to provide a far more complete description of reactive interfaces than can be accessed using conventional surface-science probes under vacuum conditions. 18 Nevertheless, the approach we outline here is not without limitations when considering actual operating conditions, motivating further developments in the community. For instance, additional physical factors should be considered to improve the realism of the simulation and characterization environments. For instance, in the case of electrochemical interfaces, voltage and photobias play a critical role under oper-

ACS Paragon Plus Environment

10

Page 11 of 17 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

include interface structure, 21,23,24,29,30,80–83 mass transport, 16 and key reaction pathways, 24,84–86 all of which are needed to develop comprehensive and rational approaches for controlling (photo)electrochemical processes through interface engineering.

ating conditions, 72 yet these are neglected in the roadmap we introduce here. On the simulation side, developments are underway to address this shortcoming, 74–76 but the methods are currently in their infancy. In addition, our investigation considered only the role of liquid water, yet actual operating electrolytes contain salt ions that can have a significant effect on performance. 77 Other shortcomings may be unique to either simulations or to experiments. For instance, first-principles simulations are typically limited to the immediate interfacial region operating under equilibrium conditions. However, many devices have behaviors that are connected to operation under non-equilibrium conditions, where dynamical changes in factors such as local concentration or activity must be considered. This requires integrating first-principles chemical simulations with broader mesoscale or device-level simulations. On the experimental side, APXPS does not provide any notable spatial resolution, and is also limited in temporal resolution. In this regard, it may be helpful to combine this technique with other advanced probes that can offer higher spatiotemporal fidelity. 78,79 The potential importance of these capabilities becomes evident when considering the possible role of local singularities and dynamical processes in promoting chemical processes. Although such features can be readily introduced in models, validation becomes difficult without appropriate experimental probes. For this reason, the roadmap introduced here avoids making specific connections between the surface chemistry and interfacial catalytic processes, focusing instead on properties such as band edges, which can be obtained via temporal and spatial averaging (i.e., equilibrium properties). Further developments in techniques with high spatiotemporal resolution could strengthen the connection between interface composition and chemical kinetics. These are a high priority for future research efforts as such approaches will offer ways to validate, and if necessary revise, the atomistic insights obtained from first-principles simulations of electrochemical interfaces that are particularly relevant for (photo)catalytic reactions. Examples

Acknowledgement Part of this work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344. T.A.P. acknowledges support from the Lawrence Fellowship. T.O and B.C.W. are supported by the U.S. Deparment of Energy, Office of Energy Efficiency and Renewable Energy, Fuel Cell Technologies Office. X.Z. and S.P. acknowledge the U.S. Department of Energy Office of Science, Office of Basic Energy Sciences under Award Number DE-FC02-04ER15533 (NDRL no: 5152). Computational support was from the LLNL Grand Challenge Program. We thank Eric Schwegler and Giulia Galli for fruitful discussions.

Biographies Tuan Anh Pham is a Lawrence Postdoc Fellow in the Quantum Simulations Group and Materials Science Division at the Lawrence Livermore National Laboratory (LLNL). He received a Ph.D. in Chemistry at the University of California, Davis in 2014. He is also a recipient of the 2011 LLNL Graduate Scholar, 2014 UC Davis Outstanding Chemistry Dissertation Award, and 2017 PCTC Postdoctoral Fellow Award. His research interest includes the development and applications of first-principles techniques to validate, understand and predict materials properties for energy and enviromental applications. Xueqiang Zhang is currently a postdoctoral research associate in the Department of Chemistry at University of Illinois at UrbanaChampaign. He received a Ph.D in Physical Chemistry from the University of Notre Dame

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

in 2016. His research interests include the design and mechanistic understanding of heterogeneous catalysis and photocatalysis, and the study of fundamental surface chemistry at vapor/solid interfaces using in-situ/operando spectroscopic techniques.

Page 12 of 17

phasis on energy-related materials, using advanced in situ techniques. Tadashi Ogitsu is a staff scientist in the Quantum Simulations Group and Materials Science Division at LLNL. He received his Ph.D. in Materials Science in 1994 from the University of Tsukuba in Japan. His research interests include high pressure physics, ultrafast non-equilibrium physics, and materials and interfaces for energy applications such as carbon aerogels for supercapacitor applications and photoelectrochemical hydrogen production. He is the LLNL Laboratory Lead for the HydroGEN Advanced Water Splitting Materials consortium, funded by the Fuel Cell Technologies Office within the DOE Office of Energy Efficiency and Renewable Energy.

Brandon C. Wood is a staff scientist in the Quantum Simulations Group and Materials Science Division at LLNL. He received a Ph.D. in Materials Science and Engineering from MIT in 2007. His primary research interests lie in the application of high-performance computing to ab initio and mesoscale simulations of complex interfaces for energy storage and conversion. He is the LLNL Laboratory Lead for the Hydrogen Storage Materials—Advanced Research Consortium (HyMARC), funded by the Fuel Cell Technologies Office within the DOE Office of Energy Efficiency and Renewable Energy.

References

David Prendergast is a staff scientist at The Molecular Foundry, Lawrence Berkeley National Laboratory. He has been the Facility Director for the Theory of Nanostructured Materials Facility at The Molecular Foundry since 2012. He obtained a Ph.D. in Physics from University College Cork, Ireland (2002). His current research focuses on developments and applications of first-principles electronic structure theory in the context of energy-relevant phenomena in materials and chemistry. He has particular expertise in the simulation and interpretation of synchrotron X-ray spectroscopy measurements as a means of connecting characterization of complex interfacial systems to atomistic or molecular scale structural and dynamical models.

(1) McKone, J. R.; Lewis, N. S.; Gray, H. B. Will solar-driven water-splitting devices see the light of day? Chem. Mat. 2013, 26, 407–414. (2) Lewis, N. S.; Nocera, D. G. Powering the planet: Chemical challenges in solar energy utilization. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 15729–15735. (3) Lewis, N. S. Research opportunities to advance solar energy utilization. Science 2016, 351, aad1920. (4) Khaselev, O.; Turner, J. A. A monolithic photovoltaic-photoelectrochemical device for hydrogen production via water splitting. Science 1998, 280, 425–427.

Sylwia Ptasinska is an Associate Professor of Physics at the University of Notre Dame. She received a Ph.D. in Physics at Leopold-Franzens-University, Innsbruck, Austria in 2004. She leads an interdisciplinary research group at the Notre Dame Radiation Laboratory which focuses on fundamental studies driven by applications in areas such as energy, medicine, and industry. One of her experimental efforts is to develop realistic physicochemical mechanisms occurring at the interphase of water/semiconductor surfaces, with specific em-

(5) Hu, S.; Shaner, M. R.; Beardslee, J. A.; Lichterman, M.; Brunschwig, B. S.; Lewis, N. S. Amorphous TiO2 coatings stabilize Si, GaAs, and GaP photoanodes for efficient water oxidation. Science 2014, 344, 1005–1009. (6) Ager, J. W.; Shaner, M. R.; Walczak, K. A.; Sharp, I. D.; Ardo, S. Experimental demonstrations of spontaneous, solar-driven photoelectrochemical water

ACS Paragon Plus Environment

12

Page 13 of 17 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

splitting. Energy Environ. Sci. 2015, 8, 2811–2824.

O 2 (g) in 1.0 M KOH (aq) using conformal coatings of amorphous TiO 2. Energy Environ. Sci. 2015, 8, 203–207.

(7) May, M. M.; Lewerenz, H.-J.; Lackner, D.; Dimroth, F.; Hannappel, T. Efficient direct solar-to-hydrogen conversion by in situ interface transformation of a tandem structure. Nat. Com. 2015, 6, 8286.

(15) Kaiser, B.; Fertig, D.; Ziegler, J.; Klett, J.; Hoch, S.; Jaegermann, W. Solar Hydrogen Generation with Wide-BandGap Semiconductors: GaP (100) Photoelectrodes and Surface Modification. ChemPhysChem 2012, 13, 3053–3060.

(8) Lichterman, M. F.; Sun, K.; Hu, S.; Zhou, X.; McDowell, M. T.; Shaner, M. R.; Richter, M. H.; Crumlin, E. J.; Carim, A. I.; Saadi, F. H. et al. Protection of inorganic semiconductors for sustained, efficient photoelectrochemical water oxidation. Catalysis Today 2016, 262, 11–23.

(16) Wood, B. C.; Schwegler, E.; Choi, W. I.; Ogitsu, T. Hydrogen-Bond Dynamics of Water at the Interface with InP/GaP (001) and the Implications for Photoelectrochemistry. J. Am. Chem. Soc. 2013, 135, 15774–15783.

(9) Wood, B. C.; Ogitsu, T.; Schwegler, E. Local structural models of complex oxygenand hydroxyl-rich GaP/InP (001) surfaces. J. Chem. Phys. 2012, 136, 064705.

(17) Smith, W. A.; Sharp, I. D.; Strandwitz, N. C.; Bisquert, J. Interfacial bandedge energetics for solar fuels production. Energy Environ. Sci. 2015, 8, 2851–2862.

(10) Szklarczyk, M.; Bockris, J. O. Photoelectrochemical evolution of hydrogen on pindium phosphide. J. Phys. Chem. 1984, 88, 5241–5245.

(18) Zaera, F. Probing liquid/solid interfaces at the molecular level. Chem. Rev. 2012, 112, 2920–2986. (19) Salmeron, M.; Schl¨ogl, R. Ambient pressure photoelectron spectroscopy: A new tool for surface science and nanotechnology. Surf. Sci. Rep. 2008, 63, 169–199.

(11) Esposito, D. V.; Levin, I.; Moffat, T. P.; Talin, A. A. H2 evolution at Si-based metal–insulator–semiconductor photoelectrodes enhanced by inversion channel charge collection and H spillover. Nat. Mater. 2013, 12, 562–568.

(20) Starr, D.; Liu, Z.; H¨avecker, M.; KnopGericke, A.; Bluhm, H. Investigation of solid/vapor interfaces using ambient pressure X-ray photoelectron spectroscopy. Chem. Soc. Rev. 2013, 42, 5833–5857.

(12) Labrador, N. Y.; Li, X.; Liu, Y.; Tan, H.; Wang, R.; Koberstein, J. T.; Moffat, T. P.; Esposito, D. V. Enhanced Performance of Si MIS Photocathodes Containing Oxide-Coated Nanoparticle Electrocatalysts. Nano Letters 2016,

(21) Gayathri, N.; Izvekov, S.; Voth, G. A. Ab initio molecular dynamics simulation of the H/InP (100)–water interface. J. Chem. Phys. 2002, 117, 872–884.

(13) Hu, S.; Lewis, N. S.; Ager, J. W.; Yang, J.; McKone, J. R.; Strandwitz, N. C. ThinFilm Materials for the Protection of Semiconducting Photoelectrodes in Solar-Fuel Generators. J. Phys. Chem. C 2015, 119, 24201–24228.

(22) Cicero, G.; Grossman, J. C.; Catellani, A.; Galli, G. Water at a hydrophilic solid surface probed by ab initio molecular dynamics: inhomogeneous thin layers of dense fluid. J. Am. Chem. Soc 2005, 127, 6830– 6835.

(14) Shaner, M. R.; Hu, S.; Sun, K.; Lewis, N. S. Stabilization of Si microwire arrays for solar-driven H 2 O oxidation to

(23) Taylor, C. D.; Wasileski, S. A.; Filhol, J.S.; Neurock, M. First principles reaction

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

modeling of the electrochemical interface: Consideration and calculation of a tunable surface potential from atomic and electronic structure. Phys. Rev. B 2006, 73, 165402.

Page 14 of 17

semiconductor interfaces. Phys. Rev. Lett. 2014, 113, 176802. (33) Pham, T. A.; Lee, D.; Schwegler, E.; Galli, G. Interfacial Effects on the Band Edges of Functionalized Si Surfaces in Liquid Water. J. Am. Chem. Soc. 2014, 136, 17071–17077.

(24) Otani, M.; Hamada, I.; Sugino, O.; Morikawa, Y.; Okamoto, Y.; Ikeshoji, T. Structure of the water/platinum interface—-a first principles simulation under bias potential. Phys. Chem. Chem. Phys. 2008, 10, 3609–3612.

(34) De Angelis, F.; Di Valentin, C.; Fantacci, S.; Vittadini, A.; Selloni, A. Theoretical studies on anatase and less common TiO2 phases: bulk, surfaces, and nanomaterials. Chem. Rev. 2014, 114, 9708–9753.

(25) Schnur, S.; Groß, A. Properties of metal– water interfaces studied from first principles. New J. Phys. 2009, 11, 125003.

(35) Velasco-Velez, J.-J.; Wu, C. H.; Pascal, T. A.; Wan, L. F.; Guo, J.; Prendergast, D.; Salmeron, M. The structure of interfacial water on gold electrodes studied by x-ray absorption spectroscopy. Science 2014, 1259437.

(26) Cheng, J.; Sprik, M. Aligning electronic energy levels at the TiO 2/H 2 O interface. Phys. Rev. B 2010, 82, 081406. (27) Carrasco, J.; Hodgson, A.; Michaelides, A. A molecular perspective of water at metal interfaces. Nat. Mater. 2012, 11, 667–674.

(36) Wood, B. C.; Schwegler, E.; Choi, W. I.; Ogitsu, T. Surface chemistry of GaP (001) and InP (001) in contact with water. J. Phys. Chem. C 2014, 118, 1062–1070.

(28) Wang, J.; Pedroza, L. S.; Poissier, A.; Fern´andez-Serra, M. Water dissociation at the GaN (1010) surface: structure, dynamics and surface acidity. J. Phys. Chem. C 2012, 116, 14382–14389.

(37) Nielsen, M.; Bj¨orketun, M. E.; Hansen, M. H.; Rossmeisl, J. Towards first principles modeling of electrochemical electrode–electrolyte interfaces. Sur. Sci. 2015, 631, 2–7.

(29) Ando, Y.; Gohda, Y.; Tsuneyuki, S. Ab initio molecular dynamics study of the Helmholtz layer formed on solid–liquid interfaces and its capacitance. Chem. Phys. Lett. 2013, 556, 9–12.

(38) Pham, T. A.; Ping, Y.; Galli, G. Modelling heterogeneous interfaces for solar water splitting. Nat. Mater. 2017, 16, 401–408. (39) Car, R.; Parrinello, M. Unified approach for molecular dynamics and densityfunctional theory. Phys. Rev. Lett. 1985, 55, 2471.

(30) Farnesi Camellone, M.; Marx, D. On the impact of solvation on a Au/TiO2 nanocatalyst in contact with water. J. Phys. Chem. Lett. 2013, 4, 514–518.

(40) Hedin, L. New method for calculating the one-particle Green’s function with application to the electron-gas problem. Phys. Rev. 1965, 139, A796.

(31) Huang, P.; Pham, T. A.; Galli, G.; Schwegler, E. Alumina (0001)/water interface: structural properties and infrared spectra from first-principles molecular dynamics simulations. J. Phys. Chem. C 2014, 118, 8944–8951.

(41) Nguyen, H.-V.; Pham, T. A.; Rocca, D.; Galli, G. Improving accuracy and efficiency of calculations of photoemission spectra within the many-body perturbation theory. Phys. Rev. B 2012, 85, 081101.

(32) Kharche, N.; Muckerman, J. T.; Hybertsen, M. S. First-principles approach to calculating energy level alignment at aqueous

ACS Paragon Plus Environment

14

Page 15 of 17 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(42) Pham, T. A.; Nguyen, H.-V.; Rocca, D.; Galli, G. GW calculations using the spectral decomposition of the dielectric matrix: Verification, validation, and comparison of methods. Phys. Rev. B 2013, 87, 155148.

(51) Wood, B. C.; Ogitsu, T.; Schwegler, E. Ab initio modeling of water–semiconductor interfaces for photocatalytic water splitting: role of surface oxygen and hydroxyl. Journal of Photonics for Energy 2011, 1, 016002–016002.

(43) Pham, T. A.; Li, T.; Nguyen, H.-V.; Shankar, S.; Gygi, F.; Galli, G. Band offsets and dielectric properties of the amorphous Si3N4/Si (100) interface: A firstprinciples study. Appl. Phys. Lett. 2013, 102, 241603.

(52) Van Laar, J.; Huijser, A.; Van Rooy, T. Electronic surface properties of Ga and In containing III–V compounds. J. Vac. Sci. Technol. 1977, 14, 894–898. (53) Weiss, W.; Wiemh¨ofer, H.; G¨opel, W. Role of ionic defects at semiconductorinsulator interfaces: Spectroscopic results on CaF 2/InP (001) and SrF 2/InP (001) structures. Phys. Rev. B 1992, 45, 8478.

(44) Pham, T. A.; Zhang, C.; Schwegler, E.; Galli, G. Probing the electronic structure of liquid water with many-body perturbation theory. Phys. Rev. B 2014, 89, 060202.

(54) Cheng, J.; Sprik, M. Alignment of electronic energy levels at electrochemical interfaces. Phys. Chem. Chem. Phys. 2012, 14, 11245–11267.

(45) Opalka, D.; Pham, T. A.; Sprik, M.; Galli, G. The ionization potential of aqueous hydroxide computed using manybody perturbation theory. J. Chem. Phys. 2014, 141, 034501.

(55) Hybertsen, M. S.; Louie, S. G. FirstPrinciples Theory of Quasiparticles: Calculation of Band Gaps in Semiconductors and Insulators. Phys. Rev. Lett. 1985, 55, 1418–1421.

(46) Opalka, D.; Pham, T. A.; Galli, G.; Sprik, M. Electronic Energy Levels and Band Alignment for Aqueous Phenol and Phenolate from First Principles. J. Phys. Chem. B 2015, 119, 9651–9660.

(56) Jiang, H.; Shen, Y.-C. Ionization potentials of semiconductors from firstprinciples. J. Chem. Phys. 2013, 139, 164114.

(47) Govoni, M.; Galli, G. Large Scale GW Calculations. J. Chem. Theory. Comput. 2015, 11, 2680–2696.

(57) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865.

(48) Gaiduk, A. P.; Govoni, M.; Seidel, R.; Skone, J. H.; Winter, B.; Galli, G. Photoelectron spectra of aqueous solutions from first principles. J. Am. Chem. Soc 2016, 138, 6912–6915.

(58) Wrighton, M. S.; Bolts, J. M.; Bocarsly, A. B.; Palazzotto, M. C.; Walton, E. G. Stabilization of n-type semiconductors to photoanodic dissolution: II–VI and III–V compound semiconductors and recent results for n-type silicon. J. Vac. Sci. Technol. 1978, 15, 1429–1435.

(49) Pham, T. A.; Govoni, M.; Seidel, R.; Bradforth, S. E.; Schwegler, E.; Galli, G. Electronic structure of aqueous solutions: Bridging the gap between theory and experiments. Sci. Adv. 2017, 3, e1603210.

(59) Van Wezemael, A.-M.; Laflere, W.; Cardon, F.; Gomes, W. Study of the schottky barrier at the n-and p-type indium phosphide electrode/electrolyte interface. J. Electroanal. Chem. Interfacial. Electrochem. 1978, 87, 105–109.

(50) Schmidt, W. III-V compound semiconductor (001) surfaces. Appl. Phys. A 2002, 75, 89–99.

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(60) Henrion, O.; Klein, A.; Jaegermann, W. Water adsorption on UHV cleaved InP (110) surfaces. Surf. Sci 2000, 457, L337– L341.

Page 16 of 17

hydroxyl chemistry on metal surfaces: the Cu (110) case. J. Phys. Chem. C 2007, 111, 14493–14499. (68) Zhang, X.; Ptasinska, S. Dissociative adsorption of water on an H2O/GaAs (100) interface: In situ near-ambient pressure XPS studies. J. Phys. Chem. C 2014, 118, 4259–4266.

(61) May, M. M.; Supplie, O.; H¨ohn, C.; van de Krol, R.; Lewerenz, H.-J.; Hannappel, T. The interface of GaP (100) and H2O studied by photoemission and reflection anisotropy spectroscopy. New J. Phys. 2013, 15, 103003.

(69) Zhang, X.; Lamere, E.; Liu, X.; Furdyna, J. K.; Ptasinska, S. Interface chemistry of H 2 O on GaAs nanowires probed by near ambient pressure X-ray photoelectron spectroscopy. Chem. Phys. Lett. 2014, 605, 51–55.

(62) May, M. M.; Lewerenz, H.-J.; Hannappel, T. Optical in situ Study of InP (100) Surface Chemistry: Dissociative Adsorption of Water and Oxygen. J. Phys. Chem. C 2014, 118, 19032–19041.

(70) Zhang, X.; Ptasinska, S. Distinct and dramatic water dissociation on GaP (111) tracked by near-ambient pressure X-ray photoelectron spectroscopy. Phys. Chem. Chem. Phys. 2015, 17, 3909–3918.

(63) MacLeod, B. A.; Steirer, K. X.; Young, J. L.; Koldemir, U.; Sellinger, A.; Turner, J. A.; Deutsch, T. G.; Olson, D. C. Phosphonic Acid Modification of GaInP2 Photocathodes Towards Unbiased Photoelectrochemical Water splitting. Acs. Appl. Mater. Inter. 2015, 7, 11346–11350.

(71) Zhang, X.; Ptasinska, S. Electronic and chemical structure of the H2O/GaN (0001) interface under ambient conditions. Sci. Rep. 2016, 6, 24848.

(64) Chandra, N.; Wheeler, B. L.; Bard, A. J. Semiconductor electrodes. 59. Photocurrent efficiencies at p-indium phosphide electrodes in aqueous solutions. J. Phys. Chem. 1985, 89, 5037–5040.

(72) Lichterman, M. F.; Hu, S.; Richter, M. H.; Crumlin, E. J.; Axnanda, S.; Favaro, M.; Drisdell, W.; Hussain, Z.; Mayer, T.; Brunschwig, B. S. et al. Direct observation of the energetics at a semiconductor/liquid junction by operando X-ray photoelectron spectroscopy. Energy Environ. Sci. 2015, 8, 2409–2416.

(65) Iranzo-Marin, F.; Debiemme-Chouvy, C.; Herlem, M.; Sculfort, J.-L.; Etcheberry, A. Electrochemical techniques for the elucidation of the interface structure of the nInP/aqueous electrolyte junction. J. Electroanal. Chem 1994, 365, 283–287.

(73) Bjorneholm, O.; Hansen, M. H.; Hodgson, A.; Liu, L.-M.; Limmer, D. T.; Michaelides, A.; Pedevilla, P.; Rossmeisl, J.; Shen, H.; Tocci, G. et al. Water at interfaces. Chem. Rev. 2016, 116, 7698–7726.

(66) Andersson, K.; Ketteler, G.; Bluhm, H.; Yamamoto, S.; Ogasawara, H.; Pettersson, L. G.; Salmeron, M.; Nilsson, A. Autocatalytic water dissociation on Cu (110) at near ambient conditions. J. Am. Chem. Soc 2008, 130, 2793–2797.

(74) Otani, M.; Sugino, O. First-principles calculations of charged surfaces and interfaces: A plane-wave nonrepeated slab approach. Phys. Rev. B 2006, 73, 115407.

(67) Andersson, K.; Ketteler, G.; Bluhm, H.; Yamamoto, S.; Ogasawara, H.; Pettersson, L. G.; Salmeron, M.; Nilsson, A. Bridging the pressure gap in water and

(75) Bonnet, N.; Morishita, T.; Sugino, O.; Otani, M. First-Principles Molecular Dy-

ACS Paragon Plus Environment

16

Page 17 of 17 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(85) Desai, S. K.; Neurock, M. First-principles study of the role of solvent in the dissociation of water over a Pt-Ru alloy. Phys. Rev. B 2003, 68, 075420.

namics at a Constant Electrode Potential. Phys. Rev. Lett. 2012, 109, 266101. (76) Radin, M. D.; Ogitsu, T.; Biener, J.; Otani, M.; Wood, B. C. Capacitive charge storage at an electrified interface investigated via direct first-principles simulations. Phys. Rev. B 2015, 91, 125415.

(86) Tocci, G.; Michaelides, A. Solvent-induced proton hopping at a water–oxide interface. J. Phys. Chem. Lett. 2014, 5, 474–480.

(77) Hill, J. C.; Choi, K.-S. Effect of Electrolytes on the Selectivity and Stability of n-type WO3 Photoelectrodes for Use in Solar Water Oxidation. J. Phys. Chem. C 2012, 116, 7612–7620. (78) Papp, C.; Steinr¨ uck, H.-P. In situ high-resolution X-ray photoelectron spectroscopy–Fundamental insights in surface reactions. Surf. Sci. Rep. 2013, 68, 446–487. (79) Esposito, D. V.; Baxter, J. B.; John, J.; Lewis, N. S.; Moffat, T. P.; Ogitsu, T.; O’Neil, G. D.; Pham, T. A.; Talin, A. A.; Velazquez, J. M. et al. Methods of photoelectrode characterization with high spatial and temporal resolution. Energy Environ. Sci. 2015, 8, 2863–2885. (80) Aschauer, U.; Selloni, A. Structure of the rutile TiO 2 (011) surface in an aqueous environment. Phys. Rev. Lett. 2011, 106, 166102. (81) Nadler, R.; Sanz, J. F. First-principles molecular dynamics simulations of the H2O/Cu (111) interface. J. Mol. Model. 2012, 18, 2433–2442. (82) Mosconi, E.; Selloni, A.; De Angelis, F. Solvent effects on the adsorption geometry and electronic structure of dye-sensitized TiO2: a first-principles investigation. J. Phys. Chem. C 2012, 116, 5932–5940. (83) Roman, T.; Groß, A. Structure of water layers on hydrogen-covered Pt electrodes. Catal. Today. 2013, 202, 183–190. (84) Gohda, Y.; Schnur, S.; Groß, A. Influence of water on elementary reaction steps in electrocatalysis. Farad. Discuss 2009, 140, 233–244. ACS Paragon Plus Environment

17