Topological Characterization of the Electron Density Laplacian in Crystals. The Case of the Group IV Elements

Electron Density Learning of Z-Bonds in Ionic Liquids and Its Application

Ionic liquids (ILs) exhibit fascinating properties due to special Z-bonds and have been widely used in electrochemical systems. The local Z-bond networks potentially cause a discrepancy in electrochemical properties. Understanding the correlations between the Z-bond energy (EZ-bond) and the electrochemical properties is helpful to identify appropriate ILs. It is difficult to estimate the correlations from single density functional theory calculations or molecular dynamic simulations. In this work, a machine learning model targeting the electronic density (ρBCP) of Z-bonds has been trained successfully, as expected for use in systems above the nanoscale size. The connection between the EZ-bond and the electrochemical potential window in ILs@TiO2, as well as that between the EZ-bond and the charge carrier mobility in ILs-PEDOT:Tos@SiO2, was separately investigated. This study highlights an efficient model for predicting ρBCP in nanoscale systems and anticipates exploring the connection between Z-bonds and the electrochemical properties of IL-based systems.

TopIso3D Viewer: Enhancing Topological Analysis through 3D Isosurfaces

We present TopIso3D Viewer, a software with a user-friendly graphical interface, that generates three-dimensional maps to analyze descriptors based on the Quantum Theory of Atoms in Molecules (QTAIM), applied in periodic and nonperiodic systems. The software also automates the launching of topological analysis calculations through the Topond package and generates a report that facilitates the identification of the values of the calculated descriptors, in the Bond Critical Points (BCP) and Critical Points of the Laplacian of the electron density (LCP), facilitating the classification of chemical interactions. The map projects created can be stored in the form of HTML files, for later consultation through any type of browser. For validation of the software, several systems with 0-3D dimensions were studied. In addition, the topology of urea molecular crystal and its isolated molecule were revisited.

Finding critical points and reconstruction of electron densities on grids

The quantum theory of atoms in molecules (QTAIM), developed by Bader and co-workers, is one of the most popular ways of extracting chemical insight from the results of quantum mechanical calculations. One of the basic tasks in QTAIM is to locate the critical points of the electron density and calculate various quantities (density, Laplacian, etc.) on them since these have been found to correlate with molecular properties of interest. If the electron density is given analytically, this process is relatively straightforward. However, locating the critical points is more challenging if the density is known only on a three-dimensional uniform grid. A density grid is common in periodic solids because it is the natural expression for the electron density in plane-wave calculations. In this article, we explore the reconstruction of the electron density from a grid and its use in critical point localization. The proposed reconstruction method employs polyharmonic spline interpolation combined with a smoothing function based on the promolecular density. The critical point search based on this reconstruction is accurate, trivially parallelizable, works for periodic and non-periodic systems, does not present directional lattice bias when the grid is non-orthogonal, and locates all critical points of the underlying electron density in all tests studied. The proposed method also provides an accurate reconstruction of the electron density over the space spanned by the grid, which may be useful in other contexts besides critical point localization.

Structural insights into the anti-cancer activity of quercetin on G-tetrad, mixed G-tetrad, and G-quadruplex DNA using quantum chemical and molecular dynamics simulations

Human telomerase referred as 'terminal transferase' is a nucleoprotein enzyme which inhibits the disintegration of telomere length and act as a drug target for the anticancer therapy. The tandem repeating structure of telomere sequence forms the guanine-rich quadruplex structures that stabilize stacked tetrads. In our present work, we have investigated the interaction of quercetin with DNA tetrads using DFT. Geometrical analysis revealed that the influence of quercetin drug induces the structural changes into the DNA tetrads. Among DNA tetrads, the quercetin stacked with GCGC tetrad has the highest interaction energy of -88.08 kcal/mol. The binding mode and the structural stability are verified by the absorption spectroscopy method. The longer wavelength was found at 380 nm and it exhibits bathochromic shift. The findings help us to understand the binding nature of quercetin drug with DNA tetrads and it also inhibits the telomerase activity. Further, the quercetin drug interacted with G-quadruplex DNA by using molecular dynamics (MD) simulation studies for 100 ns simulation at different temperatures and different pH levels (T = 298 K, 320 K and pH = 7.4, 5.4). The structural stability of the quercetin with G-quadruplex structure is confirmed by RMSD. For the acidic condition (pH = 5.4), the binding affinity is higher toward G-quadruplex DNA, this result resembles that the quercetin drug is well interacted with G-quadruplex DNA at acidic condition (pH = 7.4) than the neutral condition. The obtained results show that quercetin drug stabilizes the G-quadruplex DNA, which regulates telomerase enzyme and it potentially acts as a novel anti-cancer agent.Communicated by Ramaswamy H. Sarma.

Copper(i) complexes of functionalized sulfur-containing ligands: structural and theoretical insights into chalcogen bonding

Four new copper(i) thiocyanate complexes were studied using geometrical parameters and the lump–hole approach for justification of the strength and nature of chalcogen bonding.

Curly arrows, electron flow, and reaction mechanisms from the perspective of the bonding evolution theory

