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.

Nanoscale electronic properties of conjugated polymers and nanocrystals

The objective of this thesis is the exploration and characterisation of the nanoscale electronic properties of conjugated polymers and nanocrystals. In Chapter 2, the first application of conducting-probe atomic force microscopy (CP-AFM)-based displacement-voltage (z-V) spectroscopy to local measurement of electronic properties of conjugated polymer thin films is reported. Charge injection thresholds along with corresponding single particle gap and exciton binding energies are determined for a poly[2-methoxy-5-(2-ethylhexyloxy)-1,4-phenylenevinylene] thin film. By performing measurements across a grid of locations on the film, a series of exciton binding energy distributions are identified. The variation in measured exciton binding energies is in contrast to the smoothness of the film suggesting that the variation may be attributable to differences in the nano-environment of the polymer molecules within the film at each measurement location. In Chapter 3, the CP-AFM-based z-V spectroscopy method is extended for the first time to local, room temperature measurements of the Coulomb blockade voltage thresholds arising from sequential single electron charging of 28 kDa Au nanocrystal arrays. The fluid-like properties of the nanocrystal arrays enable reproducible formation of nanoscale probe-array-substrate junctions, allowing the influence of background charge on the electronic properties of the array to be identified. CP-AFM also allows complementary topography and phase data to be acquired before and after spectroscopy measurements, enabling comparison of local array morphology with local measurements of the Coulomb blockade thresholds. In Chapter 4, melt-assisted template wetting is applied for the first time to massively parallel fabrication of poly-(3-hexylthiophene) nanowires. The structural characteristics of the wires are first presented. Two-terminal electrical measurements of individual nanowires, utilising a CP-AFM tip as the source electrode, are then used to obtain the intrinsic nanowire resistivity and the total nanowire-electrode contact resistance subsequently allowing single nanowire hole mobility and mean nanowire-electrode barrier height values to be estimated. In Chapter 5, solution-assisted template wetting is used for fabrication of fluorene-dithiophene co-polymer nanowires. The structural characteristics of these wires are also presented. Two-terminal electrical measurements of individual nanowires indicate barrier formation at the nanowire-electrode interfaces and measured resistivity values suggest doping of the nanowires, possibly due to air exposure. The first report of single conjugated polymer nanowires as ultra-miniature photodetectors is presented, with single wire devices yielding external quantum efficiencies ~ 0.1 % and responsivities ~ 0.4 mA/W under monochromatic illumination.

Dilute nitride semiconductors : band structure, scattering and high field transport

The substitution of a small fraction x of nitrogen atoms, for the group V elements in conventional III-V semiconductors such as GaAs and GaSb strongly perturbs the conduction band of the host semiconductor. In this thesis we investigate the effects of nitrogen states on the band dispersion, carrier scattering and mobility of dilute nitride alloys. In the supercell model we solve the single particle Hamiltonian for a very large supercell containing randomly placed nitrogen. This model predicts a gap in the density of states of GaNxAs1−x, where this gap is filled in the Green’s function model. Therefore we develop a self-consistent Green’s function (SCGF) approach, which provides excellent agreement with supercell calculations and reveals a gap in the DOS, in contrast with the results of previous non-self-consistent Green’s function calculations. However, including the distribution of N states destroys this gap, as seen in experiment. We then examine the high field transport of carriers by solving the steadystate Boltzmann transport equation and find that it is necessary to include the full distribution of N levels in order to account for the small, low-field mobility and the absence of a negative differential velocity regime observed experimentally with increasing x. Overall the results account well for a wide range of experimental data. We also investigate the band structure, scattering and mobility of carriers by finding the poles of the SCGF, which gives lower carrier mobility for GaNxAs1−x, compared to those already calculated, in better agreement with experiments. The calculated optical absorption spectra for InyGa1−yNxAs1−x and GaNxSb1−x using the SCGF agree well with the experimental data, confirming the validity of this approach to study the band structure of these materials.

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.

Theory of the electronic and optical properties of dilute bismide alloys

