982 resultados para DENSITY-FUNCTIONAL METHODS
Resumo:
Flikkema, E., & Bromley, S. T. (2004). Dedicated global optimization search for ground state silica nanoclusters: (SiO2)(N) (N=6-12). Journal of Physical Chemistry B, 108 (28), 9638-9645. RAE2008
Resumo:
This thesis is focused on the application of numerical atomic basis sets in studies of the structural, electronic and transport properties of silicon nanowire structures from first-principles within the framework of Density Functional Theory. First we critically examine the applied methodology and then offer predictions regarding the transport properties and realisation of silicon nanowire devices. The performance of numerical atomic orbitals is benchmarked against calculations performed with plane waves basis sets. After establishing the convergence of total energy and electronic structure calculations with increasing basis size we have shown that their quality greatly improves with the optimisation of the contraction for a fixed basis size. The double zeta polarised basis offers a reasonable approximation to study structural and electronic properties and transferability exists between various nanowire structures. This is most important to reduce the computational cost. The impact of basis sets on transport properties in silicon nanowires with oxygen and dopant impurities have also been studied. It is found that whilst transmission features quantitatively converge with increasing contraction there is a weaker dependence on basis set for the mean free path; the double zeta polarised basis offers a good compromise whereas the single zeta basis set yields qualitatively reasonable results. Studying the transport properties of nanowire-based transistor setups with p+-n-p+ and p+-i-p+ doping profiles it is shown that charge self-consistency affects the I-V characteristics more significantly than the basis set choice. It is predicted that such ultrascaled (3 nm length) transistors would show degraded performance due to relatively high source-drain tunnelling currents. Finally, it is shown the hole mobility of Si nanowires nominally doped with boron decreases monotonically with decreasing width at fixed doping density and increasing dopant concentration. Significant mobility variations are identified which can explain experimental observations.
Resumo:
Carbon nanotubes (CNTs) are hollow tubes of sp2-hybridised carbon with diameters of the order of nanometres. Due to their unique physical properties, which include ballistic transport and high mechanical strength, they are of significant interest for technological applications. The electronic properties of CNTs are of particular interest for use as gas sensors, interconnect materials in the semi-conductor industry and as the channel material in CNT based field effect transistors. The primary difficulty associated with the use of CNTs in electronic applications is the inability to control electronic properties at the growth stage; as grown CNTs consist of a mixture of metallic and semi-conducting CNTs. Doping has the potential to solve this problem and is a focus of this thesis. Nitrogen-doped CNTs typically have defective structures; the usual hollow CNT structure is replaced by a series of compartments. Through density functional theory (DFT) calculations and experimental results, we propose an explanation for the defective structures obtained, based on the stronger binding of N to the growth catalyst in comparison to C. In real electronic devices, CNTs need to be contacted to metal, we generate the current-voltage (IV) characteristics of metal-contacted CNTs considering both the effect of dopants and the structure of the interface region on electronic properties. We find that substitutionally doped CNTs produce Ohmic contacts and that scattering at the interface is strongly influenced by structure. In addition, we consider the effect of the common vacancy defects on the electronic properties of large diameter CNTs. Defects increase scattering in the CNT, with the greatest scattering occurring for the largest defect (555777). We validate the independent scattering approximation for small diameter CNTs, which enables mean free paths in large diameter CNTs to be calculated, with a smaller mean free paths found for larger defects.
Resumo:
Thin film dielectrics based on titanium, zirconium or hafnium oxides are being introduced to increase the permittivity of insulating layers in transistors for micro/nanoelectronics and memory devices. Atomic layer deposition (ALD) is the process of choice for fabricating these films, as it allows for high control of composition and thickness in thin, conformal films which can be deposited on substrates with high aspect-ratio features. The success of this method depends crucially on the chemical properties of the precursor molecules. A successful ALD precursor should be volatile, stable in the gas-phase, but reactive on the substrate and growing surface, leading to inert by-products. In recent years, many different ALD precursors for metal oxides have been developed, but many of them suffer from low thermal stability. Much promise is shown by group 4 metal precursors that contain cyclopentadienyl (Cp = C5H5-xRx) ligands. One of the main advantages of Cp precursors is their thermal stability. In this work ab initio calculations were carried out at the level of density functional theory (DFT) on a range of heteroleptic metallocenes [M(Cp)4-n(L)n], M = Hf/Zr/Ti, L = Me and OMe, in order to find mechanistic reasons for their observed behaviour during ALD. Based on optimized monomer structures, reactivity is analyzed with respect to ligand elimination. The order in which different ligands are eliminated during ALD follows their energetics which was in agreement with experimental measurements. Titanocene-derived precursors, TiCp*(OMe)3, do not yield TiO2 films in atomic layer deposition (ALD) with water, while Ti(OMe)4 does. DFT was used to model the ALD reaction sequence and find the reason for the difference in growth behaviour. Both precursors adsorb initially via hydrogen-bonding. The simulations reveal that the Cp* ligand of TiCp*(OMe)3 lowers the Lewis acidity of the Ti centre and prevents its coordination to surface O (densification) during both of the ALD pulses. Blocking this step hindered further ALD reactions and for that reason no ALD growth is observed from TiCp*(OMe)3 and water. The thermal stability in the gas phase of Ti, Zr and Hf precursors that contain cyclopentadienyl ligands was also considered. The reaction that was found using DFT is an intramolecular α-H transfer that produces an alkylidene complex. The analysis shows that thermal stabilities of complexes of the type MCp2(CH3)2 increase down group 4 (M = Ti, Zr and Hf) due to an increase in the HOMO-LUMO band gap of the reactants, which itself increases with the electrophilicity of the metal. The reverse reaction of α-hydrogen abstraction in ZrCp2Me2 is 1,2-addition reaction of a C-H bond to a Zr=C bond. The same mechanism is investigated to determine if it operates for 1,2 addition of the tBu C-H across Hf=N in a corresponding Hf dimer complex. The aim of this work is to understand orbital interactions, how bonds break and how new bonds form, and in what state hydrogen is transferred during the reaction. Calculations reveal two synchronous and concerted electron transfers within a four-membered cyclic transition state in the plane between the cyclopentadienyl rings, one π(M=X)-to-σ(M-C) involving metal d orbitals and the other σ(C-H)-to-σ(X-H) mediating the transfer of neutral H, where X = C or N. The reaction of the hafnium dimer complex with CO that was studied for the purpose of understanding C-H bond activation has another interesting application, namely the cleavage of an N-N bond and resulting N-C bond formation. Analysis of the orbital plots reveals repulsion between the occupied orbitals on CO and the N-N unit where CO approaches along the N-N axis. The repulsions along the N-N axis are minimized by instead forming an asymmetrical intermediate in which CO first coordinates to one Hf and then to N. This breaks the symmetry of the N-N unit and the resultant mixing of MOs allows σ(NN) to be polarized, localizing electrons on the more distant N. This allowed σ(CO) and π(CO) donation to N and back-donation of π*(Hf2N2) to CO. Improved understanding of the chemistry of metal complexes can be gained from atomic-scale modelling and this provides valuable information for the design of new ALD precursors. The information gained from the model decomposition pathway can be additionally used to understand the chemistry of molecules in the ALD process as well as in catalytic systems.
Resumo:
Atomic layer deposition (ALD) is now used in semiconductor fabrication lines to deposit nanometre-thin oxide films, and has thus enabled the introduction of high-permittivity dielectrics into the CMOS gate stack. With interest increasing in transistors based on high mobility substrates, such as GaAs, we are investigating the surface treatments that may improve the interface characteristics. We focus on incubation periods of ALD processes on III-V substrates. We have applied first principles Density Functional Theory (DFT) to investigate detailed chemistry of these early stages of growth, specifically substrate and ALD precursor interaction. We have modelled the ‘clean-up’ effect by which organometallic precursors: trimethylaluminium (TMA) or hafnium and titanium amides clean arsenic oxides off the GaAs surface before ALD growth of dielectric commences and similar effect on Si3N4 substrate. Our simulations show that ‘clean-up’ of an oxide film strongly depends on precursor ligand, its affinity to the oxide and the redox character of the oxide. The predominant pathway for a metalloid oxide such as arsenic oxide is reduction, producing volatile molecules or gettering oxygen from less reducible oxides. An alternative pathway is non-redox ligand exchange, which allows non-reducible oxides (e.g. SiO2) to be cleaned-up. First principles study shows also that alkylamides are more susceptible to decomposition rather than migration on the oxide surface. This improved understanding of the chemical principles underlying ‘clean-up’ allows us to rationalize and predict which precursors will perform the reaction. The comparison is made between selection of metal chlorides, methyls and alkylamides precursors.
Resumo:
Metal oxide clusters of sub-nm dimensions dispersed on a metal oxide support are an important class of catalytic materials for a number of key chemical reactions, showing enhanced reactivity over the corresponding bulk oxide. In this paper we present the results of a density functional theory study of small sub-nm TiO2 clusters, Ti2O4, Ti3O6 and Ti4O8 supported on the rutile (110) surface. We find that all three clusters adsorb strongly with adsorption energies ranging from -3 eV to -4.5 eV. The more stable adsorption structures show a larger number of new Ti-O bonds formed between the cluster and the surface. These new bonds increase the coordination of cluster Ti and O as well as surface oxygen, so that each has more neighbours. The electronic structure shows that the top of the valence band is made up of cluster derived states, while the conduction band is made up of Ti 3d states from the surface, resulting in a reduction of the effective band gap and spatial separation of electrons and holes after photon absorption, which shows their potential utility in photocatalysis. To examine reactivity, we study the formation of oxygen vacancies in the cluster-support system. The most stable oxygen vacancy sites on the cluster show formation energies that are significantly lower than in bulk TiO2, demonstrating the usefulness of this composite system for redox catalysis.
Resumo:
High-permittivity ("high-k") dielectric materials are used in the transistor gate stack in integrated circuits. As the thickness of silicon oxide dielectric reduces below 2 nm with continued downscaling, the leakage current because of tunnelling increases, leading to high power consumption and reduced device reliability. Hence, research concentrates on finding materials with high dielectric constant that can be easily integrated into a manufacturing process and show the desired properties as a thin film. Atomic layer deposition (ALD) is used practically to deposit high-k materials like HfO2, ZrO2, and Al2O3 as gate oxides. ALD is a technique for producing conformal layers of material with nanometer-scale thickness, used commercially in non-planar electronics and increasingly in other areas of science and technology. ALD is a type of chemical vapor deposition that depends on self-limiting surface chemistry. In ALD, gaseous precursors are allowed individually into the reactor chamber in alternating pulses. Between each pulse, inert gas is admitted to prevent gas phase reactions. This thesis provides a profound understanding of the ALD of oxides such as HfO2, showing how the chemistry affects the properties of the deposited film. Using multi-scale modelling of ALD, the kinetics of reactions at the growing surface is connected to experimental data. In this thesis, we use density functional theory (DFT) method to simulate more realistic models for the growth of HfO2 from Hf(N(CH3)2)4/H2O and HfCl4/H2O and for Al2O3 from Al(CH3)3/H2O.Three major breakthroughs are discovered. First, a new reaction pathway, ’multiple proton diffusion’, is proposed for the growth of HfO2 from Hf(N(CH3)2)4/H2O.1 As a second major breakthrough, a ’cooperative’ action between adsorbed precursors is shown to play an important role in ALD. By this we mean that previously-inert fragments can become reactive once sufficient molecules adsorb in their neighbourhood during either precursor pulse. As a third breakthrough, the ALD of HfO2 from Hf(N(CH3)2)4 and H2O is implemented for the first time into 3D on-lattice kinetic Monte-Carlo (KMC).2 In this integrated approach (DFT+KMC), retaining the accuracy of the atomistic model in the higher-scale model leads to remarkable breakthroughs in our understanding. The resulting atomistic model allows direct comparison with experimental techniques such as X-ray photoelectron spectroscopy and quartz crystal microbalance.
Resumo:
Diminishing non-renewable energy resources and planet-wide de-pollution on our planet are among the major problems which mankind faces into the future. To solve these problems, renewable energy sources such as readily available and inexhaustible sunlight will have to be used. There are however no readily available photocatalysts that are photocatalytically active under visible light; it is well established that the band gap of the prototypical photocatalyst, titanium dioxide, is the UV region with the consequence that only 4% of sun light is utilized. For this reason, this PhD project focused on developing new materials, based on titanium dioxide, which can be used in visible light activated photocatalytic hydrogen production and destruction of pollutant molecules. The main goal of this project is to use simulations based on first principles to engineer and understand rationally, materials based on modifying TiO2 that will have the following properties: (1) a suitable band gap in order to increase the efficiency of visible light absorption, with a gap around 2 – 2.5 eV considered optimum. (2). The second key aspect in the photocatalytic process is electron and hole separation after photoexcitation, which enable oxidation/reduction reactions necessary to i.e. decompose pollutants. (3) Enhanced activity over unmodified TiO2. In this thesis I present results on new materials based on modifying TiO2 with supported metal oxide nanoclusters, from two classes, namely: transition metal oxides (Ti, Ni, Cu) and p-block metal oxides (Sn, Pb, Bi). We find that the deposited metal oxide nanoclusters are stable at rutile and anatase TiO2 surfaces and present an analysis of changes to the band gap of TiO2, identifying those modifiers that can change the band gap to the desirable range and the origin of this. A successful collaboration with experimental researchers in Japan confirms many of the simulation results where the origin of improved visible light photocatalytic activity of oxide nanocluster-modified TiO2 is now well understood. The work presented in this thesis, creates a road map for the design of materials with desired photocatalytic properties and contributes to better understanding these properties which are of great application in renewable energy utilization.
Resumo:
This thesis describes modelling, synthesis, spectroscopic and physical characterisation, as well as application of Magnesium, Calcium and Copper β-diketonate, β-ketoiminate, β-diiminate, Schiff base, amide and fluorenyl compounds. The selected compounds could potentially find application in materials deposition using Atomic Layer Deposition (ALD), MOCVD, CVD and Sol-Gel techniques. Quantum chemical modelling was used as a tool to perform the comprehensive and rapid study of magnesium and calcium precursor molecules in order to predict which of them would be more successful in ALD of metal oxides. Precursor chemistry plays a key role in ALD, since precursors must be volatile, thermally stable, chemisorb on the surface and react rapidly with existing surface groups. This Thesis describes one aspect of this, surface reactivity between ligands and hydroxyl groups, via a gas-phase model with energetics computed at the level of Density Functional Theory (DFT). A number of different synthetic strategies, both aerobic and anaerobic, were investigated for the synthesis of the described metal complexes. These included the use of different metal starting reagents such as, anhydrous and hydrated inorganic metal salts, metal alkyls and Grignard reagents. Some of previously unreported metal complexes of homoleptic and heteroleptic magnesium, calcium and copper β-diketonates, β-ketoiminates, β-diiminates, amides and Schiff base type were synthesised and characterised: [Mg(hfpd)2(DipPa)], [Mg(hfpd)2(MapH)2], [Mg(hf-ebp)(THF)2], [Mg(tf-Pap)Cl(THF)2], [Ca(PhNacnac)2], [Cu(tf-Pap)2], [Cu(PhNacnac)2], [Cu(hf-ebp)], [Cu(DipPa)] and [Cu(DipPa)2(4,4’-bypy)]. A comprehensive study on the thermal properties of magnesium, calcium and copper β-diketonates, β-ketoiminates, β-diiminates, Schiff base, amide and fluorenyl complexes was performed using TGA and sublimation of selected compounds. Atomic Layer Deposition of MgO using magnesium β-ketoiminate – [bis{(4-N-phenyl)-2-pentonato} magnesium] and β-diketonate - [bis(1,1,1,5,5,5-hexafluoropentane-2,4-dionato)(THF)magnesium hydrate] was performed on Si(100) substrates at 180°C and 0.2 Torr using O2 plasma.
Resumo:
Treatment of Zn(Si(SiMe3)3)2 with ZnX2 (X = Cl, Br, I) in tetrahydrofuran (THF) at 23 °C afforded [Zn(Si(SiMe3)3)X(THF)]2 in 83–99% yield. X-ray crystal structures revealed dimeric structures with Zn2X2 cores. Thermogravimetric analyses of [Zn(Si(SiMe3)3)X(THF)]2 demonstrated a loss of coordinated THF between 50 and 155 °C and then single-step weight losses between 200 and 275 °C. The nonvolatile residue was zinc metal in all cases. Bulk thermolyses of [Zn(Si(SiMe3)3)X(THF)]2 between 210 and 250 °C afforded zinc metal in 97–99% yield, Si(SiMe3)3X in 91–94% yield, and THF in 81–98% yield. Density functional theory calculations confirmed that zinc formation becomes energetically favorable upon THF loss. Similar reactions are likely to be general for M(SiR3)n/MXn pairs and may lead to new metal-film-growth processes for chemical vapor deposition and atomic layer deposition.
Resumo:
Surface modification of silicon with organic monolayers tethered to the surface by different linkers is an important process in realizing future (opto-)electronic devices. Understanding the role played by the nature of the linking group and the chain length on the adsorption structures and electronic properties of these assemblies is vital to advance this technology. This Thesis is a study of such properties and contributes in particular to a microscopic understanding of induced changes in the work function of experimentally studied functionalized silicon surfaces. Using first-principles density functional theory (DFT), at the first step, we provide predictions for chemical trends in the work function of hydrogenated silicon (111) surfaces modified with various terminations. For nonpolar terminating atomic species such as F, Cl, Br, and I, the change in the work function is directly proportional to the amount of charge transferred from the surface, thus relating to the difference in electronegativity of the adsorbate and silicon atoms. The change is a monotonic function of coverage in this case, and the work function increases with increasing electronegativity. Polar species such as −TeH, −SeH, −SH, −OH, −NH2, −CH3, and −BH2 do not follow this trend due to the interaction of their dipole with the induced electric field at the surface. In this case, the magnitude and sign of the surface dipole moment need to be considered in addition to the bond dipole to generally describe the change in work function. Compared to hydrogenated surfaces, there is slight increase in the work function of H:Si(111)-XH, where X = Te, Se, and S, whereas reduction is observed for surfaces covered with −OH, −CH3, and −NH2. Next, we study the hydrogen passivated Si(111) surface modified with alkyl chains of the general formula H:Si–(CH2)n–CH2 and H:Si–X–(CH2)n–CH3, where X = NH, O, S and n = (0, 1, 3, 5, 7, 9, 11), at half coverage. For (X)–Hexyl and (X)–Dodecyl functionalization, we also examined various coverages up to full monolayer grafting in order to validate the result of half covered surface and the linker effect on the coverage. We find that it is necessary to take into account the van der Waals interaction between the alkyl chains. The strongest binding is for the oxygen linker, followed by S, N, and C, irrespective of chain length. The result revealed that the sequence of the stability is independent of coverage; however, linkers other than carbon can shift the optimum coverage considerably and allow further packing density. For all linkers apart from sulfur, structural properties, in particular, surface-linker-chain angles, saturate to a single value once n > 3. For sulfur, we identify three regimes, namely, n = 0–3, n = 5–7, and n = 9–11, each with its own characteristic adsorption structures. Where possible, our computational results are shown to be consistent with the available experimental data and show how the fundamental structural properties of modified Si surfaces can be controlled by the choice of linking group and chain length. Later we continue by examining the work function tuning of H:Si(111) over a range of 1.73 eV through adsorption of alkyl monolayers with general formula -[Xhead-group]-(CnH2n)-[Xtail-group], X = O(H), S(H), NH(2). The work function is practically converged at 4 carbons (8 for oxygen), for head-group functionalization. For tail-group functionalization and with both head- and tail-groups, there is an odd-even effect in the behavior of the work function, with peak-to-peak amplitudes of up to 1.7 eV in the oscillations. This behavior is explained through the orientation of the terminal-group's dipole. The shift in the work function is largest for NH2-linked and smallest for SH-linked chains and is rationalized in terms of interface dipoles. Our study reveals that the choice of the head- and/or tail-groups effectively changes the impact of the alkyl chain length on the work function tuning using self-assembled monolayers and this is an important advance in utilizing hybrid functionalized Si surfaces. Bringing together the understanding gained from studying single type functionalization of H:Si(111) with different alkyl chains and bearing in mind how to utilize head-group, tail-group or both as well as monolayer coverage, in the final part of this Thesis we study functionalized H:Si(111) with binary SAMs. Aiming at enhancing work function adjustment together with SAM stability and coverage we choose a range of terminations and linker-chains denoted as –X–(Alkyl) with X = CH3, O(H), S(H), NH(2) and investigate the stability and work function of various binary components grafted onto H:Si(111) surface. Using binary functionalization with -[NH(2)/O(H)/S(H)]-[Hexyl/Dodecyl] we show that work function can be tuned within the interval of 3.65-4.94 eV and furthermore, enhance the SAM’s stability. Although direct Si-C grafted SAMs are less favourable compared to their counterparts with O, N or S linkage, regardless of the ratio, binary functionalized alkyl monolayers with X-alkyl (X = NH, O) is always more stable than single type alkyl functionalization with the same coverage. Our results indicate that it is possible to go beyond the optimum coverage of pure alkyl functionalized SAMs (50%) by adding a linker with the correct choice of the linker. This is very important since dense packed monolayers have fewer defects and deliver higher efficiency. Our results indicate that binary anchoring can modify the charge injection and therefore bond stability while preserving the interface electronic structure.
Resumo:
It is known that the exact density functional must give ground-state energies that are piecewise linear as a function of electron number. In this work we prove that this is also true for the lowest-energy excited states of different spin or spatial symmetry. This has three important consequences for chemical applications: the ground state of a molecule must correspond to the state with the maximum highest-occupied-molecular-orbital energy, minimum lowest-unoccupied-molecular-orbital energy, and maximum chemical hardness. The beryllium, carbon, and vanadium atoms, as well as the CH(2) and C(3)H(3) molecules are considered as illustrative examples. Our result also directly and rigorously connects the ionization potential and electron affinity to the stability of spin states.
Resumo:
Localized molecular orbitals (LMOs) are much more compact representations of electronic degrees of freedom than canonical molecular orbitals (CMOs). The most compact representation is provided by nonorthogonal localized molecular orbitals (NOLMOs), which are linearly independent but are not orthogonal. Both LMOs and NOLMOs are thus useful for linear-scaling calculations of electronic structures for large systems. Recently, NOLMOs have been successfully applied to linear-scaling calculations with density functional theory (DFT) and to reformulating time-dependent density functional theory (TDDFT) for calculations of excited states and spectroscopy. However, a challenge remains as NOLMO construction from CMOs is still inefficient for large systems. In this work, we develop an efficient method to accelerate the NOLMO construction by using predefined centroids of the NOLMO and thereby removing the nonlinear equality constraints in the original method ( J. Chem. Phys. 2004 , 120 , 9458 and J. Chem. Phys. 2000 , 112 , 4 ). Thus, NOLMO construction becomes an unconstrained optimization. Its efficiency is demonstrated for the selected saturated and conjugated molecules. Our method for fast NOLMO construction should lead to efficient DFT and NOLMO-TDDFT applications to large systems.
Resumo:
Investigations of the vibrational spectra of cyclo(Gly-Gly), cyclo(L-Ala-L-Ala) and cyclo(t-Ala-Gly) are reported. Raman scattering and Fourier transform infrared (FTIR) spectra of solid-state and aqueous protonated samples, as well as their corresponding N-deuterated isotopomers, have been examined. In addition, density functional theory (DFT) (B3-LYP/cc-pVDZ) calculations of molecular structures and their associated vibrational modes were carried out. In each case, the calculated structures of lowest energy for the isolated gas-phase molecules have boat conformations. Assignments have been made for the observed Raman and FTIR vibrational bands of the cyclic di-amino acid peptides (CDAPs) examined. Raman polarization studies of aqueous phase samples are consistent with C-2 and C-1 symmetries for the six-membered rings of cyclo(L-Ala-L-Ala) and cydo(L-Ala-Gly), respectively. There is a good correlation between experimental and calculated vibrational bands for the three CDAPs. These data are in keeping with boat conformations for cydo(L-Ala-L-Ala) and cyclo(L-Ala-Gly) molecules, predicted by the ab initio calculations, in both the solid and aqueous solution states. However, Raman spectroscopic results might infer that cyclo(L-AlaGly) deviates only slightly from planarity in the solid state. The potential energy distributions of the amide I and II modes of a cis-peptide linkage are shown to be significantly different from those of the trans-peptides. For example, deuterium shifts have shown that the cis-amide I vibrations found in cyclo(Gly-Gly), cyclo(L-Ala-L-Ala), and cyclo(L-Ala-Gly) have larger N-H contributions compared to their trans-amide counterparts. Compared to trans-amide II vibrations, cis-amide II vibrations show a considerable decrease in N-H character.
Resumo:
The study of catalytic behavior begins with one seemingly simple process, namely the hydrogenation of O to H2O on platinum. Despite the apparent simplicity its mechanism has been much debated. We have used density functional theory with,gradient corrections to examine microscopic reaction pathways for several elementary steps implicated in this fundamental catalytic process. We find that H2O formation from chemisorbed O and H atoms is a highly activated process. The largest barrier along this route, with a value of similar to1 eV, is the addition of the first H to O to produce OH. Once formed, however, OH groups are easily hydrogenated to H2O with a barrier of similar to0.2 eV. Disproportionation reactions with 1:1 and 2:1 stoichiometries of H2O and O have been examined as alternative routes for OH formation. Both stoichiometries of reaction produce OH groups with barriers that are much lower than that associated with the O + H reaction. H2O, therefore, acts as an autocatalyst in the overall H O formation process. Disproportionation with a 2:1 stoichiometry is thermodynamically and kinetically favored over disproportionation with a l:I stoichiometry. This highlights an additional (promotional) role of the second H2O molecule in this process. In support of our previous suggestion that the key intermediate in the low-temperature H2O formation reaction is a mixed OH and H2O overlayer we find that then is a very large barrier for the dissociation of the second H2O molecule in the 2:1 disproportionation process. We suggest that the proposed intermediate is then hydrogenated to H2O through a very facile proton transfer mechanism.