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.

Structural manifestation of the delocalization error of density functional approximations: C(4N+2) rings and C(20) bowl, cage, and ring isomers.

The ground state structure of C(4N+2) rings is believed to exhibit a geometric transition from angle alternation (N < or = 2) to bond alternation (N > 2). All previous density functional theory (DFT) studies on these molecules have failed to reproduce this behavior by predicting either that the transition occurs at too large a ring size, or that the transition leads to a higher symmetry cumulene. Employing the recently proposed perspective of delocalization error within DFT we rationalize this failure of common density functional approximations (DFAs) and present calculations with the rCAM-B3LYP exchange-correlation functional that show an angle-to-bond-alternation transition between C(10) and C(14). The behavior exemplified here manifests itself more generally as the well known tendency of DFAs to bias toward delocalized electron distributions as favored by Huckel aromaticity, of which the C(4N+2) rings provide a quintessential example. Additional examples are the relative energies of the C(20) bowl, cage, and ring isomers; we show that the results from functionals with minimal delocalization error are in good agreement with CCSD(T) results, in contrast to other commonly used DFAs. An unbiased DFT treatment of electron delocalization is a key for reliable prediction of relative stability and hence the structures of complex molecules where many structure stabilization mechanisms exist.

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.

Theory of the near K-edge structure in electron energy loss spectroscopy

Arguments are given that lead to a formalism for calculating near K-edge structure in electron energy loss spectroscopy (EELS). This is essentially a one electron picture, while many body effects may be introduced at different levels, such as the local density approximation to density functional theory or the GW approximation to the electron self-energy. Calculations are made within the all electron LMTO scheme in crystals with complex atomic and electronic structures, and these are compared with experiment. (c) 2004 Elsevier B.V. All rights reserved.

Escape times for rigid Brownian rotators in a bistable potential from the time evolution of the Green function and the characteristic time of the probability evolution

The greatest relaxation time for an assembly of three- dimensional rigid rotators in an axially symmetric bistable potential is obtained exactly in terms of continued fractions as a sum of the zero frequency decay functions (averages of the Legendre polynomials) of the system. This is accomplished by studying the entire time evolution of the Green function (transition probability) by expanding the time dependent distribution as a Fourier series and proceeding to the zero frequency limit of the Laplace transform of that distribution. The procedure is entirely analogous to the calculation of the characteristic time of the probability evolution (the integral of the configuration space probability density function with respect to the position co-ordinate) for a particle undergoing translational diffusion in a potential; a concept originally used by Malakhov and Pankratov (Physica A 229 (1996) 109). This procedure allowed them to obtain exact solutions of the Kramers one-dimensional translational escape rate problem for piecewise parabolic potentials. The solution was accomplished by posing the problem in terms of the appropriate Sturm-Liouville equation which could be solved in terms of the parabolic cylinder functions. The method (as applied to rotational problems and posed in terms of recurrence relations for the decay functions, i.e., the Brinkman approach c.f. Blomberg, Physica A 86 (1977) 49, as opposed to the Sturm-Liouville one) demonstrates clearly that the greatest relaxation time unlike the integral relaxation time which is governed by a single decay function (albeit coupled to all the others in non-linear fashion via the underlying recurrence relation) is governed by a sum of decay functions. The method is easily generalized to multidimensional state spaces by matrix continued fraction methods allowing one to treat non-axially symmetric potentials, where the distribution function is governed by two state variables. (C) 2001 Elsevier Science B.V. All rights reserved.

Electronic structure of the nitrogen-vacancy center in diamond from first-principles theory

The nitrogen-vacancy (NV) center is a paramagnetic defect in diamond with applications as a qubit. Here, we investigate its electronic structure by using ab initio density functional theory for five different NV center models of two different cluster sizes. We describe the symmetry and energetics of the low-lying states and compare the optical frequencies obtained to experimental results. We compute the major transition of the negatively charged NV centers to within 25–100 meV accuracy and find that it is energetically favorable for substitutional nitrogens to donate an electron to NV0. The excited state of the major transition and the NV0 state with a neutral donor nitrogen are found to be close in energy.

Electronic Structure of an XUV Photogenerated Solid-Density Aluminum Plasma

By use of high intensity XUV radiation from the FLASH free-electron laser at DESY, we have created highly excited exotic states of matter in solid-density aluminum samples. The XUV intensity is sufficiently high to excite an inner-shell electron from a large fraction of the atoms in the focal region. We show that soft-x-ray emission spectroscopy measurements reveal the electronic temperature and density of this highly excited system immediately after the excitation pulse, with detailed calculations of the electronic structure, based on finite-temperature density functional theory, in good agreement with the experimental results.

Dispersion interactions in room-temperature ionic liquids: Results from a non-empirical density functional

