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.

A density functional theory study of the a-olefin selectivity in Fischer-Tropsch synthesis

We studied the alpha-olefin selectivity in Fischer-Tropsch (FT) synthesis using density functional theory (131717) calculations. We calculated the relevant elementary steps from C-2 to C-6 species. Our results showed that the barriers of hydrogenation and dehydrogenation reactions were constant with different chain lengths, and the chemisorption energies of alpha-olefins from DFT calculations also were very similar, except for C-2 species. A simple expression of the paraffin/olefin ratio was obtained based on a kinetic model. Combining the expression of the paraffin/olefin ratio and our calculation results, experimental findings are satisfactorily explained. We found that the physical origin of the chain length dependence of paraffin/olefin ratio is the chain length dependence of both the van der Waals interaction between adsorbed alpha-olefins and metal surfaces and the entropy difference between adsorbed and gaseous alpha-olefins, and that the greater chemisorption energy of ethylene is the main reason for the abnormal ethane/ethylene ratio. (c) 2008 Elsevier Inc. All rights reserved.

Positronium-atom scattering at low energies

A pseudopotential for positronium-atom interaction, based on electron-atom and positron-atom phase shifts, is constructed, and the phase shifts for Ps-Kr and Ps-Ar scattering are calculated. This approach allows us to extend the Ps-atom cross sections, obtained previously in the impulse approximation [I. I. Fabrikant and G. F. Gribakin, Phys. Rev. Lett. 112, 243201 (2014)], to energies below the Ps ionization threshold. Although experimental data are not available in this low-energy region, our results describe well the tendency of the measured cross sections to drop with decreasing velocity at v < 1 a.u. Our results show that the effect of the Ps-atom van der Waals interaction is weak compared to the polarization interaction in electron-atom and positron-atom scattering. As a result, the Ps scattering length for both Ar and Kr is positive, and the Ramsauer-Townsend minimum is not observed for Ps scattering from these targets. This makes Ps scattering quite different from electron scattering in the low-energy region, in contrast to the intermediate energy range from the Ps ionization threshold up to v ∼ 2 a.u., where the two are similar.

Ultra-broadband Polarisers Based on Metastable Free-Standing Aligned Carbon Nanotube Membranes

A carbon nanotube free-standing linearly dichroic polariser is developed using solid-state extrusion. Membrane cohesion is experimentally and numerically demonstrated to derive from inter-tube van der Waals interactions in this family of planar metastable morphologies, controlled by the chemical vapour deposition conditions. Ultra-broadband polarisation (400 nm – 2.5 mm) is shown and corroborated by effective medium and full numerical simulations.

Estudo teórico QTAIM e DFT dos compostos de coordenação: efeito quelato, titanocenos e ligação química

This work is a study of coordination compounds by quantum theory of atoms in molecules (QTAIM), based on the topological analysis of the electron density of molecular systems, both theoretically and experimentally obtained. The coordination chemistry topics which were studied are the chelate effect, bent titanocene and chemical bond in coordination complexes. The chelate effect was investigated according to topological and thermodynamic parameters. The exchange of monodentate ligands on polydentate ligands from same transition metal increases the stability of the complex both from entropy and enthalpy contributions. In some cases, the latter had a higher contribution to the stability of the complex in comparison with entropy. This enthalpic contribution is explained according to topological analysis of the M-ligand bonds where polidentate complex had higher values of electron density of bond critical point, Laplacian of electron density of bond critical point and delocalization index (number of shared electrons between two atoms). In the second chapter, was studied bent titanocenes with bulky cyclopentadienyl derivative π-ligand. The topological study showed the presence of secondary interactions between the atoms of π-ligands or between atoms of π-ligand and -ligand. It was found that, in the case of titanocenes with small difference in point group symmetry and with bulky ligands, there was an nearly linear relationship between stability and delocalization index involving the ring carbon atoms (Cp) and the titanium. However, the titanocene stability is not only related to the interaction between Ti and C atoms of Cp ring, but secondary interactions also play important role on the stability of voluminous titanocenes. The third chapter deals with the chemical bond in coordination compounds by means of QTAIM. The quantum theory of atoms in molecules so far classifies bonds and chemical interactions in two categories: closed shell interaction (ionic bond, hydrogen bond, van der Waals interaction, etc) and shared interaction (covalent bond). Based on topological parameters such as electron density, Laplacian of electron density, delocalization index, among others, was classified the chemical bond in coordination compounds as an intermediate between closed shell and shared interactions

