Approximation to density functional theory for the calculation of band gaps of semiconductors

The local-density approximation (LDA) together with the half occupation (transitionstate) is notoriously successful in the calculation of atomic ionization potentials. When it comes to extended systems, such as a semiconductor infinite system, it has been very difficult to find a way to half ionize because the hole tends to be infinitely extended (a Bloch wave). The answer to this problem lies in the LDA formalism itself. One proves that the half occupation is equivalent to introducing the hole self-energy (electrostatic and exchange correlation) into the Schrodinger equation. The argument then becomes simple: The eigenvalue minus the self-energy has to be minimized because the atom has a minimal energy. Then one simply proves that the hole is localized, not infinitely extended, because it must have maximal self-energy. Then one also arrives at an equation similar to the self- interaction correction equation, but corrected for the removal of just 1/2 electron. Applied to the calculation of band gaps and effective masses, we use the self- energy calculated in atoms and attain a precision similar to that of GW, but with the great advantage that it requires no more computational effort than standard LDA.

Correlation Between Molecular Conformation, Packing, and Dynamics in Oligofluorenes: A Theoretical/Experimental Study

Fluorene-based systems have shown great potential as components in organic electronics and optoelectronics (organic photovoltaics, OPVs, organic light emitting diodes, OLEDs, and organic transistors, OTFTs). These systems have drawn attention primarily because they exhibit strong blue emission associated with relatively good thermal stability. It is well-known that the electronic properties of polymers are directly related to the molecular conformations and chain packing of polymers. Here, we used three oligofluorenes (trimer, pentamer, and heptamer) as model systems to theoretically investigate the conformational properties of fluorene molecules, starting with the identification of preferred conformations. The hybrid exchange correlation functional, OPBE, and ZINDO/S-CI showed that each oligomer exhibits a tendency to adopt a specific chain arrangement, which could be distinguished by comparing their UV/vis electronic absorption and C-13 NMR spectra. This feature was used to identify the preferred conformation of the oligomer chains in chloroform-cast films by comparing experimental and theoretical UV/vis and C-13 NMR spectra. Moreover, the oligomer chain packing and dynamics in the films were studied by DSC and several solid state NMR techniques, which indicated that the phase behavior of the films may be influenced by the tendency that each oligomeric chain has to adopt a given conformation.

Adiabatic and local approximations for the kohn-sham potential in the hubbard model

We obtain the exact time-dependent Kohn-Sham potentials Vks for 1D Hubbard chains, driven by a d.c. external field, using the time-dependent electron density and current density obtained from exact many-body time-evolution. The exact Vxc is compared to the adiabatically-exact Vad-xc and the “instantaneous ground state” Vigs-xc. The effectiveness of these two approximations is analyzed. Approximations for the exchange-correlation potential Vxc and its gradient, based on the local density and on the local current density, are also considered and both physical quantities are observed to be far outside the reach of any possible local approximation. Insight into the respective roles of ground-state and excited-state correlation in the time-dependent system, as reflected in the potentials, is provided by the pair correlation function.

Polarized–unpolarized ground state of small polycyclic aromatic hydrocarbons

Do polyacenes, circumacenes, periacenes, nanographenes, and graphene nanoribbons show a spin polarized ground state? In this work, we present monodeterminantal (Hartree–Fock (HF) and density functional theory (DFT) types), and multideterminantal calculations (Møller–Plesset and Coupled Cluster), for several families of unsaturated organic molecules (n-Acenes, n-Periacenes and n-Circumacenes). All HF calculations and many DFT show a spin-polarized (antiferromagnetic) ground state, in agreement with previous calculations. Nevertheless, the multideterminantal calculations, carried out with perturbative and variational wavefunctions, show that the more stable state is obtained starting from the unpolarized HF wavefunction. The trend of the stabilization of wavefunctions (polarized or unpolarized) with respect to exchange and correlation potentials, and to the number of benzene rings, has been analyzed. A study of the spin (〈Ŝ2〉) and the spin density on the carbon atoms has also been carried out.

Mn-doped cubic BN as an atomiclike memory device: A density functional study

We investigate the electronic properties of Mn(B) substitutional doping in cubic boron nitride (BN), for different charge states, using density functional theory (DFT) calculations. We show that the neutral Mn has a nonmagnetic ground state (S=0). Upon charge injection, it is unambiguously shown that the Mn(B)(-) has a high-spin configuration with a strong, localized magnetic moment of 5 mu(Bohr). We developed a simple model, parameterized by the DFT results, that allows us to interpret the rules played by the crystal-field and exchange-correlation splitting in the magnetization process.

Many-Body Effects on the rho(xx) Ringlike Structures in Two-Subband Wells

