Mechanical and chemical bonding properties of ground state BeH2

The crystal structure, mechanical properties and electronic structure of ground state BeH2 are calculated employing the first-principles methods based on the density functional theory. Our calculated structural parameters at equilibrium volume are well consistent with experimental results. Elastic constants, which well obey the mechanical stability criteria, are firstly theoretically acquired. The bulk modulus B, Shear modulus G, Young's modulus E, and Poisson's ratio upsilon are deduced from the elastic constants. The bonding nature in BeH2 is fully interpreted by combining characteristics in band structure, density of states, and charge distribution. The ionicity in the Be-H bond is mainly featured by charge transfer from Be 2s to H 1s atomic orbitals while its covalency is dominated by the hybridization of H 1s and Be 2p states. The Bader analysis of BeH2 and MgH2 are performed to describe the ionic/covalent character quantitatively and we find that about 1.61 (1.6) electrons transfer from each Be (Mg) atom to H atoms.

Strain relaxation and band-gap tunability in ternary InxGa1-xN nanowires

The alloy formation enthalpy and band structure of InGaN nanowires were studied by a combined approach of the valence-force field model, Monte Carlo simulation, and density-functional theory (DFT). For both random and ground-state structures of the coherent InGaN alloy, the nanowire configuration was found to be more favorable for the strain relaxation than the bulk alloy. We proposed an analytical formula for computing the band gap of any InGaN nanowires based on the results from the screened exchange hybrid DFT calculations, which in turn reveals a better band-gap tunability in ternary InGaN nanowires than the bulk alloy.

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.

Effects of residual C and O impurities on photoluminescence in undoped GaN epilayers

Photoluminescence (PL) was investigated in undoped GaN from 4.8 K to room temperature. The 4.8 K spectra exhibited recombinations of free exciton, donor-acceptor pair (DAP), blue and yellow bands (Ybs). The blue band (BB) was also identified to be a DAP recombination. The YB was assigned to a recombination from deep levels. The energy-dispersive X-ray spectroscopy show that C and O are the main residual impurities in undoped GaN and that C concentration is lower in the epilayers with the stronger BB. The electronic structures of native defects, C and O impurities, and their complexes were calculated using ab initio local-density-functional (LDF) methods with linear muffin-tin-orbital and 72-atomic supercell. The theoretical analyses suggest that the electron transitions from O-N states to C-N and to V-Ga states are responsible for DAP and the BB, respectively, and the electron transitions between the inner levels of the C-N-O-N complex may be responsible for the YB in our samples. (C) 2002 Elsevier Science B.V. All rights reserved.

Effects of residual C and O impurities on photoluminescence in undoped GaN epilayers

Photoluminescence (PL) was investigated in undoped GaN from 4.8 K to room temperature. The 4.8 K spectra exhibited recombinations of free exciton, donor-acceptor pair (DAP), blue and yellow bands (Ybs). The blue band (BB) was also identified to be a DAP recombination. The YB was assigned to a recombination from deep levels. The energy-dispersive X-ray spectroscopy show that C and O are the main residual impurities in undoped GaN and that C concentration is lower in the epilayers with the stronger BB. The electronic structures of native defects, C and O impurities, and their complexes were calculated using ab initio local-density-functional (LDF) methods with linear muffin-tin-orbital and 72-atomic supercell. The theoretical analyses suggest that the electron transitions from O-N states to C-N and to V-Ga states are responsible for DAP and the BB, respectively, and the electron transitions between the inner levels of the C-N-O-N complex may be responsible for the YB in our samples. (C) 2002 Elsevier Science B.V. All rights reserved.

Golden-rule treatment on the ClO/ClO+ electron-transfer system

The structures, properties and electron transfer reactivity of the ClO/ClO+ coupling system are studied in this paper at ab initio (HF and MP2) levels and the density functional theory (DFT: B3LYP, B3P86, B3PW91) levels employing 6311 + G(3df) basis set and on the basis of the golden-rule of the time-dependent perturbation theory. Investigations indicate that the results got from the B3LYP method employing 6-311 + G(3df) basis set is in excellent agreement with the experiment. The activation energies, the stabilization energies and the electronic coupling matrix elements have also been calculated by using the B3LYP/6-311 + G(3df) method, and then the electron transfer rates are determined at this level. The electronic coupling matrix element of EC.6 is very small, only 0.03 kcal/mol, while that of EC.7 is the biggest, being 12.41 kcal/mol, the corresponding electron transfer rate is also the fastest among these seven encounter complexes. The averaged electron transfer rate is about 1.672 X 10(11) M-1 s(-1). It is indicated that the structures optimized by B3LYP method are more reliable than the results got from the other four methods. It also testified that the electronic coupling matrix element is the vital factor that significantly affects the electron transfer rate. (C) 2003 Elsevier B.V. All rights reserved.

