Mechanistic and kinetic studies of the surface chemistry of iridium and ruthenium

Quantitative investigations of the mechanisms and the kinetics of the surface-catalyzed activation of C-H, N-H, C-C, and C-N bonds on the close-packed surfaces of Ir(111) and Ru(001) have been performed. The interaction of CH_3NH_2 with Ru(001) was investigated in ultrahigh vacuum with the techniques of high-resolution electron energy loss spectroscopy and thermal desorption mass spectrometry. Activation of the central C-N bond is observed, but it is less favored than the competing channel of complete dehydrogenation, by a ratio between 2:1 to 3:1. The decomposition mechanism has been characterized with several surface intermediates and gas-phase products identified. A pronounced preference for the activation of C-H over N-H and C-N bonds has been established. Additionally, the kinetics of the initial dissociation of short chain alkanes on Ir(111) has been examined, and the rate parameters of the activation of C-C bonds and primary, secondary, and tertiary C-H bonds have been determined. The formation of primary alkyl products is favored, over most of the experimental temperature range, despite the thermodynamic preference for the activation of individual secondary and tertiary C-H bonds in comparison to individual primary C-H bonds. At higher surface temperatures, the activation of C-C bonds occurs at competitive rates to the C-H reaction channel. The measured deuterium kinetic isotope effect implicates substantial deformation of the terminal methyl group in the transition state of C-C bond cleavage. Finally, the surface structure sensitivity of C-H bond cleavage has been quantified for smooth (111) and corrugated (110) surfaces of iridium and platinum, as well as for step edge defect sites on Ir(111).

Kinetic investigations of heterogeneously catalyzed reactions on the Ru(001) and Pt(110)-(1x2) surfaces

The initial probabilities of activated, dissociative chemisorption of methane and ethane on Pt(110)-(1 x 2) have been measured. The surface temperature was varied from 450 to 900 K with the reactant gas temperature constant at 300 K. Under these conditions, we probe the kinetics of dissociation via trapping-mediated (as opposed to 'direct') mechanism. It was found that the probabilities of dissociation of both methane and ethane were strong functions of the surface temperature with an apparent activation energies of 14.4 kcal/mol for methane and 2.8 kcal/mol for ethane, which implys that the methane and ethane molecules have fully accommodated to the surface temperature. Kinetic isotope effects were observed for both reactions, indicating that the C-H bond cleavage was involved in the rate-limiting step. A mechanistic model based on the trapping-mediated mechanism is used to explain the observed kinetic behavior. The activation energies for C-H bond dissociation of the thermally accommodated methane and ethane on the surface extracted from the model are 18.4 and 10.3 kcal/mol, respectively.

The studies of the catalytic decomposition of formic acid on the Ru(001) surface with thermal desorption mass spectrometry following the adsorption of DCOOH and HCOOH on the surface at 130 and 310 K are described. Formic acid (DCOOH) chemisorbs dissociatively on the surface via both the cleavage of its O-H bond to form a formate and a hydrogen adatom, and the cleavage of its C-O bond to form a carbon monoxide, a deuterium adatom and an hydroxyl (OH). The former is the predominant reaction. The rate of desorption of carbon dioxide is a direct measure of the kinetics of decomposition of the surface formate. It is characterized by a kinetic isotope effect, an increasingly narrow FWHM, and an upward shift in peak temperature with Ɵ_T, the coverage of the dissociatively adsorbed formic acid. The FWHM and the peak temperature change from 18 K and 326 K at Ɵ_T = 0.04 to 8 K and 395 K at Ɵ_T = 0.89. The increase in the apparent activation energy of the C-D bond cleavage is largely a result of self-poisoning by the formate, the presence of which on the surface alters the electronic properties of the surface such that the activation energy of the decomposition of formate is increased. The variation of the activation energy for carbon dioxide formation with Ɵ_T accounts for the observed sharp carbon dioxide peak. The coverage of surface formate can be adjusted over a relatively wide range so that the activation energy for C-D bond cleavage in the case of DCOOH can be adjusted to be below, approximately equal to, or well above the activation energy for the recombinative desorption of the deuterium adatoms. Accordingly, the desorption of deuterium was observed to be governed completely by the desorption kinetics of the deuterium adatoms at low Ɵ_T, jointly by the kinetics of deuterium desorption and C-D bond cleavage at intermediate Ɵ_T, and solely by the kinetics of C-D bond cleavage at high Ɵ_T. The overall branching ratio of the formate to carbon dioxide and carbon monoxide is approximately unity, regardless the initial coverage Ɵ_T, even though the activation energy for the production of carbon dioxide varies with Ɵ_T. The desorption of water, which implies C-O bond cleavage of the formate, appears at approximately the same temperature as that of carbon dioxide. These observations suggest that the cleavage of the C-D bond and that of the C-O bond of two surface formates are coupled, possibly via the formation of a short-lived surface complex that is the precursor to to the decomposition.

