1ST PRINCIPLE CALCULATIONS OF QUASI-PARTICLE ENERGIES FOR BAND STRUCTURES OF SEMICONDUCTORS

We successfully applied the Green function theory in GW approximation to calculate the quasiparticle energies for semiconductors Si and GaAs. Ab initio pseudopotential method was adopted to generate basis wavefunctions and charge densities for calculating dielectric matrix elements and electron self-energies. To evaluate dynamical effects of screened interaction, GPP model was utilized to extend dieletric matrix elements from static results to finite frequencies. We give a full account of the theoretical background and the technical details for the first principle pseudopotential calculations of quasiparticle energies in semiconductors and insulators. Careful analyses are given for the effective and accurate evaluations of dielectric matrix elements and quasiparticle self-energies by using the symmetry properties of basis wavefunctions and eigenenergies. Good agreements between the calculated excitation energies and fundamental energy gaps and the experimental band structures were achieved.

2-DELTA-FUNCTION GENERALIZED PLASMA POLE MODEL FOR 1ST PRINCIPLE CALCULATIONS OF QUASI-PARTICLE ENERGIES IN MANY-ELECTRON SYSTEMS

To evaluate the dynamical effects of the screened interaction in the calculations of quasiparticle energies in many-electron systems a two-delta-function generalized plasma pole model (GPP) is introduced to simulate the dynamical dielectric function. The usual single delta-function GPP model has the drawback of over simplifications and for the crystals without the center of symmetry is inappropriate to describe the finite frequency behavior for dielectric function matrices. The discrete frequency summation method requires too much computation to achieve converged results since ab initio calculations of dielectric function matrices are to be carried out for many different frequencies. The two-delta GPP model is an optimization of the two approaches. We analyze the two-delta GPP model and propose a method to determine from the first principle calculations the amplitudes and effective frequencies of these delta-functions. Analytical solutions are found for the second order equations for the parameter matrices entering the model. This enables realistic applications of the method to the first principle quasiparticle calculations and makes the calculations truly adjustable parameter free.

Theoretical studies on electronic structures, spectra and charge transporting properties of a series of Pt((CN)-N-Lambda)(2) complexes

We report a quantum-chemical study of electronic, optical and charge transporting properties of four platinum (II) complexes, pt((CN)-N-Lambda)(2) ((CN)-N-Lambda=phenylpyridine or thiophenepyridine). The lowest-lying absorptions at 442, 440, 447 and 429 nm are all attributed to the mixed transition characters of metal-to-ligand charge transfer (MLCT) and ligand-centered (LC) pi - pi(*) transition. While, unexpectedly, the lowest-lying phosphorescent emissions at 663, 660, 675 and 742 nm are mainly from metal-to-ligand charge transfer ((MLCT)-M-3) ligand-centered (LC) pi ->pi* transition. Ionization potential (IP), electron affinities (EA) and reorganization energy P (lambda(hole/electron)) were obtained to evaluate the charge transfer and balance properties between hole and electron.

Peierls-Nabarro model of interfacial misfit dislocation: An analytic solution

We propose a method to treat the interfacial misfit dislocation array following the original Peierls-Nabarro's ideas. A simple and exact analytic solution is derived in the extended Peierls-Nabarro's model, and this solution reflects the core structure and the energy of misfit dislocation, which depend on misfit and bond strength. We also find that only with beta < 0.2 the structure of interface can be represented by an array of singular Volterra dislocations, which conforms to those of atomic simulation. Interfacial energy and adhesive work can be estimated by inputting ab initio calculation data into the model, and this shows the method can provide a correlation between the ab initio calculations and elastic continuum theory.

Quantum Monte Carlo study of the CO interaction with a dimer model surface for Cr(110)

