Charge density topological study of bonding in lithium clusters: Part I: Planar Li n clusters (n=4, 5, 6)

Prediction of Covalent Metal-Metal Bonding in Cp-M-M'-Nacnac Complexes of Group 2 and 12 Metals (Be, Mg, Ca, Zn, Cd, Hg)

Bimetallic CpMM'Nacnac molecules with group 2 and 12 metals (M=Be, Mg, Ca, Zn, Cd, Hg) that contain novel metal-metal bonding have been investigated in a theoretical study of their molecular and electronic structure, thermodynamic stability, and metal-metal bonding. In all cases the metal-metal bonds are characterized as electron-sharing covalent single bonds from natural bond orbital (NBO) and energy-decomposition analysis with natural orbitals of chemical valence (EDA-NOCV) analysis. The sum of [MM'] charges is relatively constant, with all complexes exhibiting a [MM']2+ core. Quantum theory of atoms in molecules (QTAIM) analysis indicates the presence of non-nuclear attractors (NNA) in the metal-metal bonds of the BeBe, MgMg, and CaCa complexes. There is substantial electron density (0.75-1.33 e) associated with the NNAs, which indicates that these metal-metal bonds, while classified as covalent electron-sharing bonds, retain significant metallic character that can be associated with reducing reactivity of the complex. The predicted stability of these complexes, combined with their novel covalent metal-metal bonding and potential as reducing agents, make them appealing targets for the synthesis of new metal-metal bonds.

Electronic structure of Li1,2,3 +,0,- and nature of the bonding in Li2,3 +,0,

The current study of the small lithium molecules Li2 +,0,- and Li3 +,0,- focuses on the nature of the bonding in these molecules as well as their structures and energetics (bond energies, ionization energies, and electron affinities). Valence CASSCF (2s,2p) calculations incorporate nondynamical electron correlation in the calculations, while the corresponding multireference configuration interaction and coupled cluster calculations incorporate dynamical electron correlation. Treatment of nondynamical correlation is critical for properly describing the Li2,3 +,0,- molecules as well as the Li- anion with dynamical correlation, in general, only fine-tuning the predictions. All lithium molecules and ions are bound, with the Li3 + and Li2 + ions being the most strongly bound, followed by Li3 - , Li2 , Li2 - and Li3 . The minimum energy structures of Li3 +,0,- are, respectively, an equilateral triangle, an isosceles triangle, and a linear structure. The results of SCGVB calculations are analyzed to obtain insights into the nature of the bonding in these molecules. An important finding of this work is that interstitial orbitals, a concept first put forward by McAdon and Goddard in 1985, play an essential role in the bonding of all lithium molecules considered here except for Li2 . The interstitial orbitals found in the Li3 +,0 molecules likely give rise to the non-nuclear attractors/maxima observed in these molecules.

On the Electride Nature of Na-hP4

Early quantum mechanical models suggested that pressure drives solids towards free-electron metal behavior where the ions are locked into simple close-packed structures. The prediction and subsequent discovery of high-pressure electrides (HPEs), compounds assuming open structures where the valence electrons are localized in interstitial voids, required a paradigm shift. Our quantum chemical calculations on the iconic insulating Na-hP4 HPE show that increasing density causes a 3s→3pd electronic transition due to Pauli repulsion between the 1s2s and 3s states, and orthogonality of the 3pd states to the core. The large lobes of the resulting Na-pd hybrid orbitals point towards the center of an 11-membered penta-capped trigonal prism and overlap constructively, forming multicentered bonds, which are responsible for the emergence of the interstitial charge localization in Na-hP4. These multicentered bonds facilitate the increased density of this phase, which is key for its stabilization under pressure.

Bonding in the high spin lithium clusters: Non-nuclear attractors play a crucial role

The bonding in lithium high-spin clusters contradicts the usual chemical bonding concept since there are no electron pairs between the atoms, and they are bound with parallel spin electrons. Quantum theory of atoms in molecules and interacting quantum atom analysis (IQA) were used to investigate the nature of bonding in the high-spin Li n n + 1 n = 2 - 5 clusters. Our findings demonstrate that the non-nuclear attractors (NNAs) are an essential component of the high-spin lithium clusters and play a key role in keeping them stable. Based on IQA energy terms, an electrostatic destabilizing interaction between the lithium atoms works against the cluster formation. On the other hand, the interactions between lithium atoms and NNA basins are stabilizing and outweigh the lithium-lithium destabilizing effects. In fact, NNAs tend to draw lithium atoms together and stabilize the resulting cluster. The high-spin clusters of lithium can be regarded as electrostatically driven compounds since the electrostatic components are primarily responsible for the stabilizing interactions between NNAs and Li atoms. The only exception is ³ Li₂ , which lacks NNA and has a non-repellent lithium-lithium interaction. Indeed, in the ³ Li₂ , the interatomic electrostatic component is negligibly small, and the exchange-correlation term leads to a weak bonding interaction.

Electrostatic Potential Topology for Probing Molecular Structure, Bonding and Reactivity

Following the pioneering investigations of Bader on the topology of molecular electron density, the topology analysis of its sister field viz. molecular electrostatic potential (MESP) was taken up by the authors’ groups. Through these studies, MESP topology emerged as a powerful tool for exploring molecular bonding and reactivity patterns. The MESP topology features are mapped in terms of its critical points (CPs), such as bond critical points (BCPs), while the minima identify electron-rich locations, such as lone pairs and π-bonds. The gradient paths of MESP vividly bring out the atoms-in-molecule picture of neutral molecules and anions. The MESP-based characterization of a molecule in terms of electron-rich and -deficient regions provides a robust prediction about its interaction with other molecules. This leads to a clear picture of molecular aggregation, hydrogen bonding, lone pair–π interactions, π-conjugation, aromaticity and reaction mechanisms. This review summarizes the contributions of the authors’ groups over the last three decades and those of the other active groups towards understanding chemical bonding, molecular recognition, and reactivity through topology analysis of MESP.

New insights in chemical reactivity from quantum chemical topology

Based on the quantum chemical topology of the modified electron localization function ELF_x , an efficient and robust mechanistic methodology designed to identify the favorable reaction pathway between two reactants is proposed. We first recall and reshape how the supermolecular interaction energy can be evaluated from only three distinct terms, namely the intermolecular coulomb energy, the intermolecular exchange-correlation energy and the intramolecular energies of reactants. Thereafter, we show that the reactivity between the reactants is driven by the first-order variation in the coulomb intermolecular energy defined in terms of the response to changes in the number of electrons. Illustrative examples with the formation of the dative bond B-N involved in the BH₃ NH₃ molecule and the typical formation of the hydrogen bond in the canonical water dimer are presented. For these selected systems, our approach unveils a noticeable mimicking of E_dual onto the DFT intermolecular interaction energy surface calculated between the both reactants. An automated reaction-path algorithm aimed to determine the most favorable relative orientations when the two molecules approach each other is also outlined.

Seven confluence principles: a case study of standardized statistical analysis for 26 methods that assign net atomic charges in molecules

This article studies two kinds of information extracted from statistical correlations between methods for assigning net atomic charges (NACs) in molecules. First, relative charge transfer magnitudes are quantified by performing instant least squares fitting (ILSF) on the NACs reported by Cho et al. (ChemPhysChem, 2020, 21, 688-696) across 26 methods applied to ∼2000 molecules. The Hirshfeld and Voronoi deformation density (VDD) methods had the smallest charge transfer magnitudes, while the quantum theory of atoms in molecules (QTAIM) method had the largest charge transfer magnitude. Methods optimized to reproduce the molecular dipole moment (e.g., ACP, ADCH, CM5) have smaller charge transfer magnitudes than methods optimized to reproduce the molecular electrostatic potential (e.g., CHELPG, HLY, MK, RESP). Several methods had charge transfer magnitudes even larger than the electrostatic potential fitting group. Second, confluence between different charge assignment methods is quantified to identify which charge assignment method produces the best NAC values for predicting via linear correlations the results of 20 charge assignment methods having a complete basis set limit across the dataset of ∼2000 molecules. The DDEC6 NACs were the best such predictor of the entire dataset. Seven confluence principles are introduced explaining why confluent quantitative descriptors offer predictive advantages for modeling a broad range of physical properties and target applications. These confluence principles can be applied in various fields of scientific inquiry. A theory is derived showing confluence is better revealed by standardized statistical analysis (e.g., principal components analysis of the correlation matrix and standardized reversible linear regression) than by unstandardized statistical analysis. These confluence principles were used together with other key principles and the scientific method to make assigning atom-in-material properties non-arbitrary. The N@C60 system provides an unambiguous and non-arbitrary falsifiable test of atomic population analysis methods. The HLY, ISA, MK, and RESP methods failed for this material.

Application of a quantum genetic algorithm and QTAIM analysis in the study of structural and electronic properties of neutral bimetallic clusters NaxLiy (4 ≤ x + y ≤ 10)

Alloy clusters of Na_xLi_y (4 ≤ x + y ≤ 10) are studied by exploring the potential energy surface in the ab initio MP2 level with the support of a quantum genetic algorithm (QGA). In some cases, the structures have been also refined with DFT and coupled-cluster methods. The general trends of sodium-lithium structures are in line with previous studies. The ionization potentials and polarizabilities to all structures were calculated with MP2 method and the average error between these two properties compared with experimental data was 6% and 13%, respectively. The topological analysis based on quantum theory of atoms in molecules (QTAIM) showed that by increasing the cluster size of the diatomic system there was a decrease of atomic interaction energies. The degree of degeneracy from D3BIA aromaticity index and the analysis of the atomic charges showed the influence (by charge transfer) of the chemical element in lower quantity in the cluster with respect to the other atoms. Our achievements of comparing our theoretical results with available experimental data have demonstrated that our approach can also predict satisfactorily quantum atomic and alloy clusters properties, at least, for low nuclearities.

Exploring the Structure, Energetic, and Magnetic Properties of Neutral Small Lithium Clusters Doped with Yttrium: Supermagnetic Atom Research

Density functional theory calculations based on magnetic and energetic stability criteria were performed to study a series of yttrium-doped lithium neutral clusters. A relativistic approximation was employed to properly describe the energy and multiplicity of the given clusters' fundamental states. The interaction of the 4d-Y atomic orbitals with the sp-Li states had an important role in the magnetic and energetic behavior of the selected systems. The spin density was concentrated over the yttrium atom regardless of the size of the cluster. Li7Y is a new stable superatom due to its enhanced magnetic properties.

Decay rate of real space delocalization measures: a comparison between analytical and test systems

We examine in this contribution the possible relation between the spatial decay rate of real space delocalization measures and the insulating- or metallic-like character of molecular and extended systems. We first show that in simple one-electron models, like the Hückel or tight binding approximations, delocalization indices (DIs) are intimately linked to the first-order reduced density matrix (1RDM), whose decay rate is known to be exponential in gapped systems and algebraic in gapless ones. DIs are shown to behave equivalently, with wild oscillations in gapless 1D, 2D and 3D models that do only persist in one-dimensional real cases, as computed at the Hartree-Fock or Kohn-Sham levels. Oscillations are shown to be directly related to Pauling resonant structures and chemical mesomerism. DIs in insulating-like moieties decay extremely fast. We propose that examining the decay of DIs along different directions in real materials may be used to detect facile and non-facile conductivity channels.

Topology vs. thermodynamics in chemical reactions: the instability of PH5

The topological approach, based on Bader theory, is compared to the common thermodynamical methodology to study chemical reactivity.

Source function and plane waves: Toward complete bader analysis

The source function (SF) is a topological descriptor that was introduced and developed by C. Gatti and R.W. Bader in 1998. The SF describes the contribution of each atom to the total electron density at a given point. To date, this descriptor has only been calculable from electron densities generated by all-electron (AE) methods for the investigation of single molecules or periodic systems. This study broadens the accessibility of the SF, offering its calculation from electron densities generated by plane wave (PW) methods. The new algorithm has been implemented in the open source code, CRITIC2. Our novel approach has been validated on a series of test systems, comparing the results obtained at PW level with those previously obtained through AE methods. © 2016 Wiley Periodicals, Inc.

On the existence and characterization of molecular electrides

Molecular electrides exist and can be designed using some guidelines and a set of computational criteria.

Rapid topography mapping of scalar fields: large molecular clusters

An efficient and rapid algorithm for topography mapping of scalar fields, molecular electron density (MED) and molecular electrostatic potential (MESP) is presented. The highlight of the work is the use of fast function evaluation by Deformed-atoms-in-molecules (DAM) method. The DAM method provides very rapid as well as sufficiently accurate function and gradient evaluation. For mapping the topography of large systems, the molecular tailoring approach (MTA) is invoked. This new code is tested out for mapping the MED and MESP critical points (CP's) of small systems. It is further applied to large molecular clusters viz. (H(2)O)(25), (C(6)H(6))(8) and also to a unit cell of valine crystal at MP2/6-31+G(d) level of theory. The completeness of the topography is checked by extensive search as well as applying the Poincaré-Hopf relation. The results obtained show that the DAM method in combination with MTA provides a rapid and efficient route for mapping the topography of large molecular systems.

Relativistic-consistent electron densities of the coinage metal clusters M2, M4, M4(2-), and M4Na2 (M = Cu, Ag, Au): a QTAIM study

We employ second-order Møller-Plesset perturbation theory level in combination with recently developed pseudopotential-based correlation consistent basis sets to obtain accurate relativistic-consistent electron densities for small coinage metal clusters. Using calculated electron densities, we employ Bader's quantum theory of atoms in molecules (QTAIM) to gain insights into the nature of metal-metal bonding in the clusters M(2), M(4), M(4)(2-), and M(4)Na(2) (M = Cu, Ag, Au). For the simplest case of the metal dimer, M(2), we correlate the strength of the metal-metal bond with the value of the electron density at the bond critical point, the total energy density at the bond critical point, the sharing (delocalization) index, and the values of the two principle negative curvatures. We then consider changes to the metal-metal bonding and charge density distribution upon the addition of two metal atoms to form the metal tetramer, M(4), and then followed by the addition of an electron pair to form M(4)(2-) and finally followed by the addition of two alkali metal (sodium) ions to form M(4)Na(2). Using topological properties of the electron density, we present evidence for the existence of σ-aromaticity in Au(4)(2-). We also report the existence of two non-nuclear attractors in the molecular graph of Cu(4)(2-) and large negative charge accumulation in the nonbonded Cu basins of this cluster.

A pseudoatom in a cage: trimetallofullerene Y(3)@C(80) mimics y(3)n@c(80) with nitrogen substituted by a pseudoatom

Y(3)C(80) obtained in the synthesis of nitride clusterfullerenes Y(3)N@C(2n) (2n = 80-88) by the reactive atmosphere method is found to be a genuine trimetallofullerene, Y(3)@C(80), with low ionization potential and divalent state of yttrium atoms. DFT studies of the electronic structure of Y(3)@C(80) show that this molecule mimics Y(3)N@C(80) with the pseudoatom (PA) instead of the nitrogen atom. Topology analysis of the electron density and electron localization function show that yttrium atoms form Y-PA bonds rather than direct Y-Y bonds. Molecular dynamics simulations show that the Y(3)PA cluster is as rigid as Y(3)N and rotates inside the fullerene cage as a single entity.

Electron localizability indicators ELI-D and ELIA for highly correlated wavefunctions of homonuclear dimers. I. Li2, Be2, B2, and C2

Electron localizability indicators based on the parallel-spin electron pair density (ELI-D) and the antiparallel-spin electron pair density (ELIA) are studied for the correlated ground-state wavefunctions of Li(2), Be(2), B(2), and C(2) diatomic molecules. Different basis sets and reference spaces are used for the multireference configuration interaction method following the complete active space calculations to investigate the local effect of electron correlation on the extent of electron localizability in position space determined by the two functionals. The results are complemented by calculations of effective bond order, vibrational frequency, and Laplacian of the electron density at the bond midpoint. It turns out that for Li(2), B(2), and C(2) the reliable topology of ELI-D is obtained only at the correlated level of theory.