Screened Hybrid Exact Exchange Correction Scheme for Adsorption Energies on Perovskite Oxides

The bond formation between an oxide surface and oxygen, which is of importance for numerous surface reactions including catalytic reactions, is investigated within the framework of hybrid density functional theory that includes nonlocal Fock exchange. We show that there exists a linear correlation between the adsorption energies of oxygen on LaMO3 (M = Sc–Cu) surfaces obtained using a hybrid functional (e.g., Heyd–Scuseria–Ernzerhof) and those obtained using a semilocal density functional (e.g., Perdew–Burke–Ernzerhof) through the magnetic properties of the bulk phase as determined with a hybrid functional. The energetics of the spin-polarized surfaces follows the same trend as corresponding bulk systems, which can be treated at a much lower computational cost. The difference in adsorption energy due to magnetism is linearly correlated to the magnetization energy of bulk, that is, the energy difference between the spin-polarized and the non-spin-polarized solutions. Hence, one can estimate the correction ...

Exchange potential from the common energy denominator approximation for the Kohn-Sham Green's function: Application to (hyper)polarizabilities of molecular chains

An approximate Kohn-Sham (KS) exchange potential v(xsigma)(CEDA) is developed, based on the common energy denominator approximation (CEDA) for the static orbital Green's function, which preserves the essential structure of the density response function. v(xsigma)(CEDA) is an explicit functional of the occupied KS orbitals, which has the Slater v(Ssigma) and response v(respsigma)(CEDA) potentials as its components. The latter exhibits the characteristic step structure with "diagonal" contributions from the orbital densities \psi(isigma)\(2), as well as "off-diagonal" ones from the occupied-occupied orbital products psi(isigma)psi(j(not equal1)sigma). Comparison of the results of atomic and molecular ground-state CEDA calculations with those of the Krieger-Li-Iafrate (KLI), exact exchange (EXX), and Hartree-Fock (HF) methods show, that both KLI and CEDA potentials can be considered as very good analytical "closure approximations" to the exact KS exchange potential. The total CEDA and KLI energies nearly coincide with the EXX ones and the corresponding orbital energies epsilon(isigma) are rather close to each other for the light atoms and small molecules considered. The CEDA, KLI, EXX-epsilon(isigma) values provide the qualitatively correct order of ionizations and they give an estimate of VIPs comparable to that of the HF Koopmans' theorem. However, the additional off-diagonal orbital structure of v(xsigma)(CEDA) appears to be essential for the calculated response properties of molecular chains. KLI already considerably improves the calculated (hyper)polarizabilities of the prototype hydrogen chains H-n over local density approximation (LDA) and standard generalized gradient approximations (GGAs), while the CEDA results are definitely an improvement over the KLI ones. The reasons of this success are the specific orbital structures of the CEDA and KLI response potentials, which produce in an external field an ultranonlocal field-counteracting exchange potential. (C) 2002 American Institute of Physics.

Systematic construction of models for the exchange hole of density functional theory.

Dans ce travail, nous étendons le nombre de conditions physiques actuellement con- nues du trou d’échange exact avec la dérivation de l’expansion de quatrième ordre du trou d’échange sphérique moyenne exacte. Nous comparons les expansions de deux- ième et de quatrième ordre avec le trou d’échange exact pour des systèmes atomiques et moléculaires. Nous avons constaté que, en général, l’expansion du quatrième ordre reproduit plus fidèlement le trou d’échange exact pour les petites valeurs de la distance interélectronique. Nous démontrons que les ensembles de base de type gaussiennes ont une influence significative sur les termes de cette nouvelle condition, en étudiant com- ment les oscillations causées par ces ensembles de bases affectent son premier terme. Aussi, nous proposons quatre modèles de trous d’échange analytiques auxquels nous imposons toutes les conditions actuellement connues du trou d’échange exact et la nou- velle présentée dans ce travail. Nous évaluons la performance des modèles en calculant des énergies d’échange et ses contributions à des énergies d’atomisation. On constate que les oscillations causeés par les bases de type gaussiennes peuvent compromettre la précision et la solution des modèles.

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.

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 B3LYP

Theoretical investigation of LaC2 and LaC2+ clusters