The longitudinal resistivity rho(xx) of two-dimensional electron gases formed in wells with two subbands displays ringlike structures when plotted in a density-magnetic-field diagram, due to the crossings of spin-split Landau levels (LLs) from distinct subbands. Using spin density functional theory and linear response, we investigate the shape and spin polarization of these structures as a function of temperature and magnetic-field tilt angle. We find that (i) some of the rings ""break'' at sufficiently low temperatures due to a quantum Hall ferromagnetic phase transition, thus exhibiting a high degree of spin polarization (similar to 50%) within, consistent with the NMR data of Zhang et al. [Phys. Rev. Lett. 98, 246802 (2007)], and (ii) for increasing tilting angles the interplay between the anticrossings due to inter-LL couplings and the exchange-correlation effects leads to a collapse of the rings at some critical angle theta(c), in agreement with the data of Guo et al. [Phys. Rev. B 78, 233305 (2008)].

Nonempirical hyper-generalized-gradient functionals constructed from the Lieb-Oxford bound

A simple and completely general representation of the exact exchange-correlation functional of density-functional theory is derived from the universal Lieb-Oxford bound, which holds for any Coulomb-interacting system. This representation leads to an alternative point of view on popular hybrid functionals, providing a rationale for why they work and how they can be constructed. A similar representation of the exact correlation functional allows to construct fully nonempirical hyper-generalized-gradient approximations (HGGAs), radically departing from established paradigms of functional construction. Numerical tests of these HGGAs for atomic and molecular correlation energies and molecular atomization energies show that even simple HGGAs match or outperform state-of-the-art correlation functionals currently used in solid-state physics and quantum chemistry.

Transition-metal 13-atom clusters assessed with solid and surface-biased functionals

First-principles density-functional theory studies have reported open structures based on the formation of double simple-cubic (DSC) arrangements for Ru(13), Rh(13), Os(13), and Ir(13), which can be considered an unexpected result as those elements crystallize in compact bulk structures such as the face-centered cubic and hexagonal close-packed lattices. In this work, we investigated with the projected augmented wave method the dependence of the lowest-energy structure on the local and semilocal exchange-correlation (xc) energy functionals employed in density-functional theory. We found that the local-density approximation (LDA) and generalized-gradient formulations with different treatment of the electronic inhomogeneities (PBE, PBEsol, and AM05) confirm the DSC configuration as the lowest-energy structure for the studied TM(13) clusters. A good agreement in the relative total energies are obtained even for structures with small energy differences, e. g., 0.10 eV. The employed xc functionals yield the same total magnetic moment for a given structure, i.e., the differences in the bond lengths do not affect the moments, which can be attributed to the atomic character of those clusters. Thus, at least for those systems, the differences among the LDA, PBE, PBEsol, and AM05 functionals are not large enough to yield qualitatively different results. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3577999]

The role of electron localization in the atomic structure of transition-metal 13-atom clusters: the example of Co(13), Rh(13), and Hf(13)

The crystalline structure of transition-metals (TM) has been widely known for several decades, however, our knowledge on the atomic structure of TM clusters is still far from satisfactory, which compromises an atomistic understanding of the reactivity of TM clusters. For example, almost all density functional theory (DFT) calculations for TM clusters have been based on local (local density approximation-LDA) and semilocal (generalized gradient approximation-GGA) exchange-correlation functionals, however, it is well known that plain DFT fails to correct the self-interaction error, which affects the properties of several systems. To improve our basic understanding of the atomic and electronic properties of TM clusters, we report a DFT study within two nonlocal functionals, namely, the hybrid HSE (Heyd, Scuseria, and Ernzerhof) and GGA + U functionals, of the structural and electronic properties of the Co(13), Rh(13), and Hf(13) clusters. For Co(13) and Rh(13), we found that improved exchange-correlation functionals decrease the stability of open structures such as the hexagonal bilayer (HBL) and double simple-cubic (DSC) compared with the compact icosahedron (ICO) structure, however, DFT-GGA, DFT-GGA + U, and DFT-HSE yield very similar results for Hf(13). Thus, our results suggest that the DSC structure obtained by several plain DFT calculations for Rh(13) can be improved by the use of improved functionals. Using the sd hybridization analysis, we found that a strong hybridization favors compact structures, and hence, a correct description of the sd hybridization is crucial for the relative energy stability. For example, the sd hybridization decreases for HBL and DSC and increases for ICO in the case of Co(13) and Rh(13), while for Hf(13), the sd hybridization decreases for all configurations, and hence, it does not affect the relative stability among open and compact configurations.

Density functional energy decomposition into one- and two-atom contributions

