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 PRINCIPLES PSEUDOPOTENTIAL CALCULATIONS OF ZONE BOUNDARY PHONON FREQUENCIES OF III-V-ZINCBLENDE SEMICONDUCTORS

Longitudinal zone boundary X phonon frequencies have been calculated by a first principles pseudopotential method for III-V zincblende semiconductors AlP, AlAs, AlSb, GaP, GaAs, GaSb, InP, InAs and InSb. The phonon frequencies have been evaluated from total energy calculations in the frozen phonon approximation. The calculated phonon frequencies agree very well with the experimental values.

Electronic structure of the nitrogen-vacancy center in diamond from first-principles theory

The nitrogen-vacancy (NV) center is a paramagnetic defect in diamond with applications as a qubit. Here, we investigate its electronic structure by using ab initio density functional theory for five different NV center models of two different cluster sizes. We describe the symmetry and energetics of the low-lying states and compare the optical frequencies obtained to experimental results. We compute the major transition of the negatively charged NV centers to within 25–100 meV accuracy and find that it is energetically favorable for substitutional nitrogens to donate an electron to NV0. The excited state of the major transition and the NV0 state with a neutral donor nitrogen are found to be close in energy.

GaMnAs: Position of Mn-d levels and majority spin band gap predicted from GGA-1/2 calculations

Among all magnetic semiconductors, GaMnAs seems to be the most important one. In this work, we present accurate first-principles calculations of GaMnAs within the GGA-1/2 approach: We concentrate our efforts in obtaining the position of the peak of Mn-d levels in the valence band and also the majority spin band gap. For the position of the Mn-d peak, we find a value of 3.3 eV below the Fermi level, in good agreement with the most recent experimental results of 3.5 and 3.7 eV. An analytical expression that fits the calculated E-g(x) for majority spin is derived in order to provide ready access to the band gap for the composition range from 0 to 0.25. We found a value of 3.9 eV for the gap bowing parameter. The results agree well with the most recent experimental data. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4718602]

Digital magnetic heterostructures based on GaN using GGA-1/2 approach

We present ab-initio calculations of seven digital magnetic heterostructures, GaN delta-doped with V, Cr, Mn, Fe, Co, Ni, and Cu, forming two-dimensional systems. Only GaN delta-doped with V or Cr present a ferromagnetic ground state with high Curie temperatures. For both, to better describe the electronic properties, we used the GGA-1/2 approach. The ground state of GaN/Cr resulted in a two dimensional half-metal, with 100% spin polarization. For GaN/V, we obtained an insulating state: integer magnetic moment of 2.0 mu(B), a minority spin gap of 3.0 eV close to the gap of GaN, but a majority spin gap of 0.34 eV. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4751285]

Magnetic coupling properties of mn-doped ZnO nanowires: First-principles calculations

Based on the density functional theory, we study the magnetic coupling properties of Mn-doped ZnO nanowires. For the nanowires with passivated surfaces, the antiferromagnetic state is found and the Mn atoms have a clustering tendency. When the distance between two Mn atoms is large, the system energetically favors the paramagnetic or spin-glass state. For the nanowires with unpassivated surfaces, the ferromagnetic (FM) coupling states appear between the two nearest Mn atoms, and the zinc vacancies can further stabilize the FM states between them. The electrons with enough concentration possibly mediate the FM coupling due to the negative exchange splitting of conduction band minimum induced by the s-d coupling, which could be useful in nanomaterial design for spintronics. (C) 2008 American Institute of Physics.

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 structural, mechanical and electronic properties of OsC2

The structural, mechanical and electronic properties Of OsC2 were investigated by use of the density functional theory. Seven structures were considered, i.e., orthorhombic Cmca (No. 12, OsSi2), Pmmn (No. 59, 002) and Pnnm (No. 58, OsN2); tetragonal P4(2)/mnm (No. 136, OsO2) and 14/mmm (No. 139, CaC2); cubic Fm-3m (No. 225, CaF2) and Pa-3 (No. 205, PtN2). The results indicate that Cmca in OsSi2 type structure is energetically the most stable phase among the considered structures. It is also stable mechanically. OsC2 in Pmmn phase has the largest bulk modulus 319 GPa and shear modulus 194 GPa. The elastic anisotropy is discussed. (C) 2009 Elsevier B.V. All rights reserved.

Activity and coke formation of nickel and nickel carbide in dry reforming: A deactivation scheme from density functional theory

