First-principles study of transition metal impurities in Si

By using ab initio electronic structure calculations within density functional theory, we study the structural, electronic, and magnetic properties of Si doped with a transition metal impurity. We consider the transition metals of the 3d series V, Cr, Mn, Fe, Co, and Ni. To get insight into the level filling mechanism and the magnetization saturation, we first investigate the transition metal-Si alloys in the zinc-blende structure. Next, we investigate the doping of bulk Si with a transition metal atom, in which it occupies the substitutional site, the interstitial site with tetrahedral symmetry, and the interstitial site with hexagonal symmetry. It is found that all of these transition metal impurities prefer an interstitial position in Si. Furthermore, we show that it is possible to interpret the electronic and magnetic properties by using a simple level filling picture and a comparison is made to Ge doped with the same transition metal atoms. In order to get insight into the effect of a strained environment, we calculate the formation energy as a function of an applied homogeneous pressure and we show that an applied pressure can stabilize the substitutional position of transition metal impurities in Si. Finally, the energies of the ferromagnetic states are compared to those of the antiferromagnetic states. It is shown that the interstitial site of the Mn dopant helps us to stabilize the nearest neighbor substitutional site to realize the ferromagnetic state. For doping of Si with Cr, a ferrimagnetic behavior is predicted.

First-principles study of the electronic structures and magnetic properties of 3d transition metal-doped anatase TiO2

We study the electronic structures and magnetic properties of the anatase TiO2 doped with 3d transition metals (V, Cr, Mn, Fe, Co, Ni), using first-principles total energy calculations based on density functional theory (DFT). Using a molecular-orbital bonding model, the electronic structures of the doped anatase TiO2 are well understood. A band coupling model based on d-d level repulsions between the dopant ions is proposed to understand the chemical trend of the magnetic ordering. Ferromagnetism is found to be stabilized in the V-, Cr-, and Co-doped samples if there are no other carrier native defects or dopants. The ferromagnetism in the Cr- and Co-doped samples may be weakened by the donor defects. In the Mn-, and Fe-doped samples, the ferromagnetism can be enhanced by the acceptor and donor defects, respectively.

First-principles study of native defects in rutile TiO2

Native point defects in the rutile TiO2 are studied via first-principles pseudopotential calculations. Except for the two antisite defects, all the native point defects have low formation energies. Under the Ti-rich growth condition, high concentrations of titanium interstitials and oxygen vacancies would form spontaneously in p-type samples; whereas high concentrations of titanium vacancies would form spontaneously in n-type samples regardless of the oxygen partial pressure. (c) 2007 Elsevier B.V. All rights reserved.

Nitrogen defects and ferromagnetism in Cr-doped dilute magnetic semiconductor AlN from first principles

High Curie temperature of 900 K has been reported in Cr-doped AlN diluted magnetic semiconductors prepared by various methods, which is exciting for spintronic applications. It is believed that N defects play important roles in achieving the high-temperature ferromagnetism in good samples. Motivated by these experimental advances, we use a full-potential density-functional-theory method and supercell approach to investigate N defects and their effects on ferromagnetism of (Al,Cr)N with N vacancies (V-N). We investigate the structural and electronic properties of V-N, single Cr atom, Cr-Cr atom pairs, Cr-V-N pairs, and so on. In each case, the most stable structure is obtained by comparing different atomic configurations optimized in terms of the total energy and the force on every atom, and then it is used to calculate the defect formation energy and study the electronic structures. Our total-energy calculations show that the nearest substitutional Cr-Cr pair with the two spins in parallel is the most favorable and the nearest Cr-V-N pair makes a stable complex. Our formation energies indicate that V-N regions can be formed spontaneously under N-poor condition because the minimal V-N formation energy equals -0.23 eV or Cr-doped regions with high enough concentrations can be formed under N-rich condition because the Cr formation energy equals 0.04 eV, and hence real Cr-doped AlN samples are formed by forming some Cr-doped regions and separated V-N regions and through subsequent atomic relaxation during annealing. Both of the single Cr atom and the N vacancy create filled electronic states in the semiconductor gap of AlN. N vacancies enhance the ferromagnetism by adding mu(B) to the Cr moment each but reduce the ferromagnetic exchange constants between the spins in the nearest Cr-Cr pairs. These calculated results are in agreement with experimental observations and facts of real Cr-doped AlN samples and their synthesis. Our first-principles results are useful to elucidate the mechanism for the ferromagnetism and to explore high-performance Cr-doped AlN diluted magnetic semiconductors.