Despite the usefulness of curly arrows in chemistry, their relationship with real electron density flows is still imprecise, and even their direct connection to quantum chemistry is still controversial. The paradigmatic description - from first principles - of the mechanistic aspects of a given chemical process is based mainly on the relative energies and geometrical changes at the stationary points of the potential energy surface along the reaction pathway; however, it is not sufficient to describe chemical systems in terms of bonding aspects. Probing the electron density distribution during a chemical reaction can provide important insights, enabling us to understand and control chemical reactions. This aim has required an extension of the relationships between the concepts of traditional chemistry and those of quantum mechanics. Bonding evolution theory (BET), which combines the topological analysis of the electron localization function (ELF) and Thom's catastrophe theory (CT), provides a powerful method that offers insight into the molecular mechanism of chemical rearrangements. In agreement with the laws of physical and aspects of quantum theory, BET can be considered an appropriate tool to tackle chemical reactivity with a wide range of possible applications. In this work, BET is applied to address a long-standing problem: the ability to monitor the flow of electron density. BET analysis shows a connection between quantum mechanics and bond making/forming processes. Likewise, the present approach retrieves the classical curly arrows used to describe the rearrangements of chemical bonds and provides detailed physical grounds for this type of representation. We demonstrate this procedure using the test set of prototypical examples of thermal ring apertures, and the degenerated Cope rearrangement of semibullvalene.

A topological assessment of the electronic structure of mesoionic compounds

Mesoionic compounds belonging to the 1,3-oxazol-5-one, 1,3-diazole-4-thione and 1,3-thiazole-5-thione rings have been evaluated by a combination of Density Functional Theory, Quantum Theory of Atoms in Molecules, Electron Localization Function, Natural Bond Orbitals and Geodesic Electrostatic Potential Charge calculations. Atomic, bond, and ring properties have been considered to describe the electronic structure of mesoionic compounds. The results show that not only the ring type, but also the substituent groups have great influence on these properties. In addition, there is a significant and heterogeneous π-bonding contribution throughout the mesoionic rings. Finally, we conclude that some classical conceptions of charge localization and π-bonding contribution in these compounds are misleading or incomplete. © 2015 Wiley Periodicals, Inc.

Spatial dispersion of lone electron pairs? – Experimental charge density of cubic arsenic(iii) oxide

Lone electron pair dispersion into three separate domains in space is reported and discussed for the first time.

Density functional studies on photophysical properties and chemical reactivities of the triarylboranes: effect of the constraint of planarity

The geometric and electronic structures, absorption spectra, transporting properties, chemical reactivity indices and electrostatic potentials of the planar three-coordinate organoboron compounds 1-2 and twisted reference compound Mes(3)B, have been investigated by employing density functional theory (DFT) and conceptual DFT methods to shed light on the planarity effects on the photophysical properties and the chemical reactivity. The results show that the planar compounds 1-2 exhibit significantly lower HOMO level than Mes(3)B, owing to the stronger electronic induction effect of boron centers. This feature conspicuously induces a blue shifted absorption for 1, although 1 seemingly possesses more extended conjugation framework than Mes(3)B. Importantly, the reactivity strength of the boron atoms in 1-2 is much lower than that in Mes(3)B, despite the fact that the tri-coordinate boron centers of 1-2 are completely naked. The interesting and abnormal phenomenon is caused by the strong p-π electronic interactions, that is, the empty p-orbital of boron center is partly filled by π-electron of the neighbor carbon atoms in 1-2, which are confirmed by the analysis of Laplacian of the electron density and natural bond orbitals. Furthermore, the negative electrostatic potentials of the boron centers in 1-2 also interpret that they are not the most preferred sites for incoming nucleophiles. Moreover, it is also found that the planar compounds 1-2 can act as promising electron transporting materials since the internal reorganization energies for electron are really small.

The nature of the interaction of dimethylselenide with IIIA group element compounds

The first systematic theoretical study of the nature of intermolecular bonding of dimethylselenide as donor and IIIA group element halides as acceptors was made with the help of the approach of Quantum Theory of Atoms in Molecules. Density Functional Theory with "old" Sapporo triple-ζ basis sets was used to calculate geometry, thermodynamics, and wave function of Me2Se···AX3 complexes. The analysis of the electron density distribution and the Laplacian of the electron density allowed us to reveal and explain the tendencies in the influence of the central atom (A = B, Al, Ga, In) and halogen (X = F, Cl, Br, I) on the nature of Se···A bonding. Significant changes in properties of the selenium lone pair upon complexation were described by means of the analysis of the Laplacian of the charge density. Charge transfer characteristics and the contributions to it from electron localization and delocalization were analyzed in terms of localization and delocalization indexes. Common features of the complexation and differences in the nature of bonding were revealed. Performed analysis evidenced that gallium and indium halide complexes can be attributed to charge transfer-driven complexes; aluminum halides complexes seem to be mainly of an electrostatic nature. The nature of bonding in different boron halides essentially varies; these complexes are stabilized mainly by covalent Se···B interaction. In all the complexes under study covalence of the Se···A interaction is rather high.

Revealing non-covalent interactions in solids: NCI plots revisited

In this article, the NCI method [Johnson et al., J. Am. Chem. Soc., 2010, 132, 6498] for plotting and analysing non-covalent interactions (NCI) is extended to periodic (solid-state) electron densities and implemented in the critic program. The new code uses self-consistent electron densities from a variety of electronic structure methods (pseudopotentials/plane-wave, FP-LAPW, local orbitals, etc.), and it can also build the promolecular density from the crystal geometry alone. As an example of the new code, intermolecular interactions in several molecular crystals are presented and analyzed. The connection with QTAIM studies is established and a reinterpretation of the NCI domains is given regarding the current knowledge of the field. The connection between NCI domains and intermolecular vibrations is made apparent, as well as the ability of the method to reveal the locality of bonding.

Topological partition of the elastic constants of crystals

We present a partitioning of the elastic constants of a crystal into atomic contributions by using the atomic basin concept inherent to Bader's Quantum Theory of Atoms in Molecules. The partition is made by following the evolution of the cell volume and the atomic basin volumes under appropriately defined cell deformations. The method is carefully examined, including internal consistency checks. The transferability of atomic contributions between different crystals is determined by obtaining and comparing the oxygen contribution to the elastic constants of a selection of cubic oxides that includes the rock-salt, perovskite, antifluorite, and cuprite crystal families.