Theoretical studies on the reaction mechanism of palladium(0)-catalyzed addition of thiocyanates to alkynes

The reaction mechanism of the Pd(0)-catalyzed alkyne cyanothiolation reaction is investigated by MP2, CCSD(T) and the density functional method B3LYP. The overall reaction mechanism is examined. The B3LYP results are consistent with the results of CCSD(T) and MP2 methods for the isomerization, acetylene insertion and reductive elimination steps, but not for the oxidative addition step. For the oxidative addition, the bisphosphine and monophosphine pathways are competitive in B3LYP, while the bisphosphine one is preferred for CCSD(T) and MP2 methods.

Injection, Transport, Absorption and Phosphorescence Properties of a Series of Blue-Emitting Ir(III) Emitters in OLEDs: a DFT and Time-Dependent DFT Study

Quantum-chemistry methods were explored to investigate the electronic structures, injection and transport properties, absorption and phosphorescence mechanism of a series of blue-emitting Ir(III) complexes {[(F-2-ppy)(2)Ir(pta -X/pyN4)], where F-2-ppy = (2,4-difluoro)phenylpyridine; pta = pyridine-1,2,4-triazole; X = phenyl(1); p-tolyl (2); 2,6-difluororophenyl (3); -CF3 (4), and pyN4 = pyridine-1,2,4-tetrazolate (5)}, which are used as emitters in organic light-emitting diodes (OLEDs). The mobility of hole and electron were studied computationally based on the Marcus theory. Calculations of Ionization potentials (IPs) and electron affinities (EAs) were used to evaluate the injection abilities of holes and electrons into these complexes.

Structures and vibrational spectra of C-2 and LaC2+ clusters

C-2 and LaC2+ were studied using Hartree-Fock(HF), B3LYP (Becke 3-paremeter-Lee-Yang-Parr) density functional method, second-order Moller-Plesset perturbation (MP2) and coupled cluster singles and doubles with non-iterative triples(CCSD(T)) methods. The basis set employed was LANL1DZ. Geometries, vibrational frequencies and other quantities were reported. The results showed that for C-2, all the methods performed well for low spin state (singlet), while only HF and B3LYP remained so for high spin state (triplet). For LaC2+, four isomers were presented and fully optimized. The results suggested that linear isomers with C-infinity v and D-infinity h symmetries were predicted to be saddle points on the energy surface for all the methods, while for isomers with C-2 upsilon and C-s symmetries, they were local minima except C-2 upsilon at B3LYP level, and were isoenergetic at HF, MP2 and CCSD(T) levels, near isoenergetic at B3LYP level. From the differences between HOMO and LUMO, it is also known that the isomers with C-2 upsilon and C-s symmetries offer the largest values and therefore correspond to the most stable structure. For La-C bond lengths, B3LYP gives the shortest, the order is B3LYPmethods.

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.

Addition reaction of nitrones on the reconstructed Si(100)-2 x 1 surface

The reaction of nitrone, N-methyl nitrone, and their hydroxylamine tautomers (vinyl-hydroxylamine and N-methyl vinyl-hydroxylamine) on the reconstructed Si(100)-2 x 1 surface has been investigated by means of hybrid density functional theory (B3LYP) and Moller-Plesset second-order perturbation (MP2) methods. The calculations predicted that both of the nitrones should react with the surface dimer via facile concerted 1,3-dipolar cycloaddition leading to 5-member-ring compounds. The reaction of hydroxylamine tautomers on the Si(100) surface follows pi-complex (intermediate) mechanism. For the reaction of N-methyl vinyl-hydroxylamine, the pi-complex intermediate undergoes [2+2] cycloaddition leading to a 4-member-ring compound. But in the reaction of vinyl-hydroxylamine, the intermediate undergoes H-migration reaction ("ene" reaction) resulting in the oxime-terminated Si surface. All the surface reactions result in the hydroxyl-terminated silicon surfaces, which are very useful for the further modification of the semiconductor.

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.

