Coulomb-hole summations and energies for GW calculations with limited number of empty orbitals: a modified static remainder approach

Ab initio GW calculations are a standard method for computing the spectroscopic properties of many materials. The most computationally expensive part in conventional implementations of the method is the generation and summation over the large number of empty orbitals required to converge the electron self-energy. We propose a scheme to reduce the summation over empty states by the use of a modified static remainder approximation, which is simple to implement and yields accurate self-energies for both bulk and molecular systems requiring a small fraction of the typical number of empty orbitals.

Improved quasiparticle wave functions and mean field for G(0)W(0) calculations: Initialization with the COHSEX operator

The GW approximation to the electron self-energy has become a standard method for ab initio calculation of excited-state properties of condensed-matter systems. In many calculations, the G W self-energy operator, E, is taken to be diagonal in the density functional theory (DFT) Kohn-Sham basis within the G0 W0 scheme. However, there are known situations in which this diagonal Go Wo approximation starting from DFT is inadequate. We present two schemes to resolve such problems. The first, which we called sc-COHSEX-PG W, involves construction of an improved mean field using the static limit of GW, known as COHSEX (Coulomb hole and screened exchange), which is significantly simpler to treat than GW W. In this scheme, frequency-dependent self energy E(N), is constructed and taken to be diagonal in the COHSEX orbitals after the system is solved self-consistently within this formalism. The second method is called off diagonal-COHSEX G W (od-COHSEX-PG W). In this method, one does not self-consistently change the mean-field starting point but diagonalizes the COHSEX Hamiltonian within the Kohn-Sham basis to obtain quasiparticle wave functions and uses the resulting orbitals to construct the G W E in the diagonal form. We apply both methods to a molecular system, silane, and to two bulk systems, Si and Ge under pressure. For silane, both methods give good quasiparticle wave functions and energies. Both methods give good band gaps for bulk silicon and maintain good agreement with experiment. Further, the sc-COHSEX-PGW method solves the qualitatively incorrect DFT mean-field starting point (having a band overlap) in bulk Ge under pressure.

Defects in gallium nitride nanowires: First principles calculations

Atomic configurations and formation energies of native defects in an unsaturated GaN nanowire grown along the [001] direction and with (100) lateral facets are studied using large-scale ab initio calculation. Cation and anion vacancies, antisites, and interstitials in the neutral charge state are all considered. The configurations of these defects in the core region and outermost surface region of the nanowire are different. The atomic configurations of the defects in the core region are same as those in the bulk GaN, and the formation energy is large. The defects at the surface show different atomic configurations with low formation energy. Starting from a Ga vacancy at the edge of the side plane of the nanowire, a N-N split interstitial is formed after relaxation. As a N site is replaced by a Ga atom in the suboutermost layer, the Ga atom will be expelled out of the outermost layers and leaves a vacancy at the original N site. The Ga interstitial at the outmost surface will diffuse out by interstitialcy mechanism. For all the tested cases N-N split interstitials are easily formed with low formation energy in the nanowires, indicating N-2 molecular will appear in the GaN nanowire, which agrees well with experimental findings.

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.

1ST PRINCIPLE CALCULATIONS OF QUASI-PARTICLE ENERGIES FOR BAND STRUCTURES OF SI AND GAAS

We have applied the Green function theory in GW approximation to calculate the quasiparticle energies for semiconductors Si and GaAs. Good agreements of the calculated excitation energies and fundamental energy gaps with the experimental band structures were achieved. We obtained the calculated fundamental gaps of Si and GaAs to be 1.22 and 1.42 eV in comparison to the experimental values of 1.17 and 1.52 eV, respectively. Ab initio pseudopotential method has been used to generate basis wavefunctions and charge densities for calculating dielectric matrix elements and electron self-energies.

1ST-PRINCIPLES CALCULATIONS FOR QUASI-PARTICLE ENERGIES OF GAP AND GAAS