Structural properties and energetics of diffuse Rb-87 clusters in three-dimension

A correlated two-body basis function is used to describe the three-dimensional bosonic clusters interacting via two-body van der Waals potential. We calculate the ground state and the zero orbital angular momentum excited states for Rb-N clusters with up to N = 40. We solve the many-particle Schrodinger equation by potential harmonics expansion method, which keeps all possible two-body correlations in the calculation and determines the lowest effective many-body potential. We study energetics and structural properties for such diffuse clusters both at dimer and tuned scattering length. The motivation of the present study is to investigate the possibility of formation of N-body clusters interacting through the van der Waals interaction. We also compare the system with the well studied He, Ne, and Ar clusters. We also calculate correlation properties and observe the generalised Tjon line for large cluster. We test the validity of the shape-independent potential in the calculation of the ground state energy of such diffuse cluster. These are the first such calculations reported for Rb clusters. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4730972]

Modeling and simulation of fibrinogen and its adsorption behavior

In this thesis different approaches for the modeling and simulation of the blood protein fibrinogen are presented. The approaches are meant to systematically connect the multiple time and length scales involved in the dynamics of fibrinogen in solution and at inorganic surfaces. The first part of the thesis will cover simulations of fibrinogen on an all atom level. Simulations of the fibrinogen protomer and dimer are performed in explicit solvent to characterize the dynamics of fibrinogen in solution. These simulations reveal an unexpectedly large and fast bending motion that is facilitated by molecular hinges located in the coiled-coil region of fibrinogen. This behavior is characterized by a bending and a dihedral angle and the distribution of these angles is measured. As a consequence of the atomistic detail of the simulations it is possible to illuminate small scale behavior in the binding pockets of fibrinogen that hints at a previously unknown allosteric effect. In a second step atomistic simulations of the fibrinogen protomer are performed at graphite and mica surfaces to investigate initial adsorption stages. These simulations highlight the different adsorption mechanisms at the hydrophobic graphite surface and the charged, hydrophilic mica surface. It is found that the initial adsorption happens in a preferred orientation on mica. Many effects of practical interest involve aggregates of many fibrinogen molecules. To investigate such systems, time and length scales need to be simulated that are not attainable in atomistic simulations. It is therefore necessary to develop lower resolution models of fibrinogen. This is done in the second part of the thesis. First a systematically coarse grained model is derived and parametrized based on the atomistic simulations of the first part. In this model the fibrinogen molecule is represented by 45 beads instead of nearly 31,000 atoms. The intra-molecular interactions of the beads are modeled as a heterogeneous elastic network while inter-molecular interactions are assumed to be a combination of electrostatic and van der Waals interaction. A method is presented that determines the charges assigned to beads by matching the electrostatic potential in the atomistic simulation. Lastly a phenomenological model is developed that represents fibrinogen by five beads connected by rigid rods with two hinges. This model only captures the large scale dynamics in the atomistic simulations but can shed light on experimental observations of fibrinogen conformations at inorganic surfaces.

Electrostatically mediated specific adsorption of small molecules in metallo-organic frameworks

We investigate the interaction of ethylene and ethane with a Cu-tricarboxylate complex and show that at low loadings the lighter molecule has a higher binding energy as a result of an increased interaction with the framework Cu and stronger hydrogen bonding with the basic framework oxygens. This leads to selective adsorption of ethylene by a factor of about 2 at low pressure, which is overcome by the stronger van der Waals interaction of ethane at high loadings, explaining recent literature data. The results suggest the possibility of separation of light hydrocarbons at low pressures or in trace amounts.

Structural and Electrical Investigation of C 60 –Graphene Vertical Heterostructures

