Adsorption and desorption of hydrogen at nonpolar GaN(1-100) surfaces: Kinetics and impact on surface vibrational and electronic properties

The adsorption of hydrogen at nonpolar GaN(1-100) surfaces and its impact on the electronic and vibrational properties is investigated using surface electron spectroscopy in combination with density functional theory (DFT) calculations. For the surface mediated dissociation of H2 and the subsequent adsorption of H, an energy barrier of 0.55 eV has to be overcome. The calculated kinetic surface phase diagram indicates that the reaction is kinetically hindered at low pressures and low temperatures. At higher temperatures ab-initio thermodynamics show, that the H-free surface is energetically favored. To validate these theoretical predictions experiments at room temperature and under ultrahigh vacuum conditions were performed. They reveal that molecular hydrogen does not dissociatively adsorb at the GaN(1-100) surface. Only activated atomic hydrogen atoms attach to the surface. At temperatures above 820 K, the attached hydrogen gets desorbed. The adsorbed hydrogen atoms saturate the dangling bonds of the gallium and nitrogen surface atoms and result in an inversion of the Ga-N surface dimer buckling. The signatures of the Ga-H and N-H vibrational modes on the H-covered surface have experimentally been identified and are in good agreement with the DFT calculations of the surface phonon modes. Both theory and experiment show that H adsorption results in a removal of occupied and unoccupied intragap electron states of the clean GaN(1-100) surface and a reduction of the surface upward band bending by 0.4 eV. The latter mechanism largely reduces surface electron depletion.


I. INTRODUCTION
Gallium nitride (GaN) based optoelectronic devices are well established in solid state lighting 1-3 and power electronics. [4][5][6] High electron mobility transistor (HEMT) structures based on AlGaN/GaN have for example been modified in order to use these polar thin film devices as sensors with open or functionalized gate 7,8 , which essentially consist of the bare surface in polar orientation interacting with the surrounding species. In recent years, three-dimensional GaN nanowire structures have also attracted much attention for sensing applications 9-12 , light emission and detection [13][14][15] as well as solar water splitting [16][17][18] and photovoltaics. 19,20 They exhibit a high surface to volume ratio and are mainly composed of side facets consisting of the nonpolar m-plane (1100) surface. Such nanowire structures are typically grown by catalyst-free molecular beam epitaxy 21,22 and exhibit superior structural quality being almost free of defects and strain 23,24 with the capability to integrate vertical coreshell [25][26][27] and embedded lateral heterostructures 28,29 and to intentionally dope the material. 13,30 For sensor applications as well as to identify optimum growth conditions a detailed understanding of the interaction and adsorption of molecules and atoms in the gas phase with the surface is crucial. Of special importance is also the existence of free or saturated bonds and related surface, adsorbate or trap states with well defined electronic structure. 31 They can induce charge transfer processes, band bending including accumulation or depletion of electrons and/or form a dipole at the surface or interface in focus. 32,33 For example, the characteristics of GaN single-nanowire transistors have been found to be dependent on the valence band bending at the mplane side facets which is directly influenced by surface adsorbates. 34 Hydrogen is a simple model adsorbate system and is known to affect doping in GaN. [35][36][37] Understanding its influence is of great technological relevance, since GaN bulk and thin film growth techniques involve hydrogen directly or indirectly as possible dissociation product from ammonia or metalorganic precursors with impact on the growth and properties of the resulting material. [38][39][40] For polar GaN surfaces, the interaction with hydrogen has been studied both by experimental methods [41][42][43][44][45][46][47] and theory. [48][49][50][51][52][53][54][55] The nonpolar m-plane GaN(1100) orientation is a low energy surface of wurtzite GaN 56 and consists of buckled Ga-N dimers in the outermost layer. 57,58 The dangling bonds at the nitrogen and gallium surface atoms give rise to occupied and unoccupied surface states that have been predicted by density functional theory (DFT) calculations [59][60][61] and experimentally verified. 61,62 The chemically clean surface has been experimentally observed to exhibit a distinct electron depletion layer with a surface upward band bending of ∼ 0.6 eV. The position changes in the presence of gas molecules that adsorb on the surface. 62 The kinetics and thermodynamics of adsorption and desorption, the resulting surface/adatom geometric structures and their influence on the electronic properties of the surface and subsurface region are crucial to understand charge transfer processes across the semiconductor/adlayer interface as well as ionized adsorbate induced formation of surface dipoles.
DFT studies have focused on the stability of hydrogen species at nonpolar m-plane 63 and a-plane 64 GaN surfaces and found that under hydrogen-rich conditions hydrogen adsorbates attach at the free surface dangling bonds of both the Ga and N dimer atoms. It was further predicted that water molecules interacting with this surface spontaneously dissociate and form H and OH that bond with the surface. 65 In this study we combine first-principles calculations with surface adsorption/desorption experiments to clarify the mechanisms by which H and H 2 adsorb onto this surface. Based on this insight we study how the adsorbed H atoms modify the vibrational and electronic properties of this surface. We show that these aspects have consequences on growing GaN in hydrogen-rich environments and derive consequences when using nonpolar surfaces in electronic and chemical sensing devices.