The measurement of steady-state rate is demonstrated here to be valuable in determining kinetics associated with short-lived, molecularly adsorbed precursor to further reactions on the surface, by determining the kinetic parameters of the molecular precursor of formaldehyde to its dissociation on the Pt(110)-(1 x 2) surface.

Overlayers of nitrogen adatoms on Ru(001) have been characterized both by thermal desorption mass spectrometry and low-energy electron diffraction, as well as chemically via the postadsorption and desorption of ammonia and carbon monoxide.

The nitrogen-adatom overlayer was prepared by decomposing ammonia thermally on the surface at a pressure of 2.8 x 10^(-6) Torr and a temperature of 480 K. The saturated overlayer prepared under these conditions has associated with it a (√247/10 x √247/10)R22.7° LEED pattern, has two peaks in its thermal desorption spectrum, and has a fractional surface coverage of 0.40. Annealing the overlayer to approximately 535 K results in a rather sharp (√3 x √3)R30° LEED pattern with an associated fractional surface coverage of one-third. Annealing the overlayer further to 620 K results in the disappearance of the low-temperature thermal desorption peak and the appearance of a rather fuzzy p(2x2) LEED pattern with an associated fractional surface coverage of approximately one-fourth. In the low coverage limit, the presence of the (√3 x √3)R30° N overlayer alters the surface in such a way that the binding energy of ammonia is increased by 20% relative to the clean surface, whereas that of carbon monoxide is reduced by 15%.

A general methodology for the indirect relative determination of the absolute fractional surface coverages has been developed and was utilized to determine the saturation fractional coverage of hydrogen on Ru(001). Formaldehyde was employed as a bridge to lead us from the known reference point of the saturation fractional coverage of carbon monoxide to unknown reference point of the fractional coverage of hydrogen on Ru(001), which is then used to determine accurately the saturation fractional coverage of hydrogen. We find that ƟSAT/H = 1.02 (±0.05), i.e., the surface stoichiometry is Ru : H = 1 : 1. The relative nature of the method, which cancels systematic errors, together with the utilization of a glass envelope around the mass spectrometer, which reduces spurious contributions in the thermal desorption spectra, results in high accuracy in the determination of absolute fractional coverages.

Vibrational pooling and constrained equilibration on surfaces

In this thesis, we provide a statistical theory for the vibrational pooling and fluorescence time dependence observed in infrared laser excitation of CO on an NaCl surface. The pooling is seen in experiment and in computer simulations. In the theory, we assume a rapid equilibration of the quanta in the substrate and minimize the free energy subject to the constraint at any time t of a fixed number of vibrational quanta N(t). At low incident intensity, the distribution is limited to one- quantum exchanges with the solid and so the Debye frequency of the solid plays a key role in limiting the range of this one-quantum domain. The resulting inverted vibrational equilibrium population depends only on fundamental parameters of the oscillator (ω_e and ω_eχ_e) and the surface (ω_D and T). Possible applications and relation to the Treanor gas phase treatment are discussed. Unlike the solid phase system, the gas phase system has no Debye-constraining maximum. We discuss the possible distributions for arbitrary N-conserving diatom-surface pairs, and include application to H:Si(111) as an example.

Computations are presented to describe and analyze the high levels of infrared laser-induced vibrational excitation of a monolayer of absorbed ¹³CO on a NaCl(100) surface. The calculations confirm that, for situations where the Debye frequency limited n domain restriction approximately holds, the vibrational state population deviates from a Boltzmann population linearly in n. Nonetheless, the full kinetic calculation is necessary to capture the result in detail.

We discuss the one-to-one relationship between N and γ and the examine the state space of the new distribution function for varied γ. We derive the Free Energy, F = NγkT − kTln(∑P_n), and effective chemical potential, μn ≈ γkT, for the vibrational pool. We also find the anti correlation of neighbor vibrations leads to an emergent correlation that appears to extend further than nearest neighbor.

Structure, chemistry, and energetics of organic and inorganic adsorbates on Ga-rich GaAs and GaP(001) surfaces