The chemisorption of CO on a Cr( 110) surface is investigated using the quantum Monte Carlo method in the diffusion Monte Carlo (DMC) variant and a model Cr2CO cluster. The present results are consistent with the earlier ab initio HF study with this model that showed the tilted/ near-parallel orientation as energetically favoured over the perpendicular arrangement. The DMC energy difference between the two orientations is larger (1.9 eV) than that computed in the previous study. The distribution and reorganization of electrons during CO adsorption on the model surface are analysed using the topological electron localization function method that yields electron populations, charge transfer and clear insight on the chemical bonding that occurs with CO adsorption and dissociation on the model surface.

Hybrid Qm/Mm Simulation Of The Hydration Phenomena Of Dipalmitoylphosphatidylcholine Headgroup

The polar headgroup of dipalmitoylphosphatidylcholine (DPPC) molecule both in gas phase and aqueous Solution is investigated by the hybrid quantum mechanical/molecular mechanical (QM/MM) method, in which the polar head of DPPC molecule and the bound water molecules are treated with density functional theory (DFT), while the apolar hydrocarbon chain of DPPC molecule is treated with MM method. It is demonstrated that the hybrid QM/MM method is both accurate and efficient to describe the conformations of DPPC headgroup. Folded structures of headgroup are found in gas phase calculations. In this work, both monohydration and polyhydration phenomena are investigated. In monohydration, different water association sites are studied. Both the hydration energy and the quantum properties of DPPC and water molecules are calculated at the DFT level of theory after geometry optimization. The binding force of monohydration is estimated by using the scan method. In polyhydration, more extended conformations are found and hydration energies in different polyhydration styles are estimated. (C) 2008 Elsevier Inc. All rights reserved.

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.

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.

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.

Energy bands and acceptor binding energies of GaN

The energy bands of zinc-blende and wurtzite GaN are calculated with the empirical pseudopotential method, and the pseudopotential parameters for Ga and N atoms are-given. The calculated energy bands are in agreement with those obtained by the ab initio method. The effective-mass theory for the semiconductors of wurtzite structure is established, and the effective-mass parameters of GaN for both structures are given The binding energies of acceptor states are calculated by solving strictly the effective-mass equations. The binding energies of donor and acceptor are 24 and 142 meV for the zinc-blende structure, 20 and 131, and 97 meV for the wurtzite structure, respectively, which are consistent with recent experimental results. It is proposed that there are two kinds of acceptor in wurtzite GaN. One kind is the general acceptor such as C, which substitutes N, which satisfies the effective-mass theory. The other kind of acceptor includes Mg, Zn, Cd, etc., the binding energy of these accepters is deviated from that given by the effective mass theory. In this report, wurtzite GaN is grown by the molecular-beam epitaxy method, and the photoluminescence spectra were measured. Three main peaks are assigned to the donor-acceptor transitions from two kinds of accepters. Some of the transitions were identified as coming from the cubic phase of GaN, which appears randomly within the predominantly hexagonal material. [S0163-1829(99)15915-0].

INTERVALLEY GAMMA-CHI DEFORMATION POTENTIALS IN III-V ZINCBLENDE SEMICONDUCTORS BY ABINITIO PSEUDOPOTENTIAL CALCULATIONS

