Multiwfn is a multifunctional program for wavefunction analysis. Its main functions are: (1) Calculating and visualizing real space function, such as electrostatic potential and electron localization function at point, in a line, in a plane or in a spatial scope. (2) Population analysis. (3) Bond order analysis. (4) Orbital composition analysis. (5) Plot density-of-states and spectrum. (6) Topology analysis for electron density. Some other useful utilities involved in quantum chemistry studies are also provided. The built-in graph module enables the results of wavefunction analysis to be plotted directly or exported to high-quality graphic file. The program interface is very user-friendly and suitable for both research and teaching purpose. The code of Multiwfn is substantially optimized and parallelized. Its efficiency is demonstrated to be significantly higher than related programs with the same functions. Five practical examples involving a wide variety of systems and analysis methods are given to illustrate the usefulness of Multiwfn. The program is free of charge and open-source. Its precompiled file and source codes are available from http://multiwfn.codeplex.com.

Multiwfn: a multifunctional wavefunction analyzer.

: The unpolarized absorption and circular dichroism spectra of the fundamental vibrational transitions of the chiral molecule, 4-methyl-2-oxetanone, are calculated ab initio. Harmonic force fields are obtained using Density Functional Theory (DFT), MP2, and SCF methodologies and a 5S4P2D/3S2P (TZ2P) basis set. DFT calculations use the Local Spin Density Approximation (LSDA), BLYP, and Becke3LYP (B3LYP) density functionals. Mid-IR spectra predicted using LSDA, BLYP, and B3LYP force fields are of significantly different quality, the B3LYP force field yielding spectra in clearly superior, and overall excellent, agreement with experiment. The MP2 force field yields spectra in slightly worse agreement with experiment than the B3LYP force field. The SCF force field yields spectra in poor agreement with experiment.The basis set dependence of B3LYP force fields is also explored: the 6-31G* and TZ2P basis sets give very similar results while the 3-21G basis set yields spectra in substantially worse agreements with experiment. jg

Ab initio calculation of vibrational absorption and circular dichroism spectra using density functional force fields

A method of description of the chemical bonding combining the compactness and intuitive simplicity of Lewis theory with the flexibility and generality of canonical molecular orbital theory is presented, which is called adaptive natural density partitioning. The objects of chemical bonding in this method are n-center 2-electron bonds, where n goes from one (lone-pair) to the maximum number of atoms in the system (completely delocalized bonding). The algorithm is a generalization of the natural bonding orbital analysis and is based on the diagonalization of the blocks of the first-order density matrix in the basis of natural atomic orbitals. The results obtained by the application of the algorithm to the systems with non-classical bonding can be readily interpreted from the point of view of aromaticity/antiaromaticity concepts. The considered examples include Li4 cluster and a family of planar boron clusters observed in molecular beams.

Developing paradigms of chemical bonding: adaptive natural density partitioning.

Benzene is the emblematic example of an aromatic molecule, and the problem of its structure has given rise to a chemical serial story running over several decades. The epistemological digest of this story written by Stephen G. Brush1,2 shows how this problem has been at the root of important concepts such as those of mesomery and resonance. Before the advent of quantum mechanics, chemists had thought of the benzene structures in terms of two center bonds attempting to preserve the valence of the carbon atom and to explain its chemical properties. Kekulé’s theory of the structure of the benzene molecule3 invokes the oscillatory hypothesis in which “the fourth valence of each carbon oscillates between its neighbors, synchronously with all the other fourth valences, so that the structure switches rapidly between the two structures”,1 whereas Claus proposed a diagonal hypothesis4 in which the fourth valence of each carbon is directed toward the carbon in the para position. The latter hypothesis has been rejected because it enables only two derivatives, and it has been revised to remove this inconsistency: instead of forming a bond, the fourth valence stops near the center of the ring in the Armstrong-Baeyer formula,5 or there is only one bridging bond as in the Dewar’s bridged benzene formula.6 In Thiele’s partial valence model,7 the adjacent carbon-carbon bonds are considered as intermediate between single and double bonds. These formulas were later considered by K. C. Ingold to set up his intra-annular tautomerism,8 which appears to be the generalization of Kekulé’s oscillatory hypothesis. Ingold’s tautomerism hypothesis was later called mesomerism.9 The mesomery is an important concept in chemistry, which implicitly introduces the electron delocalization in the context of the prequantum electronic theory. The first applications of quantum chemistry to the benzene problem led on the molecular orbital (MO) side Erich Hückel to propose his famous 4n + 2 rule10 and on the valence bond (VB) side Pauling and Wheland to identify resonance with Ingold’s mesomerism.11,12 An enlighting contribution of modern VB theory on the benzene structure has been brought by Shaik et al.,13 who have shown that the hexagonal symmetry of benzene is due to the σ-system because the π component is distortive along a Kekulean distortion. * Authors to whom correspondence should be addressed (telephone +34-972-418912; fax +34-972-418356; e-mail miquel.sola@ udg.es or silvi@lct.jussieu.fr). ‡ Universitat de Girona. § Université Pierre et Marie Curie. 3911 Chem. Rev. 2005, 105, 3911−3947

Theoretical evaluation of electron delocalization in aromatic molecules by means of atoms in molecules (AIM) and electron localization function (ELF) topological approaches.