The work described in this dissertation includes fundamental investigations into three surface processes, namely inorganic film growth, water-induced oxidation, and organic functionalization/passivation, on the GaP and GaAs(001) surfaces. The techniques used to carry out this work include scanning tunneling microscopy (STM), X-ray photoelectron spectroscopy (XPS), and density functional theory (DFT) calculations. Atomic structure, electronic structure, reaction mechanisms, and energetics related to these surface processes are discussed at atomic or molecular levels.

First, we investigate epitaxial Zn₃P₂ films grown on the Ga-rich GaAs(001)(6×6) surface. The film growth mechanism, electronic properties, and atomic structure of the Zn₃P₂/GaAs(001) system are discussed based on experimental and theoretical observations. We discover that a P-rich amorphous layer covers the crystalline Zn3P2 film during and after growth. We also propose more accurate picture of the GaP interfacial layer between Zn₃P₂ and GaAs, based on the atomic structure, chemical bonding, band diagram, and P-replacement energetics, than was previously anticipated.

Second, DFT calculations are carried out in order to understand water-induced oxidation mechanisms on the Ga-rich GaP(001)(2×4) surface. Structural and energetic information of every step in the gaseous water-induced GaP oxidation reactions are elucidated at the atomic level in great detail. We explore all reasonable ground states involved in most of the possible adsorption and decomposition pathways. We also investigate structures and energies of the transition states in the first hydrogen dissociation of a water molecule on the (2×4) surface.

Finally, adsorption structures and thermal decomposition reactions of 1-propanethiol on the Ga-rich GaP(001)(2×4) surface are investigated using high resolution STM, XPS, and DFT simulations. We elucidate adsorption locations and their associated atomic structures of a single 1-propanethiol molecule on the (2×4) surface as a function of annealing temperature. DFT calculations are carried out to optimize ground state structures and search transition states. XPS is used to investigate variations of the chemical bonding nature and coverage of the adsorbate species.

Functionalization of Si(111) surfaces and the formation of mixed monolayers for the covalent attachment of molecular catalysts in photoelectrochemical devices

The functionalization of silicon surfaces with molecular catalysts for proton reduction is an important part of the development of a solar-powered, water-splitting device for solar fuel formation. The covalent attachment of these catalysts to silicon without damaging the underlying electronic properties of silicon that make it a good photocathode has proven difficult. We report the formation of mixed monolayer-functionalized surfaces that incor- porate both methyl and vinylferrocenyl or vinylbipyridyl (vbpy) moieties. The silicon was functionalized using reaction conditions analogous to those of hydrosilylation, but instead of a H-terminated Si surface, a chlorine-terminated Si precursor surface was used to produce the linked vinyl-modified functional group. The functionalized surfaces were characterized by time-resolved photoconductivity decay, X-ray photoelectron spectroscopy (XPS), electro- chemical, and photoelectrochemical measurements. The functionalized Si surfaces were well passivated, exhibited high surface coverage and few remaining reactive Si atop sites, had a very low surface recombination velocity, and displayed little initial surface oxidation. The surfaces were stable toward atmospheric and electrochemical oxidation. The surface coverage of ferrocene or bipyridine was controllably varied from 0 up to 30% of a monolayer without loss of the underlying electronic properties of the silicon. Interfacial charge transfer to the attached ferrocene group was relatively rapid, and a photovoltage of 0.4 V was generated upon illumination of functionalized n-type silicon surfaces in CH₃CN. The immobilized bipyridine ligands bound transition metal ions, and thus enabled the assembly of metal complexes on the silicon surface. XPS studies demonstrated that [Cp∗Rh(vbpy)Cl]Cl, [Cp∗Ir(vbpy)Cl]Cl, and Ru(acac)₂vbpy were assembled on the surface. For the surface prepared with iridium, x-ray absorption spectroscopy at the Ir L_III edge showed an edge energy and post-edge features virtually identical to a powder sample of [Cp∗Ir(bipy)Cl]Cl (bipy is 2,2 ́-bipyridyl). Electrochemical studies on these surfaces confirmed that the assembled complexes were electrochemically active.

Logarithmic Potential Theory on Riemann Surfaces