II. EXPERIMENT AND THEORY
Clean GaN(1100) surfaces were prepared by homoepitaxial overgrowth on bulk GaN substrates from Kyma Technologies produced by hydride vapor phase epitaxy using molecular beam epitaxy (MBE). The used substrates were unintentionally n-doped crystals with a carrier concentration of ∼ 5×10 16 cm -3 and had an epi-ready surface finish achieved by a final chemo-mechanical polishing step. 66 Atomic force microscopy measurements identify atomically flat surfaces exhibiting a terrace width in correspondence with the sample miscut and a root-mean-square roughness below 0.5 nm. Growth of an a few hundred nanometer thick GaN epilayer was performed using a Knudsen cell for Ga evaporation and a SVTA RF 4.5 plasma source (13.56 MHz) for the generation of reactive nitrogen species. The growth parameters Ga flux and substrate temperature were optimized at a constant nitrogen flux (p N ∼ 5×10 -8 bar, plasma power 450 W) using reflection high energy electron diffraction (RHEED) during, and photoelectron spectroscopy (PES) after growth, in order to obtain stoichiometric surfaces which are free of excess Ga or surface defects. The properties of the GaN(1100) samples after growth have been reported and discussed earlier in Ref. 62. After epitaxy and cooling down, the samples were directly transferred under ultrahigh vacuum (UHV) conditions (base pressure < 2×10 -13 bar) to the respective position for in-situ surface analysis by ultraviolet and X-ray photoelectron spectroscopy (UPS, XPS). These measurements were performed in normal emission using a hemispherical electron analyzer and monochromated AlKα (1486.7 eV) or He I (21.2 eV) radiation for electron excitation. The de-scription of the employed experimental conditions and parameters for PES can be found in Ref. 67. To investigate their interaction with hydrogen, the asgrown samples were exposed at room temperature to hydrogen (purity 99.999 %) by backfilling the analysis chamber (p H2 = 2.0 × 10 −11 bar) for up to 55 min. Prior to each adsorption experiment, the gas supply lines were thoroughly evacuated to a pressure below 1 × 10 −10 bar and subsequently filled with 1.5 bar H 2 . The molecular hydrogen was optionally activated by a hot filament placed close to the sample front side to obtain atomic hydrogen by partial dissociation of the H 2 molecules in front of the GaN surface. Due to a limited cracking efficiency of the hot filament, the actual amount of produced atomic H species is below the calculated total H 2 exposure. During exposure, the residual gas was monitored by quadrupole mass spectrometry to control gas purity and absence of impurities. The pressure was measured with a Bayard Alpert ionization gauge and used without any further correction of specific gas sensitivity to calculate the exposure in Langmuir (1L = 1.33× 10 −9 bar·s). The changes of the surface properties were examined in-situ by UPS and XPS.
In an additional experiment, the hydrogen-covered mplane GaN sample was transferred to a second UHV recipient using a vacuum transfer chamber with a base pressure < 1×10 -11 bar. While clean surfaces were found to be extremely reactive to molecules from the residual gas, hydrogen adsorption at the GaN(1100) surface resulted in a passivation of the surface and relatively stable conditions for vacuum transfer. The second UHV recipient is optimized for electron energy loss spectroscopy (EELS) using an Ibach spectrometer. 68 Such experiments were performed on the H-covered GaN(1100) surface in specular scattering geometry with energies of monochromatic electron beams varying between 5 and 80 eV and were repeated after desorption of the hydrogen adsorbates by annealing the sample in UHV at 820 ± 50 K.
The H adsorption energies have been calculated employing DFT, the generalized gradient approximation (GGA) and the projector augmented-wave (PAW) method. 69,70 The Ga 3d electrons are included as valence states. The surfaces are modeled using a slab geometry consisting of 12 Ga-N monolayers (MLs) and a vacuum region of 20Å. A plane-wave energy cutoff of 450 eV was used and the Brillouin zone (BZ) was sampled using an equivalent 4 × 4 × 1 Monkhorst-Pack k-point mesh for the 1 × 1 surface unit cell. The lowermost Ga and N atoms were passivated with pseudohydrogen having a fractional charge of 0.75 and 1.25, respectively. Convergence with respect to k-point sampling, energy cutoff, vacuum, and slab thickness were explicitly tested and found to provide surface energies with an accuracy better than 3 meV/1 × 1.
In order to investigate the vibrational properties and the vibrational entropic contributions we have calculated the dynamical matrix of the free and adsorbate covered surfaces. The force constant matrix and in turn the dy-namical matrix have been calculated for the top four surface atomic layers and the H atoms at the surface using slabs of 8 layer thickness, 4 × 4 surface supercells and a displacement of 0.01Å in both directions. To evaluate the H 2 chemical potential we have included the translational, rotational, and vibrational contributions. More details regarding the approach and the convergence criteria can be found in Ref. 71.
The kinetics of dissociative adsorption of molecular H 2 are addressed with the harmonic transition state theory. 72 The transition states have been identified by climbing image nudged elastic band (NEB) calculations. 73 For the NEB calculations a 2 × 2 surface slab with a thickness of 8 MLs has been implemented and in total 6 images including the two stable/metastable end states, i.e. H 2 in the vacuum and H 2 bound to a surface dimer, have been used to identify the transition points. The electronic structure of clean and H-covered m-plane surfaces have been computed with the Heyd, Scuseria, and Ernzerhof hybrid functional. 74 This functional gives a bulk bandgap of E g = 3.116 eV in agreement with previous calculations. 75

