A perturbation theory study of H2 on LiF(0 0 1)

In response to recent helium atom scattering (HAS) results, Monte Carlo (MC) simulations and perturbation theory have been performed for H₂ on LiF(0 0 1). MC simulations predict that H₂ molecules form a series of structures, p(2 × 2), p(8 × 2), p(4 × 2) with coverages Θ = 0.5, 0.625 and 0.75, respectively, that are stable up to 8 K. This is in partial agreement with the HAS results that report c(2 × 2) and c(8 × 2) structures; they agree in terms of coverage and stability, but disagree in terms of symmetry. To reconcile the results of the simulations and experiments, the orientation of the adsorbed H₂ molecules was studied using perturbation theory. These calculations show that the adsorbed H₂ molecules are azimuthally delocalized and that the structures are c-type rather then p-type. The calculations also indicate that p-H₂ and helicoptering o-H₂ prefer cationic sites, while cartwheeling o-H₂ prefers anionic sites.     

A Monte Carlo simulation study of H2 layers on NaCl(0 0 1)

Monte Carlo simulations show that, at one monolayer coverage, H₂ molecules adsorbed on a NaCl(0 0 1) surface occupy all Na⁺ sites and form a commensurate c(2 × 2) structure. If the Cl⁻ sites are occupied as well, a bi-layer p(2 × 1) structure forms. An examination of the H₂ molecules’ rotational motion shows the molecular axes are azimuthally delocalized and so both of the structures acquire (1 × 1) symmetry in accord with experimental observations. These calculations also show that helicoptering o-H₂ (J = 1, m = ±1) prefer to sit on top of Na⁺ sites, while cartwheeling o-H₂ (J = 1, m = 0) prefers to locate over Cl⁻ sites, in agreement with other work.     

A valence band photoemission study of Pb adsorption on Rh(1 0 0) and Rh(1 1 0)

We have studied the adsorption of Pb on the Rh(1 0 0) and (1 1 0) surfaces by photoemission and low energy electron diffraction (LEED), and tested the chemical properties by adsorption of CO. Pb forms two distinct c(2 × 2) phases on Rh(1 0 0), according to the temperature of the substrate. The phase formed below about 570-620 K, denoted α-c(2 × 2), reduces the coverage of adsorbed CO but does not affect the valence band spectrum of the molecule. The phase formed above this temperature, denoted β-c(2 × 2), also reduces the coverage of adsorbed CO but the valence band spectrum of the adsorbed CO is strongly affected. The two phases are also characterised by a slightly different binding energy of the Pb 5d_5/2 level, 17.54 eV for the α phase and 17.70 for the β phase. The Pb/Rh(1 1 0) surface shows two ordered Pb induced phases, c(2 × 2) and p(3 × 1). CO adsorbs on the first with reduced heat of adsorption and with a valence band spectrum that is strongly altered with respect to CO adsorbed on clean Rh(1 1 0), but does not adsorb on the p(3 × 1) structure at 300 K. We compare the present results with previous results from related systems.     

Coverage temperature phase diagram of LEED spot and arced streak intensities for Li on Cu(0 0 1) surface using lattice gas model

It has been established that the arced streaks connecting four spots observed in LEED for a Li system adsorbed on a Cu(0 0 1) surface originate from the Bragg reflection from parallel adatomic lines on a c(2 × 2) lattice site. For example, one streak at about k_y = π/a originates from the parallel atomic lines including two atoms separated at a distance of d_y = 2a, which is the second-neighbor distance in a c(2 × 2) lattice.The c(2 × 2) structure sites consist of two sublattices with y = 2na and y = (2n + 1)a. Here, the difference in the number of adatoms on the two sublattices is the cause of the intensity of the midpoint of the streak, where the differences depend on the coverage of adatoms, Θ.In this study, using a lattice gas model on the substrate lattice with Monte Carlo simulation, we obtain the coverage and temperature dependence of intensities of the spots for the c(2 × 2) structure and the streaks.We found that the intensity of the streaks increase and decrease within the coverage range 0 < Θ < 0.5. That of the spots increases monotonically in this coverage range. These theoretical findings are similar to the experimental results.On the other hand, as temperature is increased, the intensity of the streaks increases and becomes saturated. We found a similar phenomenon using analytical calculation by statistical mechanics. In addition, the intensity of the spots decreased with the second-order transition.     