Dilute bismide alloys, containing small fractions of bismuth (Bi), have recently attracted interest due to their potential for applications in a range of semiconductor devices. Experiments have revealed that dilute bismide alloys such as GaBixAs1−x, in which a small fraction x of the atoms in the III-V semiconductor GaAs are replaced by Bi, exhibit a number of unusual and unique properties. For example, the band gap energy (E g) decreases rapidly with increasing Bi composition x, by up to 90 meV per % Bi replacing As in the alloy. This band gap reduction is accompanied by a strong increase in the spin-orbit-splitting energy (ΔSO) with increasing x, and both E g and ΔSO are characterised by strong, composition-dependent bowing. The existence of a ΔSO > E g regime in the GaBixAs1−x alloy has been demonstrated for x ≳10%, a band structure condition which is promising for the development of highly efficient, temperature stable semiconductor lasers that could lead to large energy savings in future optical communication networks. In addition to their potential for specific applications, dilute bismide alloys have also attracted interest from a fundamental perspective due to their unique properties. In this thesis we develop the theory of the electronic and optical properties of dilute bismide alloys. By adopting a multi-scale approach encompassing atomistic calculations of the electronic structure using the semi-empirical tight-binding method, as well as continuum calculations based on the k•p method, we develop a fundamental understanding of this unusual class of semiconductor alloys and identify general material properties which are promising for applications in semiconductor optoelectronic and photovoltaic devices. By performing detailed supercell calculations on both ordered and disordered alloys we explicitly demonstrate that Bi atoms act as isovalent impurities when incorporated in dilute quantities in III-V (In)GaAs(P) materials, strongly perturbing the electronic structure of the valence band. We identify and quantify the causes and consequences of the unusual electronic properties of GaBixAs1−x and related alloys, and our analysis is reinforced throughout by a series of detailed comparisons to the results of experimental measurements. Our k•p models of the band structure of GaBixAs1−x and related alloys, which we derive directly from detailed atomistic calculations, are ideally suited to the study of dilute bismide-based devices. We focus in the latter part of the thesis on calculations of the electronic and optical properties of dilute bismide quantum well lasers. In addition to developing an understanding of the effects of Bi incorporation on the operational characteristics of semiconductor lasers, we also present calculations which have been used explicitly in designing and optimising the first generation of GaBixAs1−x-based devices.

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.

Predicting the frequency dispersion of electronic hyperpolarizabilities on the basis of absorption data and thomas-kuhn sum rules

Successfully predicting the frequency dispersion of electronic hyperpolarizabilities is an unresolved challenge in materials science and electronic structure theory. We show that the generalized Thomas-Kuhn sum rules, combined with linear absorption data and measured hyperpolarizability at one or two frequencies, may be used to predict the entire frequency-dependent electronic hyperpolarizability spectrum. This treatment includes two- and three-level contributions that arise from the lowest two or three excited electronic state manifolds, enabling us to describe the unusual observed frequency dispersion of the dynamic hyperpolarizability in high oscillator strength M-PZn chromophores, where (porphinato)zinc(II) (PZn) and metal(II)polypyridyl (M) units are connected via an ethyne unit that aligns the high oscillator strength transition dipoles of these components in a head-to-tail arrangement. We show that some of these structures can possess very similar linear absorption spectra yet manifest dramatically different frequency dependent hyperpolarizabilities, because of three-level contributions that result from excited state-to excited state transition dipoles among charge polarized states. Importantly, this approach provides a quantitative scheme to use linear optical absorption spectra and very limited individual hyperpolarizability measurements to predict the entire frequency-dependent nonlinear optical response. Copyright © 2010 American Chemical Society.

12C2H2-Ar van der Waals complex

New theoretical and experimental results on the acetylene-Ar van der Waals complex are presented and the literature is reviewed. New ab initio calculations at the MP2 level were performed using large basis sets with diffuse functions and taking into account the basis set superposition error. It was found that the structure of acetylene is not significantly altered by the complexation and that its vibrational frequencies are only slightly lowered. Finally, it was observed that the calculated properties of the complex (structure, vibrational spectrum, bond dissociation energy) are not sensitive to the structure imposed on acetylene. Experimentally, acetylene-Ar was produced in a supersonic expansion under experimental conditions corresponding to 9 K rotational temperature. Thanks to the performances of CW-CRDS detection, the Ka = 0 ← 1, 1 ← 0, and 2 ← 1 sub- bands of the v1 + v3 band could be recorded and resolved and most of their lines assigned. Upper-state rotational constants were fitted, however not including the upper Ka = 2 state, which shows K-doubling the opposite of the expected. The Lorentzian width of most line profiles sets the mean lifetime to some 7.5 ns. Local perturbations affecting line positions and/or line widths are demonstrated. Additional series of lines tentatively attributed to acetylene-Ar are discussed.© 2009 American Chemical Society.

Accurate ground-state potential energy surfaces of the C2H2–Kr and C2H2–Xe van der Waals complexes

Accurate ab initio intermolecular potential energy surfaces (IPES) have been obtained for the first time for the ground electronic state of the C 2H2-Kr and C2H2-Xe van der Waals complexes. Extensive tests, including complete basis set and all-electron scalar relativistic results, support their calculation at the CCSD(T) level of theory, using small-core relativistic pseudopotentials for the rare-gas atoms and aug-cc-pVQZ basis sets extended with a set of 3s3p2d1f1g mid-bond functions. All results are corrected for the basis set superposition error. The importance of the scalar relativistic and rare-gas outer-core (n.1)d correlation effects is investigated. The calculated IPES, adjusted to analytical functions, are characterized by global minima corresponding to skew T-shaped geometries, in which the Jacobi vector positioning the rare-gas atom with respect to the center of mass of the C2H2 moiety corresponds to distances of 4.064 and 4.229Å, and angles of 65.22° and 68.67° for C 2H2-Kr and C2H2-Xe, respectively. The interaction energy of both complexes is estimated to be -151.88 (1.817 kJ mol-1) and -182.76 cm-1 (2.186 kJ mol-1), respectively. The evolution of the topology of the IPES as a function of the rare-gas atom, from He to Xe, is also discussed. © 2012 Taylor and Francis.

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.

