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.

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.

First principles calculation of electron-phonon and alloy scattering in strained SiGe

First-principles electronic structure methods are used to predict the mobility of n-type carrier scattering in strained SiGe. We consider the effects of strain on the electron-phonon deformation potentials and the alloy scattering parameters. We calculate the electron-phonon matrix elements and fit them up to second order in strain. We find, as expected, that the main effect of strain on mobility comes from the breaking of the degeneracy of the six Δ and L valleys, and the choice of transport direction. The non-linear effects on the electron-phonon coupling of the Δ valley due to shear strain are found to reduce the mobility of Si-like SiGe by 50% per % strain. We find increases in mobility between 2 and 11 times that of unstrained SiGe for certain fixed Ge compositions, which should enhance the thermoelectric figure of merit in the same order, and could be important for piezoresistive applications.

First-principles calculation of carrier-phonon scattering in n-type Si1−xGex alloys

First-principles electronic structure methods are used to find the rates of inelastic intravalley and intervalley n-type carrier scattering in Si1-xGex alloys. Scattering parameters for all relevant Delta and L intra- and intervalley scattering are calculated. The short-wavelength acoustic and the optical phonon modes in the alloy are computed using the random mass approximation, with interatomic forces calculated in the virtual crystal approximation using density functional perturbation theory. Optical phonon and intervalley scattering matrix elements are calculated from these modes of the disordered alloy. It is found that alloy disorder has only a small effect on the overall inelastic intervalley scattering rate at room temperature. Intravalley acoustic scattering rates are calculated within the deformation potential approximation. The acoustic deformation potentials are found directly and the range of validity of the deformation potential approximation verified in long-wavelength frozen phonon calculations. Details of the calculation of elastic alloy scattering rates presented in an earlier paper are also given. Elastic alloy disorder scattering is found to dominate over inelastic scattering, except for almost pure silicon (x approximate to 0) or almost pure germanium (x approximate to 1), where acoustic phonon scattering is predominant. The n-type carrier mobility, calculated from the total (elastic plus inelastic) scattering rate, using the Boltzmann transport equation in the relaxation time approximation, is in excellent agreement with experiments on bulk, unstrained alloys..

First-principles calculation of P-type alloy scattering in Si1-xGex

The p-type carrier scattering rate due to alloy disorder in Si1-xGex alloys is obtained from first principles. The required alloy scattering matrix elements are calculated from the energy splitting of the valence bands, which arise when one average host atom is replaced by a Ge or Si atom in supercells containing up to 128 atoms. Alloy scattering within the valence bands is found to be characterized by a single scattering parameter. The hole mobility is calculated from the scattering rate using the Boltzmann transport equation in the relaxation time approximation. The results are in good agreement with experiments on bulk, unstrained alloys..

First-principles calculation of alloy scattering in GexSi1-x

First-principles electronic structure methods are used to find the rates of intravalley and intervalley n-type carrier scattering due to alloy disorder in Si1-xGex alloys. The required alloy scattering matrix elements are calculated from the energy splitting of nearly degenerate Bloch states which arises when one average host atom is replaced by a Ge or Si atom in supercells containing up to 128 atoms. Scattering parameters for all relevant Delta and L intravalley and intervalley alloy scattering are calculated. Atomic relaxation is found to have a substantial effect on the scattering parameters. f-type intervalley scattering between Delta valleys is found to be comparable to other scattering channels. The n-type carrier mobility, calculated from the scattering rate using the Boltzmann transport equation in the relaxation time approximation, is in excellent agreement with experiments on bulk, unstrained alloys.