Cyclic cluster study of water adsorption structures on the MgO(1 0 0) surface

Cyclic cluster calculations were performed with the quantum chemical method MSINDO to elucidate the relative stabilities of c(4 × 2), p(3 × 2) and (1 × 1) overlayer structures of water molecules on the MgO(1 0 0) surface. For the c(4 × 2) and p(3 × 2) structures both molecular adsorption and partially dissociated adsorption were considered. In agreement with earlier theoretical studies partial dissociation was found to be more stable than molecular adsorption. For the c(4 × 2) structure both monolayer and double layer coverage were studied. Adsorption was found to be more stabilized with increasing degree of dissociation until 50% of the water molecules were dissociated. In the case of 50% dissociated water molecules we found that one quarter of the Mg atoms were pulled out of the MgO surface when surface relaxation was taken into account. A new structure for the fully dissociated (1 × 1) water monolayer was found which is considerably more stable than previously studied arrangements. In all cases surface relaxation was found to be important. The most stable structures of c(4 × 2), p(3 × 2) and (1 × 1) symmetry have adsorption energies which differ by no more than 13 kJ/mol. This offers an explanation of phase transitions of overlayer structures found by experiments between 180 and 300 K.     

Global SO(3) × SO(3) × U(1) symmetry of the Hubbard model on bipartite lattices

In this paper the global symmetry of the Hubbard model on a bipartite lattice is found to be larger than SO(4). The model is one of the most studied many-particle quantum problems, yet except in one dimension it has no exact solution, so that there remain many open questions about its properties. Symmetry plays an important role in physics and often can be used to extract useful information on unsolved non-perturbative quantum problems. Specifically, here it is found that for on-site interaction U ≠ 0 the local SU(2) × SU(2) × U(1) gauge symmetry of the Hubbard model on a bipartite lattice with N_a^D sites and vanishing transfer integral t = 0 can be lifted to a global [SU(2) × SU(2) × U(1)]/Z₂² = SO(3) × SO(3) × U(1) symmetry in the presence of the kinetic-energy hopping term of the Hamiltonian with t > 0. (Examples of a bipartite lattice are the D-dimensional cubic lattices of lattice constant a and edge length L = N_aa for which D = 1, 2, 3,... in the number N_a^D of sites.) The generator of the new found hidden independent charge global U(1) symmetry, which is not related to the ordinary U(1) gauge subgroup of electromagnetism, is one half the rotated-electron number of singly occupied sites operator. Although addition of chemical-potential and magnetic-field operator terms to the model Hamiltonian lowers its symmetry, such terms commute with it. Therefore, its 4^N_aD energy eigenstates refer to representations of the new found global [SU(2) × SU(2) × U(1)]/Z₂² = SO(3) × SO(3) × U(1) symmetry. Consistently, we find that for the Hubbard model on a bipartite lattice the number of independent representations of the group SO(3) × SO(3) × U(1) equals the Hilbert-space dimension 4^N_aD. It is confirmed elsewhere that the new found symmetry has important physical consequences.     

Structural fluctuation of dimers on Ge(0 0 1) surface near transition temperature

The structural fluctuation of the orientational arrangement of buckled dimers on a Ge(0 0 1) surface near the transition temperature of the order-disorder phase transition is investigated by time-resolving dynamical Monte Carlo simulations. STM images averaged in a finite period are derived from the simulation. The coexistence of the c(4 × 2) and the apparent (2 × 1) domains in the STM images observed by experiments is reproduced in the simulated STM images. We show that the coexistence on the Ge(0 0 1) surface can be attributed to the critical slowing down near the transition temperature.     