We have applied the Green-function method in the GW approximation to calculate quasiparticle energies for the semiconductors GaP and GaAs. Good agreement between the calculated excitation energies and the experimental results was achieved. We obtained calculated direct band gaps of GaP and GaAs of 2.93 and 1.42 eV, respectively, in comparison with the experimental values of 2.90 and 1.52 eV, respectively. An ab initio pseudopotential method has been used to generate basis wave functions and charge densities for calculating the dielectric matrix elements and self-enegies. To evaluate the dynamical effects of the screened interaction, the generalized-plasma-pole model has been utilized to extend the dielectric matrix elements from static results to finite frequencies. We presen the calculated quasiparticle energies at various high-symmetry points of the Brillouin zone and compare them with the experimental results and other calculations.

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.

Trends in elasticity and electronic structure of 5d transition metal diborides: first-principles calculations

We investigate the cohesive energy, heat of formation, elastic constant and electronic band structure of transition metal diborides TMB2 (TM = Hf, Ta, W, Re, Os and Ir, Pt) in the Pmmn space group using the ab initio pseudopotential total energy method. Our calculations indicate that there is a relationship between elastic constant and valence electron concentration (VEC): the bulk modulus and shear modulus achieve their maximum when the VEC is in the range of 6.8-7.2. In addition, trends in the elastic constant are well explained in terms of electronic band structure analysis, e.g., occupation of valence electrons in states near the Fermi level, which determines the cohesive energy and elastic properties. The maximum in bulk modulus and shear modulus is attributed to the nearly complete filling of TM d-B p bonding states without filling the antibonding states. On the basis of the observed relationship, we predict that alloying W and Re in the orthorhombic structure OsB2 might be harder than alloying the Ir element. Indeed, the further calculations confirmed this expectation.

GRADIENT CORRECTIONS IN DENSITY-FUNCTIONAL THEORY CALCULATIONS FOR SURFACES - CO ON PD(110)

Ab initio total energy calculations have been performed for CO chemisorption on Pd(110). Local density approximation (LDA) calculations yield chemisorption energies which are significantly higher than experimental values but inclusion of the generalised gradient approximation (GGA) gives better agreement. In general, sites with higher coordination of the adsorbate to surface atoms lead to a larger degree of overbinding with LDA, and give larger corrections with GGA. The reason is discussed using a first-order perturbation approximation. It is concluded that this may be a general failure of LDA for chemisorption energy calculations. This conclusion may be extended to many surface calculations, such as potential energy surfaces for diffusion.

Towards a spectroscopically accurate set of potentials for heavy hydride laser cooling candidates: effective core potential calculations of BaH

BaH (and its isotopomers) is an attractive molecular candidate for laser cooling to ultracold temperatures and a potential precursor for the production of ultracold gases of hydrogen and deuterium. The theoretical challenge is to simulate the laser cooling cycle as reliably as possible and this paper addresses the generation of a highly accurate ab initio $^{2}\Sigma^+$ potential for such studies. The performance of various basis sets within the multi-reference configuration-interaction (MRCI) approximation with the Davidson correction (MRCI+Q)is tested and taken to the Complete Basis Set (CBS) limit. It is shown that the calculated molecular constants using a 46 electron Effective Core-Potential (ECP) and even-tempered augmented polarized core-valence basis sets (aug-pCV$n$Z-PP, n= 4 and 5) but only including three active electrons in the MRCI calculation are in excellent agreement with the available experimental values. The predicted dissociation energy De for the X$^2\Sigma^+$ state (extrapolated to the CBS limit) is 16895.12 cm$^{-1}$ (2.094 eV), which agrees within 0.1$\%$ of a revised experimental value of <16910.6 cm$^{-1}$, while the calculated re is within 0.03 pm of the experimental result.

Electronic structure calculations of small AI_n (n = 2 - 8) clusters

The electronic states of small AI_n (n = 2 - 8) clusters have been calculated with a relativistic ab-initio MOLCAO Dirac-Fock-Slater method using numerical atomic DFS wave-functions. The excitation energies were obtained from a ground state calculation of neutral clusters, and in addition from negative clusters charged by half an electron in order to account for part of the relaxation. These energies are compared with experimental photoelectron spectra.

New calculations of P(b) curves for 1s_\sigma excitation in low-Z (Z\le10) ion-atom scattering

Ab initio fully relativistic SCF molecular calculations of energy eigenvalues as well as coupling-matrix elements are used to calculate the 1s_\sigma excitation differential cross section for Ne-Ne and Ne-O in ion-atom collisions. A relativistic perturbation treatment which allows a direct comparison with analogous non-relativistic calculations is also performed.

