Electronic structures of 5d transition metal monoxides by density functional theory

Bond distances, vibrational frequencies, electron affinities, ionization potentials, and dissociation energies of the diatomic 5d transition metal (except La) monoxides and their positively and negatively charged ions were studied by use of density functional methods B3LYP, BLYP, B3PW91, BPW91, B3P86, BP86, MPW1PW91, PBE1PBE, and SVWN. Our calculation shows that for each individual species, the calculated properties are quite sensitive to the method used. Compared with hybrid density functional method B3PW91 (B3P86), pure density functional method BPW91 (BP86) gives longer bond distance (lower vibrational frequency) from HfO to PtO for neutral species, HfO+ to IrO+ for cationic species, and HfO- to AuO- for anionic species. While for B3LYP and BLYP, the trend was observed for cationic species from HfO+ to IrO+ and anionic species from HfO- to AuO- (except TaO-), but not for neutrals. Pure density function methods BLYP, BPW91, and BP86 give larger dissociation energy compared with hybrid density functional methods B3LYP, B3PW91, and B3P86. SVWN in most cases gives the smallest bond distance, while BLYP gives the largest value. MPW1PW91 and PBE1PBE show the same performance in predicting the spectroscopic constants.

Density-functional study of lanthanum, ytterbium, and lutetium dimers

La-2, Yb-2, and Lu-2 have been studied by use of the density-functional methods B3LYP, BLYP, B3PW91, BHLYP, BP86, B3P86, MPW1PW91, and PBE1PBE. In these density-functional methods, the exchange functional is from either Becke's three-parameter HF-DFT hybrid exchange functional (B3), pure DFT exchange functional of 1988 (B), a modification of the half-and-half HF/DFT hybrid method (BH), Perdew-Wang 1991 (PW91), or Barone's modified PW91 (MPW1), while the correlation functional is from either Lee, Yang, and Parr (LYP), Perdew-Wang 1991 (PW91), or Perdew 86 (P86). PBE1PBE is the generalized-gradient-approximation exchange-correlation functional of Perdew, Burke, and Ernzerhof. For La-2, the calculated bond distance is in reasonable agreement with the experiment, but the calculated vibrational frequency is underestimated significantly compared with the experiment. Only BP86 and B3P86 have the best performance in reproducing the experimental dissociation energy for La-2. For the van der Waals dimer Yb-2, three functionals, B3LYP, BLYP, and BHLYP have excellent performance in reproducing the spectroscopic constants compared with both the experiment and previous theoretical studies.

Addition reactions of nitrones on the reconstructed C(100)-2 x 1 surfaces

In this paper, the reactions of nitrone, N-methyl nitrone, N-phenyl nitrone and their hydroxylamine tautomers (vinyl-hydroxylamine, N-methyl-vinyl-hydroxylamine and N-phenyl-vinyl-hydroxylamine) on the reconstructed C(100)-2 x 1 surface have been investigated using hybrid density functional theory (B3LYP), Moller-Plesset second-order perturbation (MP2) and multi-configuration complete-active-space self-consistent-field (CASSCF) methods. The calculations showed that all the nitrones can react with the surface "dimer" via facile 1.3-dipolar cycloaddition with small activation barriers (less than 12.0 kJ/mol at B3LYP/6-31g(d) level). The [2+2] cycloaddition of hydroxylamine tautomers on the C(100) surface follows a diradical mechanism. Hydroxylamine tautomers first form diradical intermediates with the reconstructed C(I 00)-2 x I surface by overcoming a large activation barrier of 50-60 kJ/mol (B3LYP), then generate [2+2] cycloaddition products via diradical transition states with negligible activation barriers. The surface reactions result in hydroxyl or amino-terminated diamond surfaces, which offers new opportunity for further modifications. (C) 2007 Elsevier B.V. All rights reserved.

Theoretical study on adsorption of Na+ and Na+(H2O)(n) (n=1-6) on a clean Si(111) surface

To model the adsorption of Na+ in aqueous solution on the semiconductor surface, the interactions of Na+ and Na+(H2O)(n) (n = 1-6) with a clean Si(111) surface were investigated by using hybrid density functional theory (B3LYP) and Moller-Plesset second-order perturbation (MP2) methods. The Si(111) surface was described with Si8H12, Si16H20, and Si22H21 Cluster models. The effect of the basis set superposition error (BSSE) was taken into account by applying the counterpoise (CP) correction. The calculated results indicated that the interactions between the Na+ cation and the dangling bonds of the Si(111) surface are primarily electrostatic with partial orbital interactions. The magnitude of the binding energies depends weakly on the adsorption sites and the size of the clusters. When water molecules are present, the interaction between the Nal and Si(I 11) surfaces weakens and the binding energy has the tendency to saturate. On a Si22H21 cluster described surface, the optimized Na+-surface distance for Na+(H2O)(5) adsorbed at on-top site is 4.16 angstrom and the CP-corrected binding energy (MP2) is -35.4 kJ/mol, implying a weakly adsorption of hydrated Na+ cation on clean Si(111) surface.

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.

Structure maps for hcp metals from first-principles calculations