Electronic Properties of Nonstoichiometric PbSe Quantum Dots from First Principles

One-Hundred Talent Plan of the Chinese Academy of Sciences；National Science Fund for Distinguished Young Scholars 60925016；National High Technology Research and Development program of China 2009AA034101

First principles based predictions of the toughness of a metal/oxide interface

We describe a first-principles-based strategy to predict the macroscopic toughness of a gamma-Ni(Al)/alpha-Al2O3 interface. Density functional theory calculations are used to ascertain energy changes upon displacing the two materials adjacent to the interface, with relaxation conducted over all atoms located within adjoining rows. Traction/displacernent curves are obtained from derivatives of the energy. Calculations are performed in mode I (opening), mode II (shear) and at a phase angle of 45 degrees. The shear calculations are conducted for displacements along < 110 > and < 112 > of the Ni lattice. A generalized interface potential function is used to characterize the results. Initial fitting to both the shear and normal stress results is required to calibrate the unknowns. Thereafter, consistency is established by using the potential to predict other traction quantities. The potential is incorporated as a traction/displacement function within a cohesive zone model and used to predict the steady-state toughness of the interface. For this purpose, the plasticity of the Ni alloy must be known, including the plasticity length scale. Measurements obtained for a gamma-Ni superalloy are used and the toughness predicted over the full range of mode mixity. Additional results for a range of alloys are used to demonstrate the influences of yield strength and length scale.

First principles based predictions of the toughness of a metal/oxide interface

We describe a first-principles-based strategy to predict the macroscopic toughness of a gamma-Ni(Al)/alpha-Al2O3 interface. Density functional theory calculations are used to ascertain energy changes upon displacing the two materials adjacent to the interface, with relaxation conducted over all atoms located within adjoining rows. Traction/displacernent curves are obtained from derivatives of the energy. Calculations are performed in mode I (opening), mode II (shear) and at a phase angle of 45 degrees. The shear calculations are conducted for displacements along < 110 > and < 112 > of the Ni lattice. A generalized interface potential function is used to characterize the results. Initial fitting to both the shear and normal stress results is required to calibrate the unknowns. Thereafter, consistency is established by using the potential to predict other traction quantities. The potential is incorporated as a traction/displacement function within a cohesive zone model and used to predict the steady-state toughness of the interface. For this purpose, the plasticity of the Ni alloy must be known, including the plasticity length scale. Measurements obtained for a gamma-Ni superalloy are used and the toughness predicted over the full range of mode mixity. Additional results for a range of alloys are used to demonstrate the influences of yield strength and length scale.

Electronic and mechanical properties of 5d transition metal mononitrides via first principles

The electronic and mechanical properties of 5d transition metal mononitrides from LaN to AuN are systematically investigated by use of the density-functional theory. For each nitride, Six Structures are considered, i.e., rocksalt, zinc blende, CsCl, wurtzite, NiAs and WC structures. Among the considered structures, rocksalt structure is the most stable for LaN, HfN and ALIN, WC structure for TaN, NiAs structure for WN, wurtzite structure for ReN, OsN, IrN and PtN. The most stable Structure for each nitride is mechanically stable. The formation enthalpy increases from LaN to AuN.

First Principles Investigation on the Ultra-Incompressible and Hard TaN