Recent Advances in Wave Function-Based Methods of Molecular-Property Calculations Trygve Helgaker,* Sonia Coriani, Poul Jørgensen, Kasper Kristensen, Jeppe Olsen, and Kenneth Ruud Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway Dipartimento di Scienze Chimiche e Farmaceutiche, Universit a degli Studi di Trieste, Via Giorgieri 1, I-34127 Trieste, Italy Lundbeck Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø, N-9037 Tromsø, Norway

Recent advances in wave function-based methods of molecular-property calculations.

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.

/pdf/the-aromatic-fluctuation-index-flu-a-new-aromaticity-index-33v1gz4z7j.pdf

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

Aromaticity cannot be measured directly by any physical or chemical experiment because it is not a well-defined magnitude. Its quantification is done indirectly from the measure of different properties that are usually found in aromatic compounds such as bond length equalisation, energetic stabilisation, and particular magnetic behaviour associated with induced ring currents. These properties have been used to set up the myriad of structural-, energetic-, and magnetic-based indices of aromaticity known to date. The cyclic delocalisation of mobile electrons in two or three dimensions is probably one of the key aspects that characterise aromatic compounds. However, it has not been until the last decade that electron delocalisation measures have been widely employed to quantify aromaticity. Some of these new indicators of aromaticity such as the PDI, FLU, ING, and INB were defined in our group. In this paper, we review the different existing descriptors of aromaticity that are based on electron delocalisation properties, we compare their performance with indices based on other properties, and we summarise a number of applications of electronic-based indices for the analysis of aromaticity in interesting chemical problems.

Quantifying aromaticity with electron delocalisation measures

Aromaticity is a central chemical concept widely used in modern chemistry for the interpretation of molecular structure, stability, reactivity, and magnetic properties of many compounds. As such, its reliable prediction is an important task of computational chemistry. In recent years, many methods to quantify aromaticity based on different physicochemical properties of molecules have been proposed. However, the nonobservable nature of aromaticity makes difficult to assess the performance of the numerous existing indices. In the present work, we introduce a series of fifteen aromaticity tests that can be used to analyze the advantages and drawbacks of a group of aromaticity descriptors. On the basis of the results obtained for a set of ten indicators of aromaticity, we conclude that indices based on the study of electron delocalization in aromatic species are the most accurate among those examined in this work.

/pdf/on-the-performance-of-some-aromaticity-indices-a-critical-12oytlpyf9.pdf

On the performance of some aromaticity indices: a critical assessment using a test set.

Electron sharing indexes (ESI) have been applied to numerous bonding situations to provide an insight into the nature of the molecular electronic structures. Some of the most popular ESI given in the literature, namely, the delocalization index (DI), defined in the context of the quantum theory of atoms in molecules (QTAIM), and the Fuzzy-Atom bond order (FBO), are here calculated at a correlated level for a wide set of molecules. Both approaches are based on the same quantity, the exchange–correlation density, to recover the electron sharing extent, and their differences lie in the definition of an atom in a molecule. In addition, while FBO atomic regions enable accurate and fast integrations, QTAIM definition of an atom leads to atomic domains that occasionally make the integration over these ones rather cumbersome. Besides, when working with a many-body wavefunction one can decide whether to calculate the ESI from first-order density matrices, or from second-order ones. The former way is usually preferred, since it avoids the calculation of the second-order density matrix, which is difficult to handle. Results from both definitions are discussed. Although these indexes are quite similar in their definition and give similar descriptions, when analyzed in greater detail, they reproduce different features of the bonding. In this manuscript DI is shown to explain certain bonding situations that FBO fails to cope with. Finally, these indexes are applied to the description of the aromaticity, through the aromatic fluctuation (FLU) and the para-DI (PDI) indexes. FLU and PDI indexes have been successfully applied using the DI measures, but other ESI based on other partitions such as Fuzzy-Atom can be used. The results provided in this manuscript for carbon skeleton molecules encourage the use of FBO for FLU and PDI indexes even at the correlated level.

/pdf/electron-sharing-indexes-at-the-correlated-level-application-235m3ilzb8.pdf

Electron sharing indexes at the correlated level. Application to aromaticity calculations.

In this paper the behavior of three well-known electron-sharing indexes, namely, the AIM delocalization index and the Mayer and fuzzy atom bond orders are studied at the Hartree-Fock level. A large number of five-membered ring molecules, containing several types of bonding, constitute the training set chosen for such purpose. A detailed analysis of the results obtained shows that the three indexes studied exhibit strong correlations, especially for homonuclear bonds. The correlation is somewhat poorer but still significant for polar bonds. In this case, the bond orders obtained with the Mayer and fuzzy atom approaches are normally closer to the formally predicted bond orders than those given by the AIM delocalization indexes, which are usually smaller than those expected from chemical intuition. In some particular cases, the use of diffuse functions in the calculation of Mayer bond orders leads to unrealistic results. In particular, noticeable trends are found for C-C bonds, encouraging the substitution of the delocalization index by the cheaper fuzzy atom or even the Mayer bond orders in the calculation of aromaticity indexes based on the delocalization index such as the para-delocalization index and the aromatic fluctuation index.

/pdf/comparison-of-the-aim-delocalization-index-and-the-mayer-and-1m6d8idz3o.pdf

Eduard Matito

Papers

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

Quantifying aromaticity with electron delocalisation measures

On the performance of some aromaticity indices: a critical assessment using a test set.

Electron sharing indexes at the correlated level. Application to aromaticity calculations.

Comparison of the AIM Delocalization Index and the Mayer and Fuzzy Atom Bond Orders