The role of dispersion or van de Waals (VDW) interactions in imidazolium-based room-temperature ionic liquids is studied within the framework of density functional theory, using a recently developed non-empirical functional [M. Dion, H. Rydberg, E. Schroder, D. C. Langreth, and B. I. Lundqvist, Phys. Rev. Lett. 92, 246401 (2004)], as efficiently implemented in the SIESTA code [G. Roman-Perez and J. M. Soler, Phys. Rev. Lett. 103, 096102 (2009)]. We present results for the equilibrium structure and lattice parameters of several crystalline phases, finding a general improvement with respect to both the local density (LDA) and the generalized gradient approximations (GGA). Similar to other systems characterized by VDW bonding, such as rare gas and benzene dimers as well as solid argon, equilibrium distances and volumes are consistently overestimated by approximate to 7%, compared to -11% within LDA and 11% within GGA. The intramolecular geometries are retained, while the intermolecular distances and orientations are significantly improved relative to LDA and GGA. The quality is superior to that achieved with tailor-made empirical VDW corrections ad hoc [M. G. Del Popolo, C. Pinilla, and P. Ballone, J. Chem. Phys. 126, 144705 (2007)]. We also analyse the performance of an optimized version of this non-empirical functional, where the screening properties of the exchange have been tuned to reproduce high-level quantum chemical calculations [J. Klimes, D. Bowler, and A. Michaelides, J. Phys.: Condens. Matter 22, 074203 (2010)]. The results for solids are even better with volumes and geometries reproduced within 2% of experimental data. We provide some insight into the issue of polymorphism of [bmim][Cl] crystals, and we present results for the geometry and energetics of [bmim][Tf] and [mmim][Cl] neutral and charged clusters, which validate the use of empirical force fields. (C) 2011 American Institute of Physics. [doi:10.1063/1.3652897]

Electronic properties of zircon and hafnon from many-body perturbation theory

The electronic properties of zircon and hafnon, two wide-gap high-kappa materials, are investigated using many-body perturbation theory (MBPT) combined with the Wannier interpolation technique. For both materials, the calculated band structures differ from those obtained within density-functional theory and MBPT by (i) a slight displacement of the highest valence-band maximum from the Gamma point and (ii) an opening of the indirect band gap to 7.6 and 8.0 eV for zircon and hafnon, respectively. The introduction of vertex corrections in the many-body self-energy does not modify the results except for a global rigid shift of the many-body corrections.

On the required shape corrections to the local density and generalized gradient approximations to the Kohn-Sham potentials for molecular response calculations of (hyper)polarizabilities and excitation energies

It is well known that shape corrections have to be applied to the local-density (LDA) and generalized gradient (GGA) approximations to the Kohn-Sham exchange-correlation potential in order to obtain reliable response properties in time dependent density functional theory calculations. Here we demonstrate that it is an oversimplified view that these shape corrections concern primarily the asymptotic part of the potential, and that they affect only Rydberg type transitions. The performance is assessed of two shape-corrected Kohn-Sham potentials, the gradient-regulated asymptotic connection procedure applied to the Becke-Perdew potential (BP-GRAC) and the statistical averaging of (model) orbital potentials (SAOP), versus LDA and GGA potentials, in molecular response calculations of the static average polarizability alpha, the Cauchy coefficient S-4, and the static average hyperpolarizability beta. The nature of the distortions of the LDA/GGA potentials is highlighted and it is shown that they introduce many spurious excited states at too low energy which may mix with valence excited states, resulting in wrong excited state compositions. They also lead to wrong oscillator strengths and thus to a wrong spectral structure of properties like the polarizability. LDA, Becke-Lee-Yang-Parr (BLYP), and Becke-Perdew (BP) characteristically underestimate contributions to alpha and S-4 from bound Rydberg-type states and overestimate those from the continuum. Cancellation of the errors in these contributions occasionally produces fortuitously good results. The distortions of the LDA, BLYP, and BP spectra are related to the deficiencies of the LDA/GGA potentials in both the bulk and outer molecular regions. In contrast, both SAOP and BP-GRAC potentials produce high quality polarizabilities for 21 molecules and also reliable Cauchy moments and hyperpolarizabilities for the selected molecules. The analysis for the N-2 molecule shows, that both SAOP and BP-GRAC yield reliable energies omega(i) and oscillator strengths f(i) of individual excitations, so that they reproduce well the spectral structure of alpha and S-4.(C) 2002 American Institute of Physics.

Density Functional study of the photoactive yellow protein's chromophore

We have investigated the structural and electronic properties of p-coumaric acid, the chromophore of the photoactive yellow protein (PYP), by means of first-principles molecular dynamics based on density functional theory (DFT). We have studied the chromophore both in the vacuum and in an extended model which includes the nearest residues in the binding pocket of PYP, as derived from crystallographic data. We have characterized the ground state of the isolated chromophore in its protonated and deprotonated forms and computed the energy barrier involved in the trans to cis isomerization process around the carbon-carbon double bond. A comparison of the optimized structures of the chromophore in the vacuum and in the extended protein model, both in the trans (ground state of PYP in the dark) and cis (first light-activated intermediate) configuration, shows how the protein environment affects the chromophore in the first step of the photocycle. Our model gives an energy storage of 25 kcal/mol associated with the trans-to-cia photoisomerization. Finally, we have elucidated the nature of the electronic excitation relevant for the photochemistry of PYP by means of time-dependent DFT calculations.