地幔压力下CaSiO3 Perovskite状态方程和弹性常数的第一性原理研究

In the past decade density functional theory (DFT) has made its way from a peripheral position in quantum chemistry to center. Of course the often excellent accuracy of the DFT based methods has provided the primary driving force of this development. This dissertation is devoted to the study of physical and chemical properties of planetary materials by first-principle calculation. The concerned properties include the geometry, elastic constants and anisotropy. In the first chapter, we give a systematic introduction to theoretical background and review its progress. Development of quantum chemistry promotes the establishment of DFT. Theorem of Hohenberg-Kohn is the fundament of DFT and is developed to Kohn-Sham equation, which can be used to perform real calculations. Now, new corrections and extensions, together with developed exchange-correlation, have made DFT more accurate and suitable for larger systems. In the second chapter, we focus on the calculational methods and technical aspects of DFT. Although it is important to develop methods and program, external package are still often used. At the end of this chapter, we briefly some widely used simulation package and the application of DFT. In the third chapter, we begin to focus on properties of real materials by first principles calculation. We study a kind of minerals named Ca perovskite, investigate its possible structure and anisotropy at Earth’s mental condition. By understanding and predicting geo-physically important materials properties at extreme conditions, we can get the most accurate information to interpret seismic data in the context of likely geophysical processes.

压力和铁含量对铁方镁石(Mg1-xFex)O(x<=0.25)弹性常数影响的第一性原理研究

By first principle methods based on density functional theory (DFT),the equation of state(EOS) and elastic constants of both periclase and ferropericlase are calculated. The pressure and iron doping effects on the elastic constants of ferropericlase are investigated systematically. Firstly, we calculate the elastic constants of periclase and compare the obtained results with experimental data and other theoretical calculations, which shows a encouraging consistence and demonstrates the practicability of first-principle methods. Secondly, by adding iron into periclase crystal model, we build up ferropericlase with iron contents ranging from 0% to 25% mole percent. The corresponding elastic constants are calculated in a large pressure range(0~120GPa). Emphatically, the strong correlation of 3d electrons in transitional elements, such as iron, is difficult to treat in first-principle methods for a long time. The current solution is to make additional correction. During the initial stage of this study, the strong correlation of 3d electrons in iron is not considered, and we observed that addition of iron decreases the volume of ferropericlase, which is totally contradictory to the experimental data. By applying LDA+U approximation in order to solve the strongly correlated 3d electron of iron, we observed the expansion of volume by iron as expected. On the basis of the LDA+U approximation, the elastic constants of ferropericlase are calculated. After a detailed analysis of data obtained from theoretical calculations, we have reached the following conclusions:(1)pressure imposes positive effects on all elastic constants, and the degree of effects is C11>C12>C44. (2) Iron has no distinctive effects on C11 and C12, although some fluctuations are observed around 60GPa. However, iron has obvious softening effects on C44 The softening effects on C44 are intensified as pressure increases. Above the 100GPa, the effects increase greatly, even surpasses the pressure's positive effects in ferropericlase crystal models with iron mole percent of having 12.5%, 18.75% and 25% iron content. (3)As to the modulus deprived from elastic constants, iron has no effect on the adiabatic bulk module BS, only a little fluctuation around 60GPa. We find iron's softening effects on shear modulus G. (4)We find out that, compared with low iron content, elastic constants with iron content approaching 25mole% is consistently fluctuated，which may be caused by the limitations of the LDA+U approximation method itself. (5)We investigate the pressure and Fe doping effects on elastic anisotropy factor（A=(2C44+C12-C11)/C11） of ferropericlase and find out that iron contents will lower the critical isotropic pressure. At the same pressure, when the pressure is below the isotropic pressure, iron softens the anisotropy factor ; when pressure surpasses the isotropic pressure, iron increases the anisotropy factor.

A new interatomic potential for nanoscale silica

Flikkema, Edwin; Bromley, S.T., (2003) 'A new interatomic potential for nanoscale silica', Chemical Physics Letters 378(5-6) pp.622-629 RAE2008