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.

Computational study of the growth of copper thin films by atomic layer deposition

Copper is the main interconnect material in microelectronic devices, and a 2 nm-thick continuous Cu film seed layer needs to be deposited to produce microelectronic devices with the smallest features and more functionality. Atomic layer deposition (ALD) is the most suitable method to deposit such thin films. However, the reaction mechanism and the surface chemistry of copper ALD remain unclear, which is deterring the development of better precursors and design of new ALD processes. In this thesis, we study the surface chemistries during ALD of copper by means of density functional theory (DFT). To understand the effect of temperature and pressure on the composition of copper with substrates, we used ab initio atomistic thermodynamics to obtain phase diagram of the Cu(111)/SiO2(0001) interface. We found that the interfacial oxide Cu2O phases prefer high oxygen pressure and low temperature while the silicide phases are stable at low oxygen pressure and high temperature for Cu/SiO2 interface, which is in good agreement with experimental observations. Understanding the precursor adsorption on surfaces is important for understanding the surface chemistry and reaction mechanism of the Cu ALD process. Focusing on two common Cu ALD precursors, Cu(dmap)2 and Cu(acac)2, we studied the precursor adsorption on Cu surfaces by means of van der Waals (vdW) inclusive DFT methods. We found that the adsorption energies and adsorption geometries are dependent on the adsorption sites and on the method used to include vdW in the DFT calculation. Both precursor molecules are partially decomposed and the Cu cations are partially reduced in their chemisorbed structure. It is found that clean cleavage of the ligand−metal bond is one of the requirements for selecting precursors for ALD of metals. 2 Bonding between surface and an atom in the ligand which is not coordinated with the Cu may result in impurities in the thin film. To have insight into the reaction mechanism of a full ALD cycle of Cu ALD, we proposed reaction pathways based on activation energies and reaction energies for a range of surface reactions between Cu(dmap)2 and Et2Zn. The butane formation and desorption steps are found to be extremely exothermic, explaining the ALD reaction scheme of original experimental work. Endothermic ligand diffusion and re-ordering steps may result in residual dmap ligands blocking surface sites at the end of the Et2Zn pulse, and in residual Zn being reduced and incorporated as an impurity. This may lead to very slow growth rate, as was the case in the experimental work. By investigating the reduction of CuO to metallic Cu, we elucidated the role of the reducing agent in indirect ALD of Cu. We found that CuO bulk is protected from reduction during vacuum annealing by the CuO surface and that H2 is required in order to reduce that surface, which shows that the strength of reducing agent is important to obtain fully reduced metal thin films during indirect ALD processes. Overall, in this thesis, we studied the surface chemistries and reaction mechanisms of Cu ALD processes and the nucleation of Cu to form a thin film.

Time-dependent transport through molecular junctions.

We investigate transport properties of molecular junctions under two types of bias--a short time pulse or an ac bias--by combining a solution for Green's functions in the time domain with electronic structure information coming from ab initio density functional calculations. We find that the short time response depends on lead structure, bias voltage, and barrier heights both at the molecule-lead contacts and within molecules. Under a low frequency ac bias, the electron flow either tracks or leads the bias signal (resistive or capacitive response) depending on whether the junction is perfectly conducting or not. For high frequency, the current lags the bias signal due to the kinetic inductance. The transition frequency is an intrinsic property of the junctions.

Efficient construction of nonorthogonal localized molecular orbitals in large systems.

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.

Vibrational spectroscopy and DFT calculations of di-amino acid cyclic peptides. Part I: cyclo(Gly-Gly), cyclo(L-Ala-L-Ala) and cyclo(L-Ala-Gly) in the solid state and in aqueous solution

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.

Catalytic Water Formation on Platinum: A First-Principles Study

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.

DFT and in situ EXAFS investigation of gold/ceria-zirconia low-temperature water gas shift catalysts: Identification of the nature of the active form of gold

A combined experimental and theoretical investigation of the nature of the active form of gold in oxide-supported gold catalysts for the water gas shift reaction has been performed. In situ extended X-ray absorption fine structure (EXAFS) and X-ray absorption near-edge structure (XANES) experiments have shown that in the fresh catalysts the gold is in the form of highly dispersed gold ions. However, under water gas shift reaction conditions, even at temperatures as low as 100 degrees C, the evidence from EXAFS and XANES is only 14 consistent with rapid, and essentially complete, reduction of the gold to form metallic clusters containing about 50 atoms. The presence of Au-Ce distances in the EXAFS spectra, and the fact that about 15% of the gold atoms can be reoxidized after exposure to air at 150 degrees C, is indicative of a close interaction between a fraction (ca. 15%) of the gold atoms and the oxide support. Density functional theory (DFT) calculations are entirely consistent with this model and suggest that an important aspect of the active and stable form of gold under water gas shift reaction conditions is the location of a partially oxidized gold (Audelta+) species at a cerium cation vacancy in the surface of the oxide support. It is found that even with a low loading gold catalysts (0.2%) the fraction of ionic gold under water gas shift conditions is below the limit of detection by XANES (<5%). It is concluded that under water gas shift reaction conditions the active form of gold comprises small metallic gold clusters in intimate contact with the oxide support.

On the solvation of L-aspartic acid

We use molecular statics and dynamics to study the stability of L-aspartic acid both in vacuo and solvated by polar and non-polar molecules using density functional theory in the generalized gradient approximation. We find that structures stable in vacuo are unstable in aqueous solution and vice versa. From our simulations we are able to come to some conclusions about the mechanism of stabilisation of zwitterions by polar protic solvents, water and methanol.