Inclusive probability calculations for the K-vacancy transfer in collisions of S{^15+} on Ar

Using the single-particle amplitudes from a 20-level coupled-channel calculation with ab initio relativistic self consistent LCAO-MO Dirac-Fock-Slater energy eigenvalues and matrix elements we calculate within the frame of the inclusive probability formalism impact-parameter-dependent K-hole transfer probabilities. As an example we show results for the heavy asymmetric collision system S{^15+} on Ar for impact energies from 4.7 to 16 MeV. The inclusive probability formalism which reinstates the many-particle aspect of the collision system permits a qualitative and quantitative agreement with the experiment which is not achieved by the single-particle picture.

An investigation of the germylene addition reaction, GeH2+C2H2: Time-resolved gas-phase kinetic studies and quantum chemical calculations of the reaction energy surface

Time resolved studies of germylene, GeH2, generated by laser flash photolysis of 3,4-dimethylgermacyclopentene-3, have been carried out to obtain rate constants for its bimolecular reaction with acetylene, C2H2. The reaction was studied in the gas-phase over the pressure range 1-100 Tort, with SF6 as bath gas, at 5 temperatures in the range 297-553 K. The reaction showed a very slight pressure dependence at higher temperatures. The high pressure rate constants (obtained by extrapolation at the three higher temperatures) gave the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1)) (-10.94 +/- 0.05) + (6.10 +/- 0.36 kJ mol(-1))/RTln10. These Arrhenius parameters are consistent with a fast reaction occurring at approximately 30% of the collision rate at 298 K. Quantum chemical calculations (both DFT and ab initio G2//B3LYP and G2//QCISD) of the GeC2H4 potential energy surface (PES), show that GeH2 + C2H2 react initially to form germirene which can isomerise to vinylgermylene with a relatively low barrier. RRKM modelling, based on a loose association transition state, but assuming vinylgermylene is the end product (used in combination with a weak collisional deactivation model) predicts a strong pressure dependence using the calculated energies, in conflict with the experimental evidence. The detailed GeC2H4 PES shows considerable complexity with ten other accessible stable minima (B3LYP level), the three most stable of which are all germylenes. Routes through this complex surface were examined in detail. The only product combination which appears capable of satisfying the (P-3) + C2H4.C2H4 was confirmed as a product by GC observed lack of a strong pressure dependence is Ge(P-3) + C2H4. C2H4 was confirmed as a product by GC analysis. Although the formation of these products are shown to be possible by singlet-triplet curve crossing during dissociation of 1-germiranylidene (1-germacyclopropylidene), it seems more likely (on thermochernical grounds) that the triplet biradical, (GeCH2CH2.)-Ge-., is the immediate product precursor. Comparisons are made with the reaction of SiH2 with C2H2.

Investigation of the prototype silylene reaction, SiH2+H2O(and D2O): time-resolved gas-phase kinetic studies, isotope effects, RRKM calculations, and quantum chemical calculations of the reaction energy surface

Time-resolved kinetic studies of the reaction of silylene, SiH2, with H2O and with D2O have been carried out in the gas phase at 296 and at 339 K, using laser flash photolysis to generate and monitor SiH2. The reaction was studied over the pressure range 10-200 Torr with SF6 as bath gas. The second-order rate constants obtained were pressure dependent, indicating that the reaction is a third-body assisted association process. Rate constants at 339 K were about half those at 296 K. Isotope effects, k(H)/k(D), were small averaging 1.076 0.080, suggesting no involvement of H- (or D-) atom transfer in the rate determining step. RRKM modeling was undertaken based on a transition state appropriate to formation of the expected zwitterionic donoracceptor complex, H2Si...OH2. Because the reaction is close to the low pressure (third order) region, it is difficult to be definitive about the activated complex structure. Various structures were tried, both with and without the incorporation of rotational modes, leading to values for the high-pressure limiting (i.e., true secondorder) rate constant in the range 9.5 x 10(-11) to 5 x 10(-10) cm(3) molecule' s(-1). The RRKM modeling and mechanistic interpretation is supported by ab initio quantum calculations carried out at the G2 and G3 levels. The results are compared and contrasted with the previous studies.