LaC2 (with doublet and quartet states) and LaC2+ (with singlet and triplet states) cluster have been studied by using the B3LYP (Becke three-parameter/Lee-Yang-Parr) density functional method and the HF (Hartree-Fock) method with LANLIDZ basis set. For each cluster, four possible isomers in C-2v, C-s, C-proportional to v and D-proportional to h symmetries have been investigated. The results indicate that structures in C-s symmetry are local minima in all cases and, in most cases (particularly for high spin states), our initial guess in C-s symmetry converges to structures in C-2v symmetry. For the isomers in C-2v, C-proportional to v and D-proportional to h symmetries, local minima were found to be dependent on the method and spin state. The two clusters may also exist as linear chains. The ordering of the binding energies for the isomers in all spin states is C-s similar to C-2v < C-proportional to v < D-proportional to h. The ionization potential of LaC2 is reported as well. (C) 1998 Elsevier Science B.V.

Theoretical investigation of LaC3n+ (n = 0, 1, 2) clusters by density functional theory

LaC3n+ (n = 0, 1, 2) clusters have been studied using B3LYP (Becke 3-parameter-Lee-Yang-Parr) density functional method. The basis set is Dunning/ Huzinaga valence double zeta for carbon and [2s2p2d] for lanthanum, denoted LANL1DZ. Four isomers are presented for each cluster; two of them are edge binding isomers with C-2 upsilon symmetry, the other two are Linear chains with C-infinity upsilon symmetry. Meanwhile, two spin states for each isomer, that is, singlet and triplet for LaC3+, doublet and quartet for LaC3 and LaC32+, respectively, are also considered. Geometries, vibrational frequencies, infrared intensities, and other quantities are reported and discussed. The results indicate that at some spin states; the C-2 upsilon symmetry isomers are the dominant structures, while for the other spin states, linear isomers are energetically favored. (C) 1998 John Wiley & Sons, Inc.

Theoretical investigation of LaCn+ (n=2-8) clusters

LaCn+ (n = 2-8) have been studied using Hartree-Fock (HF) and B3LYP density functional method. The results indicated that at both levels, isomers with C-2v, C-s symmetry for n = 2, and edge insertion isomer for n = 4, 6, 8, as well as edge binding isomer for n = 3, 5, 7 were found as ground states. This is in good agreement with experimental results. The exceptional case is for n = 6 at B3LYP level, in which edge insertion and edge binding isomers were computed to be near isoenergetic. (C) 1997 Elsevier Science B.V.

Ab initio study of C-9 and LaC9+ clusters

The energies and geometries of C-9 and LaC9+ clusters were calculated at HF, MP2 and DFT levels. For C-9, all theoretical levels show that the linear chain is the most stable structure. For LaC9+, two isomers were considered. In the first case La has two single bonds (A), while it forms a double bond in the second (B). Our results showed that in HF calculation, B is marginally more stable than A, while for MP2 and DFT, A is favored. Our results also revealed that there is not enough space for C-9 ring to accommodate lanthanum. Our conclusion agrees well with experiment.

Density functionals from many-body perturbation theory: The band gap for semiconductors and insulators

Theoretically the Kohn-Sham band gap differs from the exact quasiparticle energy gap by the derivative discontinuity of the exchange-correlation functional. In practice for semiconductors and insulators the band gap calculated within any local or semilocal density approximations underestimates severely the experimental energy gap. On the other hand, calculations with an "exact" exchange potential derived from many-body perturbation theory via the optimized effective potential suggest that improving the exchange-correlation potential approximation can yield a reasonable agreement between the Kohn-Sham band gap and the experimental gap. The results in this work show that this is not the case. In fact, we add to the exact exchange the correlation that corresponds to the dynamical (random phase approximation) screening in the GW approximation. This accurate exchange-correlation potential provides band structures similar to the local density approximation with the corresponding derivative discontinuity that contributes 30%-50% to the energy gap. Our self-consistent results confirm substantially the results for Si and other semiconductors obtained perturbatively [R. W. Godby , Phys. Rev. B 36, 6497 (1987)] and extend the conclusion to LiF and Ar, a wide-gap insulator and a noble-gas solid. (c) 2006 American Institute of Physics.

