1
|
Castillo-Orellana C, Vöhringer-Martinez E, Villegas-Escobar N. Non-covalent interactions and charge transfer in the CO 2 activation by low-valent group 14 complexes. J Mol Model 2024; 30:365. [PMID: 39365341 DOI: 10.1007/s00894-024-06150-5] [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: 06/13/2024] [Accepted: 09/16/2024] [Indexed: 10/05/2024]
Abstract
CONTEXT The CO2 activation by low-valent group 14 catalysts encompasses the rupture of varied covalent bonds in a single transition state through a concerted pathway. The bond between the central main group atom and the hydride in the complex is elongated to trigger the formation of the C-H bond with CO2 accompanied by the concomitant formation of the E-O bond between the complex and CO2 to lead the corresponding formate product. Prior studies have established that besides the apolar nature of CO2 , its initial interaction with the complex is primarily governed by electrostatic interactions. Notably, other stabilizing interactions and the transfer of charge between catalysts and CO2 during the initial phases of the reaction have been ignored. In this study, we have quantified the non-covalent interactions and charge transfer that facilitate the activation of CO2 by group 14 main group complex. Our findings indicate that electrostatic interactions predominantly stabilize the complex and CO2 in the reactant region. However, induction energy becomes the main stabilizing force as the reaction progresses towards the transition state, surpassing electrostatics. Induction contributes about 50% to the stabilization at the transition state, followed by electrostatics (40%) and dispersion interactions (10%). Atomic charges calculated with the minimal basis iterative stockholder (MBIS) method reveal larger charge transfer for the back-side reaction path in which CO2 approaches the catalysts in contrast to the front-side approach. Notably, it was discovered that a minor initial bending of CO2 to approximately 176 ∘ initiates the charge transfer process for all systems. Furthermore, our investigation of group 14 elements demonstrates a systematic reduction in both activation energies and charge transfer to CO2 while descending in group 14. METHODS All studied reactions were characterized along the reaction coordinate obtained with the intrinsic reaction coordinate (IRC) methodology at the M06-2X/6-31 g(d,p) level of theory. Gibbs free energy in toluene was computed using electronic energies at the DLPNO-CCSD(T)/cc-pVTZ-SSD(E) level of theory. Vibrational and translational entropy corrections were applied to provide a more accurate description of the obtained Gibbs free energies. To better characterize the changes in the reaction coordinate for all reactions, the reaction force analysis (RFA) has been employed. It enables the partition of the reaction coordinate into the reactant, transition state, and product regions where different stages of the mechanism occur. A detailed characterization of the main non-covalent driving forces in the initial stages of the activation of CO2 by low-valent group 14 complexes was performed using symmetry-adapted perturbation theory (SAPT). The SAPT0-CT/def2-SVP method was employed for these computations. Charge transfer descriptors based on atomic population using the MBIS scheme were also obtained to complement the SAPT analyses.
Collapse
Affiliation(s)
- Carlos Castillo-Orellana
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Concepción, 4070139, Chile
| | - Esteban Vöhringer-Martinez
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Concepción, 4070139, Chile.
| | - Nery Villegas-Escobar
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Concepción, 4070139, Chile.
| |
Collapse
|
2
|
Villegas-Escobar N. Insights into the variations of kinetic and potential energies in a multi-bond reaction: the reaction electronic flux perspective. J Mol Model 2024; 30:262. [PMID: 38990414 DOI: 10.1007/s00894-024-06024-w] [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: 01/24/2024] [Accepted: 06/11/2024] [Indexed: 07/12/2024]
Abstract
CONTEXT The debate over whether kinetic energy (KE) or potential energy (PE) are the fundamental energy components that contribute to forming covalent bonds has been enduring and stimulating over time. However, the supremacy of these energy components in reactions where multiple bonds are simultaneously formed or broken has yet to be explored. In this study, we use the reaction electronic flux (REF), an effective tool for investigating changes in driving electronic activity when bond formation or dissociation occurs in a chemical reaction, to examine the fluctuations in the KE and PE in a multi-bond reaction. To that end, the activation of CO2 by low-valent group 14 catalysts through a concerted σ -bond metathesis mechanism is analyzed. The findings of this preliminary study suggest that the REF can be utilized as a tool to rationalize alterations in the KE and PE in a multi-bond reaction. Specifically, analyses across the reaction coordinate reveal that changes in the KE and PE precede activation in the REF, stimulating the electronic activity where bond formation or dissociation processes dominate. METHODS The activation of CO2 by the low-valent LEH catalysts (L = N,N'-bis(2,6-diisopropyl phenyl)- β -diketiminate; E = Si, Ge, Sn, and Pb) was studied along the reaction coordinate at the M06-2X/6-31 G(d,p)-LANL2DZ(E) level of theory. The respective minimum energy path calculations were obtained using the intrinsic reaction coordinate (IRC) procedure. The reaction electronic flux (REF) was calculated through the computation of the electronic chemical potential using the frontier molecular orbital approximation. Mayer bond orders along the reaction coordinate have been determined using the NBO 3.1 program in Gaussian16. Most of the reaction coordinate quantities reported in this study (REF, KE, PE, among others) have been determined using the Kudi program and custom Python scripts.
Collapse
Affiliation(s)
- Nery Villegas-Escobar
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Concepción, 4070139, Chile.
| |
Collapse
|
3
|
Bovolenta GM, Silva-Vera G, Bovino S, Molpeceres G, Kästner J, Vogt-Geisse S. In-depth exploration of catalytic sites on amorphous solid water: I. The astrosynthesis of aminomethanol. Phys Chem Chem Phys 2024; 26:18692-18706. [PMID: 38922674 DOI: 10.1039/d4cp01865f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/27/2024]
Abstract
Chemical processes taking place on ice-grain mantles are pivotal to the complex chemistry of interstellar environments. In this study, we conducted a comprehensive analysis of the catalytic effects of an amorphous solid water (ASW) surface on the reaction between ammonia (NH3) and formaldehyde (H2CO) to form aminomethanol (NH2CH2OH) using density functional theory. We identified potential catalytic sites based on the binding energy distribution of NH3 and H2CO reactants, on a set-of-clusters surface model composed of 22 water molecules and found a total of 14 reaction paths. Our results indicate that the catalytic sites can be categorized into four groups, depending on the interactions of the carbonyl oxygen and the amino group with the ice surface in the reactant complex. A detailed analysis of the reaction mechanism using Intrinsic Reaction Coordinate and reaction force analysis, revealed three distinct chemical events for this reaction: formation of the C-N bond, breaking of the N-H bond, and formation of the O-H hydroxyl bond. Depending on the type of catalytic site, these events can occur within a single, concerted, albeit asynchronous, step, or can be isolated in a step-wise mechanism, with the lowest overall transition state energy observed at 1.3 kcal mol-1. A key requirement for the low-energy mechanism is the presence of a pair of dangling OH bonds on the surface, found at 5% of the potential catalytic sites on an ASW porous surface.
Collapse
Affiliation(s)
- Giulia M Bovolenta
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Concepción, Chile.
- Atomistic Simulations, Italian Institute of Technology, 16152 Genova, Italy
| | - Gabriela Silva-Vera
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Concepción, Chile.
| | - Stefano Bovino
- Chemistry Department, Sapienza University of Rome, P.le A. Moro, 00185 Rome, Italy
- INAF - Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125 Firenze, Italy
- Departamento de Astronomía, Facultad de Ciencias Físicas y Matemáticas, Universidad de Concepción, Av. Esteban Iturra s/n Barrio Universitario, Concepción, Chile
| | - German Molpeceres
- Departamento de Astrofísica Molecular Instituto de Física Fundamental (IFF-CSIC), Madrid, Spain
| | - Johannes Kästner
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Stefan Vogt-Geisse
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Concepción, Chile.
| |
Collapse
|
4
|
Ryzhkov FV, Ryzhkova YE, Elinson MN. Python tools for structural tasks in chemistry. Mol Divers 2024:10.1007/s11030-024-10889-7. [PMID: 38744790 DOI: 10.1007/s11030-024-10889-7] [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: 03/14/2024] [Accepted: 04/27/2024] [Indexed: 05/16/2024]
Abstract
In recent decades, the use of computational approaches and artificial intelligence in the scientific environment has become more widespread. In this regard, the popular and versatile programming language Python has attracted considerable attention from scientists in the field of chemistry. It is used to solve a variety of chemical and structural problems, including calculating descriptors, molecular fingerprints, graph construction, and computing chemical reaction networks. Python offers high-quality visualization tools for analyzing chemical spaces and compound libraries. This review is a list of tools for the above tasks, including scripts, libraries, ready-made programs, and web interfaces. Inevitably this manuscript does not claim to be an all-encompassing handbook including all the existing Python-based structural chemistry codes. The review serves as a starting point for scientists wishing to apply automatization or optimization to routine chemistry problems.
Collapse
Affiliation(s)
- Fedor V Ryzhkov
- N. D. Zelinsky Institute of Organic Chemistry Russian Academy of Sciences, 47 Leninsky Prospekt, Moscow, 119991, Russia.
| | - Yuliya E Ryzhkova
- N. D. Zelinsky Institute of Organic Chemistry Russian Academy of Sciences, 47 Leninsky Prospekt, Moscow, 119991, Russia
| | - Michail N Elinson
- N. D. Zelinsky Institute of Organic Chemistry Russian Academy of Sciences, 47 Leninsky Prospekt, Moscow, 119991, Russia
| |
Collapse
|
5
|
Barrales-Martínez C, Gutiérrez-Oliva S, Toro-Labbé A, Pendás ÁM. Interacting Quantum Atoms Analysis of the Reaction Force: A Tool to Analyze Driving and Retarding Forces in Chemical Reactions. Chemphyschem 2021; 22:1976-1988. [PMID: 34293240 DOI: 10.1002/cphc.202100428] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2021] [Revised: 07/21/2021] [Indexed: 12/22/2022]
Abstract
The analysis of the reaction force and its topology has provided a wide range of fruitful concepts in the theory of chemical reactivity over the years, allowing to identify chemically relevant regions along a reaction profile. The reaction force (RF), a projection of the Hellmann-Feynman forces acting on the nuclei of a molecular system onto a suitable reaction coordinate, is partitioned using the interacting quantum atoms approach (IQA). The exact IQA molecular energy decomposition is now shown to open a unique window to identify and quantify the chemical entities that drive or retard a chemical reaction. The RF/IQA coupling offers an extraordinarily detailed view of the type and number of elementary processes that take reactants into products, as tested on two sets of simple reactions.
Collapse
Affiliation(s)
- César Barrales-Martínez
- Departamento de Química Orgánica y Fisicoquímica, Facultad de Ciencias Químicas y Farmacéuticas, Universidad de Chile, Sergio Livingstone 1007, Independencia, Santiago, Chile
| | - Soledad Gutiérrez-Oliva
- Laboratorio de Química Teórica Computacional (QTC), Departamento de Química-Física, Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, Macul, Santiago, Chile
| | - Alejandro Toro-Labbé
- Laboratorio de Química Teórica Computacional (QTC), Departamento de Química-Física, Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, Macul, Santiago, Chile
| | - Ángel Martín Pendás
- Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, 33006, Oviedo, Spain
| |
Collapse
|
6
|
Laloo JZA, Savoo N, Laloo N, Rhyman L, Ramasami P. ExcelAutomat 1.3: Fragment analysis based on the distortion/interaction-activation strain model. J Comput Chem 2018; 40:619-624. [DOI: 10.1002/jcc.25546] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2018] [Revised: 06/27/2018] [Accepted: 07/07/2018] [Indexed: 12/15/2022]
Affiliation(s)
- Jalal Z. A. Laloo
- Computational Chemistry Group, Department of Chemistry, Faculty of Science; University of Mauritius; Réduit 80837 Mauritius
| | - Nandini Savoo
- Computational Chemistry Group, Department of Chemistry, Faculty of Science; University of Mauritius; Réduit 80837 Mauritius
| | - Nassirah Laloo
- School of Innovative Technologies and Engineering, Department of Creative Arts, Film and Media Technologies, University of Technology; Mauritius, La Tour Koenig, Pointe-aux-Sables 11129 Mauritius
| | - Lydia Rhyman
- Computational Chemistry Group, Department of Chemistry, Faculty of Science; University of Mauritius; Réduit 80837 Mauritius
- Department of Applied Chemistry; University of Johannesburg, Doornfontein Campus; Johannesburg 2028 South Africa
| | - Ponnadurai Ramasami
- Computational Chemistry Group, Department of Chemistry, Faculty of Science; University of Mauritius; Réduit 80837 Mauritius
- Department of Applied Chemistry; University of Johannesburg, Doornfontein Campus; Johannesburg 2028 South Africa
| |
Collapse
|
7
|
Villegas-Escobar N, Larsen née Vilhelmsen MH, Gutiérrez-Oliva S, Hashmi ASK, Toro-Labbé A. Double Gold Activation of 1-Ethynyl-2-(Phenylethynyl)Benzene Toward 5-exo
-dig and 6-endo
-dig Cyclization Reactions. Chemistry 2017; 23:13360-13368. [DOI: 10.1002/chem.201701595] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2017] [Indexed: 12/15/2022]
Affiliation(s)
- Nery Villegas-Escobar
- Laboratorio de Química Teórica Computacional (QTC), Facultad de Química; Pontificia Universidad Católica de Chile; Santiago Chile
| | | | - Soledad Gutiérrez-Oliva
- Laboratorio de Química Teórica Computacional (QTC), Facultad de Química; Pontificia Universidad Católica de Chile; Santiago Chile
| | - A. Stephen K. Hashmi
- Institute of Organic Chemistry; Heidelberg University; Heidelberg Germany
- Chemistry Department, Faculty of Science; King Abdulaziz University; Jeddah 21589 Saudi Arabia
| | - Alejandro Toro-Labbé
- Laboratorio de Química Teórica Computacional (QTC), Facultad de Química; Pontificia Universidad Católica de Chile; Santiago Chile
- Freiburg Institute for Advanced Studies (FRIAS); Albert-Ludwigs-Universität Freiburg; Freiburg Germany
| |
Collapse
|
8
|
Vogt-Geisse S, Mata RA, Toro-Labbé A. High level potential energy surface and mechanism of Al(CH3)2OCH3-promoted lactone polymerization: initiation and propagation. Phys Chem Chem Phys 2017; 19:8989-8999. [DOI: 10.1039/c7cp00809k] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A hitherto unreported, second transition state (TS2) is the stationary state with the highest relative energy of the Al(CH3)2OCH3 + glycolide initiation reaction.
Collapse
Affiliation(s)
- Stefan Vogt-Geisse
- Departamento de Físico-Química
- Facultad de Ciencias Químicas
- Universidad de Concepción
- Millenium Nucleus Chemical Processes and Catalysis (CPC)
- Concepción
| | - Ricardo A. Mata
- Institut für Physikalische Chemie
- Georg-August-Universität Göttingen
- D-37077 Göttingen
- Germany
| | - Alejandro Toro-Labbé
- Nucleus Millennium Chemical Processes and Catalysis
- Laboratorio de Química Teórica Computacional (QTC)
- Departamento de Química-Física
- Facultad de Química
- Pontificia Universidad Católica de Chile
| |
Collapse
|