Atomic layer deposition of copper – study through density functional theory

The wonder of the last century has been the rapid development in technology. One of the sectors that it has touched immensely is the electronic industry. There has been exponential development in the field and scientists are pushing new horizons. There is an increased dependence in technology for every individual from different strata in the society. Atomic Layer Deposition (ALD) is a unique technique for growing thin films. It is widely used in the semiconductor industry. Films as thin as few nanometers can be deposited using this technique. Although this process has been explored for a variety of oxides, sulphides and nitrides, a proper method for deposition of many metals is missing. Metals are often used in the semiconductor industry and hence are of significant importance. A deficiency in understanding the basic chemistry at the nanoscale for possible reactions has delayed the improvement in metal ALD. In this thesis, we study the intrinsic chemistry involved for Cu ALD. This work reports computational study using Density Functional Theory as implemented in TURBOMOLE program. Both the gas phase and surface reactions are studied in most of the cases. The merits and demerits of a promising transmetallation reaction have been evaluated at the beginning of the study. Further improvements in the structure of precursors and coreagent have been proposed. This has led to the proposal of metallocenes as co-reagents and Cu(I) carbene compounds as new set of precursors. A three step process for Cu ALD that generates ligand free Cu layer after every ALD pulse has also been studied. Although the chemistry has been studied under the umbrella of Cu ALD the basic principles hold true for ALD of other metals (e.g. Co, Ni, Fe ) and also for other branches of science like thin film deposition other than ALD, electrochemical reactions, etc.

Hybrid density functional theory description of N- and C-doping of NiO

The large intrinsic bandgap of NiO hinders its potential application as a photocatalyst under visible-light irradiation. In this study, we have performed first-principles screened exchange hybrid density functional theory with the HSE06 functional calculations of N- and C-doped NiO to investigate the effect of doping on the electronic structure of NiO. C-doping at an oxygen site induces gap states due to the dopant, the positions of which suggest that the top of the valence band is made up primarily of C 2p-derived states with some Ni 3d contributions, and the lowest-energy empty state is in the middle of the gap. This leads to an effective bandgap of 1.7 eV, which is of potential interest for photocatalytic applications. N-doping induces comparatively little dopant-Ni 3d interactions, but results in similar positions of dopant-induced states, i.e., the top of the valence band is made up of dopant 2p states and the lowest unoccupied state is the empty gap state derived from the dopant, leading to bandgap narrowing. With the hybrid density functional theory (DFT) results available, we discuss issues with the DFT corrected for on-site Coulomb description of these systems.

Theory of elasticity and electric polarization effects in the group-III nitrides

In this work, the properties of strained tetrahedrally bonded materials are explored theoretically, with special focus on group-III nitrides. In order to do so, a multiscale approach is taken: accurate quantitative calculations of material properties are carried out in a quantum first-principles frame, for small systems. These properties are then extrapolated and empirical methods are employed to make predictions for larger systems, such as alloys or nanostructures. We focus our attention on elasticity and electric polarization in semiconductors. These quantities serve as input for the calculation of the optoelectronic properties of these systems. Regarding the methods employed, our first-principles calculations use highly- accurate density functional theory (DFT) within both standard Kohn-Sham and generalized (hybrid functional) Kohn-Sham approaches. We have developed our own empirical methods, including valence force field (VFF) and a point-dipole model for the calculation of local polarization and local polarization potential. Our local polarization model gives insight for the first time to local fluctuations of the electric polarization at an atomistic level. At the continuum level, we have studied composition-engineering optimization of nitride nanostructures for built-in electrostatic field reduction, and have developed a highly efficient hybrid analytical-numerical staggered-grid computational implementation of continuum elasticity theory, that is used to treat larger systems, such as quantum dots.

Precursor adsorption on copper surfaces as the first step during the deposition of copper: a density functional study with van der Waals correction