The present work provides a generalization of Mayer's energy decomposition for the density-functional theory (DFT) case. It is shown that one- and two-atom Hartree-Fock energy components in Mayer's approach can be represented as an action of a one-atom potential VA on a one-atom density ρ A or ρ B. To treat the exchange-correlation term in the DFT energy expression in a similar way, the exchange-correlation energy density per electron is expanded into a linear combination of basis functions. Calculations carried out for a number of density functionals demonstrate that the DFT and Hartree-Fock two-atom energies agree to a reasonable extent with each other. The two-atom energies for strong covalent bonds are within the range of typical bond dissociation energies and are therefore a convenient computational tool for assessment of individual bond strength in polyatomic molecules. For nonspecific nonbonding interactions, the two-atom energies are low. They can be either repulsive or slightly attractive, but the DFT results more frequently yield small attractive values compared to the Hartree-Fock case. The hydrogen bond in the water dimer is calculated to be between the strong covalent and nonbonding interactions on the energy scale

Extent and limitations of density functional theory in describing magnetic systems

The performance of density-functional theory to solve the exact, nonrelativistic, many-electron problem for magnetic systems has been explored in a new implementation imposing space and spin symmetry constraints, as in ab initio wave function theory. Calculations on selected systems representative of organic diradicals, molecular magnets and antiferromagnetic solids carried out with and without these constraints lead to contradictory results, which provide numerical illustration on this usually obviated problem. It is concluded that the present exchange-correlation functionals provide reasonable numerical results although for the wrong physical reasons, thus evidencing the need for continued search for more accurate expressions.

Extent and limitations of density functional theory in describing magnetic systems

The performance of density-functional theory to solve the exact, nonrelativistic, many-electron problem for magnetic systems has been explored in a new implementation imposing space and spin symmetry constraints, as in ab initio wave function theory. Calculations on selected systems representative of organic diradicals, molecular magnets and antiferromagnetic solids carried out with and without these constraints lead to contradictory results, which provide numerical illustration on this usually obviated problem. It is concluded that the present exchange-correlation functionals provide reasonable numerical results although for the wrong physical reasons, thus evidencing the need for continued search for more accurate expressions.

Validação computacional de métodos compostos no estudo de propriedades moleculares

Composite methods using ONIOM and different basis sets have been used to calculate proton and electron affinities for a set of alcohols at QCISD(T)/6-311++G(2df,p) level of theory. The study was carried out considering HF, MP2 and DFT (25 exchange correlation functional) methods. The calculation performed at ONIOM2(QCISD(T)/6-311++G(2df,p):HF/6-31G(d))//ONIOM2(O3LYP/6-31G(d):HF/6-31G(d)) resulted in the smallest average absolute deviation for AP and AE, 4,75 kJ/mol e 0,43 eV, respectively.

The role of van der Waals interactions in the case of H2 adsorption on ruthenium (0001) surface

Computational material science with the Density Functional Theory (DFT) has recently gained a method for describing, for the first time the non local bonding i.e., van der Waals (vdW) bonding. The newly proposed van der Waals-Density Functional (vdW-DF) is employed here to address the role of non local interactions in the case of H2 adsorption on Ru(0001) surface. The later vdW-DF2 implementation with the DFT code VASP (Vienna Ab-initio Simulation Package) is used in this study. The motivation for studying H2 adsorption on ruthenium surface arose from the interest to hydrogenation processes. Potential energy surface (PES) plots are created for adsorption sites top, bridge, fcc and hcp, employing the vdW-DF2 functional. The vdW-DF yields 0.1 eV - 0.2 eV higher barriers for the dissociation of the H2 molecule; the vdW-DF seems to bind the H2 molecule more tightly together. Furthermore, at the top site, which is found to be the most reactive, the vdW functional suggests no entrance barrier or in any case smaller than 0.05 eV, whereas the corresponding calculation without the vdW-DF does. Ruthenium and H2 are found to have the opposite behaviors with the vdW-DF; Ru lattice constants are overestimated while H2 bond length is shorter. Also evaluation of the CPU time demand of the vdW-DF2 is done from the PES data. From top to fcc sites the vdW-DF computational time demand is larger by 4.77 % to 20.09 %, while at the hcp site it is slightly smaller. Also the behavior of a few exchange correlation functionals is investigated along addressing the role of vdW-DF. Behavior of the different functionals is not consistent between the Ru lattice constants and H2 bond lengths. It is thus difficult to determine the quality of a particular exchange correlation functional by comparing equilibrium separations of the different elements. By comparing PESs it would be computationally highly consuming.

Extent and limitations of density functional theory in describing magnetic systems

The performance of density-functional theory to solve the exact, nonrelativistic, many-electron problem for magnetic systems has been explored in a new implementation imposing space and spin symmetry constraints, as in ab initio wave function theory. Calculations on selected systems representative of organic diradicals, molecular magnets and antiferromagnetic solids carried out with and without these constraints lead to contradictory results, which provide numerical illustration on this usually obviated problem. It is concluded that the present exchange-correlation functionals provide reasonable numerical results although for the wrong physical reasons, thus evidencing the need for continued search for more accurate expressions.