Graphene, with its unique electronic and structural qualities, has become an important playground for studying adsorption and assembly of various materials including organic molecules. Moreover, organic/graphene vertical structures assembled by van der Waals interaction have potential for multifunctional device applications. Here, we investigate structural and electrical properties of vertical heterostructures composed of C60 thin film on graphene. The assembled film structure of C60 on graphene is investigated using transmission electron microscopy, which reveals a uniform morphology of C60 film on graphene with a grain size as large as 500 nm. The strong epitaxial relations between C60 crystal and graphene lattice directions are found, and van der Waals ab initio calculations support the observed phenomena. Moreover, using C60-graphene heterostructures, we fabricate vertical graphene transistors incorporating n-type organic semiconducting materials with an on/off ratio above 3 × 10(3). Our work demonstrates that graphene can serve as an excellent substrate for assembly of molecules, and attained organic/graphene heterostructures have great potential for electronics applications.

Transmission Spectra of Rb 87 atoms near an Optical Nanofiber.

We present measurements of the transmission spectra of 87Rb atoms at 780 nm in the vicinity of a nanofiber. A uniform distribution of fixed atoms around a nanofiber should produce a spectrum that is broadened towards the red due to shifts from the van der Waals potential. If the atoms are free, this also produces an attractive force that accelerates them until they collide with the fiber which depletes the steady-state density of near-surface atoms. It is for this reason that measurements of the van der Waals interaction are sparse. We confirm this by measuring the spectrum cold atoms from a magneto-optical trap around the fiber, revealing a symmetric line shape with nearly the natural linewidth of the transition. When we use an auxiliary 750 nm laser we are able to controllably desorb a steady flux of atoms from the fiber that reside near the surface (less than 50 nm) long enough to feel the van der Walls interaction and produce an asymmetric spectrum. We quantify the spectral asymmetry as a function of 750 nm laser power and find a maximum. Our model, which that takes into account the change in the density distribution, qualitatively explains the observations. In the future this can be used as a tool to more comprehensively study atom-surface interactions.

Interaction of nucleic acids with carbon nanotubes and dendrimers

Nucleic acid interaction with nanoscale objects like carbon nanotubes (CNTs) and dendrimers is of fundamental interest because of their potential application in CNT separation, gene therapy and antisense therapy. Combining nucleic acids with CNTs and dendrimers also opens the door towards controllable self-assembly to generate various supra-molecular and nano-structures with desired morphologies. The interaction between these nanoscale objects also serve as a model system for studying DNA compaction, which is a fundamental process in chromatin organization. By using fully atomistic simulations, here we report various aspects of the interactions and binding modes of DNA and small interfering RNA (siRNA) with CNTs, graphene and dendrimers. Our results give a microscopic picture and mechanism of the adsorption of single- and double-strand DNA (ssDNA and dsDNA) on CNT and graphene. The nucleic acid-CNT interaction is dominated by the dispersive van der Waals (vdW) interaction. In contrast, the complexation of DNA (both ssDNA and dsDNA) and siRNA with various generations of poly-amido-amine (PAMAM) dendrimers is governed by electrostatic interactions. Our results reveal that both the DNA and siRNA form stable complex with the PAMAM dendrimer at a physiological pH when the dendrimer is positively charged due to the protonation of the primary amines. The size and binding energy of the complex increase with increase in dendrimer generation. We also give a summary of the current status in these fields and discuss future prospects.

Non-Covalent Interaction: Revealed by Rotational Spectroscopy