The structural, electronic, and mechanical properties of TaN were investigated by use of the density functional theory (DFT). Eight structures were considered, i.e.. hexagonal WC TaN, NiAs, wurtzite, and CoSn structures. cubic NaCl. zinc-blende and CsCl structures. The results indicate that TaN in TaN-type structure is the most stable at ambient conditions among the considered structures. Above 5 GPa, TaN in WC-type structure becomes energetically the most stable phase. They are also stable both thermodynamically and mechanically. TaN in WC-type has the largest shear Modulus 243 GPa and large bulk modulus 337 GPa among the considered structures. The Volume compressibility is slightly larger than diamond, but smaller than c-BN at pressures from 0 to 100 GPa. The compressibility along the c axis is smaller than the linear compressibility of both diamond and c-BN.

Charge ordering induced metal-semiconductor transition in Ag2BiO3: A first-principles study

Accurate ab initio density-functional calculations are performed to investigate the relationship of the ground-state crystal structures and electronic properties of Ag2BiO3 compound. The results indicate that Ag2BiO3 in Pnna phase, in which the bismuth atoms occupy the same Wyckoff positions, exhibits metallic conductivity, while in Pnn2 and Pn phases, Ag2BiO3 exhibits semiconducting character, which is in agreement with the experimental results. Charge ordering is indeed induced by the crystal inversion twin in the Pnn2 phase compared with the Pnna phase. In the low temperature phase Pn, the charge ordering is similar to that of Pnn2 phase although it is more distorted in Pn phase. In addition, the calculation indicates that the charge ordering is caused in the 6s electron rearrangement.

A first-principles study of the different magnetoresistance mechanisms in CaCu3Mn4O12 and LaCu3Mn4O12

The electronic structure of CaCu3Mn4O12 and LaCu3Mn4O12 was investigated using a full-potential linearized augmented plane wave method within the Generalized Gradient Approximation (GGA). The ferrimagnetic and ferromagnetic states in these two compounds were investigated and the calculated spin magnetic moments were found to be close to the available experimental values. Calculations of spin polarization for these two oxides show that the ferrimagnetic configurations are the energetically favored ground state, which is consistent with experimental observation. The calculations predict that CaCu3Mn4O12 is a semiconductor and that LaCu3Mn4O12 is a half-metallic material. Furthermore, the relevance of these different electronic structures to the magnetoresistance is discussed.

Low-compressibility and hard materials ReB2 and WB2: Prediction from first-principles study

First-principle calculations are performed to investigate the structural, elastic, and electronic properties of ReB2 and WB2. The calculated equilibrium structural parameters of ReB2 are consistent with the available experimental data. The calculations indicate that WB2 in the P6(3)/mmc space group is more energetically stable under the ambient condition than in the P6/mmm. Based on the calculated bulk modulus, shear modulus of polycrystalline aggregate, ReB2 and WB2 can be regarded as potential candidates of ultra-incompressible and hard materials. Furthermore, the elastic anisotropy is discussed by investigating the elastic stiffness constants. Density of states and electron density analysis unravel the covalent bonding between the transition metal atoms and the boron atoms as the driving force of the high bulk modulus and high shear modulus as well as small Poisson's ratio.

First principles simulation of amorphous silicon bulk, interfaces, and nanowires for photovoltaics