Magnetic properties of copper hexadecaphthalocyanine (F16CuPc) thin films and powders

The structural and magnetic properties of F16CuPc thin films and powder, including x-ray diffraction (XRD), superconducting quantum interference device (SQUID) magnetometry, and theoretical modelling of exchange interactions are reported. Analysis of XRD from films, with thickness ranging between 100 and 160 nm, deposited onto Kapton and a perylene-3,4,9,10-tetracarboxylic-3,4,9,10-dianhydride (PTCDA) interlayer shows that the stacking angle (defined in the text) of the film is independent of the thickness, but that the texture is modified by both film thickness and substrate chemistry. The SQUID measurements suggest that all samples are paramagnetic, a result that is confirmed by our theoretical modelling including density functional theory calculations of one-dimensional molecular chains and Green's function perturbation theory calculations for a molecular dimer. By investigating theoretically a range of different geometries, we predict that the maximum possible exchange interaction between F16CuPc molecules is twice as large as that in unfluorinated copper-phthalocyanine (CuPc). This difference arises from the smaller intermolecular spacing in F16CuPc. Our density functional theory calculation for isolated F16CuPc molecule also shows that the energy levels of Kohn-Sham orbitals are rigidly shifted similar to 1 eV lower in F16CuPc compared to CuPc without a significant modification of the intramolecular spin physics, and that therefore the two molecules provide a suitable platform for independently varying magnetism and charge transport. 

The Role of Arg228 in the phosphorylation of galactokinase: The mechanism of GHMP kinases by QM/MM studies

GHMP kinases are a group of structurally-related small molecule kinases. They have been found in all kingdoms of life and are mostly responsible for catalysing the ATP-dependent phosphorylation of intermediary metabolites. Although the GHMP kinases are of clinical, pharmaceutical and biotechnological importance, the mechanism of GHMP-kinases is controversial. A catalytic base mechanism was suggested for mevalonate kinase that has a structural feature of the ?-phosphate of ATP close to an aspartate residue; however, for one GHMP member, homoserine kinase, where the residue acting as general base is absent, a direct phosphorylation mechanism was suggested. Furthermore, it has been proposed by some authors that all the GHMP kinases function via the direct phosphorylation mechanism. This controversy in mechanism has limited our ability to exploit these enzymes as drug targets and in biotechnology. Here the phosphorylation reaction mechanism of the human galactokinase, a member of GHMP kinase was investigated using molecular dynamics simulations and density functional theory-based QM/MM calculations (B3LYP-D/AMBER99). The reaction coordinates were localized by potential energy scan using adiabatic mapping method. Our results indicate that a highly conserved Glu174 captures Arg105 to the proximity of the a-phosphate of ATP forming a H-bond network, therefore the mobility of ATP in the large oxyanion hole is restricted. Arg228 functions to stabilize the negative charge developed at the ß,?-bridging oxygen of the ATP during bond cleavage. The reaction occurs via direct phosphorylation mechanism and the Asp186 in proximity of ATP does not directly participate in the reaction pathway. Since Arg228 is not conserved among GHMP kinases, reagents which form interactions with Arg228, and therefore can interrupt its function in phosphorylation may be developed into potential selective inhibitors for galactokinase.

Mechanistic Study of 1,3-Butadiene Formation in Acetylene Hydrogenation over the Pd-Based Catalysts Using Density Functional Calculations

Green oil, which leads to the deactivation of the catalysts used for the selective hydrogenation of acetylene, has long been observed but its formation mechanism is not fully understood. In this work, the formation of 1,3-butadiene, known to be the precursor of green oil, on both Pd(111) and Pd(211) surfaces is examined using density functional theory calculations. The pathways containing C-2 + C-2 coupling reactions as well as the corresponding hydrogenation reactions are studied in detail. Three pathways for 1,3-butadiene production, namely coupling plus hydrogenation and further hydrogenation, hydrogenation plus coupling plus hydrogenation, and a two step hydrogenation followed by coupling, are determined. By comparing the effective barriers, we identify the favored pathway on both surfaces. A general understanding toward the deactivation process of the industrial catalysts is also provided. In addition, the effects of the formation of subsurface carbon atoms as well as the Ag alloying on the 1,3-butadiene formation on Pd-based catalysts are also investigated and compared with experimental results.

Density Functional Study on Reaction Activity of Steam Decomposition Products Adsorbed on Carbonaceous Surface

A systematic theoretical study on the adsorption of steam and its thermal decomposition products on carbon both zigzag and armchair surface was performed to provide molecular-level understanding of the reaction activity of all these reactants in biomass steam gasification process. All the calculations were carried out using density functional theory (DFT) at the B3LYP/6-31+g(d,p) level. The structures of carbonaceous surfaces, all reactants and surface complexes were optimized and characterized. Based on the value of adsorption heat been obtained from the calculation, the activity of all reactants can be ordered as: O > O2 >H2 >H >OH >H2O for both zigzag and armchair surface, and the adsorption style is physisorption to water molecule and chemisorption to the other dissociated components.