Effect of relaxation on the oxygen K-edge electron energy-loss near-edge structure in yttria-stabilized zirconia

The electron energy-loss near-edge structure (ELNES) at the oxygen K-edge has been investigated in a range of yttria-stabilized zirconia (YSZ) materials. The electronic structure of the three polymorphs of pure ZrO2 and of the doped YSZ structure close to the 33 mol %Y2O3 composition have been calculated using a full-potential linear muffin-tin orbital method (NFP-LMTO) as well as a pseudopotential based technique. Calculations of the ELNES dipole transition matrix elements in the framework of the NFP-LMTO scheme and inclusion of core hole screening within Slater's transition state theory enable the ELNES to be computed. Good agreement between the experimental and calculated ELNES is obtained for pure monoclinic ZrO2. The agreement is less good with the ideal tetragonal and cubic structures. This is because the inclusion of defects is essential in the calculation of the YSZ ELNES. If the model used contains ordered defects such as vacancies and metal Y planes, agreement between the calculated and experimental O K-edges is significantly improved. The calculations show how the five different O environments of Zr,Y,O, are connected with the features observed in the experimental spectra and demonstrate clearly the power of using ELNES to probe the stabilization mechanism in doped metal oxides.

Molecular electrostatic properties of ions in an ionic liquid

We have analysed the electronic wave functions from an ab initio simulation of the ionic liquid (room temperature molten salt) dimethyl imidazolium chloride ([dmim][Cl] or [C1mim][Cl]) using localized Wannier orbitals. This allows us to assign electron density to individual ions. The probability distributions of the ionic dipole moments for an isolated ion and for ions in solution are compared. The liquid environment is found to polarize the cation by about 0.7 D and to increase the amplitude of the fluctuations in the dipole moments of both cation and anion. The relative changes in nuclear and electronic contributions are shown. The implications for classical force fields are discussed.

Probing warm dense lithium by inelastic X-ray scattering

One of the grand challenges of contemporary physics is understanding strongly interacting quantum systems comprising such diverse examples as ultracold atoms in traps, electrons in high-temperature superconductors and nuclear matter. Warm dense matter, defined by temperatures of a few electron volts and densities comparable with solids, is a complex state of such interacting matter. Moreover, the study of warm dense matter states has practical applications for controlled thermonuclear fusion, where it is encountered during the implosion phase, and it also represents laboratory analogues of astrophysical environments found in the core of planets and the crusts of old stars, Here we demonstrate how warm dense matter states can be diagnosed and structural properties can be obtained by inelastic X-ray scattering measurements on a compressed lithium sample. Combining experiments and ab initio simulations enables us to determine its microscopic state and to evaluate more approximate theoretical models for the ionic structure.

Geometric bonding effects in the X2A1, A2∑+u, and B2IIg states of Li2F

Published ab-initio and pseudopotential calculations for the dialkali halide systems suggest that the preferred co-linear geometry is for the metal to approach the metal end of the alkali halide. Here, ab-initio calculations on the Li2F system reveal that the well depth on the halide side in this radical is much deeper and is a local saddle-point associated with the ionic non-linear global minima. Although many features of the pseudopotential surfaces are confirmed, significant differences are apparent including the existence of a linear excited state instead of a triangular one, a considerably deeper global minimum some 50% lower in energy and a close approach between the X2A1 and the states, with the minimum 87 kJ mol-1 below the ground state asymptote. All the results can be rationalised as the avoided crossings between a long range, covalent potential dominant within the LiLiF geometry and an ionic state that forms the global minimum. Calculations on the 3rd 2A' potential indicate that even for Li + LiF collisions at ultracold temperatures the collision dynamics could involve as many as three electronic states.

Ultracold fluorine production via Doppler cooled BeF

Large parts of the periodic table cannot be cooled by current laser-based methods. We investigate whether zero energy fragmentation of laser cooled fluorides is a potential source of ultracold fluorine atoms. We report new ab initio calculations on the lowest electronic states of the BeF diatomic molecule including spin-orbit coupling, the calculated minima for the valence electronic states being within 1 pm of the spectroscopic values. A four colour cooling scheme based on the A2? ? X2S+ transition is shown to be feasible for this molecule. Multi-Reference Configuration Interaction (MRCI) potentials of the lowest energy Rydberg states are reported for the first time and found to be in good agreement with experimental data. A series of multi-pulse excitation schemes from a single rovibrational level of the cooled molecule are proposed to produce cold fluorine atoms.