Amorphous silicon has become the material of choice for many technologies, with major applications in large area electronics: displays, image sensing and thin film photovoltaic cells. This technology development has occurred because amorphous silicon is a thin film semiconductor that can be deposited on large, low cost substrates using low temperature. In this thesis, classical molecular dynamics and first principles DFT calculations have been performed to generate structural models of amorphous and hydrogenated amorphous silicon and interfaces of amorphous and crystalline silicon, with the ultimate aim of understanding the photovoltaic properties of core-shell crystalline amorphous Si nanowire structures. We have shown, unexpectedly, from the simulations, that our understanding of hydrogenated bulk a-Si needs to be revisited, with our robust finding that when fully saturated with hydrogen, bulk a-Si exhibits a constant optical energy gap, irrespective of the hydrogen concentration in the sample. Unsaturated a-Si:H, with a lower than optimum hydrogen content, shows a smaller optical gap, that increases with hydrogen content until saturation is reached. The mobility gaps obtained from an analysis of the electronic states show similar behavior. We also obtained that the optical and mobility gaps show a volcano curve as the H content is varied from 7% (undersaturation) to 18% (mild oversaturation). In the case of mild over saturation, the mid-gap states arise exclusively from an increase in the density of strained Si-Si bonds. Analysis of our structures shows the extra H atoms in this case form a bridge between neighboring silicon atoms which increases the corresponding Si-Si distance and promotes bond length disorder in the sample. That has the potential to enhance the Staebler-Wronski effect. Planar interface models of amorphous-crystalline silicon have been generated in Si (100), (110) and (111) surfaces. The interface models are characterized by structure, RDF, electronic density of states and optical absorption spectrum. We find that the least stable (100) surface will result in the formation of the thickest amorphous silicon layer, while the most stable (110) surface forms the smallest amorphous region. We calculated for the first time band offsets of a-Si:H/c-Si heterojunctions from first principles and examined the influence of different surface orientations and amorphous layer thickness on the offsets and implications for device performance. The band offsets depend on the amorphous layer thickness and increase with thickness. By controlling the amorphous layer thickness we can potentially optimise the solar cell parameters. Finally, we have successfully generated different amorphous layer thickness of the a-Si/c-Si and a-Si:H/c-Si 5 nm nanowires from heat and quench. We perform structural analysis of the a-Si-/c-Si nanowires. The RDF, Si-Si bond length distributions, and the coordination number distributions of amorphous regions of the nanowires reproduce similar behaviour compared to bulk amorphous silicon. In the final part of this thesis we examine different surface terminating chemical groups, -H, - OH and –NH2 in (001) GeNW. Our work shows that the diameter of Ge nanowires and the nature of surface terminating groups both play a significant role in both the magnitude and the nature of the nanowire band gaps, allowing tuning of the band gap by up to 1.1 eV. We also show for the first time how the nanowire diameter and surface termination shifts the absorption edge in the Ge nanowires to longer wavelengths. Thus, the combination of nanowire diameter and surface chemistry can be effectively utilised to tune the band gaps and thus light absorption properties of small diameter Ge nanowires.

Catalytic Role of Metal Oxides in Gold-Based Catalysts: A First Principles Study of CO Oxidation on TiO2 Supported Au

CO oxidation on TiO2 supported Au has been studied using density functional theory calculations. Important catalytic roles of the oxide have been identified: (i) CO oxidation occurs at the interface between Au and the oxide with a very small barrier; and (ii) O-2 adsorption at the interface is the key step in the reaction. The physical origin of the oxide promotion effect has been further investigated: The oxide enhances electron transfer from the Au to the antibonding states of O-2, giving rise to (i) strong ionic bonding between the adsorbed O-2, Au, and the Ti cation; and (ii) a significant activation of O-2 towards CO oxidation.

First-principles study of oxygenates on Co surfaces in Fischer-Tropsch synthesis

Extensive density function theory calculations are performed to study the mechanism of the formation of aldehyde and alcohol on Co surfaces in Fischer-Tropsch synthesis, a challenging issue in heterogeneous catalysis. Three possible pathways for the production of formaldehyde and methanol on flat and stepped Co(0001) surfaces are investigated: (i) CO + 4H -> CHO + 3H -> CH2O + 2H -> CH3O + H -> CH3OH; (ii) CO + 4H -> COH + 3H -> CHOH + 2H -> CH2OH + H -> CH3OH; and (iii) the coupling reactions of CH2 + O -> CH2O and CH3 + OH -> CH3OH. It is found that these pathways are generally favored at step sites, and the preferred mechanism is pathway (i) via CHO. Furthermore, the three traditional chain growth mechanisms in Fischer-Tropsch synthesis are semi quantitatively compared and discussed. Our results suggest that the two mechanisms involving oxygenate intermediates (the CO-insertion and hydroxycarbene mechanisms) are less important than the carbene mechanism in the production of long chain hydrocarbons. However, the CO-insertion mechanism may be responsible for the production of long-chain oxygenates.