Combined STM, LEED and DFT study of Ag(1 0 0) exposed to oxygen near atmospheric pressures

We have investigated the interaction of molecular oxygen with the Ag(1 0 0) surface in a temperature range from 130 K to 470 K and an oxygen partial pressure ranging up to 10 mbar by scanning tunneling microscopy, low electron energy diffraction, Auger electron spectroscopy and ab initio density functional calculations. We find that at 130 K, following oxygen exposures of 6000 Langmuirs O₂, the individual oxygen atoms are randomly distributed on the surface. When the sample is exposed to 10 mbar O₂ at room temperature, small, p(2 × 2) reconstructed patches are formed on the surface. After oxidation at ≈470 K and 10 mbar O₂ pressure the surface undergoes a c(4 × 6) reconstruction coexisting with a (6 × 6) superstructure. By ab initio thermodynamic calculations it is shown that the c(4 × 6) reconstruction is an oxygen adsorption induced superstructure which is thermodynamically stable for an intermediate range of oxygen chemical potential.     

Study of c(2 × 2)-MnAu(0 0 1) layers on Mn(0 0 1) by means of scanning tunneling microscopy/spectroscopy

Intermixing, growth, geometric and electronic structures of gold films grown on antiferromagnetic stacking body-centered-tetragonal manganese (0 0 1) films were studied by means of scanning tunneling microscopy/spectroscopy at room temperature in ultra-high vacuum. We found stable ordered c(2 × 2)-MnAu(0 0 1) alloy layers after depositing Au on pure Mn layers. Since at the fourth layer (5 × 23)-like Au reconstruction appears instead of the c(2 × 2) structure and local density of states peaks obtained on the c(2 × 2)-MnAu surface disappear, pure Au layers likely grow from the fourth layer.     

CO2 adsorption on Cr(1 1 0) and Cr2O3(0 0 0 1)/Cr(1 1 0)

We attempt to correlate qualitatively the surface structure with the chemical activity for a metal surface, Cr(1 1 0), and one of its surface oxides, Cr₂O₃(0 0 0 1)/Cr(1 1 0). The kinetics and dynamics of CO₂ adsorption have been studied by low energy electron diffraction (LEED), Aug er electron spectroscopy (AES), and thermal desorption spectroscopy (TDS), as well as adsorption probability measurements conducted for impact energies of E_i = 0.1-1.1 eV and adsorption temperatures of T_s = 92-135 K. The Cr(1 1 0) surface is characterized by a square shaped LEED pattern, contamination free Cr AES, and a single dominant TDS peak (binding energy E_d = 33.3 kJ/mol, first order pre-exponential 1 × 10¹³ s⁻¹). The oxide exhibits a hexagonal shaped LEED pattern, Cr AES with an additional O-line, and two TDS peaks (E_d = 39.5 and 30.5 kJ/mol). The initial adsorption probability, S₀, is independent of T_s for both systems and decreases exponentially from 0.69 to 0.22 for Cr(1 1 0) with increasing E_i, with S₀ smaller by ∼0.15 for the surface oxide. The coverage dependence of the adsorption probability, S(Θ), at low E_i is approx. independent of coverage (Kisliuk-shape) and increases initially at large E_i with coverage (adsorbate-assisted adsorption). CO₂ physisorbs on both systems and the adsorption is non-activated and precursor mediated. Monte Carlo simulations (MCS) have been used to parameterize the beam scattering data. The coverage dependence of E_d has been obtained by means of a Redhead analysis of the TDS curves.     

Oxygen-induced p(3 × 1) reconstruction of the W(1 0 0) surface

