1
|
Ćeranić K, Milovanović B, Petković M. Density functional theory study of crown ether-magnesium complexes: from a solvated ion to an ion trap. Phys Chem Chem Phys 2023; 25:32656-32665. [PMID: 38010878 DOI: 10.1039/d3cp03991a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2023]
Abstract
Metal ion detection rests on host-guest recognition. We propose a theoretical protocol for designing an optimal trap for a desired metal cation. A host for magnesium ions was sought for among derivatives of crown ethers 12-crown-4, 15-crown-5, and 18-crown-6. Mg-crown complexes and their hydrated counterparts with water molecules bound to the cation were optimized using density functional theory. Based on specific geometric criteria, Interacting quantum atoms analysis and density functional theory-based molecular dynamics of Mg-crown complexes immersed in water, crown ethers for optimal accommodation of Mg2+ in aqueous solution were identified. Selectivity of the chosen crowns towards Na+, K+, and Ca2+ ions is addressed.
Collapse
Affiliation(s)
- Katarina Ćeranić
- Innovative Centre of the Faculty of Chemistry, Studentski trg 12-16, 11158 Belgrade, Serbia
- University of Belgrade - Faculty of Physical Chemistry, Studentski trg 12-16, 11158 Belgrade, Serbia.
| | - Branislav Milovanović
- University of Belgrade - Faculty of Physical Chemistry, Studentski trg 12-16, 11158 Belgrade, Serbia.
| | - Milena Petković
- University of Belgrade - Faculty of Physical Chemistry, Studentski trg 12-16, 11158 Belgrade, Serbia.
| |
Collapse
|
2
|
Jabłoński M. Bader's Topological Bond Path Does Not Necessarily Indicate Stabilizing Interaction-Proof Studies Based on the Ng@[3 n]cyclophane Endohedral Complexes. Molecules 2023; 28:6353. [PMID: 37687183 PMCID: PMC10490063 DOI: 10.3390/molecules28176353] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2023] [Revised: 08/25/2023] [Accepted: 08/28/2023] [Indexed: 09/10/2023] Open
Abstract
According to Bader's quantum theory of atoms in molecules (QTAIM), the simultaneous presence of a bond path and the corresponding bond critical point between any two atoms is both a necessary and sufficient condition for the atoms to be bonded to one another. In principle, this means that this pair of atoms should make a stabilizing contribution to the molecular system. However, the multitude of so-called counterintuitive bond paths strongly suggests that this statement is not necessarily true. Particularly 'troublesome' are endohedral complexes, in which encapsulation-enforced proximity between the trapped guest (e.g., an atom) and the host's cage system usually 'produces' many counterintuitive bond paths. In the author's opinion, the best evidence to demonstrate the repulsive nature of the intra-cage guest⋯host interaction is the use of some trapping systems containing small escape channels and then showing that the initially trapped entity spontaneously escapes outside the host's cage during geometry optimization of the initially built guest@host endohedral complex. For this purpose, a group of 24 Ng@[3n]cyclophane (3≤n≤6) endohedral complexes is used. As a result, arguments are presented showing that Bader's topological bond path does not necessarily indicate a stabilizing interaction.
Collapse
Affiliation(s)
- Mirosław Jabłoński
- Faculty of Chemistry, Nicolaus Copernicus University, Gagarina 7, 87-100 Toruń, Poland
| |
Collapse
|
3
|
The Ultrashort Spike-Ring Interaction in Substituted Iron Maiden Molecules. Molecules 2023; 28:molecules28052244. [PMID: 36903489 PMCID: PMC10004599 DOI: 10.3390/molecules28052244] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2023] [Revised: 02/23/2023] [Accepted: 02/24/2023] [Indexed: 03/06/2023] Open
Abstract
The in forms of molecular iron maidens are known for their unique ultrashort interaction between the apical hydrogen atom or its small substituent and the surface of the benzene ring. It is generally believed that this forced ultrashort X⋯π contact is associated with high steric hindrance, which is responsible for specific properties of iron maiden molecules. The main aim of this article is to investigate the influence of significant charge enrichment or depletion of the benzene ring on the characteristics of the ultrashort C-X⋯π contact in iron maiden molecules. For this purpose, three strongly electron-donating (-NH2) or strongly electron-withdrawing (-CN) groups were inserted into the benzene ring of in-[34,10][7]metacyclophane and its halogenated (X = F, Cl, Br) derivatives. It is shown that, despite such extremely electron-donating or electron-accepting properties, the considered iron maiden molecules surprisingly reveal quite high resistance to changes in electronic properties.
Collapse
|
4
|
Jabłoński M. The physical nature of the ultrashort spike-ring interaction in iron maiden molecules. J Comput Chem 2022; 43:1206-1220. [PMID: 35593685 DOI: 10.1002/jcc.26879] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2022] [Revised: 04/11/2022] [Accepted: 04/13/2022] [Indexed: 11/08/2022]
Abstract
The so-called 'iron maiden' molecules belong to one of the most interesting subgroups of cyclophanes due to the presence of the ultrashort interaction between the CX apical bond and the benzene ring. This article presents an in-depth theoretical study of 16 'iron maiden' molecules, in which X = H, F, Cl or Br and the side chains are of various lengths and types: CSC, CSCC, CCC, and CCCC. It is shown that the H → F → Cl → Br substitution leads to a significant expansion of the 'iron maiden' molecule. Shorter chains lead to more pronounced effects, while insertion of sulfur atoms into the side chains lowers them. Structural changes are associated with an increase in energetic destabilization of X. Moreover, unlike for H, in the case of X = halogen, the out → in isomerization is energetically disadvantageous. The 'iron maiden' molecules are characterized by the presence of only three X⋯CAr bond paths. Particularly noteworthy are unusually large (even up to 32) values of the X⋯CAr bond ellipticity, which results from flat electron density distribution. The X⋯π interaction in each of the investigated 'iron maiden' molecule turned out to be multi-center, stabilizing and almost purely covalent in nature as indicated by the definitely dominant percentage (94.8%-101.6%) of the exchange-correlation energy. The spatial hindrance within the 'iron maiden' molecules appears to be not so much due to the X⋯π repulsion, but due to unfavorable steric interactions between X and the CC side bonds. It is also confirmed that some CH⋯HC interactions in aliphatic chains can be very weakly stabilizing.
Collapse
Affiliation(s)
- Mirosław Jabłoński
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, Toruń, Poland
| |
Collapse
|
5
|
Fast Quantum Approach for Evaluating the Energy of Non-Covalent Interactions in Molecular Crystals: The Case Study of Intermolecular H-Bonds in Crystalline Peroxosolvates. Molecules 2022; 27:molecules27134082. [PMID: 35807323 PMCID: PMC9268483 DOI: 10.3390/molecules27134082] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2022] [Revised: 06/20/2022] [Accepted: 06/21/2022] [Indexed: 11/17/2022] Open
Abstract
Energy/enthalpy of intermolecular hydrogen bonds (H-bonds) in crystals have been calculated in many papers. Most of the theoretical works used non-periodic models. Their applicability for describing intermolecular H-bonds in solids is not obvious since the crystal environment can strongly change H-bond geometry and energy in comparison with non-periodic models. Periodic DFT computations provide a reasonable description of a number of relevant properties of molecular crystals. However, these methods are quite cumbersome and time-consuming compared to non-periodic calculations. Here, we present a fast quantum approach for estimating the energy/enthalpy of intermolecular H-bonds in crystals. It has been tested on a family of crystalline peroxosolvates in which the H∙∙∙O bond set fills evenly (i.e., without significant gaps) the range of H∙∙∙O distances from ~1.5 to ~2.1 Å typical for strong, moderate, and weak H-bonds. Four of these two-component crystals (peroxosolvates of macrocyclic ethers and creatine) were obtained and structurally characterized for the first time. A critical comparison of the approaches for estimating the energy of intermolecular H-bonds in organic crystals is carried out, and various sources of errors are clarified.
Collapse
|
6
|
Determining Repulsion in Cyclophane Cages. Molecules 2022; 27:molecules27133969. [PMID: 35807214 PMCID: PMC9268502 DOI: 10.3390/molecules27133969] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2022] [Revised: 06/15/2022] [Accepted: 06/18/2022] [Indexed: 02/01/2023] Open
Abstract
Superphane, i.e., [2.2.2.2.2.2](1,2,3,4,5,6)cyclophane, is a very convenient molecule in studying the nature of guest⋯host interactions in endohedral complexes. Nevertheless, the presence of as many as six ethylene bridges in the superphane molecule makes it practically impossible for the trapped entity to escape out of the superphane cage. Thus, in this article, I have implemented the idea of using the superphane derivatives with a reduced number of ethylene linkers, which leads to the [2n] cyclophanes where n<6. Seven such cyclophanes are then allowed to form endohedral complexes with noble gas (Ng) atoms (He, Ne, Ar, Kr). It is shown that in the vast majority of cases, the initially trapped Ng atom spontaneously escapes from the cyclophane cage, creating an exohedral complex. This is the best proof that the Ng⋯cyclophane interaction in endohedral complexes is indeed highly repulsive, i.e., destabilizing. Apart from the ‘sealed’ superphane molecule, endohedral complexes are only formed in the case of the smallest He atom. However, it has been shown that in these cases, the Ng⋯cyclophane interaction inside the cyclophane cage is nonbonding, i.e., repulsive. This highly energetically unfavorable effect causes the cyclophane molecule to ‘swell’.
Collapse
|
7
|
Does the Presence of a Bond Path Really Mean Interatomic Stabilization? The Case of the Ng@Superphane (Ng = He, Ne, Ar, and Kr) Endohedral Complexes. Symmetry (Basel) 2021. [DOI: 10.3390/sym13122241] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
Using a fairly structurally flexible and, therefore, very suitable for this type of research, superphane molecule, we demonstrate that the inclusion of a noble gas atom (Ng = He, Ne, Ar, and Kr) inside it and, thus, the formation of the Ng@superphane endohedral complex, leads to its ‘swelling’. Positive values of both the binding and strain energies prove that encapsulation and in turn ‘swelling’ of the superphane molecule is energetically unfavorable and that the Ng⋯C interactions in the interior of the cage are destabilizing, i.e., repulsive. Additionally, negative Mayer Bond Orders indicate the antibonding nature of Ng⋯C contacts. This result in combination with the observed Ng⋯C bond paths shows that the presence of a bond path in the molecular graph does not necessarily prove interatomic stabilization. It is shown that the obtained conclusions do not depend on the computational methodology, i.e., the method and the basis set used. However, on the contrary, the number of bond paths may depend on the methodology. This is yet another disadvantageous finding that does not favor the treatment of bond paths on molecular graphs as indicators of chemical bonds. The Kr@superphane endohedral complex features one of the longest C–C bonds ever reported (1.753 Å).
Collapse
|
8
|
Cukrowski I. A unified molecular‐wide and electron density based concept of chemical bonding. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1579] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Ignacy Cukrowski
- Department of Chemistry, Faculty of Natural and Agricultural Sciences University of Pretoria Pretoria South Africa
| |
Collapse
|
9
|
Stationary states of systems with intermolecular interactions dominated by electrostatics: Structure of trimethylammonium and tetramethylammonium chlorides and bromides in the gas phase, monomers and dimers. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2021.138809] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
10
|
de Sousa DWO, Nascimento MAC. Substituent Effects on the Quantum Interference of Two-Center One-Electron Bonds: [B 2X 6] - (X = H, F, Cl, CN, OH, CH 3, and OCH 3). J Phys Chem A 2021; 125:4558-4564. [PMID: 34014679 DOI: 10.1021/acs.jpca.1c02771] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The interference energy analysis (IEA) provided by the generalized product function energy partitioning (GPF-EP) method was applied to investigate the influence of the neighboring atoms on the nature of the two-center one-electron (2c1e) bonds in the anion dimers of BX3 species (X = H, F, Cl, CN, OH, CH3, and OCH3). The species were studied at the GVB-PP(6/12).SC(1,2)/6-31**G++ level of calculation. The IEA has revealed that there is a balance between two main factors determining the chemical stability of the species. Quantum interference acts as the sole stabilizing effect in the formation of the chemical bonds, particularly as the result of the drop in kinetic energy, and the electronegativity of the substituent has a direct influence on the magnitude of this effect. The quasi-classical energy is responsible for the destabilizing factors, mainly the group bulkiness, and the "electron-withdrawing" effect in the case of the cyano group.
Collapse
|
11
|
Molecular Tailoring Approach for the Estimation of Intramolecular Hydrogen Bond Energy. Molecules 2021; 26:molecules26102928. [PMID: 34069140 PMCID: PMC8155843 DOI: 10.3390/molecules26102928] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2021] [Revised: 05/11/2021] [Accepted: 05/12/2021] [Indexed: 11/29/2022] Open
Abstract
Hydrogen bonds (HBs) play a crucial role in many physicochemical and biological processes. Theoretical methods can reliably estimate the intermolecular HB energies. However, the methods for the quantification of intramolecular HB (IHB) energy available in the literature are mostly empirical or indirect and limited only to evaluating the energy of a single HB. During the past decade, the authors have developed a direct procedure for the IHB energy estimation based on the molecular tailoring approach (MTA), a fragmentation method. This MTA-based method can yield a reliable estimate of individual IHB energy in a system containing multiple H-bonds. After explaining and illustrating the methodology of MTA, we present its use for the IHB energy estimation in molecules and clusters. We also discuss the use of this method by other researchers as a standard, state-of-the-art method for estimating IHB energy as well as those of other noncovalent interactions.
Collapse
|
12
|
Two Faces of Water in the Formation and Stabilization of Multicomponent Crystals of Zwitterionic Drug-Like Compounds. Symmetry (Basel) 2021. [DOI: 10.3390/sym13030425] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
Two new hydrated multicomponent crystals of zwitterionic 2-aminonicotinic acid with maleic and fumaric acids have been obtained and thoroughly characterized by a variety of experimental (X-ray analysis and terahertz Raman spectroscopy) and theoretical periodic density functional theory calculations, followed by Bader analysis of the crystalline electron density) techniques. It has been found that the Raman-active band in the region of 300 cm−1 is due to the vibrations of the intramolecular O-H...O bond in the maleate anion. The energy/enthalpy of the intermolecular hydrogen bonds was estimated by several empirical approaches. An analysis of the interaction networks reflects the structure-directing role of the water molecule in the examined multicomponent crystals. A general scheme has been proposed to explain the proton transfer between the components during the formation of multicomponent crystals in water. Water molecules were found to play the key role in this process, forming a “water wire” between the COOH group of the dicarboxylic acid and the COO– group of the zwitterion and the rendering crystal lattice of the considered multicomponent crystals.
Collapse
|
13
|
Medvedev AG, Churakov AV, Prikhodchenko PV, Lev O, Vener MV. Crystalline Peroxosolvates: Nature of the Coformer, Hydrogen-Bonded Networks and Clusters, Intermolecular Interactions. Molecules 2020; 26:E26. [PMID: 33374602 PMCID: PMC7793138 DOI: 10.3390/molecules26010026] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2020] [Revised: 12/19/2020] [Accepted: 12/19/2020] [Indexed: 01/22/2023] Open
Abstract
Despite the technological importance of urea perhydrate (percarbamide) and sodium percarbonate, and the growing technological attention to solid forms of peroxide, fewer than 45 peroxosolvates were known by 2000. However, recent advances in X-ray diffractometers more than tripled the number of structurally characterized peroxosolvates over the last 20 years, and even more so, allowed energetic interpretation and gleaning deeper insight into peroxosolvate stability. To date, 134 crystalline peroxosolvates have been structurally resolved providing sufficient insight to justify a first review article on the subject. In the first chapter of the review, a comprehensive analysis of the structural databases is carried out revealing the nature of the co-former in crystalline peroxosolvates. In the majority of cases, the coformers can be classified into three groups: (1) salts of inorganic and carboxylic acids; (2) amino acids, peptides, and related zwitterions; and (3) molecular compounds with a lone electron pair on nitrogen and/or oxygen atoms. The second chapter of the review is devoted to H-bonding in peroxosolvates. The database search and energy statistics revealed the importance of intermolecular hydrogen bonds (H-bonds) which play a structure-directing role in the considered crystals. H2O2 always forms two H-bonds as a proton donor, the energy of which is higher than the energy of analogous H-bonds existing in isostructural crystalline hydrates. This phenomenon is due to the higher acidity of H2O2 compared to water and the conformational mobility of H2O2. The dihedral angle H-O-O-H varies from 20 to 180° in crystalline peroxosolvates. As a result, infinite H-bonded 1D chain clusters are formed, consisting of H2O2 molecules, H2O2 and water molecules, and H2O2 and halogen anions. H2O2 can form up to four H-bonds as a proton acceptor. The third chapter of the review is devoted to energetic computations and in particular density functional theory with periodic boundary conditions. The approaches are considered in detail, allowing one to obtain the H-bond energies in crystals. DFT computations provide deeper insight into the stability of peroxosolvates and explain why percarbamide and sodium percarbonate are stable to H2O2/H2O isomorphic transformations. The review ends with a description of the main modern trends in the synthesis of crystalline peroxosolvates, in particular, the production of peroxosolvates of high-energy compounds and mixed pharmaceutical forms with antiseptic and analgesic effects.
Collapse
Affiliation(s)
- Alexander G. Medvedev
- Kurnakov Institute of General and Inorganic Chemistry, Russian Academy of Sciences, Leninskii Prosp. 31, 119991 Moscow, Russia; (A.G.M.); (A.V.C.); (P.V.P.)
| | - Andrei V. Churakov
- Kurnakov Institute of General and Inorganic Chemistry, Russian Academy of Sciences, Leninskii Prosp. 31, 119991 Moscow, Russia; (A.G.M.); (A.V.C.); (P.V.P.)
| | - Petr V. Prikhodchenko
- Kurnakov Institute of General and Inorganic Chemistry, Russian Academy of Sciences, Leninskii Prosp. 31, 119991 Moscow, Russia; (A.G.M.); (A.V.C.); (P.V.P.)
| | - Ovadia Lev
- The Casali Center of Applied Chemistry, The Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Mikhail V. Vener
- Kurnakov Institute of General and Inorganic Chemistry, Russian Academy of Sciences, Leninskii Prosp. 31, 119991 Moscow, Russia; (A.G.M.); (A.V.C.); (P.V.P.)
- Department of Quantum Chemistry, Mendeleev University of Chemical Technology, Miusskaya Square 9, 125047 Moscow, Russia
| |
Collapse
|
14
|
Jabłoński M. A Critical Overview of Current Theoretical Methods of Estimating the Energy of Intramolecular Interactions. Molecules 2020; 25:molecules25235512. [PMID: 33255559 PMCID: PMC7728086 DOI: 10.3390/molecules25235512] [Citation(s) in RCA: 44] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2020] [Revised: 11/19/2020] [Accepted: 11/23/2020] [Indexed: 11/16/2022] Open
Abstract
This article is probably the first such comprehensive review of theoretical methods for estimating the energy of intramolecular hydrogen bonds or other interactions that are frequently the subject of scientific research. Rather than on a plethora of numerical data, the main focus is on discussing the theoretical rationale of each method. Additionally, attention is paid to the fact that it is very often possible to use several variants of a particular method. Both of the methods themselves and their variants often give wide ranges of the obtained estimates. Attention is drawn to the fact that the applicability of a particular method may be significantly limited by various factors that disturb the reliability of the estimation, such as considerable structural changes or new important interactions in the reference system.
Collapse
Affiliation(s)
- Mirosław Jabłoński
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, 87-100 Toruń, Poland
| |
Collapse
|
15
|
|
16
|
Guevara-Vela JM, Francisco E, Rocha-Rinza T, Martín Pendás Á. Interacting Quantum Atoms-A Review. Molecules 2020; 25:E4028. [PMID: 32899346 PMCID: PMC7504790 DOI: 10.3390/molecules25174028] [Citation(s) in RCA: 55] [Impact Index Per Article: 13.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2020] [Revised: 08/24/2020] [Accepted: 08/26/2020] [Indexed: 12/16/2022] Open
Abstract
The aim of this review is threefold. On the one hand, we intend it to serve as a gentle introduction to the Interacting Quantum Atoms (IQA) methodology for those unfamiliar with it. Second, we expect it to act as an up-to-date reference of recent developments related to IQA. Finally, we want it to highlight a non-exhaustive, yet representative set of showcase examples about how to use IQA to shed light in different chemical problems. To accomplish this, we start by providing a brief context to justify the development of IQA as a real space alternative to other existent energy partition schemes of the non-relativistic energy of molecules. We then introduce a self-contained algebraic derivation of the methodological IQA ecosystem as well as an overview of how these formulations vary with the level of theory employed to obtain the molecular wavefunction upon which the IQA procedure relies. Finally, we review the several applications of IQA as examined by different research groups worldwide to investigate a wide variety of chemical problems.
Collapse
Affiliation(s)
- José Manuel Guevara-Vela
- Institute of Chemistry, National Autonomous University of Mexico, Circuito Exterior, Ciudad Universitaria, Delegación Coyoacán C.P., Mexico City 04510, Mexico; (J.M.G.-V.); (T.R.-R.)
| | - Evelio Francisco
- Department of Analytical and Physical Chemistry, University of Oviedo, E-33006 Oviedo, Spain;
| | - Tomás Rocha-Rinza
- Institute of Chemistry, National Autonomous University of Mexico, Circuito Exterior, Ciudad Universitaria, Delegación Coyoacán C.P., Mexico City 04510, Mexico; (J.M.G.-V.); (T.R.-R.)
| | - Ángel Martín Pendás
- Department of Analytical and Physical Chemistry, University of Oviedo, E-33006 Oviedo, Spain;
| |
Collapse
|
17
|
Anisimov AA, Ananyev IV. Interatomic exchange-correlation interaction energy from a measure of quantum theory of atoms in molecules topological bonding: A diatomic case. J Comput Chem 2020; 41:2213-2222. [PMID: 32731310 DOI: 10.1002/jcc.26390] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2020] [Revised: 07/05/2020] [Accepted: 07/06/2020] [Indexed: 01/01/2023]
Abstract
The potential relations between the measure of topological interatomic bonding-integrals of electron density with respect to internuclear axis over the corresponding quantum theory of atoms in molecules (QTAIM)-defined interatomic surface (IAS)-and interatomic exchange-correlation contributions from the interacting quantum atoms approach are discussed. The quantum chemical computations of 38 equilibrium diatomic systems at different levels of theory (HF, MP2, MP4SDQ, and CCSD) are invoked to support abstract considerations. Parameters of excellent correlations between IAS integrals and interatomic exchange-correlation energy are found by the optimization. The performance of these trends depends on the accuracy of the electronic correlation treatment. The resulting trends are a unique feature of equilibrium states, whereas more complicated dependencies are explored for several systems at non-equilibrium conditions. The relations of established trends with other IAS-based estimations of strength of bonding interactions between topological atoms and issues explored for multiatomic systems are briefly discussed.
Collapse
Affiliation(s)
- Aleksei A Anisimov
- A.N. Nesmeyanov Institute of Organoelement Compounds, Russian Academy of Sciences, Vavilova str. 28, Moscow, 119991, GSP-1, Russia
| | - Ivan V Ananyev
- A.N. Nesmeyanov Institute of Organoelement Compounds, Russian Academy of Sciences, Vavilova str. 28, Moscow, 119991, GSP-1, Russia.,National Research University Higher School of Economics, Miasnitskaya Str. 20, Moscow, 101000, Russia
| |
Collapse
|
18
|
Milovanović B, Stanojević A, Etinski M, Petković M. Intriguing Intermolecular Interplay in Guanine Quartet Complexes with Alkali and Alkaline Earth Cations. J Phys Chem B 2020; 124:3002-3014. [DOI: 10.1021/acs.jpcb.0c01165] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Branislav Milovanović
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11158 Belgrade, Serbia
| | - Ana Stanojević
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11158 Belgrade, Serbia
| | - Mihajlo Etinski
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11158 Belgrade, Serbia
| | - Milena Petković
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11158 Belgrade, Serbia
| |
Collapse
|
19
|
Hayashi S, Nishide T, Nakanishi W. Nature of intramolecular O-H⋯π interactions as elucidated by QTAIM dual functional analysis with QC calculations. RSC Adv 2019; 9:15521-15530. [PMID: 35514837 PMCID: PMC9064312 DOI: 10.1039/c9ra01788g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2019] [Accepted: 05/09/2019] [Indexed: 11/21/2022] Open
Abstract
The intrinsic dynamic and static nature of intramolecular OH-*-π interactions is elucidated using a QTAIM dual functional analysis (QTAIM-DFA) after clarifying the structural features. Asterisks (*) are employed to emphasize the presence of bond critical points (BCPs) on the bond paths (BPs), which correspond to the interactions in question. Data from the fully optimized structures correspond to the static nature of the interactions. In our treatment, data from the perturbed structures, which are based around the fully optimized structure, are employed for the analysis in addition to those from the fully optimized structure, which represent the dynamic nature of the interaction. Seven intramolecular OH-*-C(π) interactions were detected in six-membered rings, with six BPs and BCPs for each, among the 72 conformers of the species examined here (1-15). The interactions are predicted to have a vdW or t-HBnc (typical hydrogen bonds with no covalency) nature, which appeared in the pure closed shell region. They appear to be stronger than the corresponding intermolecular interactions. Nine BPs with BCPs were also detected for the intramolecular O-*-X interactions (X = C(π) and H(π), joined to C(π)) in the 5-7-membered rings. The E(2) values of the interactions, as obtained by NBO, are discussed in relation to the stabilities of the conformers and the BPs with BCPs.
Collapse
Affiliation(s)
- Satoko Hayashi
- Faculty of Systems Engineering, Wakayama University 930 Sakaedani Wakayama 640-8510 Japan +81 73 457 8253 +81 73 457 8252
| | - Taro Nishide
- Faculty of Systems Engineering, Wakayama University 930 Sakaedani Wakayama 640-8510 Japan +81 73 457 8253 +81 73 457 8252
| | - Waro Nakanishi
- Faculty of Systems Engineering, Wakayama University 930 Sakaedani Wakayama 640-8510 Japan +81 73 457 8253 +81 73 457 8252
| |
Collapse
|
20
|
de Kock S, Dillen J, Esterhuysen C. Steric and Electronic Effects in Gold N-Heterocyclic Carbene Complexes Revealed by Computational Analysis. ChemistryOpen 2019; 8:539-550. [PMID: 31061779 PMCID: PMC6488199 DOI: 10.1002/open.201900076] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2019] [Indexed: 11/12/2022] Open
Abstract
A computational analysis of a series of cationic and neutral gold imidazolylidene and benzimidizolylidene complexes is reported. The Bond Dissociation Energies of the various ligands in the complexes calculated at the PBE0-D3/def2-TZVP level of theory increase with increasing ligand volume, except for those of complexes containing t-butyl-substituted ligands, which are anomalously low particularly for the benzimidazolylidene species. Atoms in Molecules studies show the presence of a variety of weak intramolecular interactions, characterised by the presence of bond critical points with a range of different properties. Energy Decomposition Analysis and calculation of Electrostatic Surface Potentials indicate that some interactions are weakly attractive dispersion-type interactions, while others are repulsive. The octanol/water partition coefficients (log P values) were calculated as a measure of the lipophilicities of the complexes and were found to increase with increasing volume.
Collapse
Affiliation(s)
- Sunel de Kock
- Department of Chemistry and Polymer ScienceStellenbosch UniversityPrivate Bag X1, MatielandStellenbosch7602South Africa
| | - Jan Dillen
- Department of Chemistry and Polymer ScienceStellenbosch UniversityPrivate Bag X1, MatielandStellenbosch7602South Africa
| | - Catharine Esterhuysen
- Department of Chemistry and Polymer ScienceStellenbosch UniversityPrivate Bag X1, MatielandStellenbosch7602South Africa
| |
Collapse
|
21
|
Tsirelson VG, Stash AI, Tokatly IV. Quantum pressure focusing in solids: a reconstruction from experimental electron density. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2019; 75:201-209. [PMID: 32830745 DOI: 10.1107/s2052520619001471] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/26/2018] [Accepted: 01/25/2019] [Indexed: 06/11/2023]
Abstract
Here an approach is presented for reconstructing the distribution of electronic internal quantum pressure in the electronic continuum of solids from the experimental electron density. Using the formalism of the density functional theory, the spatial inner-crystal map of the quantum pressure is obtained. The results are visualized via the indicator of quantum pressure focusing (IQPF) which reveals the regions where the pressure is concentrated or depleted due to quantum effects. IQPF contains all quantum electron-shell structure-forming contributions resulting from kinetic, exchange and correlation effects, and presents a clear picture of the chemical bond features in crystals with different type of bonding mechanisms.
Collapse
Affiliation(s)
- Vladimir G Tsirelson
- Quantum Chemistry Department, Mendeleev University of Chemical Technology, Moscow, 125047, Russian Federation
| | - Adam I Stash
- Quantum Chemistry Department, Mendeleev University of Chemical Technology, Moscow, 125047, Russian Federation
| | - Ilya V Tokatly
- Departamento de Fisica de Materials, Universidad del Pais Vasco UPV/EHU, San Sebastian, 20018, Spain
| |
Collapse
|
22
|
Jabłoński M. On the Uselessness of Bond Paths Linking Distant Atoms and on the Violation of the Concept of Privileged Exchange Channels. ChemistryOpen 2019; 8:497-507. [PMID: 31019875 PMCID: PMC6470636 DOI: 10.1002/open.201900109] [Citation(s) in RCA: 47] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2019] [Indexed: 11/09/2022] Open
Abstract
We refer to frequently used determinants suggesting dominant interactions between distant atoms in various dimers. First of all, we show, against the still-prevailling opinion, that, in general, bond paths have nothing in common with dominant intermolecular interactions and therefore they are useless in such cases. Quite the contrary, reliable information about dominant intermolecular interactions can be obtained by means of electrostatic potential maps, which very convincingly explain mutual orientation of molecules in a dimer. For the first time, numerous examples of interactions that violate both the concept of privileged exchange channels proposed by Pendás and his collaborators as well as inequalities obtained by Tognetti and Joubert for the β parameter related to secondary interactions are presented. The possible cause of this violation is suggested. We also show that the so-called counterintuitive bond paths result from quite natural behavior of the electron density gradient vector, i. e. searching for those areas of space that are characterized by large values of electron density or the most expanded its distributions.
Collapse
Affiliation(s)
- Mirosław Jabłoński
- Department of Quantum Chemistry, Faculty of ChemistryNicolaus Copernicus University in Toruń7-Gagarina St.87-100ToruńPoland
| |
Collapse
|
23
|
σ-Holes and Si···N intramolecular interactions. J Mol Model 2019; 25:101. [DOI: 10.1007/s00894-019-3962-2] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2018] [Accepted: 02/07/2019] [Indexed: 10/27/2022]
|
24
|
Jabłoński M. Bond paths between distant atoms do not necessarily indicate dominant interactions. J Comput Chem 2018; 39:2183-2195. [PMID: 30298926 DOI: 10.1002/jcc.25532] [Citation(s) in RCA: 39] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2018] [Revised: 06/27/2018] [Accepted: 06/27/2018] [Indexed: 12/16/2022]
Abstract
The goal of the article is to revive discussion on the interpretation of bond paths linking distant atoms, particularly tracing weak interactions in dimers. According to the Pendás' concept of privileged exchange channel, a bond path is formed between this pair of competing atoms, which is associated with larger value of the exchange energy. We point out that, due to the short-range nature of the exchange energy, bond paths linking distant atoms clearly become doubtful indicators of dominant intermolecular interactions, particularly if some other characteristics (geometric, spectroscopic, based on electrostatic parameters, etc.) indicate other intermolecular interactions as dominant. Several such cases are thoroughly investigated. We show that electrostatic parameters are much more reliable indicators of dominant intermolecular interactions than bond paths. Then, we pay attention that the presence of ("unexpected", i.e., not necessarily indicating dominant intermolecular interactions) bond paths between pairs of atoms featuring highly expanded charge distributions can be easily explained by visual exploration of isodensity contour plots. As always pointing in the direction of the steepest increase, the gradient vector of the electron density favors areas of its high values gaining higher exchange energy, yet being blind to highly electron deficient areas nearby, which, however, can quite often be involved in dominant intermolecular interactions as strongly suggested by many other bonding analysis. We also suggest that an interatomic component of Hellmann-Feynman force would most likely be the most reliable indicator of attractive or repulsive character of individual interatomic interaction. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Mirosław Jabłoński
- Department of Quantum Chemistry, Faculty of Chemistry, Nicolaus Copernicus University in Toruń, 7-Gagarina Street, 87-100, Toruń, Poland
| |
Collapse
|
25
|
van Niekerk DME, Gerber WJ. A DFT Mechanistic Study of the trans-[Os VIO 2(OH) 4] 2– and [Os VIIIO 4(OH) n] n− ( n = 1, 2 cis) Comproportionation Proton-Coupled Electron Transfer Reaction. Inorg Chem 2018; 57:8909-8922. [DOI: 10.1021/acs.inorgchem.8b00840] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Daniël M. E. van Niekerk
- Department of Chemistry and Polymer Science, Stellenbosch University, Private Bag X1, Stellenbosch, 7602 Western Cape, South Africa
| | - Wilhelmus J. Gerber
- Department of Chemistry and Polymer Science, Stellenbosch University, Private Bag X1, Stellenbosch, 7602 Western Cape, South Africa
| |
Collapse
|
26
|
Wick CR, Clark T. On bond-critical points in QTAIM and weak interactions. J Mol Model 2018; 24:142. [DOI: 10.1007/s00894-018-3684-x] [Citation(s) in RCA: 73] [Impact Index Per Article: 12.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2018] [Accepted: 05/21/2018] [Indexed: 11/29/2022]
|
27
|
Hayashi S, Sugibayashi Y, Nakanishi W. Behaviour of the XH-*-π and YX-*-π interactions (X, Y = F, Cl, Br and I) in the coronene π-system, as elucidated by QTAIM dual functional analysis with QC calculations. RSC Adv 2018; 8:16349-16361. [PMID: 35542236 PMCID: PMC9080335 DOI: 10.1039/c8ra01862f] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2018] [Accepted: 04/24/2018] [Indexed: 11/24/2022] Open
Abstract
The dynamic and static nature of XH-*-π and YX-*-π in the coronene π-system (π(C24H12)) is elucidated by QTAIM dual functional analysis, where * emphasizes the presence of bond critical points (BCPs) in the interactions. The nature of the interactions is elucidated by analysing the plots of the total electron energy densities H b(r c) versus H b(r c) - V b(r c)/2 [=(ħ 2/8m)∇2 ρ b(r c)] for the interactions at BCPs, where V b(r c) are the potential energy densities at the BCPs. The data for the perturbed structures around the fully optimized structures are employed for the plots in addition to those of the fully optimized structures. The plots are analysed using the polar coordinate of (R, θ) for the data of the fully optimized structures, while those containing the perturbed structures are analysed using (θ p, κ p), where θ p corresponds to the tangent line of each plot and κ p is the curvature. Whereas (R, θ) show the static nature, (θ p, κ p) represent the dynamic nature of the interactions. All interactions in X-H-*-π(C24H12) (X = F, Cl, Br and I) and Y-X-*-π(C24H12) (Y-X = F-F, Cl-Cl, Br-Br, I-I, F-Cl, F-Br and F-I) are classified by pure CS (closed shell) interactions and are characterized as having the vdW nature, except for X-H = F-H and Y-X = F-Cl, F-Br and F-I, which show the typical-HB nature without covalency. The structural features of the complexes are also discussed.
Collapse
Affiliation(s)
- Satoko Hayashi
- Faculty of Systems Engineering, Wakayama University 930 Sakaedani Wakayama 640-8510 Japan +81 73 457 8253 +81 73 457 8252
| | - Yuji Sugibayashi
- Faculty of Systems Engineering, Wakayama University 930 Sakaedani Wakayama 640-8510 Japan +81 73 457 8253 +81 73 457 8252
| | - Waro Nakanishi
- Faculty of Systems Engineering, Wakayama University 930 Sakaedani Wakayama 640-8510 Japan +81 73 457 8253 +81 73 457 8252
| |
Collapse
|
28
|
Tsubomoto Y, Hayashi S, Nakanishi W, Mapp LK, Coles SJ. High-resolution X-ray diffraction determination of the electron density of 1-(8-PhSC 10H 6)SS(C 10H 6SPh-8')-1' with the QTAIM approach: evidence for S 4 σ(4c-6e) at the naphthalene peri-positions. RSC Adv 2018; 8:9651-9660. [PMID: 35540809 PMCID: PMC9078645 DOI: 10.1039/c7ra13636f] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/25/2017] [Accepted: 02/20/2018] [Indexed: 12/18/2022] Open
Abstract
An extended hypervalent S4 σ(4c-6e) system was confirmed for the linear BS-∗-AS-∗-AS-∗-BS interaction in 1-(8-PhBSC10H6)AS-AS(C10H6 BSPh-8')-1' (1) via high-resolution X-ray diffraction determination of electron densities. The presence of bond critical points (BCPs; ∗) on the bond paths confirms the nature and extent of this interaction. The recently developed QTAIM dual functional analysis (QTAIM-DFA) approach was also applied to elucidate the nature of the interaction. Total electron energy densities H b( r c) were plotted versus H b( r c) - V b( r c)/2 for the interaction at the BCPs, where V b( r c) represents the potential energy densities at the BCP. The results indicate that although the data for an interaction in the fully optimized structure corresponds to a static nature, the data obtained for the perturbed structures around it represent the dynamic nature of the interaction in QTAIM-DFA. The former classifies the interaction and the latter characterises it. Although AS-∗-AS in 1 is classified by a shared shell interaction and exhibits weak covalent character, AS-∗-BS is characterized as having typical hydrogen-bond nature with covalent properties in the region of the regular closed shell interactions. The experimental results are supported by matching theoretical calculations throughout, particularly for the extended hypervalent E4 σ(4c-6e) (E = S) interaction.
Collapse
Affiliation(s)
- Yutaka Tsubomoto
- Faculty of Systems Engineering, Wakayama University 930 Sakaedani Wakayama 640-8510 Japan +81 73 457 8253 +81 73 457 8252
| | - Satoko Hayashi
- Faculty of Systems Engineering, Wakayama University 930 Sakaedani Wakayama 640-8510 Japan +81 73 457 8253 +81 73 457 8252
| | - Waro Nakanishi
- Faculty of Systems Engineering, Wakayama University 930 Sakaedani Wakayama 640-8510 Japan +81 73 457 8253 +81 73 457 8252
| | - Lucy K Mapp
- Chemistry, Faculty of Natural and Environmental Sciences, University of Southampton Southampton SO17 1BJ UK +44 (0)2380596721
| | - Simon J Coles
- Chemistry, Faculty of Natural and Environmental Sciences, University of Southampton Southampton SO17 1BJ UK +44 (0)2380596721
| |
Collapse
|
29
|
de Lange JH, van Niekerk DME, Cukrowski I. FALDI-based decomposition of an atomic interaction line leads to 3D representation of the multicenter nature of interactions. J Comput Chem 2018; 39:973-985. [DOI: 10.1002/jcc.25175] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2017] [Revised: 12/07/2017] [Accepted: 01/14/2018] [Indexed: 11/09/2022]
Affiliation(s)
- Jurgens H. de Lange
- Department of Chemistry; Faculty of Natural and Agricultural Sciences, University of Pretoria, Lynnwood Road; Hatfield Pretoria 0002 South Africa
| | - Daniel M. E. van Niekerk
- Department of Chemistry; Faculty of Natural and Agricultural Sciences, University of Pretoria, Lynnwood Road; Hatfield Pretoria 0002 South Africa
| | - Ignacy Cukrowski
- Department of Chemistry; Faculty of Natural and Agricultural Sciences, University of Pretoria, Lynnwood Road; Hatfield Pretoria 0002 South Africa
| |
Collapse
|
30
|
Ferreira LM, Eaby A, Dillen J. The topology of the Coulomb potential density. A comparison with the electron density, the virial energy density, and the Ehrenfest force density. J Comput Chem 2017; 38:2784-2790. [PMID: 28963855 DOI: 10.1002/jcc.25071] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2017] [Revised: 09/05/2017] [Accepted: 09/06/2017] [Indexed: 01/25/2023]
Abstract
The topology of the Coulomb potential density has been studied within the context of the theory of Atoms in Molecules and has been compared with the topologies of the electron density, the virial energy density and the Ehrenfest force density. The Coulomb potential density is found to be mainly structurally homeomorphic with the electron density. The Coulomb potential density reproduces the non-nuclear attractor which is observed experimentally in the molecular graph of the electron density of a Mg dimer, thus, for the first time ever providing an alternative and energetic foundation for the existence of this critical point. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Lizé-Mari Ferreira
- Department of Chemistry and Polymer Science, Stellenbosch University, Private Bag ×1, Matieland, 7602, South Africa
| | - Alan Eaby
- Department of Chemistry and Polymer Science, Stellenbosch University, Private Bag ×1, Matieland, 7602, South Africa
| | - Jan Dillen
- Department of Chemistry and Polymer Science, Stellenbosch University, Private Bag ×1, Matieland, 7602, South Africa
| |
Collapse
|
31
|
Dem'yanov PI, Polestshuk PM, Kostin VV. The Nature of Metal-Metal Interactions in Dimeric Hydrides and Halides of Group 11 Elements in the Light of High Level Relativistic Calculations. Chemistry 2017; 23:3257-3261. [PMID: 28093856 DOI: 10.1002/chem.201605519] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2016] [Indexed: 11/08/2022]
Abstract
The titular calculations show that charges at metal atoms M are apparently the main factor governing the nature of M⋅⋅⋅M interactions in two-nuclear coinage-metal complexes, and there are certain critical values of positive charges on M atoms, on exceeding which the pair-wise M⋅⋅⋅M interactions and/or the binding between M atoms in such complexes become repulsive despite negative formation energies of such complexes, short M-M internuclear distances, and the existence of a bond critical point (BCP) between M atoms.
Collapse
Affiliation(s)
- Piotr I Dem'yanov
- Chemistry Department, M. V. Lomonosov Moscow State University, Lenin Hills 1, Building 3, 119991, Moscow, Russia
| | - Pavel M Polestshuk
- Chemistry Department, M. V. Lomonosov Moscow State University, Lenin Hills 1, Building 3, 119991, Moscow, Russia
| | - Vladimir V Kostin
- Chemistry Department, M. V. Lomonosov Moscow State University, Lenin Hills 1, Building 3, 119991, Moscow, Russia
| |
Collapse
|
32
|
de Lange JH, Cukrowski I. Toward deformation densities for intramolecular interactions without radical reference states using the fragment, atom, localized, delocalized, and interatomic (FALDI) charge density decomposition scheme. J Comput Chem 2017; 38:981-997. [DOI: 10.1002/jcc.24772] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2016] [Revised: 01/09/2017] [Accepted: 02/02/2017] [Indexed: 01/30/2023]
Affiliation(s)
- Jurgens Hendrik de Lange
- Department of Chemistry, Faculty of Natural and Agricultural Sciences; University of Pretoria; Lynnwood Road Hatfield Pretoria 0002 South Africa
| | - Ignacy Cukrowski
- Department of Chemistry, Faculty of Natural and Agricultural Sciences; University of Pretoria; Lynnwood Road Hatfield Pretoria 0002 South Africa
| |
Collapse
|
33
|
Hayashi S, Sugibayashi Y, Nakanishi W. Behavior of interactions between hydrogen chalcogenides and an anthracene π-system elucidated by QTAIM dual functional analysis with QC calculations. RSC Adv 2017. [DOI: 10.1039/c7ra04224h] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
The nature of the interactions between chalcogenides and the anthracene p-system, EH2-*-p(C14H10), is predicted to be close to that of EH2-*-p(C10H8), although the partial structures around the central rings can be found in EH2-*-p(C6H6).
Collapse
Affiliation(s)
- Satoko Hayashi
- Department of Material Science and Chemistry
- Faculty of Systems Engineering
- Wakayama University
- Wakayama
- 640-8510 Japan
| | - Yuji Sugibayashi
- Department of Material Science and Chemistry
- Faculty of Systems Engineering
- Wakayama University
- Wakayama
- 640-8510 Japan
| | - Waro Nakanishi
- Department of Material Science and Chemistry
- Faculty of Systems Engineering
- Wakayama University
- Wakayama
- 640-8510 Japan
| |
Collapse
|
34
|
Non-Covalent Interactions in Hydrogen Storage Materials LiN(CH3)2BH3 and KN(CH3)2BH3. CRYSTALS 2016. [DOI: 10.3390/cryst6030028] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
35
|
Cukrowski I, Mangondo P. Interacting quantum fragments-rooted preorganized-interacting fragments attributed relative molecular stability of the Be(II) complexes of nitrilotriacetic acid and nitrilotri-3-propionic acid. J Comput Chem 2016; 37:1373-87. [PMID: 26993356 DOI: 10.1002/jcc.24346] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2015] [Revised: 01/26/2016] [Accepted: 01/26/2016] [Indexed: 11/11/2022]
Abstract
A method designed to investigate, on a fundamental level, the origin of relative stability of molecular systems using Be(II) complexes with nitrilotriacetic acid (NTA) and nitrilotri-3-propionic acid (NTPA) is described. It makes use of the primary and molecular fragment energy terms as defined in the IQA/F (Interacting Quantum Atoms/Fragments) framework. An extensive classical-type investigation, focused on single descriptors (bond length, density at critical point, the size of metal ion or coordination ring, interaction energy between Be(II) and a donor atom, etc.) showed that it is not possible to explain the experimental trend. The proposed methodology is fundamentally different in that it accounts for the total energy contributions coming from all atoms of selected molecular fragments, and monitors changes in defined energy terms (e.g., fragment deformation, inter- and intra-fragment interaction) on complex formation. By decomposing combined energy terms we identified the origin of relative stability of Be(II) (NTA) and Be(II) (NTPA) complexes. We found that the sum of coordination bonds' strength, as measured by interaction energies between Be(II) ion and donor atoms, favours Be(II) (NTA) but the binding energy of Be(II) ion to the entire ligand correlates well with experimental trend. Surprisingly, the origin of Be(II) (NTPA) being more stable is due to less severe repulsive interactions with the backbone of NTPA (C and H-atoms). This general purpose protocol can be employed not only to investigate the origin of relative stability of any molecular system (e.g., metal complexes) but, in principle, can be used as a predictive tool for, e.g., explaining reaction mechanism. © 2016 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Ignacy Cukrowski
- Department of Chemistry, Faculty of Natural and Agricultural Sciences, University of Pretoria, Lynnwood Road, Hatfield, Pretoria, 0002, South Africa
| | - Paidamwoyo Mangondo
- Department of Chemistry, Faculty of Natural and Agricultural Sciences, University of Pretoria, Lynnwood Road, Hatfield, Pretoria, 0002, South Africa
| |
Collapse
|
36
|
Remya PR, Suresh CH. Planar tetracoordinate carbon in tungstenacyclobutadiene from alkyne metathesis and expanded structures. Dalton Trans 2016; 45:1769-78. [DOI: 10.1039/c5dt03922c] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Establishing the Cβ of tungstenacyclobutadiene (WCBD) as a ptC center paves the way for a new strategy to make novel materials containing multiple ptC centers. The 1-, 2- and 3-dimensional expansion of the WCBD motifs provides access to ptC-incorporated new metal–organic frameworks.
Collapse
Affiliation(s)
- Premaja R. Remya
- Chemical Sciences and Technology Division
- CSIR-National Institute for Interdisciplinary Science and Technology
- Trivandrum 695 019
- India
- Academy of Scientific & Innovative Research
| | - Cherumuttathu H. Suresh
- Chemical Sciences and Technology Division
- CSIR-National Institute for Interdisciplinary Science and Technology
- Trivandrum 695 019
- India
- Academy of Scientific & Innovative Research
| |
Collapse
|
37
|
|
38
|
Tsirelson VG, Stash AI, Tokatly IV. Bonding in molecular crystals from the local electronic pressure viewpoint. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1101173] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Affiliation(s)
- Vladimir G. Tsirelson
- Quantum Chemistry Department, Mendeleev University of Chemical Technology, Moscow, Russia
| | - Adam I. Stash
- Laboratory of Oxide Materials, Karpov Research Institute of Physical Chemistry, Moscow, Russia
| | - Ilya V. Tokatly
- Departamento de Fisica de Materials, Universidad del Pais Vasco UPV/EHU, San Sebastian, Spain
- IKERBASQUE, Basgue Foundation for Science, Bilbao, Spain
| |
Collapse
|
39
|
IQA-embedded fragment attributed molecular system energy change in exploring intramolecular interactions. COMPUT THEOR CHEM 2015. [DOI: 10.1016/j.comptc.2015.04.018] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
40
|
Wu LC, Lee GH. Experimental and Theoretical Charge Density Study on Di-2-pyrazylamine (Hdpza) Molecule in Crystal. J CHIN CHEM SOC-TAIP 2015. [DOI: 10.1002/jccs.201500110] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
|
41
|
Adeyinka AS, Cukrowski I. Structural-topological preferences and protonation sequence of aliphatic polyamines: a theoretical case study of tetramine trien. J Mol Model 2015; 21:162. [PMID: 26037149 DOI: 10.1007/s00894-015-2709-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2015] [Accepted: 05/19/2015] [Indexed: 11/28/2022]
Abstract
A large set of lowest and medium energy conformers of aliphatic tetramine trien was used to uncover structural-topological preferences of poliamines. Numerous common structural features among HL and H 2 L tautomers were identified, e.g., H-atoms of protonated functional groups are always involved in intramolecular NH•••N interactions and they result in as large and as many as possible rings in lowest energy conformers. Largest, 11-membered, molecular rings stabilize a molecule most and they appeared to be strain free whereas 5-memebred intramolecular rings were most strained (all formed due to NH•••N interactions). The CH•••HC interactions with QTAIM-defined atomic interaction lines were also found but, surprisingly, mainly in the lowest energy conformers of HL tautomers. According to the non-covalent interaction-based (NCI) analysis, 5-memebered rings formed by CH•••HC interactions are not strained and, in general, 3D NCI isosurfaces mimic those obtained for weaker NH•••N interactions. Also, 3D NCI isosurfaces found for NH•••N and CH•••HC interactions, regardless whether linked or not by an atomic interaction line, appeared to be indistinguishable. Using lowest energy conformers, theoretically predicted mixture of primary (HL p ) and secondary (HL s ) forms of trien was found to be in accord with the literature reports; using linear conformers resulted in predicting HL s as the only tautomer formed. In contrast to HF, the overall performance of B3LYP was found satisfactory for the purpose of the study.
Collapse
Affiliation(s)
- Adedapo S Adeyinka
- Department of Chemistry, Faculty of Natural and Agricultural Sciences, University of Pretoria, Lynnwood Road, Hatfield, 0002, Pretoria, South Africa
| | | |
Collapse
|
42
|
Bartashevich EV, Pendás ÁM, Tsirelson VG. An anatomy of intramolecular atomic interactions in halogen-substituted trinitromethanes. Phys Chem Chem Phys 2015; 16:16780-9. [PMID: 25001471 DOI: 10.1039/c4cp01257g] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
The intramolecular interactions in substituted trinitromethanes, XC(NO2)3 (X = F, Cl, I, H) are studied and clarified by using a combination of the Quantum Theory of Atoms in Molecules (QTAIM), the non-covalent interaction analysis and the Interacting Quantum Atoms (IQA) methods. The stretching vibration modes are formed by the concerted displacements of atoms involved in the covalent bonds showing the significant multiatomic influence in substituted trinitromethanes. In agreement with that, the arrangement of the local reduced density gradient minima indicates that the electron density favors the non-covalent intramolecular interactions X···O and N···O. However, the corresponding QTAIM bond paths are not formed; instead, contacts, which we call uncompleted links in this context, are accompanied by "quasi-bonding channels" corresponding to the λ2() ≤ 0 regions on the sign[λ2(r)]ρ(r) contour maps. The intramolecular IQA energy contributions signal the appreciable electron exchange between the pairs of atoms associated with potential atomic interactions or the bond-path-free non-covalent links. The IQA analysis shows that the electrostatic term destabilizes FC(NO2)3 and distinctly stabilizes IC(NO2)3, whereas it is close to neutral in ClC(NO2)3. The exchange energy between the X atom and the NO2 groups, in contrast, stabilizes all the molecules.
Collapse
Affiliation(s)
- Ekaterina V Bartashevich
- Department of Chemistry, South Ural State University (National Research University), 454080 Chelyabinsk, Russia.
| | | | | |
Collapse
|
43
|
Dillen J. The topology of the Ehrenfest force density revisited. A different perspective based on Slater-type orbitals. J Comput Chem 2015; 36:883-90. [DOI: 10.1002/jcc.23869] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2014] [Revised: 01/27/2015] [Accepted: 01/29/2015] [Indexed: 01/20/2023]
Affiliation(s)
- Jan Dillen
- Department of Chemistry and Polymer Science; Private Bag X1; Matieland 7602 South Africa
| |
Collapse
|
44
|
Cukrowski I, de Lange JH, Adeyinka AS, Mangondo P. Evaluating common QTAIM and NCI interpretations of the electron density concentration through IQA interaction energies and 1D cross-sections of the electron and deformation density distributions. COMPUT THEOR CHEM 2015. [DOI: 10.1016/j.comptc.2014.10.005] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
45
|
Safin DA, Babashkina MG, Kubisiak P, Mitoraj MP, Le Duff CS, Robeyns K, Garcia Y. Supramolecular Coordination Complexes of theN-Thiophosphorylated 2,5-Dithiobiurea [NHC(S)NHP(S)(OiPr)2]2with ZnIIand CdIIIons - Cation-Induced Dinuclear Mesocate Structure versus Tetranuclear Nanoscaled Aggregate. Eur J Inorg Chem 2014. [DOI: 10.1002/ejic.201402544] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
46
|
Foroutan-Nejad C, Shahbazian S, Marek R. Toward a Consistent Interpretation of the QTAIM: Tortuous Link between Chemical Bonds, Interactions, and Bond/Line Paths. Chemistry 2014; 20:10140-52. [DOI: 10.1002/chem.201402177] [Citation(s) in RCA: 166] [Impact Index Per Article: 16.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2014] [Indexed: 11/10/2022]
|
47
|
Galembeck SE, Bickelhaupt FM, Fonseca Guerra C, Galembeck E. Effects of the protonation state in the interaction of an HIV-1 reverse transcriptase (RT) amino acid, Lys101, and a non nucleoside RT inhibitor, GW420867X. J Mol Model 2014; 20:2332. [DOI: 10.1007/s00894-014-2332-3] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2014] [Accepted: 05/30/2014] [Indexed: 01/15/2023]
|
48
|
Firouzi R, Shahbazian S. Seeking for ultrashort “non-bonded” hydrogen–hydrogen contacts in some rigid hydrocarbons and their chlorinated derivatives. Struct Chem 2014. [DOI: 10.1007/s11224-014-0411-9] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
49
|
Liu JB, Schwarz WHE, Li J. On Two Different Objectives of the Concepts of Ionic Radii. Chemistry 2013; 19:14758-67. [DOI: 10.1002/chem.201300917] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
50
|
Guevara-Vela JM, Chávez-Calvillo R, García-Revilla M, Hernández-Trujillo J, Christiansen O, Francisco E, Martín Pendás Á, Rocha-Rinza T. Hydrogen-Bond Cooperative Effects in Small Cyclic Water Clusters as Revealed by the Interacting Quantum Atoms Approach. Chemistry 2013; 19:14304-15. [DOI: 10.1002/chem.201300656] [Citation(s) in RCA: 73] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2013] [Revised: 05/11/2013] [Indexed: 11/08/2022]
|