A. Monitoring of hydrogen adsorption/desorption
After growth the clean GaN(1100) surface exhibits electron depletion with strong upward band bending. A distinct occupied surface state is found close to the valence band edge at 3.1 eV below the Fermi energy. This feature was identified as emission from an occupied surface state related to the filled dangling bond states located at the N atoms of the GaN(1100) surface dimer structure. For a detailed discussion on the electronic properties after growth we refer to our earlier study. 62 Since the focus of the present study is on hydrogen adsorption the variation of the valence band (VB) spectra has been monitored using UPS in continuing hydrogen interaction experiments increasing the exposure up to 50 Langmuir (L). In a first experiment molecular H 2 was offered to the clean GaN(1100) surface at room temperature (RT) resulting in negligible changes of the valence band features and surface band bending (not shown). Consequently, at RT no significant H 2 dissociation and H adsorption is observed for the used H 2 partial pressure of 2.0 × 10 −11 bar and the chosen reaction time.
As a consequence, in a following experiment we have initiated the H adsorption process by implementing a hot filament for partial generation of thermally activated atomic hydrogen close to the sample surface. Fig. 1 (a) shows a series of UPS (He I) valence band spectra during ongoing H and H 2 exposure up to 50 L in total. A gradual shift of the occupied states away from the Fermi level E F at 0 eV is observed. This effect is directly linked to a reduction of the surface band bending, which initially amounts to 0.6 eV for the as grown m-plane GaN surface. 62 In addition, the onset of low-energy secondary electron emission (not shown) was shifted, pointing to a reduction of the work function Φ. The temporal variation of Φ as well as the determined change in surface band bending ∆V bb are plotted in Fig. 1 (b) and (c), respectively. Both values decrease monotonically with H exposure, with the tendency of saturation at the end of the experiment. In parallel, the signal intensity of the filled N dangling bond state, initially found at 3.1 eV binding energy, is significantly reduced upon the interaction process [Figs. 1 (a) and 2 (d)]. These aspects provide indirect evidence for H adsorption at the GaN(1100) surface for this second experiment in which the sample was exposed to activated H species.
To prove the observed shift in surface band bending, we have also characterized the core level binding energies using AlKα X-ray excitation. The corresponding spectra of the Ga2p 3/2 , N1s and Ga3d states are shown in Fig. 2 (a) -(c). Their respective binding energies for the as-grown surface are 1118.0, 397.7 and 20.1 eV. We emphasize that no surface impurities were detected by XPS after the MBE growth or the subsequent H 2 exposure, and therefore the observed changes are not induced by surface impurities, while unfortunately a direct detection and analysis of hydrogen surface species is not possible by XPS. After the performed H exposure, in all cases a shift of the core level binding energy by 0.4 eV towards higher values is observed, consistent with the changes observed in the He I spectra during the adsorption process. These observations indicate that the reaction with the activated hydrogen species results in a saturation of the free dangling bonds at the Ga-N surface. The impact of the saturation will be discussed in detail in comparison with theoretical predictions below.
In order to analyze the stability or reversibility of the H adsorption process, we have afterwards annealed the sample in ultrahigh vacuum. While an annealing temperature of 520 ± 50 K did not substantially change the surface electronic properties in terms of band bending and work function indicating a certain stability of the adsorbate structure, heating the substrate up to 820 ± 50 K almost recovered the characteristics after MBE growth. Figs. 2 (a) -(d) include the photoelectron spectra obtained after annealing the hydrogenated surface at 820 ± 50 K for 10 min (black). Obviously, the occupied states shift back towards the Fermi level and most noticeably, the near VB edge emission from the occupied surface state also recovers in intensity to almost the signal strength after growth [ Fig. 2 (d)]. Consequently, the hydrogen atoms adsorbed with an effective coverage in the submonolayer regime during the reaction of activated H species and the clean GaN(1100) surface, can be reversed, i.e. desorption initiated, if sufficient energy is introduced (a) Specular EEL spectra (circles) of clean GaN(1100). The energy of incident electrons was set to 20 eV. The full line depicts calculated results. The prominent loss features are due to single and multiple electron scattering from the Fuchs-Kliewer phonon at 88 meV. (b) Surface band alignment and electron density profile resulting from modeling the experimental EEL spectra measured using different primary electron energies.
into the system, e.g. thermally induced as examined in this experiment. Fig. 2 (e) schematically summarizes the experimentally determined differences in electronic properties for the clean and H-covered GaN(1100) surface including variation in surface band bending and work function.
The aforementioned results indicate a complex interplay between thermodynamics and kinetics in hydrogen adsorption/desorption processes: they imply that it is energetically favorable for atomic hydrogen to adsorb at the surface, passivate the Ga-and N-dangling bonds of the clean GaN(1100) surface and saturate the surface states. This results in a reduction of the upward band bending by 0.4 eV and reduction of the work function by 0.2 eV as illustrated in the surface band diagram in Fig. 2 (e). However, hydrogen desorption can be mediated at elevated temperatures which can be attributed to the existence of a kinetic barrier. These aspects will be further addressed below.