The pulsed jet Fourier transform microwave spectroscopy have been applied to several molecular complexes involving H2O, freons, methane, carboxylic acids, and rare gas. The obtained results showcase the suitability of this technique for studying the intermolecular interactions. The rotational spectra of three water adducts of halogenated organic molecules, i.e. chlorotrifluoroethylene, isoflurane and alfa,alfa,alfa,-trifluoroanisole, have been investigated. It has been found that, the halogenation of the partner molecules definitely changes the way in which water will link to the partner molecule. Quadrupole hyperfine structures and/or the tunneling splittings have been observed in the rotational spectra of difluoromethane-dichloromethane, chlorotrifluorometane-fluoromethane, difluoromethane-formaldehyde and trifluoromethane-benzene. These features have been useful to describe their intermolecular interactions (weak hydrogen bonds or halogen bonds), and to size the potential energy surfaces of their internal motions. The rotational spectrum of pyridine-methane pointed out that methane prefers to locate above the ring and link to pyridine through a C-H•••π weak hydrogen bond, rather than the C-H•••n interaction. This behavior, typical of complexes of pyridine with rare gases, suggests classifying CH4, in relation to its ability to form molecular complexes with aromatic molecules, as a pseudo rare gas. The conformational equilibria of three bi-molecules of carboxylic acids, acrylic acid-trifluoroacetic acid, difluoroacetic acid-formic acid and acrylic acid-fluoroacetic acid have been studied. The increase of the hydrogen bond length upon H→D isotopic substitution (Ubbelohde effect) has been deduced from the elongation of the carboxylic carbons C•••C distance. The van der Waals complex tetrahydrofuran-krypton shows that the systematic doubling of the rotational lines has been attributed to the residual pseudo-rotation of tetrahydrofuran in the complex, based on the values of the Coriolis coupling constants, and on the type (mu_b) of the interstate transitions.

Pentafluorophenyl phenyl interaction in biphenyl DNA

We prepared and investigated oligonucleotide duplexes of the sequence d(GATGAC(X)(n)GCTAG)d(CTAGC(Y)(n)GTCATC), in which X and Y designate biphenyl- (bph) and pentafluorobiphenyl- ((5F)bph) C-nucleotides, respectively, and n varies from 0-4. These hydrophobic base substitutes are expected to adopt a zipperlike, interstrand stacking motif, in which not only bph/bph or (5F)bph/(5F)bph homo pairs, but also (5F)bph/bph mixed pairs can be formed. By performing UV-melting curve analysis we found that incorporation of a single (5F)bph/(5F)bph pair leads to a duplex that is essentially as stable as the unmodified duplex (n=0), and 2.4 K more stable than the duplex with the nonfluorinated bph/bph pair. The T(m) of the mixed bph/(5F)bph pair was in between the T(m) values of the respective homo pairs. Additional, unnatural aromatic pairs increased the T(m) by +3.0-4.4 K/couple, irrespective of the nature of the aromatic residue. A thermodynamic analysis using isothermal titration calorimetry (ITC) of a series of duplexes with n=3 revealed lower (less negative) duplex formation enthalpies (DeltaH) in the (5F)bph/(5F)bph case than in the bph/bph case, and confirmed the higher thermodynamic stability (DeltaG) of the fluorinated duplex, suggesting it to be of entropic origin. Our data are compatible with a model in which the stacking of (5F)bph versus bph is dominated by dehydration of the aromatic units upon duplex formation. They do not support a model in which van der Waals dispersive forces (induced dipoles) or electrostatic (quadrupole) interactions play a dominant role

Structural details of an interaction between cardiolipin and an integral membrane protein

Anionic lipids play a variety of key roles in biomembrane function, including providing the immediate environment for the integral membrane proteins that catalyze photosynthetic and respiratory energy transduction. Little is known about the molecular basis of these lipid–protein interactions. In this study, x-ray crystallography has been used to examine the structural details of an interaction between cardiolipin and the photoreaction center, a key light-driven electron transfer protein complex found in the cytoplasmic membrane of photosynthetic bacteria. X-ray diffraction data collected over the resolution range 30.0–2.1 Å show that binding of the lipid to the protein involves a combination of ionic interactions between the protein and the lipid headgroup and van der Waals interactions between the lipid tails and the electroneutral intramembrane surface of the protein. In the headgroup region, ionic interactions involve polar groups of a number of residues, the protein backbone, and bound water molecules. The lipid tails sit along largely hydrophobic grooves in the irregular surface of the protein. In addition to providing new information on the immediate lipid environment of a key integral membrane protein, this study provides the first, to our knowledge, high-resolution x-ray crystal structure for cardiolipin. The possible significance of this interaction between an integral membrane protein and cardiolipin is considered.