Carlos Castillo-Orellana, Esteban Vöhringer-Martinez, Nery Villegas-Escobar
{"title":"Non-covalent interactions and charge transfer in the CO\\(_{2}\\) activation by low-valent group 14 complexes","authors":"Carlos Castillo-Orellana, Esteban Vöhringer-Martinez, Nery Villegas-Escobar","doi":"10.1007/s00894-024-06150-5","DOIUrl":null,"url":null,"abstract":"<div><h3>Context</h3><p>The CO<span>\\(_{2}\\)</span> 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 CO<span>\\(_{2}\\)</span> accompanied by the concomitant formation of the E–O bond between the complex and CO<span>\\(_{2}\\)</span> to lead the corresponding formate product. Prior studies have established that besides the apolar nature of CO<span>\\(_{2}\\)</span>, its initial interaction with the complex is primarily governed by electrostatic interactions. Notably, other stabilizing interactions and the transfer of charge between catalysts and CO<span>\\(_{2}\\)</span> 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 CO<span>\\(_{2}\\)</span> by group 14 main group complex. Our findings indicate that electrostatic interactions predominantly stabilize the complex and CO<span>\\(_{2}\\)</span> 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 CO<span>\\(_{2}\\)</span> approaches the catalysts in contrast to the front-side approach. Notably, it was discovered that a minor initial bending of CO<span>\\(_{2}\\)</span> to approximately <span>\\(176^\\circ \\)</span> 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 CO<span>\\(_{2}\\)</span> while descending in group 14.</p><h3>Methods</h3><p>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 CO<span>\\(_{2}\\)</span> 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.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":null,"pages":null},"PeriodicalIF":2.1000,"publicationDate":"2024-10-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Molecular Modeling","FirstCategoryId":"92","ListUrlMain":"https://link.springer.com/article/10.1007/s00894-024-06150-5","RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q4","JCRName":"BIOCHEMISTRY & MOLECULAR BIOLOGY","Score":null,"Total":0}
引用次数: 0
Abstract
Context
The CO\(_{2}\) 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 CO\(_{2}\) accompanied by the concomitant formation of the E–O bond between the complex and CO\(_{2}\) to lead the corresponding formate product. Prior studies have established that besides the apolar nature of CO\(_{2}\), its initial interaction with the complex is primarily governed by electrostatic interactions. Notably, other stabilizing interactions and the transfer of charge between catalysts and CO\(_{2}\) 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 CO\(_{2}\) by group 14 main group complex. Our findings indicate that electrostatic interactions predominantly stabilize the complex and CO\(_{2}\) 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 CO\(_{2}\) approaches the catalysts in contrast to the front-side approach. Notably, it was discovered that a minor initial bending of CO\(_{2}\) to approximately \(176^\circ \) 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 CO\(_{2}\) 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 CO\(_{2}\) 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.
期刊介绍:
The Journal of Molecular Modeling focuses on "hardcore" modeling, publishing high-quality research and reports. Founded in 1995 as a purely electronic journal, it has adapted its format to include a full-color print edition, and adjusted its aims and scope fit the fast-changing field of molecular modeling, with a particular focus on three-dimensional modeling.
Today, the journal covers all aspects of molecular modeling including life science modeling; materials modeling; new methods; and computational chemistry.
Topics include computer-aided molecular design; rational drug design, de novo ligand design, receptor modeling and docking; cheminformatics, data analysis, visualization and mining; computational medicinal chemistry; homology modeling; simulation of peptides, DNA and other biopolymers; quantitative structure-activity relationships (QSAR) and ADME-modeling; modeling of biological reaction mechanisms; and combined experimental and computational studies in which calculations play a major role.