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 ...

Physical signatures of discontinuities of the time-dependent exchange-correlation potential

The exact exchange-correlation (XC) potential in time-dependent density-functional theory (TDDFT) is known to develop steps and discontinuities upon change of the particle number in spatially confined regions or isolated subsystems. We demonstrate that the self-interaction corrected adiabatic local-density approximation for the XC potential has this property, using the example of electron loss of a model quantum well system. We then study the influence of the XC potential discontinuity in a real-time simulation of a dissociation process of an asymmetric double quantum well system, and show that it dramatically affects the population of the resulting isolated single quantum wells. This indicates the importance of a proper account of the discontinuities in TDDFT descriptions of ionization, dissociation or charge transfer processes.

Antiferromagnetic exchange interactions from hybrid density functional theory

A hybrid theory which combines the full nonlocal ¿exact¿ exchange interaction with the local spin-density approximation of density-functional theory is shown to lead to marked improvement in the description of antiferromagnetically coupled systems. Semiquantitative agreement with experiment is found for the magnitude of the coupling constant in La2CuO4, KNiF3, and K2NiF4. The magnitude of the unpaired spin population on the metal site is in excellent agreement with experiment for La2CuO4.

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.

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.

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.

É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.

Multiscale modeling of light-harvesting systems based on rhodamine B aggregates

Rhodamine B (RB) has been successfully exploited in the synthesis of light harvesting systems, but since RB is prone to form dimers acting as quenchers for the fluorescence, high energy transfer efficiencies can be reached only when using bulky and hydrophobic counterions acting as spacers between RBs. In this PhD thesis, a multiscale theoretical study aimed at providing insights into the structural, photophysical and optical properties of RB and its aggregates is presented. At the macroscopic level (no atomistic details) a phenomenological model describing the fluorescence decay of RB networks in presence of both quenching from dimers and exciton-exciton annihiliation is presented and analysed, showing that the quenching from dimers affects the decay only at long times, a feature that can be exploited in global fitting analysis to determine relevant chemical and photophysical information. At the mesoscopic level (atomistic details but no electronic structure) the RB aggregation in water in presence of different counterions is studied with molecular dynamics (MD) simulations. A new force field has been parametrized for describing the RB flexibility and the RB-RB interaction driving the dimerization. Simulations correctly predict the RB/counterion aggregation only in presence of bulky and hydrophobic counterion and its ability to prevent the dimerization. Finally, at the microscopic level, DFT calculations are performed to demonstrate the spacing action of bulky counterions, but standard TDDFT calculations are showed to fail in correctly describing the excited states of RB and its dimers. Moreover, also standard procedures proposed in literature for obtaining ad hoc functionals are showed to not work properly. A detailed analysis on the effect of the exact exchange shows that its short-range contribution is the crucial quantity for ameliorating results, and a new functional containing a proper amount of such an exchange is proposed and successfully tested.

EXACT AND APPROXIMATE RELATIONS FOR THE SPIN-DEPENDENCE OF THE EXCHANGE ENERGY IN HIGH MAGNETIC FIELDS

The exchange energy of an arbitrary collinear-spin many-body system in an external magnetic field is a functional of the spin-resolved charge and current densities, E(x)[n(up arrow), n(down arrow), j(up arrow), j(down arrow)]. Within the framework of density-functional theory (DFT), we show that the dependence of this functional on the four densities can be fully reconstructed from either of two extreme limits: a fully polarized system or a completely unpolarized system. Reconstruction from the limit of an unpolarized system yields a generalization of the Oliver-Perdew spin scaling relations from spin-DFT to current-DFT. Reconstruction from the limit of a fully polarized system is used to derive the high-field form of the local-spin-density approximation to current-DFT and to magnetic-field DFT.

The capacity exchange protocol

This paper proposes a new strategy to integrate shared resources and precedence constraints among real-time tasks, assuming no precise information on critical sections and computation times is available. The concept of bandwidth inheritance is combined with a capacity sharing and stealing mechanism to efficiently exchange bandwidth among tasks to minimise the degree of deviation from the ideal system’s behaviour caused by inter-application blocking. The proposed Capacity Exchange Protocol (CXP) is simpler than other proposed solutions for sharing resources in open real-time systems since it does not attempt to return the inherited capacity in the same exact amount to blocked servers. This loss of optimality is worth the reduced complexity as the protocol’s behaviour nevertheless tends to be fair and outperforms the previous solutions in highly dynamic scenarios as demonstrated by extensive simulations. A formal analysis of CXP is presented and the conditions under which it is possible to guarantee hard real-time tasks are discussed.

Exact Skewness-Kurtosis Tests for Multivariate Normality and Goodness-of-fit in Multivariate Regressions with Application to Asset Pricing Models

We study the problem of testing the error distribution in a multivariate linear regression (MLR) model. The tests are functions of appropriately standardized multivariate least squares residuals whose distribution is invariant to the unknown cross-equation error covariance matrix. Empirical multivariate skewness and kurtosis criteria are then compared to simulation-based estimate of their expected value under the hypothesized distribution. Special cases considered include testing multivariate normal, Student t; normal mixtures and stable error models. In the Gaussian case, finite-sample versions of the standard multivariate skewness and kurtosis tests are derived. To do this, we exploit simple, double and multi-stage Monte Carlo test methods. For non-Gaussian distribution families involving nuisance parameters, confidence sets are derived for the the nuisance parameters and the error distribution. The procedures considered are evaluated in a small simulation experi-ment. Finally, the tests are applied to an asset pricing model with observable risk-free rates, using monthly returns on New York Stock Exchange (NYSE) portfolios over five-year subperiods from 1926-1995.

Exact Multivariate Tests of Asset Pricing Models with Stable Asymmetric Distributions

In this paper, we propose exact inference procedures for asset pricing models that can be formulated in the framework of a multivariate linear regression (CAPM), allowing for stable error distributions. The normality assumption on the distribution of stock returns is usually rejected in empirical studies, due to excess kurtosis and asymmetry. To model such data, we propose a comprehensive statistical approach which allows for alternative - possibly asymmetric - heavy tailed distributions without the use of large-sample approximations. The methods suggested are based on Monte Carlo test techniques. Goodness-of-fit tests are formally incorporated to ensure that the error distributions considered are empirically sustainable, from which exact confidence sets for the unknown tail area and asymmetry parameters of the stable error distribution are derived. Tests for the efficiency of the market portfolio (zero intercepts) which explicitly allow for the presence of (unknown) nuisance parameter in the stable error distribution are derived. The methods proposed are applied to monthly returns on 12 portfolios of the New York Stock Exchange over the period 1926-1995 (5 year subperiods). We find that stable possibly skewed distributions provide statistically significant improvement in goodness-of-fit and lead to fewer rejections of the efficiency hypothesis.