Dry reforming is a promising reaction to utilise the greenhouse gases CO2 and CH4. Nickel-based catalysts are the most popular catalysts for the reaction, and the coke formation on the catalysts is the main obstacle to the commercialisation of dry reforming. In this study, the whole reaction network of dry reformation on both flat and stepped nickel catalysts (Ni(111) and Ni(211)) as well as nickel carbide (flat: Ni3C(001); stepped: Ni3C(111)) is investigated using density functional theory calculations. The overall reaction energy profiles in the free energy landscape are obtained, and kinetic analyses are utilised to evaluate the activity of the four surfaces. By careful examination of our results, we find the following regarding the activity: (i) flat surfaces are more active than stepped surfaces for the dry reforming and (ii) metallic nickel catalysts are more active than those of nickel carbide, and therefore, the phase transformation from nickel to nickel carbide will reduce the activity. With respect to the coke formation, the following is found: (i) the coke formation probability can be measured by the rate ratio of CH oxidation pathway to C oxidation pathway (r(CH)/r(C)) and the barrier of CO dissociation, (ii) on Ni(111), the coke is unlikely to form, and (iii) the coke formations on the stepped surfaces of both nickel and nickel carbide can readily occur. A deactivation scheme, using which experimental results can be rationalised, is proposed. 

Orbital-polarization terms: From a phenomenological to a first-principles description of orbital magnetism in density-functional theory

Phenomenological orbital-polarizition (OP) terms have been repeatedly introduced in the single-particle equations of spin-density-functional theory, in order to improve the description of orbital magnetic moments in systems containing transition metal ions. Here we show that these ad hoc corrections can be interpreted as approximations to the exchange-correlation vector potential A(xc) of current-density functional theory (CDFT). This connection provides additional information on both approaches: phenomenological OP terms are connected to first-principles theory, leading to a rationale for their empirical success and a reassessment of their limitations and the approximations made in their derivation. Conversely, the connection of OP terms with CDFT leads to a set of simple approximations to the CDFT potential A(xc), with a number of desirable features that are absent from electron-gas-based functionals. (C) 2008 Wiley Periodicals, Inc.

Encapsulation of small magnetic clusters in fullerene cages: A density functional theory investigation within van der Waals corrections

The encapsulation of magnetic transition-metal (TM) clusters inside carbon cages (fullerenes, nanotubes) has been of great interest due to the wide range of applications, which spread from medical sensors in magnetic resonance imaging to photonic crystals. Several theoretical studies have been reported; however, our atomistic understanding of the physical properties of encapsulated magnetic TM 3d clusters is far from satisfactory. In this work, we will report general trends, derived from density functional theory within the generalized gradient approximation proposed by Perdew, Burke, and Ernzerhof (PBE), for the encapsulation properties of the TMm@C-n (TM = Fe, Co, Ni; m = 2-6, n = 60,70,80,90) systems. Furthermore, to understand the role of the van der Waals corrections to the physical properties, we employed the empirical Grimme's correction (PBE + D2). We found that both PBE and PBE + D2 functionals yield almost the same geometric parameters, magnetic and electronic properties, however, PBE + D2 strongly enhances the encapsulation energy. We found that the center of mass of the TMm clusters is displaced towards the inside C-n surfaces, except for large TMm clusters (m = 5 and 6). For few cases, e. g., Co-4 and Fe-4, the encapsulation changes the putative lowest-energy structure compared to the isolated TMm clusters. We identified few physical parameters that play an important role in the sign and magnitude of the encapsulation energy, namely, cluster size, fullerene equatorial diameter, shape, curvature of the inside C-n surface, number of TM atoms that bind directly to the inside C-n surface, and the van der Waals correction. The total magnetic moment of encapsulated TMm clusters decreases compared with the isolated TMm clusters, which is expected due to the hybridization of the d-p states, and strongly depends on the size and shape of the fullerene cages.

Ethanol and Water Adsorption on Close-Packed 3d, 4d, and 5d Transition-Metal Surfaces: A Density Functional Theory Investigation with van der Waals Correction

