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

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

Computational modeling of defects in nanoscale device materials

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

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.

Accurate localized resolution of identity approach for linear-scaling hybrid density functionals and for many-body perturbation theory

© 2015 IOP Publishing Ltd and Deutsche Physikalische Gesellschaft.A key component in calculations of exchange and correlation energies is the Coulomb operator, which requires the evaluation of two-electron integrals. For localized basis sets, these four-center integrals are most efficiently evaluated with the resolution of identity (RI) technique, which expands basis-function products in an auxiliary basis. In this work we show the practical applicability of a localized RI-variant ('RI-LVL'), which expands products of basis functions only in the subset of those auxiliary basis functions which are located at the same atoms as the basis functions. We demonstrate the accuracy of RI-LVL for Hartree-Fock calculations, for the PBE0 hybrid density functional, as well as for RPA and MP2 perturbation theory. Molecular test sets used include the S22 set of weakly interacting molecules, the G3 test set, as well as the G2-1 and BH76 test sets, and heavy elements including titanium dioxide, copper and gold clusters. Our RI-LVL implementation paves the way for linear-scaling RI-based hybrid functional calculations for large systems and for all-electron many-body perturbation theory with significantly reduced computational and memory cost.

DFT studies of long-chain FAMEs: theoretical justification for determining chain length and unsaturation from experimental Raman spectra

Density functional calculations, using B3LPY/6-31G(d) methods, have been used to investigate the conformations and vibrational (Raman) spectra of a series of long-chain, saturated fatty acid methyl esters (FAMEs) with the formula CH2nO2 (n = 5-21) and two series of unsaturated FAMEs. The calculations showed that the lowest energy conformer within the saturated FAMEs is the simple (all-trans) structure and, in general, it was possible to reproduce experimental data using calculations on only the all-trans conformer. The only exception was C6H12O2, where a second low-lying conformer had to be included in order to correctly simulate the experimental Raman spectrum. The objective of the work was to provide theoretical justification for the methods that are commonly used to determine the properties of the fats and oils, such as chain length and degree of unsaturation, from experimental Raman data. Here it is shown that the calculations reproduce the trends and calibration curves that are found experimentally and also allow the reasons for the failure of what would appear to be rational measurements to be understood. This work shows that although the assumption that each FAME can simply be treated as a collection of functional groups can be justified in some cases, many of the vibrational modes are complex motions of large sections of the molecules and thus would not be expected to show simple linear trends with changes in structure, such as increasing chain length and/or unsaturation. Simple linear trends obtained from experimental data may thus arise from cancellation of opposing effects, rather than reflecting an underlying simplicity.

Conformations, vibrational frequencies and Raman intensities of short-chain fatty acid methyl esters using DFT with 6-31G(d) and Sadlej pVTZ basis sets

Density functional calculations, using B3LPY/6-31G(d) methods, have been used to investigate the conformations and vibrational (Raman) spectra of three short-chain fatty acid methyl esters (FAMEs) with the formula CnH2nO2 (n = 3-5). In all three FAMEs, the lowest energy conformer has a simple 'all-trans' structure but there are other conformers, with different torsions about the backbone, which lie reasonably close in energy to the global minimum. One result of this is that the solid samples we studied do not appear to consist entirely of the lowest energy conformer. Indeed, to account for the 'extra' bands that were observed in the Raman data but were not predicted for the all-trans conformer, it was necessary to add-in contributions from other conformers before a complete set of vibrational assignments could be made. Provided this was done, the agreement between experimental Raman frequencies and 6-31G(d) values (after scaling) was excellent, RSD = 12.6 cm(-1). However, the agreement between predicted and observed intensities was much less satisfactory. To confirm the validity of the approach followed by the 6-3 1 G(d) basis set, we used a larger basis set, Sadlej pVTZ, and found that these calculations gave accurate Raman intensities and simulated spectra (summed from two different conformers) that were in quantitative agreement with experiment. In addition, the unscaled Sadlej pVTZ, and the scaled 6-3 1 G(d) calculations gave the same vibrational mode assignments for all bands in the experimental data. This work provides the foundation for calculations on longer-chain FAMEs (which are closer to those found as triglycerides in edible fats and oils) because it shows that scaled 6-3 1 G(d) calculations give equally accurate frequency predictions, and the same vibrational mode assignments, as the much more CPU-expensive Sadlej pVTZ basis set calculations.

Boron in copper: A perfect misfit in the bulk and cohesion enhancer at a grain boundary

