{"title":"Cooperativity and halonium transfer in the ternary NCI···CH3I···−CN halogen-bonded complex: An ab initio gas phase study","authors":"Rubén D. Parra","doi":"10.1007/s00894-024-06160-3","DOIUrl":null,"url":null,"abstract":"<div><h3>Context</h3><p>The strength and nature of the two halogen bonds in the NCI···CH<sub>3</sub>I···<sup>−</sup>CN halogen-bonded ternary complex are studied in the gas phase via ab initio calculations. Different indicators of halogen bond strength were employed to examine the interactions including geometries, complexation energies, Natural Bond Order (NBO) Wiberg bond indices, and Atoms in Molecules (AIM)-based charge density topological properties. The results show that the halogen bond is strong and partly covalent in nature when CH<sub>3</sub>I donates the halogen bond, but weak and noncovalent in nature when CH<sub>3</sub>I accepts the halogen bond. Significant halogen bond cooperativity emerges in the ternary complex relative to the corresponding heterodimer complexes, NCI···CH<sub>3</sub>I and CH<sub>3</sub>I···<sup>−</sup>CN, respectively. For example, the CCSD(T) complexation energy of the ternary complex (-18.27 kcal/mol) is about twice the sum of the complexation energies of the component dimers (-9.54 kcal/mol). The halonium transfer reaction that converts the ternary complex into an equivalent one was also investigated. The electronic barrier for the halonium transfer was calculated to be 6.70 kcal/mol at the CCSD(T) level. Although the MP2 level underestimates and the MP3 overestimates the barrier, their calculated MP2.5 average barrier (6.44 kcal/mol) is close to that of the more robust CCSD(T) level. Insights on the halonium ion transfer reaction was obtained by examining the reaction energy and force profiles along the intrinsic reaction coordinate, IRC. The corresponding evolution of other properties such as bond lengths, Wiberg bond indices, and Mulliken charges provides specific insight on the extent of structural rearrangements and electronic redistribution throughout the entire IRC space.</p><h3>Methods</h3><p>The MP2 method was used for geometry optimizations. Energy calculations were performed using the CCSD(T) method. The aug-cc-pVTZ basis set was employed for all atoms other than iodine for which the aug-cc-pVTZ-PP basis set was used instead.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":null,"pages":null},"PeriodicalIF":2.1000,"publicationDate":"2024-10-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s00894-024-06160-3.pdf","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-06160-3","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 strength and nature of the two halogen bonds in the NCI···CH3I···−CN halogen-bonded ternary complex are studied in the gas phase via ab initio calculations. Different indicators of halogen bond strength were employed to examine the interactions including geometries, complexation energies, Natural Bond Order (NBO) Wiberg bond indices, and Atoms in Molecules (AIM)-based charge density topological properties. The results show that the halogen bond is strong and partly covalent in nature when CH3I donates the halogen bond, but weak and noncovalent in nature when CH3I accepts the halogen bond. Significant halogen bond cooperativity emerges in the ternary complex relative to the corresponding heterodimer complexes, NCI···CH3I and CH3I···−CN, respectively. For example, the CCSD(T) complexation energy of the ternary complex (-18.27 kcal/mol) is about twice the sum of the complexation energies of the component dimers (-9.54 kcal/mol). The halonium transfer reaction that converts the ternary complex into an equivalent one was also investigated. The electronic barrier for the halonium transfer was calculated to be 6.70 kcal/mol at the CCSD(T) level. Although the MP2 level underestimates and the MP3 overestimates the barrier, their calculated MP2.5 average barrier (6.44 kcal/mol) is close to that of the more robust CCSD(T) level. Insights on the halonium ion transfer reaction was obtained by examining the reaction energy and force profiles along the intrinsic reaction coordinate, IRC. The corresponding evolution of other properties such as bond lengths, Wiberg bond indices, and Mulliken charges provides specific insight on the extent of structural rearrangements and electronic redistribution throughout the entire IRC space.
Methods
The MP2 method was used for geometry optimizations. Energy calculations were performed using the CCSD(T) method. The aug-cc-pVTZ basis set was employed for all atoms other than iodine for which the aug-cc-pVTZ-PP basis set was used instead.
期刊介绍:
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.