Bandstructure approach to near edge structure

We review the current state of the art in EELS fingerprinting by computer simulation, focusing on the bandstructure approach to the problem. Currently calculations are made using a one electron theory, but we describe in principle the way to go beyond this to include final state effects. We include these effects within the one electron framework using the Slater transition state formula and assess the errors involved. Two examples are then given which illustrate the use of the one electron approximation within density functional theory. Our approach is to combine predicted atomic structure with predicted electronic structure to assist in fingerprinting of complex crystal structures.

The near-edge structure in energy-loss spectroscopy: many-electron and magnetic effects in transition metal nitrides and carbides

We investigate the ability of the local density approximation (LDA) in density functional theory to predict the near-edge structure in electron energy-loss spectroscopy in the dipole approximation. We include screening of the core hole within the LDA using Slater's transition state theory. We find that anion K-edge threshold energies are systematically overestimated by 4.22 +/- 0.44 eV in twelve transition metal carbides and nitrides in the rock-salt (B1) structure. When we apply this 'universal' many-electron correction to energy-loss spectra calculated within the transition state approximation to LDA, we find quantitative agreement with experiment to within one or two eV for TiC, TiN and VN. We compare our calculations to a simpler approach using a projected Mulliken density which honours the dipole selection rule, in place of the dipole matrix element itself. We find remarkably close agreement between these two approaches. Finally, we show an anomaly in the near-edge structure in CrN to be due to magnetic structure. In particular, we find that the N K edge in fact probes the magnetic moments and alignments of ther sublattice.

Ab initio molecular dynamics simulation of a room temperature ionic liquid

Ab initio molecular dynamics simulations have been performed for the first time on the room-temperature organic ionic liquid dimethyl imidazolium chloride [DMIM][Cl] using density functional theory. The aim is to compare the local liquid structure with both that obtained from two different classical force fields and from neutron scattering experiments. The local structure around the cation shows significant differences compared to both the classical calculations and the neutron results. In particular, and unlike in the gas-phase ion pair, chloride ions tend to be located near a ring C-H proton in a position suggesting hydrogen bonding. The results are used to suggest ways in which the classical potentials may be improved.

Fixed-node diffusion Monte Carlo computations for closed-shell jellium spheres

Fixed-node diffusion Monte Carlo computations are used to determine the ground state energy and electron density for jellium spheres with up to N = 106 electrons and background densities corresponding to the electron gas parameter 1 less than or equal to r(s)less than or equal to5.62. We analyze the density and size dependence of the surface energy, and we extrapolate our data to the thermodynamic limit. The results agree well with the predictions of density functional computations using the local density approximation. In the case of N = 20, we extend our computation to higher densities and identify a transition between atomic- and jelliumlike nodal structures occurring at the background density corresponding to r(s)=0.13. In this case the local density approximation is unable to reproduce the changes in the correlation energy due to the discontinuous transition in the ground state nodal structure. We discuss the relevance of our results for nonlocal approximations to density functional theory.

Thin Ta films: growth, stability, and diffusion studied by molecular-dynamics simulations

We report results of classical molecular-dynamics simulations of bcc and beta-Ta thin films. Thermal PVD film growth, surface roughness, argon ion bombardment, phase stability and transformation, vacancy and adatom diffusion, and thermal relaxation kinetics are discussed. Distinct differences between the two structures are observed, including a complex vacancy diffusion mechanism in beta-Ta. Embedded atom method potentials, which were fitted to bcc properties, have been used to model the Ta-Ta interactions. In order to verify the application of these potentials to the more complex beta-Ta structure, we have also performed density functional theory calculations. Results and implications of these calculations are discussed.

Catalytic Role of Metal Oxides in Gold-Based Catalysts: A First Principles Study of CO Oxidation on TiO2 Supported Au

CO oxidation on TiO2 supported Au has been studied using density functional theory calculations. Important catalytic roles of the oxide have been identified: (i) CO oxidation occurs at the interface between Au and the oxide with a very small barrier; and (ii) O-2 adsorption at the interface is the key step in the reaction. The physical origin of the oxide promotion effect has been further investigated: The oxide enhances electron transfer from the Au to the antibonding states of O-2, giving rise to (i) strong ionic bonding between the adsorbed O-2, Au, and the Ti cation; and (ii) a significant activation of O-2 towards CO oxidation.

Molecular effects in the ionization of N-2, O-2, and F-2 by intense laser fields

In this paper we study the response in time of N2, O2, and F2 to laser pulses having a wavelength of 390 nm. We find single-ionization suppression in O2 and its absence in F2, in accordance with experimental results at lambda= 800 nm. Within our framework of time-dependent density functional theory we are able to explain deviations from the predictions of intense-field many-body S-matrix theory (IMST). We confirm the connection of ionization suppression with destructive interference of outgoing electron waves from the ionized electron orbital. However, the prediction of ionization suppression, justified within the IMST approach through the symmetry of the highest occupied molecular orbital (HOMO), is not reliable since it turns out that—e.g., in the case of F2—the electronic response to the laser pulse is rather complicated and does not lead to dominant depletion of the HOMO. Therefore, the symmetry of the HOMO is not sufficient to predict ionization suppression. However, at least for F2, the symmetry of the dominantly ionized orbital is consistent with the nonsuppression of ionization.