Using first principles electronic structure methods, we calculate the effects of boron impurities in bulk copper and at surfaces and grain boundaries. We find that boron segregation to the Sigma5(310)[001] grain boundary should strengthen the boundary up to 1.5 ML coverage (15.24 at./nm2). The maximal effect is observed at 0.5 ML and corresponds to boron atoms filling exclusively grain boundary interstices. In copper bulk, B causes significant distortion both in interstitial and regular lattice sites, for which boron atoms are either too big or too small. The distortion is compensated to a large extent when the interstitial and substitutional boron combine together to form a strongly bound dumbbell. Our prediction is that bound boron impurities should appear in a sizable proportion if not dominate in most experimental conditions. A large discrepancy between calculated heats of solution and experimental terminal solubility of B in Cu is found, indicating either a significant failure of the density functional approach or, more likely, strongly overestimated solubility limits in the existing B-Cu phase diagram.

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.

Lowest energy potentials of Li2 + F accessible at ultracold temperatures

Ab initio calculations for the strongly exoergic Li2+F harpoon reaction are presented using density-functional theory, complete active space self-consistent field, and multireference configuration interaction methods to argue that this reaction would be an ideal candidate for investigation with ultracold molecules. The lowest six states are calculated with the aug-correlation-consistent polarized valence triple-zeta basis set and at least two can be accessed by a ground rovibronic Li2 molecule with zero collision energy at all reaction geometries. The large reactive cross section (characteristic of harpoon reactions) and chemiluminescent products are additional attractive features of these reactions.

An understanding of chemoselective hydrogenation on crotonaldehyde over Pt(111) in the free energy landscape: The microkinetics study based on first-principles calculations

Microkinetic model is developed in the free energy landscape based on density functional theory (DFT) to quantitatively investigate the reaction mechanism of chemoselective partial hydrogenation of crotonaldehyde to crotyl alcohol over Pt(1 1 1) at the temperature of 353 K. Three different methods (mobile, immobile and collision theory models) were carried out to obtain free energy barrier of adsorption/desorption processes. The results from mobile and collision theory models are similar. The calculated TOFs from both models are close to the experiment value. However, for the immobile model, in which the free energy barrier of desorption approaches the energy barrier, the calculated TOF is 2 orders of magnitude lower than the other models. The difficulty of adsorption/ desorption may be overestimated in the immobile model. In addition, detailed analyses show that for the surface hydrogenation elementary steps, the entropy and internal energy effects are small under the reaction condition, while the zero-point-energy (ZPE) correction is significant, especially for the multi-step hydrogenation reaction. The total energy with the ZPE correction approaches to the full free energy calculation for the surface reaction under the reaction condition. (c) 2011 Elsevier B.V. All rights reserved.

Ab initio potentials of F + Li-2 accessible at ultracold temperatures

Ab initio calculations for the strongly exoergic Li-2 + F harpoon reaction are presented using density-functional theory, complete active space self-consistent field, and multireference configuration interaction methods to argue that this reaction would be an ideal candidate for investigation with ultracold molecules. The lowest six states are calculated with the aug-correlation-consistent polarized valence triple-zeta basis set and at least two can be accessed by a ground rovibronic Li-2 molecule with zero collision energy at all reaction geometries. The large reactive cross section (characteristic of harpoon reactions) and chemiluminescent products are additional attractive features of these reactions.

Electronic stopping power of H and He in Al and LiF from first principles

Non-linearities in the electronic stopping power of light projectiles in bulk Al and LiF are addressed from first principles using time-evolving time-dependent density functional theory. In the case of Al, the agreement of the calculations with experiments for H and He projectiles is fair, but a recently observed transition for He from one value of the electronic friction coefficient to a higher value at v ~ 0.3 a.u. is not reproduced by the calculations. For LiF, better accuracy is obtained as compared with previously published simulations, albeit the threshold remains overestimated. © 2013 Elsevier B.V.

Implementation and testing of Lanczos-based algorithms for Random-Phase Approximation eigenproblems

The treatment of the Random-Phase Approximation Hamiltonians, encountered in different frameworks, like time-dependent density functional theory or Bethe-Salpeter equation, is complicated by their non-Hermicity. Compared to their Hermitian Hamiltonian counterparts, computational methods for the treatment of non-Hermitian Hamiltonians are often less efficient and less stable, sometimes leading to the breakdown of the method. Recently [Gruning et al. Nano Lett. 8 (2009) 28201, we have identified that such Hamiltonians are usually pseudo-Hermitian. Exploiting this property, we have implemented an algorithm of the Lanczos type for Random-Phase Approximation Hamiltonians that benefits from the same stability and computational load as its Hermitian counterpart, and applied it to the study of the optical response of carbon nanotubes. We present here the related theoretical grounds and technical details, and study the performance of the algorithm for the calculation of the optical absorption of a molecule within the Bethe-Salpeter equation framework. (C) 2011 Elsevier B.V. All rights reserved.

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.