The oxidation of the W(1 0 0) surface at elevated temperatures has been studied using room temperature STM and LEED. High exposure of the clean surface to O₂ at 1500 K followed by flash-annealing to 2300 K in UHV results in the formation of a novel p(3 × 1) reconstruction, which is imaged by STM as a missing-row structure on the surface. Upon further annealing in UHV, this surface develops a floreted LEED pattern characteristic of twinned microdomains of monoclinic WO_x, while maintaining the p(3 × 1) missing-row structure. Atomically resolved STM images of this surface show a complex domain structure with single and double W〈0 1 0〉 rows coexisting on the surface in different domains.     

Electron surface states in short-period superlattices: (GaAs)2/(AlAs)2(1 0 0)-c(4 × 4)

The electronic structure of (GaAs)₂/(AlAs)₂(1 0 0)-c(4 × 4) superlattice surfaces was studied by means of angular-resolved photoelectron spectroscopy (ARUPS) in the photon energy range 20-38 eV. Four samples with different surface termination layers were grown and As-capped by molecular beam epitaxy (MBE). ARUPS measurements were performed on decapped samples with perfect c(4 × 4) reconstructed surfaces. An intensive surface state was, for the first time, observed below the top of the valence band. This surface state was found to shift with superlattices’ different surface termination in agreement with theoretical predictions.     

SrTiO3(0 0 1) reconstructions: the (2 × 2) to c(4 × 4) transition

Scanning tunneling microscopy (STM) is used to investigate the (0 0 1) surface structure of Nb doped SrTiO₃ single crystals annealed in ultra high vacuum (UHV). Atomically resolved images of the (2 × 2) reconstructed surface are obtained after annealing a chemically etched sample. With further annealing dotted row domains appear, which coexist with the (2 × 2) reconstruction. The expansion of these domains with further annealing gives rise to the formation of a TiO₂ enriched c(4 × 4) reconstruction.     

Flip motion of heterogeneous buckled dimers on Ge(0 0 1) by electron injection from STM tip

Surface motion of a topological defect between p(2×2) and c(4×2) structures, a “kink”, across buckled Sn-Ge and Si-Ge dimers on Ge(0 0 1) surfaces was investigated using scanning tunneling microscopy. Energy thresholds of π^∗ electrons for flipping these dimers in the kink are obtained by analyzing the kink surface motion. Electronic states of these systems and energy barriers for flipping the dimers are examined by first-principles calculations for considering elementary processes of the electronically-excited flip motion of the dimers. We propose that the flip motion is caused by a resonant scattering of the π^∗ electrons with localized electronic states at the kink.     

Adsorption dynamics of CO2 on Cu(1 1 0): A molecular beam study

Molecular beam scattering measurements have been conducted to examine the adsorption dynamics of CO₂ on Cu(1 1 0). The initial adsorption probability, S₀, decreases exponentially from 0.43 ± 0.03 to a value close to the detection limit (∼0.03) within the impact energy range of E_i = (0.12-1.30) eV. S₀ is independent of the adsorption temperature, T_s, and the impact angle, α_i, i.e., the adsorption is non-activated and total energy scaling is obeyed. The coverage, Θ, dependent adsorption probability, S(Θ), agrees with precursor-assisted adsorption dynamics (Kisliuk type) above T_s ∼ 91 K. However, below that temperature adsorbate-assisted adsorption (S increases with Θ) has been observed. That effect is most distinct at large E_i and low T_s. The S(Θ) data have been modeled by Monte Carlo simulations. No indications of CO₂ dissociation were obtained from Auger Electron Spectroscopy or the molecular beam scattering data.     

Morphological change of D2O layers on Ru(0 0 0 1) probed with He atom scattering

The morphological change of D₂O layers on a Ru(0 0 0 1) surface has been investigated on the basis of He atom scattering. With the increase of D₂O exposure on Ru(0 0 0 1) at 111 K, the intensity of a specularly reflected He beam continuously decreases up to the exposure of 1.0 L (Langmuir). At the D₂O coverage of 1.0 adsorbed layer (∼1.5 L), which is characterized by temperature-programmed desorption measurements, the formation of the (√3 × √3)R30° superstructure as a result of the diffusion of D₂O on the surface was confirmed by He atom diffraction. With the further increase of D₂O exposure, at 2-3 adsorbed layers, the disordered structure was found to be on the surface at 111 K. The morphological change of the disordered layers was observed during annealing, and discussed in detail.     