Copper dimethylamino-2-propoxide [Cu(dmap)2] is used as a precursor for low-temperature atomic layer deposition (ALD) of copper thin films. Chemisorption of the precursor is the necessary first step of ALD, but it is not known in this case whether there is selectivity for adsorption sites, defects, or islands on the substrate. Therefore, we study the adsorption of the Cu(dmap)2 molecule on the different sites on flat and rough Cu surfaces using PBE, PBE-D3, optB88-vdW, and vdW-DF2 methods. We found the relative order of adsorption energies for Cu(dmap)2 on Cu surfaces is Eads (PBE-D3) > Eads (optB88-vdW) > Eads (vdW-DF2) > Eads (PBE). The PBE and vdW-DF2 methods predict one chemisorption structure, while optB88-vdW predicts three chemisorption structures for Cu(dmap)2 adsorption among four possible adsorption configurations, whereas PBE-D3 predicts a chemisorbed structure for all the adsorption sites on Cu(111). All the methods with and without van der Waals corrections yield a chemisorbed molecule on the Cu(332) step and Cu(643) kink because of less steric hindrance on the vicinal surfaces. Strong distortion of the molecule and significant elongation of Cu–N bonds are predicted in the chemisorbed structures, indicating that the ligand–Cu bonds break during the ALD of Cu from Cu(dmap)2. The molecule loses its initial square-planar structure and gains linear O–Cu–O bonding as these atoms attach to the surface. As a result, the ligands become unstable and the precursor becomes more reactive to the coreagent. Charge redistribution mainly occurs between the adsorbate O–Cu–O bond and the surface. Bader charge analysis shows that electrons are donated from the surface to the molecule in the chemisorbed structures, so that the Cu center in the molecule is partially reduced.

Investigation of numerical atomic orbitals for first-principles calculations of the electronic and transport properties of silicon nanowire structures

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.

Carbon nanotubes as materials in nanotechnology

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.

Ab initio calculations of group 4 metallocene reaction mechanisms: atomic layer deposition and bond activation catalysis

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.

First principles modelling of nucleation and growth during atomic layer deposition onto III-V substrates

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.

Reactivity of sub 1 nm supported clusters: (TiO2)(n) clusters supported on rutile TiO2 (110)

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.

Multi-scale modelling of atomic layer deposition

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.

Engineering of metal oxide interfaces for renewable energy applications

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.

Precursor chemistry for oxide related materials

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.

Computational modeling of defects in nanoscale device materials

This PhD thesis concerns the computational modeling of the electronic and atomic structure of point defects in technologically relevant materials. Identifying the atomistic origin of defects observed in the electrical characteristics of electronic devices has been a long-term goal of first-principles methods. First principles simulations are performed in this thesis, consisting of density functional theory (DFT) supplemented with many body perturbation theory (MBPT) methods, of native defects in bulk and slab models of In0.53Ga0.47As. The latter consist of (100) - oriented surfaces passivated with A12O3. Our results indicate that the experimentally extracted midgap interface state density (Dit) peaks are not the result of defects directly at the semiconductor/oxide interface, but originate from defects in a more bulk-like chemical environment. This conclusion is reached by considering the energy of charge transition levels for defects at the interface as a function of distance from the oxide. Our work provides insight into the types of defects responsible for the observed departure from ideal electrical behaviour in III-V metal-oxidesemiconductor (MOS) capacitors. In addition, the formation energetics and electron scattering properties of point defects in carbon nanotubes (CNTs) are studied using DFT in conjunction with Green’s function based techniques. The latter are applied to evaluate the low-temperature, low-bias Landauer conductance spectrum from which mesoscopic transport properties such as the elastic mean free path and localization length of technologically relevant CNT sizes can be estimated from computationally tractable CNT models. Our calculations show that at CNT diameters pertinent to interconnect applications, the 555777 divacancy defect results in increased scattering and hence higher electrical resistance for electron transport near the Fermi level.

Reductive elimination of hypersilyl halides from zinc (II) complexes. Implications for electropositive metal thin film growth

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.

A microscopic study of structural and electronic properties of functionalized silicon surfaces based on first-principles

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.