Kinetics and coverage dependent reaction mechanisms of the copper atomic layer deposition from copper dimethylamino-2-propoxide and diethylzinc

Atomic layer deposition (ALD) has been recognized as a promising method to deposit conformal and uniform thin film of copper for future electronic devices. However, many aspects of the reaction mechanism and the surface chemistry of copper ALD remain unclear. In this paper, we employ plane wave density functional theory (DFT) to study the transmetalation ALD reaction of copper dimethylamino-2-propoxide [Cu(dmap)2] and diethylzinc [Et2Zn] that was realized experimentally by Lee et al. [ Angew. Chem., Int. Ed. 2009, 48, 4536−4539]. We find that the Cu(dmap)2 molecule adsorbs and dissociates through the scission of one or two Cu–O bonds into surface-bound dmap and Cu(dmap) fragments during the copper pulse. As Et2Zn adsorbs on the surface covered with Cu(dmap) and dmap fragments, butane formation and desorption was found to be facilitated by the surrounding ligands, which leads to one reaction mechanism, while the migration of ethyl groups to the surface leads to another reaction mechanism. During both reaction mechanisms, ligand diffusion and reordering are generally endothermic processes, which may result in residual ligands blocking the surface sites at the end of the Et2Zn pulse, and in residual Zn being reduced and incorporated as an impurity. We also find that the nearby ligands play a cooperative role in lowering the activation energy for formation and desorption of byproducts, which explains the advantage of using organometallic precursors and reducing agents in Cu ALD. The ALD growth rate estimated for the mechanism is consistent with the experimental value of 0.2 Å/cycle. The proposed reaction mechanisms provide insight into ALD processes for copper and other transition metals.

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.

Hybrid density functional theory description of N- and C-doping of NiO

The large intrinsic bandgap of NiO hinders its potential application as a photocatalyst under visible-light irradiation. In this study, we have performed first-principles screened exchange hybrid density functional theory with the HSE06 functional calculations of N- and C-doped NiO to investigate the effect of doping on the electronic structure of NiO. C-doping at an oxygen site induces gap states due to the dopant, the positions of which suggest that the top of the valence band is made up primarily of C 2p-derived states with some Ni 3d contributions, and the lowest-energy empty state is in the middle of the gap. This leads to an effective bandgap of 1.7 eV, which is of potential interest for photocatalytic applications. N-doping induces comparatively little dopant-Ni 3d interactions, but results in similar positions of dopant-induced states, i.e., the top of the valence band is made up of dopant 2p states and the lowest unoccupied state is the empty gap state derived from the dopant, leading to bandgap narrowing. With the hybrid density functional theory (DFT) results available, we discuss issues with the DFT corrected for on-site Coulomb description of these systems.

A first principles analysis of the effect of hydrogen concentration in hydrogenated amorphous silicon on the formation of strained Si-Si bonds and the optical and mobility gaps

In this paper, we use a model of hydrogenated amorphous silicon generated from molecular dynamics with density functional theory calculations to examine how the atomic geometry and the optical and mobility gaps are influenced by mild hydrogen oversaturation. The optical and mobility gaps show a volcano curve as the hydrogen content varies from undersaturation to mild oversaturation, with largest gaps obtained at the saturation hydrogen concentration. At the same time, mid-gap states associated with dangling bonds and strained Si-Si bonds disappear at saturation but reappear at mild oversaturation, which is consistent with the evolution of optical gap. The distribution of Si-Si bond distances provides the key to the change in electronic properties. In the undersaturation regime, the new electronic states in the gap arise from the presence of dangling bonds and strained Si-Si bonds, which are longer than the equilibrium Si-Si distance. Increasing hydrogen concentration up to saturation reduces the strained bonds and removes dangling bonds. In the case of mild oversaturation, the mid-gap states arise exclusively from an increase in the density of strained Si-Si bonds. Analysis of our structure shows that the extra hydrogen atoms form a bridge between neighbouring silicon atoms, thus increasing the Si-Si distance and increasing disorder in the sample.