Surface phase transition and electronic structure of c(5√2 × √2)R45°-Pb/Cu(1 0 0)

The c(5√2 × √2)R45°-Pb/Cu(1 0 0) surface phase is investigated by means of angle resolved ultraviolet photoemission and low energy electron diffraction in the temperature range between 300 and 550 K. We identify and characterize a temperature-induced surface phase transition at 440 K from the room temperature c(5√2 × √2) R45° phase to a (√2 × √2)R45° structure with split superstructure spots. The phase transition is fully reversible and takes place before the two-dimensional melting of the structure at 520 K. The electronic structure of the split (√2 × √2)R45° phase is characterized by a metallic free-electron like surface band. This surface band is backfolded with c(5√2 × √2)R45° periodicity phase at room temperature, giving rise to a surface band gap at the Fermi energy. We propose that a gain in electronic energy explains in part the stability of the c(5√2 × √2)R45° phase.     

Core-level shifts of the c(8 × 2)-reconstructed InAs(1 0 0) and InSb(1 0 0) surfaces

We have studied In-stabilized c(8 × 2)-reconstructed InAs(1 0 0) and InSb(1 0 0) semiconductor surfaces, which play a key role in growing improved III–V interfaces for electronics devices, by core-level photoelectron spectroscopy and first-principles calculations. The calculated surface core-level shifts (SCLSs) for the ζ and ζa models, which have been previously established to describe the atomic structures of the III–V(1 0 0)c(8 × 2) surfaces, yield hitherto not reported interpretation for the As 3d, In 4d, and Sb 4d core-level spectra of the III–V(1 0 0)c(8 × 2) surfaces, concerning the number and origins of SCLSs. The fitting analysis of the measured spectra with the calculated ζ and ζa SCLS values shows that the InSb spectra are reproduced by the ζ SCLSs better than by the ζa SCLSs. Interestingly, the ζa fits agree better with the InAs spectra than the ζ fits do, indicating that the ζa model describes the InAs surface better than the InSb surface. These results are in agreement with previous X-ray diffraction data. Furthermore, an introduction of the complete-screening model, which includes both the initial and final state effects, does not improve the fitting of the InSb spectra, proposing the suitability of the initial-state model for the SCLSs of the III–V(1 0 0)c(8 × 2) surfaces. The found SCLSs are discussed with the ab initio on-site charges.     

Initial stages of the autocatalytic oxidation of the InAs(0 0 1)-(4 × 2)/c(8 × 2) surface by molecular oxygen

The initial stages of oxidation of the In-rich InAs(0 0 1)-(4 × 2)/c(8 × 2) surface by molecular oxygen (O₂) were studied using scanning tunneling microscopy (STM) and density functional theory (DFT). It was shown that the O₂ dissociatively chemisorbs along the rows in the [1 1 0] direction on the InAs surface either by displacing the row-edge As atoms or by inserting between In atoms on the rows. The dissociative chemisorption is consistent with being autocatalytic: there is a high tendency to form oxygen chemisorption sites which grow in length along the rows in the [1 1 0] direction at preexisting oxygen chemisorption sites. The most common site size is about 21-24 Å in length at ∼25% ML coverage, representing 2-3 unit cell lengths in the [1 1 0] direction (the length of ∼5-6 In atoms on the row). The autocatalysis was confirmed by modeling the site distribution as non-Poisson. The autocatalysis and the low sticking probability (∼10⁻⁴) of O₂ on the InAs(0 0 1)-(4 × 2)/c(8 × 2) are consistent with activated dissociative chemisorption. The results show that is it critical to protect the InAs surface from oxygen during subsequent atomic layer deposition (ALD) or molecular beam epitaxy (MBE) oxide growth since oxygen will displace As atoms.