We develop a logarithmic potential theory on Riemann surfaces which generalizes logarithmic potential theory on the complex plane. We show the existence of an equilibrium measure and examine its structure. This leads to a formula for the structure of the equilibrium measure which is new even in the plane. We then use our results to study quadrature domains, Laplacian growth, and Coulomb gas ensembles on Riemann surfaces. We prove that the complement of the support of the equilibrium measure satisfies a quadrature identity. Furthermore, our setup allows us to naturally realize weak solutions of Laplacian growth (for a general time-dependent source) as an evolution of the support of equilibrium measures. When applied to the Riemann sphere this approach unifies the known methods for generating interior and exterior Laplacian growth. We later narrow our focus to a special class of quadrature domains which we call Algebraic Quadrature Domains. We show that many of the properties of quadrature domains generalize to this setting. In particular, the boundary of an Algebraic Quadrature Domain is the inverse image of a planar algebraic curve under a meromorphic function. This makes the study of the topology of Algebraic Quadrature Domains an interesting problem. We briefly investigate this problem and then narrow our focus to the study of the topology of classical quadrature domains. We extend the results of Lee and Makarov and prove (for n ≥ 3) c ≤ 5n-5, where c and n denote the connectivity and degree of a (classical) quadrature domain. At the same time we obtain a new upper bound on the number of isolated points of the algebraic curve corresponding to the boundary and thus a new upper bound on the number of special points. In the final chapter we study Coulomb gas ensembles on Riemann surfaces.

Depth-sensitive analysis of fluorine on lunar sample surfaces and in carbonaceous chondrites by a nuclear resonant reaction technique

The nuclear resonant reaction ¹⁹F(ρ,αγ)¹⁶O has been used to perform depth-sensitive analyses of fluorine in lunar samples and carbonaceous chondrites. The resonance at 0.83 MeV (center-of-mass) in this reaction is utilized to study fluorine surface films, with particular interest paid to the outer micron of Apollo 15 green glass, Apollo 17 orange glass, and lunar vesicular basalts. These results are distinguished from terrestrial contamination, and are discussed in terms of a volcanic origin for the samples of interest. Measurements of fluorine in carbonaceous chondrites are used to better define the solar system fluorine abundance. A technique for measurement of carbon on solid surfaces with applications to direct quantitative analysis of implanted solar wind carbon in lunar samples is described.

Effects of magnetostriction and superlattice formation in ferromagnetic thin films

The contribution to the magnetic uniaxial perpendicular anisotropy which arises from substrate constraint through magnetostrictive effects has been measured in Ni-Fe and Ni-Co thin films evaporated on substrates at room temperature. This was accomplished by measuring the perpendicular anisotropy before and after removal of the film from the substrate. Data are given for the fcc crystal structure regions of both alloy systems, but data for Ni-Co include compositions with less than 60% Ni which have a small percentage of the hcp phase mixed with the fcc phase. The constraint contribution to the perpendicular anisotropy correlates well with the value of the bulk magnetostriction constant using the equation ∆K_˔=3/2λ_sσ. Measured values of isotropic stress for films thicker than 600 Å were 1.6 x 10¹⁰ dyn/cm². In films less than 600 Å thick the isotropic stress decreased with decreasing thickness. After removal of the films from the substrates, the measured perpendicular anisotropy deviated from the expected geometrical shape anisotropy near pure Ni in both alloys. This indicates that additional significant sources of anisotropy exist at these compositions.

The effect of substrate constraint on the crystalline anisotropy K₁ of Ni-Fe epitaxial films has been studied by use of a film removal technique, which involves the evaporation of an epitaxial layer of LiF on MgO, the epitaxial growth of the metallic film on the LiF, and the stripping of the film with water soluble tape. Films ranging in composition from 50% to 100% Ni have been studied. For compositions below 90% Ni the experimental values agree reasonably well with the first order theoretical prediction, ∆K₁=[-9/4(C₁₁-C₁₂)λ² 100+9/2C₄₄λ²111].

In order to compare the magnetic properties of epitaxial thin films more completely with the properties of bulk single crystals, Ni-Fe films ranging in composition from 60% to 90% Ni, which were evaporated epitaxially on (100) MgO substrates, have been subsequently annealed at 400°C in a vacuum of less than 10^-7 Torr to form the ordered Ni₃Fe structure near the 75% composition. This ordered structure has been confirmed by electron diffraction.

The saturation magnetization at Ni₃Fe increased about 6% with ordering which is in good agreement with previous bulk data. Measurements of the magnetocrystalline anisotropy energy K₁ for the epitaxial films show the same large changes with ordering as observed in bulk single crystal samples. In the (001) plane the magnetostriction constants ^λ100, ^λ111 are directly related to the induced anisotropy due to a uniform uniaxial strain in the [100] and [110] directions respectively. Assuming that the elastic constants of a film are the same as in bulk material and are unchanged by ordering, the changes in strain sensitivity with ordering for the epitaxial films are found to be in good agreement with values predicted from bulk data. The exchange constant A as measured by ferromagnetic resonance has been measured at the Ni₃Fe composition and found to increase 25% with ordering. This seems to indicate a significant increase in the Curie temperature which has only been inferred indirectly for bulk material.