Étude, par principes premiers, des effets de la corrélation entre électrons sur les propriétés électroniques et magnétiques de polymères pontés et de supraconducteurs à haute température critique

La présente thèse traite de la description de systèmes complexes, notamment des polymères et des cuprates, par la théorie de la fonctionnelle de la densité. En premier lieu, la théorie de la fonctionnelle de la densité ainsi que différentes fonctionnelles utilisées pour simuler les matériaux à l’étude sont présentées. Plus spécifiquement, les fonctionnelles LDA et GGA sont décrites et leurs limites sont exposées. De plus, le modèle de Hubbard ainsi que la fonctionnelle LDA+U qui en découle sont abordés dans ce chapitre afin de permettre la simulation des propriétés de matériaux à forte corrélation électronique. Par la suite, les résultats obtenus sur les polymères sont résumés par deux articles. Le premier traite de la variation de la bande interdite entre les polymères pontés et leurs homologues non pontés. Le second se penche sur l’étude de polymères à faible largeur de bande interdite. Dans ce dernier, il sera démontré qu’une fonctionnelle hybride, contenant de l’échange exact, est nécessaire afin de décrire les propriétés électroniques des systèmes à l’étude. Finalement, le dernier chapitre est consacré à l’étude des cuprates supraconducteurs. La LDA+U pouvant rendre compte de la forte localisation dans les orbitales 3d des atomes de cuivre, une étude de l’impact de cette fonctionnelle sur les propriétés électroniques est effectuée. Un dernier article investiguant différents ordres magnétiques dans le La2CuO4 dopé termine le dernier chapitre. On trouve aussi, en annexe, un complément d’information pour le second article et une description de la théorie de la supraconductivité de Bardeen, Cooper et Schrieffer.

Dependence of Response Functions and Orbital Functionals on Occupation Numbers

Explicitly orbital-dependent approximations to the exchange-correlation energy functional of density functional theory typically not only depend on the single-particle Kohn-Sham orbitals but also on their occupation numbers in the ground-state Slater determinant. The variational calculation of the corresponding exchange-correlation potentials with the optimized effective potential (OEP) method therefore also requires a variation of the occupation numbers with respect to a variation in the effective single-particle potential, which is usually not taken into account. Here it is shown under which circumstances this procedure is justified.

Self-consistent relativistic band structure of the CH3NH3PbI3 perovskite

The electronic structure and properties of the orthorhombic phase of the CH 3 NH 3 PbI 3 perovskite are computed with density functional theory. The structure, optimized using a van der Waals functional, reproduces closely the unit cell volume. The experimental band gap is reproduced accurately by combining spin-orbit effects and a hybrid functional in which the fraction of exact exchange is tuned self-consistently to the optical dielectric constant. Including spin-orbit coupling strongly reduces the anisotropy of the effective mass tensor, predicting a low electron effective mass in all crystal directions. The computed binding energy of the unrelaxed exciton agrees with experimental data, and the values found imply a fast exciton dissociation at ambient temperature. Also polaron masses for the separated carriers are estimated. The values of all these parameters agree with recent indications that fast dynamics and large carrier diffusion lengths are key in the high photovoltaic efficiencies shown by these materials.

Modern Computational Chemistry Methods for Prediction of Ground- and Excited-State Properties in Open-Shell Systems

The main goal of the research presented in this work is to provide some important insights about computational modeling of open-shell species. Such projects are: the investigation of the size-extensivity error in Equation-of-Motion Coupled Cluster methods, the analysis of the Long-Range corrected scheme in predicting UV-Vis spectra of Cu(II) complexes with the 4-imidazole acetate and its ethylated derivative, and the exploration of the importance of choosing a proper basis set for the description of systems such as the lithium monoxide anion. The most significant findings of this research are: (i) The contribution of the left operator to the size-extensivity error of the CR-EOMCC(2,3) approach, (ii) The cause of d-d shifts when varying the range-separation parameter and the amount of the exact exchange arising from the imbalanced treatment of localized vs. delocalized orbitals via the "tuned" CAM-B3LYP* functional, (iii) The proper acidity trend of the first-row hydrides and their lithiated analogs that may be reversed if the basis sets are not correctly selected.