Intervalley GAMMA-X deformation-potential constants (IVDP's) have been calculated by use of a first-principles pseudopotential method for the III-V zinc-blende semiconductors AlP, AlAs, AlSb, GaP, GaAs, GaSb, InP, InAs, and InSb. When the calculated IVDP's of LA phonons for GaP, InP, and InAs and of LO phonons for AlAs, AlSb, GaAs, GaSb, and InSb are compared with results of a previous calculation that used the empirical pseudopotential method (EPM) and a rigid-ion approximation, good agreement is found. However, our ab initio pseudopotential results on IVDP's of LA phonons for AlAs, AlSb, GaAs, GaSb, and InSb and of LO phonons for GaP, InP, and InAs are about one order of magnitude smaller than those obtained by use of EPM calculations, indicating that the electron redistribution accompanying crystal-lattice deformation has a significant effect on GAMMA-X intervalley scattering for these phonon modes when the anions are being displaced. In our calculations the LA- and LO-phonon modes at the X point have been evaluated in the frozen-phonon approximation. We have also obtained the LAX- and LOX-phonon frequencies for these materials from total-energy calculations, which agree very well with experimental values for these semiconductors. We have also calculated GAMMA-X hole-phonon scattering matrix elements for the top valence bands in these nine semiconductors, from which the GAMMA-X IVDP's of the top valence bands for the longitudinal phonons and transverse phonons are evaluated, respectively.

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.

Theoretical study on analytical potential function and spectroscopic parameters for CaF molecule

The equilibrium properties and potential energy curves of the ground electronic state of CaF have been calculated using the Brueckner Doubles calculation with a triples contribution added [BD(T)] and the gradient-corrected density functional theory with three-parameter exact exchange mixing (B3LY-P) method, with 6-311 + G*,6-311 + G(2df,2pd) and 6-311 + G(3df,3pd) basis sets. All the computational PECs are fitted to analytical potential energy functions using Murrell-Sorbie, Huxley and Tang-Toennies potentials. Based on this, the spectroscopic parameters are calculated, and then compared with some other theoretical and experimental data. (C) 2004 Elsevier B.V. All rights reserved.

Why [6,6]- and 1,2-Benzal-3-N-4-O-Cyclic Phenylimidate C-60 Undergo Electrochemically Induced Retro-Addition Reactions while 1,4-Dibenzyl-2,3-Cyclic Phenylimidate C-60 Does Not? C-H center dot center dot center dot X (X = N, O) Intramolecular Interactions in Organofullerenes

The electrochemical properties of a series of structurally related fullerooxazoles, [6,6] cyclic phenylimidate C-60 (1), 1,2-benzal-3-N-4-O-cyclic phenylimidate C-60 (2), and 1,4-dibenzyl-2,3-cyclic phenylimidate C-60 (3), are described, and the spectroscopic characterizations of their anionic species are reported. The results show that compounds I and 2 undergo retro-cycloaddition reactions that lead to the formation of C-60 and C61HPh, respectively, upon two-electron-transfer reduction. However, compound 3 demonstrates much more electrochemical stability as no retro-cycloaddition reaction occurs under similar conditions. Natural bond orbital (NBO) calculations on charge distribution show there is no significant difference among the dianions of 1, 2, and 3, indicating that the electrochemical stability of 3 is unlikely to be caused by the charge distribution difference of the dianions of three compounds. Examination on the crystal structure of compound 3 reveals close contacts of the C-H group with the heteroatoms (N and O) of cyclic phenylimidate, suggesting the existence of C-H center dot center dot center dot X (X = N, O) intramolecular hydrogen bonding among the addends, which is further confirmed by NBO analysis.

Electronic structures and chemical bonding in diatomic ScX to ZnX (X = S, Se, Te)

Bond distances, vibrational frequencies, electron affinities, ionization potentials, dissociation energies, and dipole moments of the title molecules in neutral, positively, and negatively charged ions were studied using density functional method. Ground electronic state was assigned for each molecule. The bonding patterns were analyzed and compared with both the available data and across the series. It was found that, besides ionic component, covalent bonds are formed between the metal s, d orbitals, and the p orbital of S, Se, and Te. For neutral and cationic molecules, the covalent character increases from ScX to CrX and from FeX to CuX with an exception of decrease at MnX and ZnX, while for anionic molecules, the trend is not obvious. For both neutral and charged molecules, the sulfides have the shortest bond distance and largest vibrational frequency, while tellurides have the largest bond distance and smallest vibrational frequency. For neutral and anionic molecules, the dissociation energy of sulfides is the largest, that of tellurides is the smallest, while this only remains true for cationic molecules from ScX+ to FeX+.