Selected configuration interaction with truncation energy error and application to the Ne atom

Selected configuration interaction (SCI) for atomic and molecular electronic structure calculations is reformulated in a general framework encompassing all CI methods. The linked cluster expansion is used as an intermediate device to approximate CI coefficients BK of disconnected configurations (those that can be expressed as products of combinations of singly and doubly excited ones) in terms of CI coefficients of lower-excited configurations where each K is a linear combination of configuration-state-functions (CSFs) over all degenerate elements of K. Disconnected configurations up to sextuply excited ones are selected by Brown's energy formula, ΔEK=(E-HKK)BK2/(1-BK2), with BK determined from coefficients of singly and doubly excited configurations. The truncation energy error from disconnected configurations, Δdis, is approximated by the sum of ΔEKS of all discarded Ks. The remaining (connected) configurations are selected by thresholds based on natural orbital concepts. Given a model CI space M, a usual upper bound ES is computed by CI in a selected space S, and EM=E S+ΔEdis+δE, where δE is a residual error which can be calculated by well-defined sensitivity analyses. An SCI calculation on Ne ground state featuring 1077 orbitals is presented. Convergence to within near spectroscopic accuracy (0.5 cm-1) is achieved in a model space M of 1.4× 109 CSFs (1.1 × 1012 determinants) containing up to quadruply excited CSFs. Accurate energy contributions of quintuples and sextuples in a model space of 6.5 × 1012 CSFs are obtained. The impact of SCI on various orbital methods is discussed. Since ΔEdis can readily be calculated for very large basis sets without the need of a CI calculation, it can be used to estimate the orbital basis incompleteness error. A method for precise and efficient evaluation of ES is taken up in a companion paper

Analysis of chemical bonding and aromaticity from electronic delocalization descriptors

Interactions between electrons determine the structure and properties of matter from molecules to solids. Therefore, the understanding of the electronic structure of molecules will enable us to extract relevant chemical information. In the first part of this thesis, we focus our attention on the analysis of chemical bonding by means of the Electron Localization Function (ELF) and the Domain-Averaged Fermi Hole analysis (DAFH). In the second part, we assess the performance of some indicators of aromaticity by analyzing their advantages and drawbacks. We propose a series of tests based on well-known aromaticity trends that can be applied to evaluate the aromaticity of current and future indicators of aromaticity in both organic and inorganic species. Moreover, we investigate the nature of electron delocalization in both aromatic and antiaromatic systems in the light of Hückel’s (4n + 2) rule. Finally, we analyze the phenomenon of multiple aromaticity in all-metal clusters.

Theoretical studies of systems of biochemical interest containing Fe and Cu transition metals

La presència de la química teòrica i computacional està augmentant en quasi tots els camps de la recerca en química. Els càlculs teòrics poden ajudar a entendre millor l'estructura, les propietats i la reactivitat de compostos metàl·lics d'àrees tan diferents com la química inorgànica, organometàl·lica i bioinorgànica. No obstant això, és imprescindible utilitzar la metodologia adequada per obtenir resultats teòrics fiables. Els estudis d'aquesta tesi es poden dividir en dos grups diferents. El primer grup inclou l'estudi teòric del mecanisme de reacció de diversos sistemes que contenen coure i tenen diferents estructures Cun-O2. Aquests estudies s'han dut a terme amb l'objectiu de profunditzar en la natura dels processos oxidants químics i biològics promoguts per sistemes que contenen coure. En la segona part de la tesi, s'estudia la fiabilitat de diferents tècniques utilitzades per estudiar l'estructura electrònica i la reactivitat de sistemes que contenen coure, ferro i altres metalls de transició.

Critical analysis and extension of the Hirshfeld atoms in molecules

The computational approach to the Hirshfeld [Theor. Chim. Acta 44, 129 (1977)] atom in a molecule is critically investigated, and several difficulties are highlighted. It is shown that these difficulties are mitigated by an alternative, iterative version, of the Hirshfeld partitioning procedure. The iterative scheme ensures that the Hirshfeld definition represents a mathematically proper information entropy, allows the Hirshfeld approach to be used for charged molecules, eliminates arbitrariness in the choice of the promolecule, and increases the magnitudes of the charges. The resulting "Hirshfeld-I charges" correlate well with electrostatic potential derived atomic charges

Erratum: "The aromatic fluctuation index (FLU): a new aromaticity index based on electron delocalization" [J. Chem Phys. 122, 014109 (2005)]

An erratum to the article "The aromatic fluctuation index (FLU): a new aromaticity index based on electron delocalization", published in The Journal Chemistry of Physics, 2005, v.122, art. no.014109. Values FLU have been corrected in the last column of Table I because they were not correct

The aromatic fluctuation index (FLU): a new aromaticity index based on electron delocalization

In this work, the aromatic fluctuation index (FLU) that describes the fluctuation of electronic charge between adjacent atoms in a given ring is introduced as a new aromaticity measure. This new electronic criterion of aromaticity is based on the fact that aromaticity is related to the cyclic delocalized circulation of π electrons. It is defined not only considering the amount of electron sharing between contiguous atoms, which should be substantial in aromatic molecules, but also taking into account the similarity of electron sharing between adjacent atoms. For a series of rings in 15 planar polycyclic aromatic hydrocarbons, we have found that, in general, FLU is strongly correlated with other widely used indicators of local aromaticity, such as the harmonic-oscillator model of aromaticity, the nucleus independent chemical shift, and the para-delocalization index (PDI). In contrast to PDI, the FLU index can be applied to study the aromaticity of rings with any number of members and it can be used to analyze both the local and global aromatic character of rings and molecules

Are the maximum hardness and minimum polarizability principles always obeyed in nontotally symmetric vibrations?

An overview is given on a study which showed that not only in chemical reactions but also in the favorable case of nontotally symmetric vibrations where the chemical and external potentials keep approximately constant, the generalized maximum hardness principle (GMHP) and generalized minimum polarizability principle (GMPP) may not be obeyed. A method that allows an accurate determination of the nontotally symmetric molecular distortions with more marked GMPP or anti-GMPP character through diagonalization of the polarizability Hessian matrix is introduced