The ability to predict the existence and crystal type of ordered structures of materials from their components is a major challenge of current materials research. Empirical methods use experimental data to construct structure maps and make predictions based on clustering of simple physical parameters. Their usefulness depends on the availability of reliable data over the entire parameter space. Recent development of high-throughput methods opens the possibility to enhance these empirical structure maps by ab initio calculations in regions of the parameter space where the experimental evidence is lacking or not well characterized. In this paper we construct enhanced maps for the binary alloys of hcp metals, where the experimental data leaves large regions of poorly characterized systems believed to be phase separating. In these enhanced maps, the clusters of noncompound-forming systems are much smaller than indicated by the empirical results alone. © 2010 The American Physical Society.

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.

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.

The RMT method for many-electron atomic systems in intense short-pulse laser light

We describe a new ab initio method for solving the time-dependent Schrödinger equation for multi-electron atomic systems exposed to intense short-pulse laser light. We call the method the R-matrix with time-dependence (RMT) method. Our starting point is a finite-difference numerical integrator (HELIUM), which has proved successful at describing few-electron atoms and atomic ions in strong laser fields with high accuracy. By exploiting the R-matrix division-of-space concept, we bring together a numerical method most appropriate to the multi-electron finite inner region (R-matrix basis set) and a different numerical method most appropriate to the one-electron outer region (finite difference). In order to exploit massively parallel supercomputers efficiently, we time-propagate the wavefunction in both regions by employing Arnoldi methods, originally developed for HELIUM.

Prospects for ultracold carbon via charge exchange reactions and laser cooled carbides

Strategies to produce an ultracold sample of carbon atoms are explored and assessed with the help of quantum chemistry. After a brief discussion of the experimental difficulties using conventional methods, two strategies are investigated. The first attempts to exploit charge exchange reactions between ultracold metal atoms and sympathetically cooled C+ ions. Ab initio calculations including electron correlation have been conducted on the molecular ions [LiC]+ and [BeC]+ to determine whether alkali or alkaline earth metals are a suitable buffer gas for the formation of C atoms but strong spontaneous radiative charge exchange ensure they are not ideal. The second technique involves the stimulated production of ultracold C atoms from a gas of laser cooled carbides. Calculations on LiC suggest that the alkali carbides are not suitable but the CH radical is a possible laser cooling candidate thanks to very favourable Frank-Condon factors. A scheme based on a four pulse STIRAP excitation pathway to a Feshbach resonance is outlined for the production of atomic fragments with near zero centre of mass velocity.

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.

Strong electronic correlation in the hydrogen chain: A variational Monte Carlo study

In this paper, we report a fully ab initio variational Monte Carlo study of the linear and periodic chain of hydrogen atoms, a prototype system providing the simplest example of strong electronic correlation in low dimensions. In particular, we prove that numerical accuracy comparable to that of benchmark density-matrix renormalization-group calculations can be achieved by using a highly correlated Jastrow-antisymmetrized geminal power variational wave function. Furthermore, by using the so-called "modern theory of polarization" and by studying the spin-spin and dimer-dimer correlations functions, we have characterized in detail the crossover between the weakly and strongly correlated regimes of this atomic chain. Our results show that variational Monte Carlo provides an accurate and flexible alternative to highly correlated methods of quantum chemistry which, at variance with these methods, can be also applied to a strongly correlated solid in low dimensions close to a crossover or a phase transition.

Universal tight binding model for chemical reactions in solution and at surfaces. II. Water

A revised water model intended for use in condensed phase simulations in the framework of the self consistent polarizable ion tight binding theory is constructed. The model is applied to water monomer, dimer, hexamers, ice, and liquid, where it demonstrates good agreement with theoretical results obtained by more accurate methods, such as DFT and CCSD(T), and with experiment. In particular, the temperature dependence of the self diffusion coefficient in liquid water predicted by the model, closely reproduces experimental curves in the temperature interval between 230 K and 350 K. In addition, and in contrast to standard DFT, the model properly orders the relative densities of liquid water and ice. A notable, but inevitable, shortcoming of the model is underestimation of the static dielectric constant by a factor of two. We demonstrate that the description of inter and intramolecular forces embodied in the tight binding approximation in quantum mechanics leads to a number of valuable insights which can be missing from ab initio quantum chemistry and classical force fields. These include a discussion of the origin of the enhanced molecular electric dipole moment in the condensed phases, and a detailed explanation for the increase of coordination number in liquid water as a function of temperature and compared with ice-leading to insights into the anomalous expansion on freezing. The theory holds out the prospect of an understanding of the currently unexplained density maximum of water near the freezing point. 

Positron scattering and annihilation on noble-gas atoms

Positron scattering and annihilation on noble-gas atoms is studied ab initio using many-body theory methods for positron energies below the positronium formation threshold. We show that in this energy range, the many-body theory yields accurate numerical results and provides a near-complete understanding of the positron–noble-gas atom system. It accounts for positron-atom and electron-positron correlations, including the polarization of the atom by the positron and the nonperturbative effect of virtual positronium formation. These correlations have a large influence on the scattering dynamics and result in a strong enhancement of the annihilation rates compared to the independent-particle mean-field description. Computed elastic scattering cross sections are found to be in good agreement with recent experimental results and Kohn variational and convergent close-coupling calculations. The calculated values of the annihilation rate parameter Zeff (effective number of electrons participating in annihilation) rise steeply along the sequence of noble-gas atoms due to the increasing strength of the correlation effects, and agree well with experimental data.