Nowadays, there is a great interest in the economic success of direct ethanol fuel cells; however, our atomistic understanding of the designing of stable and low-cost catalysts for the steam reforming of ethanol is still far from satisfactory, in particular due to the large number of undesirable intermediates. In this study, we will report a first-principles investigation of the adsorption properties of ethanol and water at low coverage on close-packed transition-metal (TM) surfaces, namely, Fe(110), Co(0001), Ni(111), Cu(111), Ru(0001), Rh(111), Pd(111), Ag(111), Os(0001), Ir(111), Pt(111), and Au(111), employing density functional theory (DFT) calculations. We employed the generalized gradient approximation with the formulation proposed by Perdew, Burke, and Erzenholf (PBE) to the exchange correlation functional and the empirical correction proposed by S. Grimme (DFT+D3) for the van der Waals correction. We found that both adsorbates binds preferentially near or on the on top sites of the TM surfaces through the 0 atoms. The PBE adsorption energies of ethanol and water decreases almost linearly with the increased occupation of the 4d and 5d d-band, while there is a deviation for the 3d systems. The van der Waals correction affects the linear behavior and increases the adsorption energy for both adsorbates, which is expected as the van der Waals energy due to the correlation effects is strongly underestimated by DFT-PBE for weak interacting systems. The geometric parameters for water/TM are not affected by the van der Waals correction, i.e., both DFT and DFT+D3 yield an almost parallel orientation for water on the TM surfaces; however, DFT+D3 changes drastically the ethanol orientation. For example, DFT yields an almost perpendicular orientation of the C-C bond to the TM surface, while the C-C bond is almost parallel to the surface using DFT +D3 for all systems, except for ethanol/Fe(110). Thus, the van der Waals correction decreases the distance of the C atoms to the TM surfaces, which might contribute to break the C-C bond. The work function decreases upon the adsorption of ethanol and water, and both follow the same trends, however, with different magnitude (larger for ethanol/TM) due to the weak binding of water to the surface. The electron density increases mainly in the region between the topmost layer and the adsorbates, which explains the reduction of the substrate work function.

A Microscopic Theory of the f.c.c.-b.c.c. Interface .

We present a first-principles theory of the equilibrium b.c.c.-f.c.c. interface at coexistence using the density functional method. We assume that the interfacial region has local body-centred tetragonal (b.c.t.) symmetry and predict typical interfacial widths to be of order 2 to 3 lattice spacings with typical energies close to 0.05 J/m2. These quantities are in good agreement with laboratory measurements on coherent interfaces.

Quantum Confinement Effects and Electronic Properties of SnO2 Quantum Wires and Dots

On the basis of the density functional theory (DFT) within local density approximations (LDA) approach, we calculate the band gaps for different size SnO2 quantum wire (QWs) and quantum dots (QDs). A model is proposed to passivate the surface atoms of SnO2 QWs and QDs. We find that the band gap increases between QWs and bulk evolve as Delta E-g(wire) = 1.74/d(1.20) as the effective diameter d decreases, while being Delta E-g(dot) = 2.84/d(1.26) for the QDs. Though the similar to d(1.2) scale is significantly different from similar to d(2) of the effective mass result, the ratio of band gap increases between SnO2 QWs and QDs is 0.609, very close to the effective mass prediction. We also confirm, although the LDS calculations underestimate the band gap, that they give the trend of band gap shift as much as that obtained by the hybrid functional (PBE0) with a rational mixing of 25% Fock exchange and 75% of the conventional Perdew-Burke-Ernzerhof (PBE) exchange functional for the SnO2 QWs and QDs. The relative deviation of the LDA calculated band gap difference Lambda E-g compared with the corresponding PBE0 results is only within 5%. Additionally, it is found the states of valence band maximum (VBM) and conduction band minimum (CBM) of SnO2 QWs or QDs have a mostly p- and s-like envelope function symmetry, respectively, from both LDA and PBE0 calculations.

Studies of tetragonal PbTiO3 subjected to uniaxial stress along the c-axis

Tetragonal PbTiO3 under uniaxial stress along the c-axis is investigated from first-principles. The structural parameters, polarization, and squares of the lowest optical phonon frequencies for E(1TO) and A(1)(1TO) modes at Gamma show abrupt changes near a stress sigma(c) of 1.04 GPa, which is related to the dramatic change of elastic constant c(33) resulting from the uniaxial stress applied along the c-axis. We also find that the uniaxial compressive stress could enhance the piezoelectric stress coefficients, whereas the uniaxial tensile stress could enhance the piezoelectric strain coefficients. It is also found that when the magnitude of uniaxial compressive stress sigma(33) is greater than 12 GPa, PbTiO3 is transformed to the paraelectric tetragonal phase.