1
|
Labet V, Geoffroy-Neveux A, Alikhani ME. How to search for and reveal a hidden intermediate? The ELF topological description of non-synchronicity in double proton transfer reactions under oriented external electric field. J Mol Model 2024; 30:367. [PMID: 39365459 DOI: 10.1007/s00894-024-06163-0] [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: 06/20/2024] [Accepted: 09/23/2024] [Indexed: 10/05/2024]
Abstract
CONTEXT The nature of double intermolecular proton transfer was studied with the ELF topological approach in two model dimers (the formic acid homodimer and the 1,2,3-triazole-guanidine heterodimer) under an oriented external electric field. It has been shown that each of the two dimers can have either a one-step (one transition state structure) or two-step (two transition state structures) reaction path, depending on the intensity and orientation of the external electric field. The presence of a singularly broad shoulder (plateau in the case of homodimer and plateau-like for heterodimer) around the formal transition state structure results from the strong asynchronicity of the reaction. A careful ELF topological analysis of the nature of protons, hydride (localized) or roaming (delocalized) proton, along the reaction path allowed us to unambiguously classify the one-step mechanisms governing the double-proton transfer reactions into three distinct classes: (1) concerted-synchronous, when two events (roaming proton regions) completely overlap, (2) concerted-asynchronous, when two events (roaming proton regions) partially overlap, and (3) two-stage one-step non-concerted, when two roaming proton regions are separated by a "hidden intermediate region". All the structures belonging to this separatrix region are of the zwitterion form. METHODS Geometry optimization of the stationary points on the potential energy surface was performed using density functional theory-wB97XD functional-in combination with the 6-311+ +G(2d, 2p) basis set for all the atoms. All first-principles calculations were performed using the Gaussian 09 quantum chemical packages. We also used the electron localization function (ELF) to reveal the nature of the proton along the reaction path: a bound proton (hydride) becomes a roaming proton (carrying a tiny negative charge ≈ 0.3 e) exchanging with two adjacent atoms via two attractors (topological critical points with (3, - 3) signature). The ELF analyses were performed using the TopMod package.
Collapse
Affiliation(s)
- Vanessa Labet
- Sorbonne Université CNRS, MONARIS, UMR8233, 75005, Paris, France.
| | | | | |
Collapse
|
2
|
Krafft MP, Riess JG. About Perfluoropolyhedranes, Their Electron-Accepting Ability and Questionable Supramolecular Hosting Capacity. Angew Chem Int Ed Engl 2023; 62:e202302942. [PMID: 37208990 DOI: 10.1002/anie.202302942] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Revised: 05/16/2023] [Accepted: 05/17/2023] [Indexed: 05/21/2023]
Abstract
Polyhedral molecules are appealing for their eye-catching architecture and distinctive chemistry. Perfluorination of such, often greatly strained, compounds is a momentous challenge. It drastically changes the electron distribution, structure and properties. Notably, small high-symmetry perfluoropolyhedranes feature a centrally located, star-shaped low-energy unoccupied molecular orbital that can host an extra electron within the polyhedral frame, thus producing a radical anion, without loss of symmetry. This predicted electron-hosting capacity was definitively established for perfluorocubane, the first perfluorinated Platonic polyhedrane to be isolated pure. Hosting atoms, molecules, or ions in such "cage" structures is, however, all but forthright, if not illusionary, offering no easy access to supramolecular constructs. While adamantane and cubane have fostered numerous applications in materials science, medicine, and biology, specific uses for their perfluorinated counterparts remain to be established. Some aspects of highly fluorinated carbon allotropes, such as fullerenes and graphite, are briefly mentioned for context.
Collapse
Affiliation(s)
- Marie Pierre Krafft
- Institut Charles Sadron (CNRS), University of Strasbourg, 23 rue du Loess., 67034, Strasbourg Cedex, France
| | - Jean G Riess
- Harangoutte Institute, 68160, Ste-Croix-aux-Mines, France
| |
Collapse
|
3
|
Sugiyama M, Akiyama M, Yonezawa Y, Komaguchi K, Higashi M, Nozaki K, Okazoe T. Electron in a cube: Synthesis and characterization of perfluorocubane as an electron acceptor. Science 2022; 377:756-759. [PMID: 35951682 DOI: 10.1126/science.abq0516] [Citation(s) in RCA: 19] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Fluorinated analogs of polyhedral hydrocarbons have been predicted to localize an electron within their cages upon reduction. Here, we report the synthesis and characterization of perfluorocubane, a stable polyhedral fluorocarbon. The key to the successful synthesis was the efficient introduction of multiple fluorine atoms to cubane by liquid-phase reaction with fluorine gas. The solid-state structure of perfluorocubane was confirmed using x-ray crystallography, and its electron-accepting character was corroborated electrochemically and spectroscopically. The radical anion of perfluorocubane was examined by matrix-isolation electron spin resonance spectroscopy, which revealed that the unpaired electron accepted by perfluorocubane is located predominantly inside the cage.
Collapse
Affiliation(s)
- Masafumi Sugiyama
- Department of Chemistry and Biotechnology, Graduate School of Engineering, The University of Tokyo, Bunkyo-ku, Tokyo 113-8656, Japan
| | - Midori Akiyama
- Department of Chemistry and Biotechnology, Graduate School of Engineering, The University of Tokyo, Bunkyo-ku, Tokyo 113-8656, Japan
| | - Yuki Yonezawa
- Department of Chemistry and Biotechnology, Graduate School of Engineering, The University of Tokyo, Bunkyo-ku, Tokyo 113-8656, Japan
| | - Kenji Komaguchi
- Department of Applied Chemistry, Graduate School of Advanced Science and Engineering, Hiroshima University, Kagamiyama, Higashi-Hiroshima, Hiroshima 739-8527, Japan
| | - Masahiro Higashi
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Nishikyo-ku, Kyoto 615-8510, Japan
| | - Kyoko Nozaki
- Department of Chemistry and Biotechnology, Graduate School of Engineering, The University of Tokyo, Bunkyo-ku, Tokyo 113-8656, Japan
| | - Takashi Okazoe
- Department of Chemistry and Biotechnology, Graduate School of Engineering, The University of Tokyo, Bunkyo-ku, Tokyo 113-8656, Japan.,AGC Inc., Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan
| |
Collapse
|
4
|
Geoffroy-Neveux A, Labet V, Alikhani ME. Influence of an Oriented External Electric Field on the Mechanism of Double Proton Transfer between Pyrazole and Guanidine: from an Asynchronous Plateau Transition State to a Synchronous or Stepwise Mechanism. J Phys Chem A 2022; 126:3057-3071. [PMID: 35544749 DOI: 10.1021/acs.jpca.1c10553] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The double proton transfer (DPT) reaction between pyrazole and guanidine, a concerted reaction but strongly asynchronous and presenting a "plateau transition region", has been theoretically reinvestigated in the presence of an external uniform electric field. First, we computed the reaction path by DFT and proposed a very detailed description of the constitutive electronic events, based on the ELF topology and the bond evolution theory. Then, we studied the effect of an oriented external electric field (OEEF) on the reaction mechanism, for an OEEF oriented along the proton transfer axis. We observe that in one direction, the DPT reaction can be transformed into a stepwise reaction, going through a stabilized single proton transferred intermediate. Contrarily, the two proton transfers occur simultaneously when the electric field is applied in the opposite direction. In the latter case, the order in which the two protons are transferred in the same elementary step can even be reversed if the OEEF is intense enough. Finally, it has been shown that the evolution of the double proton transfer reaction in the presence of an electric field could be quantitatively anticipated by analyzing the ELF value at the bifurcation point between V(A, H) proton donor and V(B) proton acceptor of the double hydrogen bonded complex in the entrance channel.
Collapse
Affiliation(s)
| | - Vanessa Labet
- MONARIS UMR 8233 CNRS, Sorbonne Université, 4 place Jussieu, 75252 Paris Cedex 05, France
| | - M Esmail Alikhani
- MONARIS UMR 8233 CNRS, Sorbonne Université, 4 place Jussieu, 75252 Paris Cedex 05, France
| |
Collapse
|
5
|
Adjieufack AI, Liégeois V, Mbouombouo Ndassa I, Champagne B. Topological investigation of the reaction mechanism of glycerol carbonate decomposition by bond evolution theory. RSC Adv 2021; 11:10083-10093. [PMID: 35423535 PMCID: PMC8695523 DOI: 10.1039/d0ra09755a] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2020] [Accepted: 02/26/2021] [Indexed: 11/30/2022] Open
Abstract
The reaction mechanisms of the decomposition of glycerol carbonate have been investigated at the density functional theory level within the bond evolution theory. The four reaction pathways yield to 3-hydroxypropanal (TS1), glycidol (TS2a and TS2b), and 4-methylene-1,3-dioxolan-2-one (TS3). The study reveals non-concerted processes with the same number (four) of structural stability domains for each reaction pathway. For the two decarboxylation mechanisms, the two first steps are similar. They correspond to the cleavage of two single CO bonds to the detriment of the increased population of the lone pairs of two O atoms. These are followed, along TS1, by the transformation of a CO single bond into a double bond together with a proton transfer to create a CH bond. For TS2a and TS2b, the last step is a cyclization by CO bond formation. For the TS3 pathway, the first stage consists in the cleavage of a CH bond and the transfer of its electron population to both a proton and a C atom, the second step corresponds to the formation of an OH bond, and the last one describes the formation of a CC double bond. Moreover, the analysis of the energies, enthalpies, and free enthalpies of reaction and of activation leads to the conclusion that 3-hydroxypropanal is both the thermodynamic and kinetic product, independent of the method of calculation. ELFs of glycerol carbonate and of its kinetic and thermodynamic decomposition product, 3-hydroxypropanal (+CO2).![]()
Collapse
Affiliation(s)
- Abel Idrice Adjieufack
- Physical and Theoretical Chemistry Laboratory
- Cameroon
- Computational Chemistry Laboratory
- High Teacher Training College
- Cameroon
| | - Vincent Liégeois
- Laboratory of Theoretical Chemistry and Namur Institute of Structured Matter (NISM)
- University of Namur
- B-5000 Namur
- Belgium
| | | | - Benoît Champagne
- Laboratory of Theoretical Chemistry and Namur Institute of Structured Matter (NISM)
- University of Namur
- B-5000 Namur
- Belgium
| |
Collapse
|
6
|
Adjieufack AI, Moto Ongagna J, Kenmogne Tchidjo JF, Mbouombouo Ndassa I. Topological unraveling of the [3+2] cycloaddition (32CA) reaction between N-methylphenylnitrone and styrene catalyzed by the chromium tricarbonyl complex using electron localization function and catastrophe theory. NEW J CHEM 2021. [DOI: 10.1039/d1nj04121e] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
We have investigated the reaction mechanisms of [3+2] cycloaddition (32CA) between N-methylphenylnitrone and styrene catalyzed by the chromium tricarbonyl complex at the MPWB1K/6-311G(d,p) level of approximation.
Collapse
Affiliation(s)
- Abel Idrice Adjieufack
- Physical and Theoretical Chemistry Laboratory, University of Yaoundé I, P.O. Box 812, Yaoundé, Cameroon
- Computational Chemistry Laboratory, High Teacher Training College, University of Yaoundé 1, Cameroon. P.O. Box 47, Yaoundé, Cameroon
| | - Jean Moto Ongagna
- Computational and Theoretical Chemistry Unit, Department of Chemistry, Faculty of Science, University of Douala, P.O. Box 24157, Douala, Cameroon
| | | | - Ibrahim Mbouombouo Ndassa
- Computational Chemistry Laboratory, High Teacher Training College, University of Yaoundé 1, Cameroon. P.O. Box 47, Yaoundé, Cameroon
| |
Collapse
|
7
|
Laplaza R, Polo V, Contreras-García J. Localizing electron density errors in density functional theory. Phys Chem Chem Phys 2019; 21:20927-20938. [PMID: 31517339 DOI: 10.1039/c9cp02806d] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The accuracy of different density functional approximations is assessed through the use of quantum chemical topology on molecular electron densities. In particular, three simple yet ever-important systems are studied: N2, CO and ethane. Our results exemplify how real-space descriptors can help understand the sources of errors in density functional theory, avoiding unwanted error compensation present in simplified statistical metrics. Errors in "well-built" functionals are shown to be concentrated in chemically meaningful regions of space, and hence they are predictable. Conversely, strongly parametrized functionals show isotropic errors that cannot be traced back to chemically transferable units. Moreover, we will show that energetic corrections are mapped back into improvements in the density in chemically meaningful regions. These results point at the relevance of real-space perspectives when parametrizing or relating energy and density errors.
Collapse
Affiliation(s)
- Rubén Laplaza
- Sorbonne Université, CNRS, Laboratoire de Chimie Thèorique, LCT, F. 75005 Paris, France and Departamento de Química Física and Instituto de Biocomputación y Física de Sistemas Complejos (BIFI), Universidad de Zaragoza, E. 50009, Zaragoza, Spain
| | - Victor Polo
- Departamento de Química Física and Instituto de Biocomputación y Física de Sistemas Complejos (BIFI), Universidad de Zaragoza, E. 50009, Zaragoza, Spain
| | - Julia Contreras-García
- Sorbonne Université, CNRS, Laboratoire de Chimie Thèorique, LCT, F. 75005 Paris, France.
| |
Collapse
|
8
|
Michalski M, Gordon AJ, Berski S. Topological analysis of the electron localisation function (ELF) applied to the electronic structure of oxaziridine: the nature of N-O bond. Struct Chem 2019. [DOI: 10.1007/s11224-019-01407-9] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
9
|
Shenderovich IG, Denisov GS. Solvent effects on acid-base complexes. What is more important: A macroscopic reaction field or solute-solvent interactions? J Chem Phys 2019; 150:204505. [PMID: 31153188 DOI: 10.1063/1.5096946] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Can the geometry of an acid-base complex in solution be reproduced in calculations using an implicit accounting for the solvent effect in the form of a macroscopic reaction field? The answer is, "Yes, it can." Is this field equal to the real electric field experienced by the complex in solution? The answer is, "No, it is not." How can the geometry be correct under wrong conditions? This question is answered using density functional theory modeling of geometric and NMR parameters of pyridine⋯HF⋯(HCF3)n adducts in the absence and presence of an external electric field. This adduct under field approach shows that the N⋯H distance is a function of the H-F distance whatever method is used to change the geometry of the latter. An explicit account for solute-solvent interactions is required to get a realistic value of the solvent reaction field. Besides that, this approach reveals how certain NMR parameters depend on the solvent reaction field, the solute-solvent interactions, and the geometry of the N⋯H-F hydrogen bond. For some of them, the obtained dependences are far from self-evident.
Collapse
Affiliation(s)
- Ilya G Shenderovich
- Institute of Organic Chemistry, University of Regensburg, Universitaetstrasse 31, 93053 Regensburg, Germany
| | - Gleb S Denisov
- Department of Physics, St.Petersburg State University, Saint Petersburg, Russia
| |
Collapse
|
10
|
|
11
|
Affiliation(s)
- Robert J. Doedens
- Department of Chemistry; University of California; 92697-2025 Irvine CA USA
| | - Philip E. Eaton
- Department of Chemistry; The University of Chicago; 60637 Chicago IL USA
| | | |
Collapse
|
12
|
Jäntschi L, Bálint D, Pruteanu LL, Bolboacă SD. Elemental factorial study on one-cage pentagonal face nanostructure congeners. ACTA ACUST UNITED AC 2016. [DOI: 10.1016/j.md.2016.12.001] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|