B. Electron density profile and surface vibrations
While PES unraveled changes in the electronic structure upon H adsorption, vibrational spectroscopy was performed to confirm adsorption of atomic H and to identify adsorption sites. The discovered possibility to saturate the surface by atomic hydrogen, forming a stable adsorbate (passivation) layer was used to transfer the Hcovered GaN samples to a separate UHV chamber. The same heating procedure as used for the PES analysis was then performed to remove the hydrogen adatoms and to analyze the vibrational characteristics of the GaN(1100) surface. Figure 3 (a) shows a representative spectrum measured after transfer of the H-covered GaN sample and annealing at 820 ± 50 K in specular scattering geometry and a primary electron energy (E p ) of 20 eV. Apart from the signature of elastically scattered electrons at 0 meV the spectra of the clean surface exhibit almost equidistantly separated loss features. These peaks are assigned to the spectroscopic signatures of the Fuchs-Kliewer (FK) phonon at (88.0 ± 0.4) meV (average value of EEL spectra measured at different E p ) and its multiple excitations. This interpretation is corroborated by previous findings for GaN(0001) and GaN(0001) surfaces 43,45,46 and the simulations to be discussed next.
EELS data were simulated using a model which is similar to the theoretical approach developed in Refs. 76 and 77. In these calculations, the surface energy-loss function is derived using the continued-fraction expansion method 77 , when the subsurface region is represented by a finite number of sublayers of certain thickness to reproduce a smooth variation of the depth dependent electron density below the surface. For this purpose, electron density depth profiles are computed by solving the Schrödinger and Poisson equations self-consistently. It should be noted that the model used for fitting the measured spectra includes only two contributions, originating in collective lattice vibrations (phonons) and oscillations of the free-electron gas in the conduction band (plasmons). Also, due to sufficiently low bulk electron density [see Fig. 3 (b)] in these samples, the plasmons can manifest themselves only as a small broadening of the elastic peak, in addition to the instrumental resolution (∼ 3 meV for the implemented experimental conditions).
The resulting electron density profile from modeling the combined EEL spectral data with varying E p and the associated band edge alignment in the near surface region are depicted in Fig. 3 (b). The calculations are in very good agreement with the qualitative model of upward band bending [ Fig. 2 (e)] with an V bb value of 0.7 eV compared to 0.6 eV extracted from the PES measurements. Furthermore some important quantitative information can be extracted: the calculations indicate a bulk electron concentration of 5 × 10 16 cm −3 combined with a depletion layer thickness of ∼ 100 nm. Consequently, these results confirm the low electron concentration in these samples and that PES can be used for the determination of the band bending values since the width of the depletion layer is much larger than the information depth of a few nanometers from the surface obtained in PES. In this region the slope of the bands is found to be negligible within the uncertainty of 0.1 eV. Figure 4 shows specular EEL spectra of clean [ Fig. 4 (a)] and H-covered [ Fig. 4  those reported for GaN(0001) and (0001) surfaces. 43,45,46 For the H-covered surface, additional signatures are observed and marked by solid lines in Fig. 4. Corresponding FK phonon combination losses (e.g. ν Ga-H + FK1) are marked by dotted lines of the same color. The feature at 233 meV is caused by Ga-H stretching vibrations 46 of H adsorbed at the Ga-dangling bond of the Ga-N dimer. The corresponding N-H stretching vibration mode is also observed at 403 meV.
In Fig. 4 (c) the projected bulk phonon density of states (DOS) (gray shaded area) as well as the phonon DOS of the clean and H-covered GaN(1100) surfaces (red and blue curves, respectively) are shown. Each DOS is the sum over all states within the respective bulk or surface Brillouin zone. Besides the acoustic (0 -40 meV) and optical (60 -90 meV) bulk phonon modes 78 , which are not directly detected in the EELS experiment, three additional features are found for the m-plane surface with H atoms adsorbed at the Ga-N dimer dangling bonds. The N-H stretching mode (ν N-H ) is identified as nondispersing state at a calculated vibrational energy of 419 meV and the corresponding Ga-H stretching mode (ν Ga-H ) exhibits slight dispersion between 229 and 232 meV within the BZ (231.4 meV at the Γ-point). Both values are in fairly good agreement with the EELS experiment. In addition, the signature between 113 meV and 123 meV can be assigned to bending vibrations of H atoms adsorbed at the surface dimer structure (δ Ga-H and δ N-H ). 63 The calculations reveal that this structure consists of two states that disperse in the BZ within the mentioned energy range with a Γ-point energy of 115.6 and 117.2 meV, respectively. These vibrational energies are nearly two times larger than the energy reported for the Ga-H bending mode at the GaAs(110) surface. 79 The deviation may be ascribed to the difference of the microscopic adsorption geometry.
A faint shoulder on the high-energy side of the first FK phonon loss feature is observed around 118 meV and is attributed to these calculated Ga-H and N-H bending modes, rather than to N-OH vibrations that were previously reported to exist at a slightly lower vibrational energy of 106 meV. 46 The additional shoulder at 206 meV is attributed to a superposition of contributions from a FK phonon combination loss (δ Ga-H,N-H + FK1) and an O-H bending vibrational mode (δ O-H ). 46 The latter aspect is corroborated by the presence of a weak loss structure at 453 meV as side feature of the FK5 multiple, which is caused by O-H stretching vibrations (ν O-H ). 45,46 Consequently, the main spectroscopic features are assigned to vibrations of atomic H adsorbed to Ga and N surface atoms. A slight uptake of hydroxides is indicated by the EEL spectra and is due to the high reactivity of unsaturated GaN surfaces 47 combined with the two orders of magnitude higher base pressure in the used vacuum transfer system compared to the recipient for in-situ PES analyses. After annealing at 820 ± 50 K signatures of adsorbate vibrational modes fall below the detection limit of the spectrometer [ Fig. 4 (a)]. As a result, the EELS experiment provided important information that atomic H saturates Ga-and N-dangling bonds of the surface dimer structure.

C. Influence of hydrogen on the structural and electronic properties of GaN(1100) surfaces
In order to investigate the effect of H adsorption on the electronic properties of GaN surfaces and to develop a microscopic model of the differences in surface geometry, we have performed DFT calculations to compute the surface crystal structures and the band structures of the clean and hydrogen-covered m-plane GaN surfaces. In Table I  2.4, 2.68 and 2.7 eV, respectively. However, it is in good agreement with previous self-energycorrected LDA-1/2 calculations which predict a surface band gap of 3.03 eV but smaller than the value of 3.31 eV calculated by HSE with 32% fraction of exact exchange and the Ga 3d electrons treated as valence states. 80 In all the aforementioned high level DFT calculations the unoccupied surface state is well below the bulk conduction band minimum (CBM). The differences in the calculated surface band gaps can be attributed to the different methods as well as to the different slab thickness employed in these calculations. However, the position of the unoccupied surface state, i.e. 2.98 eV above the bulk VBM, is in good agreement with the measured band bending of ≈ 0.6 eV as discussed in detail below.
The band structure of the hydrogen-passivated mplane GaN surface is depicted in Fig. 5(b): Passivation of the surface dangling bonds by hydrogen results in Ga-H and N-H occupied bonding and unoccupied antibonding states. The former shift from above the VB edge for the clean surface into the bulk VB region. The unoccupied states, initially also found as intragap states shift into the bulk conduction band (CB) region. Hence, hydrogen passivation shifts the surface states out of the fundamental band gap, providing a suitable explanation for the discovered changes in band bending as measured by PES. even above the CB edge as calculated for the H-saturated surface, the surface Fermi level follows resulting in a reduced upward band bending V bb or even unpins the surface Fermi level resulting in flat band conditions. For the performed experiment, ∆V bb is 0.4 eV, indicating a remaining slight upward band bending/electron depletion. However, from the experimental data it is not possible to extract an exact number for the H-coverage for this experiment and one might expect an asymptotic convergence to the situation of a fully covered surface (considered in the calculations) for higher exposure.

D. Thermodynamics and Kinetics of hydrogen adsorption
In order to address the thermodynamics of hydrogen adsorption, the free energy difference ∆F of the hydrogen-covered and clean GaN(1100) surfaces was calculated as where ∆E tot = E is the difference in vibrational contributions to the free energy, and µ H2 is the chemical potential of H 2 . In Fig. 6 the difference in the surface free energies is plotted as function of temperature and H 2 pressure. Higher temperatures and/or lower pressures favor the clean surface. This is at-tributed to the large translational entropic contributions H 2 molecules have in the gas phase at these conditions. On the other side, at low temperatures and/or high partial pressures the hydrogen-covered surface is thermodynamically favorable. More specifically, at 300 K and for H 2 pressures larger than 5×10 −19 bar it is thermodynamically favorable to adsorb hydrogen at the m-plane GaN surface. However, this is in contrast to the experimental finding that at the same temperature and at 8 orders of magnitude higher pressure (i.e. 2.0×10 −11 bar) of molecular H 2 , no significant hydrogen adsorption is observed. Furthermore, the annealing experiments indicate that dehydrogenation of the surface, within the time scale of the experiments, requires elevated temperatures as high as 820 K. This further indicates that kinetic effects rather than the thermodynamic properties control the H coverage on the surfaces. In order to identify and investigate these mechanisms we next focus on the adsorption and desorption kinetics of both atomic and molecular hydrogen.
The adsorption of atomic hydrogen is barrierless. The desorption/binding energy E des of atomic hydrogen is defined as: where E surf:H and E surf are the total energies of the surface with and without adsorbed hydrogen atom, respectively and E Hatom is the total energy of a hydrogen atom. In the calculation of the desorption energies different effects have to be considered: First the N-H bond is stronger than the Ga-H. Second, both unpassivated and doubly passivated surface dimers, i.e. both Ga and N atoms of the same dimer are passivated by hydrogen, obey the electron counting rule and do not introduce occupied states deep in the fundamental gap. On the contrary, passivation of only cation or anion dangling bonds of one dimer will result in the formation of fully or partially occupied states deep in the gap region. Furthermore, as has already been discussed, surface relaxation and re-hybridization effects result in different atomic geometries for the clean and doubly passivated surface dimers (see Table I).
As discussed above, bond enthalpies, electronic structure and surface strain are expected to strongly influence the desorption energies. Thus, they have to be explicitly considered by investigating different atomic hydrogen desorption scenarios. Hence, different desorption mechanisms have been calculated. These correspond to H desorption from cations and anions at doubly and singly passivated surface dimers in the limit of clean and fully covered surfaces. Our calculations reveal that H desorption from cations is energetically preferred to desorption from anions by at least 0.5 eV. This value is considerably higher than k B T at 300 K (0.026 eV) or even at temperatures as high as 1000 K (0.086 eV). The corresponding desorption energies are ≈2.5 eV and ≈4.8 eV in the limit of a clean or fully covered surface, respectively.
The aforementioned reaction mechanism might not be the most relevant since adsorption/desorption of hydrogen may be realized by the formation of H 2 molecules. As with atomic hydrogen desorption, different reaction mechanisms have been considered: Adsorption/desorption at (i) a surface dimer, (ii) two neighboring Ga surface atoms, and (iii) two neighboring N surface atoms considering the two limits, i.e. of a clean and fully covered surface. The corresponding energy barriers are listed in Table II. The adsorption mechanism with the lowest adsorption barrier is H 2 dissociatively binding at a single surface dimer. In Fig. 7  quired to dissociate the molecule. On the other hand the desorption energy is ≈ 2.30 eV. It should be noted that the aforementioned energy barriers depend weakly on the surface coverage. In the limit of a fully hydrogen-covered surface the adsorption and desorption energy barriers are ≈ 0.62 eV and ≈ 2.30 eV, respectively. H 2 desorption from two neighboring N or Ga surface atoms in the limit of a fully covered surface has higher kinetic barriers, 5.43 and 2.85 eV, respectively. These reaction mechanisms result in single passivated surface dimers and increase both the surface strain and the electronic contributions to the surface energy (see above). Interestingly, the desorption barrier of 0.52 eV is remarkably rather small for H 2 binding at two neighboring Ga atoms in the limit of a clean surface. This is because it is highly unfavorable for H to passivate only surface cations even at extreme H-rich conditions. In order that this mechanism can actually take place, the H atoms would already have been desorbed from the N surface atoms. However, the latter has considerably higher kinetic barriers. Hence, desorption of molecular hydrogen from neighboring surface cations or anions can be neglected and molecular hydrogen adsorption and desorption is taking place by preferentially binding to and desorbing from Ga-N dimers.
The desorption barrier for the hydrogen molecule is considerably smaller than the desorption energy of atomic hydrogen. Hence, desorption is taking place as molecular H 2 rather than as atomic hydrogen. On the other hand for the adsorption of a H 2 molecule, an energy barrier has to be overcome, while adsorption of atomic hydrogen is barrierless. Thus, if both atomic and molecular hydrogen are present in the gas phase, then surface passivation by hydrogen will preferentially take place through atomic hydrogen adsorption and the rate limiting mechanism will be the flux of incident hydrogen atoms at the surface.
The flux of incident particles at a surface depends on the temperature and the corresponding partial pressure p and is given by the Hertz-Knudsen equation: 82,83 where m is the mass of the corresponding particles. For example, for an atomic hydrogen partial pressure of 10 −11 bar at T = 300 K the flux of incident atomic hydrogen is ≈ 0.025 s −1 per 1 × 1 surface cell area. Under these conditions and assuming that hydrogen desorption is kinetically suppressed, 50% or 100% of a monolayer surface coverage of an initially clean surface will be achieved within ≈ 1 min and ≈ 10 min, respectively. In contrast to atomic hydrogen the adsorption of H 2 molecules is not barrierless and the corresponding rate is given by the following equation: where A is the area of the 1 × 1 surface unit cell and E ads is the kinetic barrier for adsorption. It has to be noted here that the sticking coefficient of adsorbing H 2 depends on the orientation as well as the impinging angle of the molecule. 84,85 This dependency is not included in the used model. Nevertheless, Eq. 4 provides an upper limit for the adsorption frequency and hence a lower limit of the time to achieve thermodynamic equilibrium. Similarly the desorption frequency reads: where the attempt frequency ν 0 is calculated within the harmonic transition state theory: 72 Here ν min i are the 3N eigenfrequencies at the minimum and ν sad i are the 3N -1 nonimaginary eigenfrequencies at the transition point. In order to estimate the attempt frequency we have calculated the dynamical matrix of the 4 topmost atomic layers of an 8 layer thick 2 × 2 slab using the small displacements method for a H 2 molecule (i) adsorbed at the surface and (ii) at the transition point. The attempt frequency calculated from Eq. 6 is ν 0 = 5.6 × 10 13 s −1 .
The temporal evolution of the surface coverage is described by the following rate equation: where c is the surface coverage and t is the time. Equation (7) was solved for I c(t = 0) = 0 (clean surface) and for II c(t = 0) = 1 (fully covered surface). The former starting condition (I) corresponds to an adsorption experiment and the latter (II) to a desorption experiment. In Fig. 8 the coverage is plotted against exposure time to molecular H 2 for the two aforementioned initial boundary conditions and for various H 2 pressures and temperatures. A striking finding is that although at H 2 pressures in the order of 2 × 10 −11 bar at RT it is thermodynamically favorable to adsorb hydrogen at the surface and the equilibrium coverage is almost 100%, the time scale to achieve 10% or 100% surface coverage is more than a century or a millennium, respectively. Hence, within the time scale of the performed adsorption experiments at RT and partial pressures as low as 10 −11 bar (section A) thermodynamic equilibrium between H 2 gas and the GaN(1100) surface is kinetically hindered and only atomic hydrogen is able to adsorb quickly and to induce changes in the structural and the electronic surface properties. On the other hand, although the desorption barrier of H 2 is larger than the H 2 adsorption barrier, desorption of molecular hydrogen can take place in considerably shorter time scales. This can be attributed to (i) the higher temperatures applied to dehydrogenate the surface and (ii) the considerably larger desorption attempt frequency ν 0 than the adsorption attempt frequency f [see Eq. (4)]. For example, at the aforementioned temperature and pressure, the adsorption attempt frequency f in Eq. (4) is 2 × 10 −7 s −1 per 1 × 1 surface cell area as opposed to the desorption attempt frequency of ν 0 = 5.6 × 10 13 Hz in Eq. (5).
These results are consistent with the observations made in the UHV adsorption/desorption experiments. At low temperatures and low partial pressures thermodynamic equilibrium of the clean m-plane GaN surface with a H 2 atmosphere is kinetically hindered (compare to the solid black line in Fig. 8 which represents the experimental conditions for the performed experiment of H 2 exposure). On the other hand, for the interaction of the clean GaN(1100) surface with activated atomic H thermodynamic equilibrium between the surface and the offered H species can be established within the time scale of a few minutes even at pressures as low as 10 −11 bar.

IV. SUMMARY
Photoelectron and electron energy loss spectroscopy experiments were combined with first-principles calculations to investigate adsorption and desorption of molecular as well as atomic hydrogen on the nonpolar GaN(1100) surfaces. Our results show that passivation of the surface cation and anion dangling bonds by hydrogen is thermodynamically favored at room temperature even at hydrogen pressures as low as 10 −16 bar. Adsorption of molecular hydrogen is associated with a barrier of 0.55 eV, which leads to unrealistically high exposure times to complete a full monolayer. By contrast, only a few minutes of exposure time are required if the clean m-plane GaN surface is brought into an atomic hydrogen atmosphere. On the other hand, hydrogen desorption requires elevated temperatures and is predominantly taking place as hydrogen molecules desorbing from surface dimers. More specifically, our adsorption/desorption experiments confirm that in UHV at room temperature, H-adsorption could only be achieved using atomic hydrogen and desorption starts above ∼ 800 K.
The differences in the surface properties of the bare and the H-covered GaN(1100) surface have been characterized by electron spectroscopy and density functional calculations. Stretching and bending vibrations of the H adatoms at the Ga-N surface dimer structure were identified experimentally and are consistent with the energies and dispersion of calculated surface phonon modes. Furthermore, H adsorption was demonstrated to strongly influence the surface electronic properties. H adatoms occupy the Ga-and N-dangling bonds and induce a shift of occupied and unoccupied surface states out of the gap region across the VBM and CBM, respectively, which influences the surface electron depletion layer. Specifically, the raise of unoccupied intragap surface states from below the CBM for the clean surface to energies above the CBM for the H-covered surface induces an unpinning of the surface Fermi level and a reduction of the surface upward band bending from 0.6 to 0.2 eV.
These insights demonstrate that chemisorption in a gas exposure experiment or furthermore the probability of impurity incorporation during crystal growth is strongly dependent on the experimental conditions as well as on the energetics and kinetics of the surface dissociation/adsorption reactions at the surface. For other reactions at GaN surfaces, one might expect comparable kinetically hindered reaction mechanisms that lead to strong deviations between the actually obtained coverage and the supplied reactant load.