1
|
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
|
2
|
Jabłoński M. Halogen Bond to Experimentally Significant N-Heterocyclic Carbenes (I, IMe 2, I iPr 2, I tBu 2, IPh 2, IMes 2, IDipp 2, IAd 2; I = Imidazol-2-ylidene). Int J Mol Sci 2023; 24:ijms24109057. [PMID: 37240403 DOI: 10.3390/ijms24109057] [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: 04/24/2023] [Revised: 05/15/2023] [Accepted: 05/18/2023] [Indexed: 05/28/2023] Open
Abstract
The subjects of the article are halogen bonds between either XCN or XCCH (X = Cl, Br, I) and the carbene carbon atom in imidazol-2-ylidene (I) or its derivatives (IR2) with experimentally significant and systematically increased R substituents at both nitrogen atoms: methyl = Me, iso-propyl = iPr, tert-butyl = tBu, phenyl = Ph, mesityl = Mes, 2,6-diisopropylphenyl = Dipp, 1-adamantyl = Ad. It is shown that the halogen bond strength increases in the order Cl < Br < I and the XCN molecule forms stronger complexes than XCCH. Of all the carbenes considered, IMes2 forms the strongest and also the shortest halogen bonds with an apogee for complex IMes2⋯ICN for which D0 = 18.71 kcal/mol and dC⋯I = 2.541 Å. In many cases, IDipp2 forms as strong halogen bonds as IMes2. Quite the opposite, although characterized by the greatest nucleophilicity, ItBu2 forms the weakest complexes (and the longest halogen bonds) if X ≠ Cl. While this finding can easily be attributed to the steric hindrance exerted by the highly branched tert-butyl groups, it appears that the presence of the four C-H⋯X hydrogen bonds may also be of importance here. Similar situation occurs in the case of complexes with IAd2.
Collapse
Affiliation(s)
- Mirosław Jabłoński
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, Gagarina 7, 87-100 Torun, Poland
| |
Collapse
|
3
|
Nature of Beryllium, Magnesium, and Zinc Bonds in Carbene⋯MX 2 (M = Be, Mg, Zn; X = H, Br) Dimers Revealed by the IQA, ETS-NOCV and LED Methods. Int J Mol Sci 2022; 23:ijms232314668. [PMID: 36498996 PMCID: PMC9738500 DOI: 10.3390/ijms232314668] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2022] [Revised: 11/21/2022] [Accepted: 11/22/2022] [Indexed: 11/25/2022] Open
Abstract
The nature of beryllium−, magnesium− and zinc−carbene bonds in the cyclopropenylidene⋯MX2 (M = Be, Mg, Zn; X = H, Br) and imidazol-2-ylidene⋯MBr2 dimers is investigated by the joint use of the topological QTAIM-based IQA decomposition scheme, the molecular orbital-based ETS-NOCV charge and energy decomposition method, and the LED energy decomposition approach based on the state-of-the-art DLPNO-CCSD(T) method. All these methods show that the C⋯M bond strengthens according to the following order: Zn < Mg << Be. Electrostatics is proved to be the dominant bond component, whereas the orbital component is far less important. It is shown that QTAIM/IQA underestimates electrostatic contribution for zinc bonds with respect to both ETS-NOCV and LED schemes. The σ carbene→MX2 donation appears to be much more important than the MX2→ carbene back-donation of π symmetry. The substitution of hydrogen atoms by bromine (X in MX2) strengthens the metal−carbene bond in all cases. The physical origin of rotational barriers has been unveiled by the ETS-NOCV approach.
Collapse
|
4
|
Jabłoński M. On the Coexistence of the Carbene⋯H-D Hydrogen Bond and Other Accompanying Interactions in Forty Dimers of N-Heterocyclic-Carbenes (I, IMe 2, I iPr 2, I tBu 2, IMes 2, IDipp 2, IAd 2; I = imidazol-2-ylidene) and Some Fundamental Proton Donors (HF, HCN, H 2O, MeOH, NH 3). Molecules 2022; 27:molecules27175712. [PMID: 36080481 PMCID: PMC9457876 DOI: 10.3390/molecules27175712] [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: 07/26/2022] [Revised: 08/25/2022] [Accepted: 08/29/2022] [Indexed: 11/16/2022] Open
Abstract
The subject of research is forty dimers formed by imidazol-2-ylidene (I) or its derivative (IR2) obtained by replacing the hydrogen atoms in both N-H bonds with larger important and popular substituents of increasing complexity (methyl = Me, iso-propyl = iPr, tert-butyl = tBu, phenyl = Ph, mesityl = Mes, 2,6-diisopropylphenyl = Dipp, 1-adamantyl = Ad) and fundamental proton donor (HD) molecules (HF, HCN, H2O, MeOH, NH3). While the main goal is to characterize the generally dominant C⋯H-D hydrogen bond engaging a carbene carbon atom, an equally important issue is the often omitted analysis of the role of accompanying secondary interactions. Despite the often completely different binding possibilities of the considered carbenes, and especially HD molecules, several general trends are found. Namely, for a given carbene, the dissociation energy values of the IR2⋯HD dimers increase in the following order: NH3< H2O < HCN ≤ MeOH ≪ HF. Importantly, it is found that, for a given HD molecule, IDipp2 forms the strongest dimers. This is attributed to the multiplicity of various interactions accompanying the dominant C⋯H-D hydrogen bond. It is shown that substitution of hydrogen atoms in both N-H bonds of the imidazol-2-ylidene molecule by the investigated groups leads to stronger dimers with HF, HCN, H2O or MeOH. The presented results should contribute to increasing the knowledge about the carbene chemistry and the role of intermolecular interactions, including secondary ones.
Collapse
Affiliation(s)
- Mirosław Jabłoński
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, ul. Gagarina 7, 87-100 Toruń, Poland
| |
Collapse
|
5
|
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
|
6
|
Jabłoński M. Theoretical Study of N-Heterocyclic-Carbene-ZnX 2 (X = H, Me, Et) Complexes. MATERIALS 2021; 14:ma14206147. [PMID: 34683739 PMCID: PMC8539904 DOI: 10.3390/ma14206147] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/23/2021] [Revised: 09/24/2021] [Accepted: 10/12/2021] [Indexed: 11/16/2022]
Abstract
This article discusses the properties of as many as 30 carbene–ZnX2 (X = H, Me, Et) complexes featuring a zinc bond C⋯Zn. The group of carbenes is represented by imidazol-2-ylidene and its nine derivatives (labeled as IR), in which both hydrogen atoms of N-H bonds have been substituted by R groups with various spatial hindrances, from the smallest Me, iPr, tBu through Ph, Tol, and Xyl to the bulkiest Mes, Dipp, and Ad. The main goal is to study the relationship between type and size of R and X and both the strength of C⋯Zn and the torsional angle of the ZnX2 plane with respect to the plane of the imidazol-2-ylidene ring. Despite the considerable diversity of R and X, the range of dC⋯Zn is quite narrow: 2.12–2.20 Å. On the contrary, D0 is characterized by a fairly wide range of 18.5–27.4 kcal/mol. For the smallest carbenes, the ZnX2 molecule is either in the plane of the carbene or is only slightly twisted with respect to it. The twist angle becomes larger and more varied with the bulkier R. However, the value of this angle is not easy to predict because it results not only from the presence of steric effects but also from the possible presence of various interatomic interactions, such as dihydrogen bonds, tetrel bonds, agostic bonds, and hydrogen bonds. It has been shown that at least some of these interactions may have a non-negligible influence on the structure of the IR–ZnX2 complex. This fact should be taken into account in addition to the commonly discussed R⋯X steric repulsion.
Collapse
Affiliation(s)
- Mirosław Jabłoński
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, Gagarina 7, 87-100 Toruń, Poland
| |
Collapse
|
7
|
Jabłoński M. Study of Beryllium, Magnesium, and Spodium Bonds to Carbenes and Carbodiphosphoranes. Molecules 2021; 26:2275. [PMID: 33920004 PMCID: PMC8071025 DOI: 10.3390/molecules26082275] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2021] [Revised: 04/09/2021] [Accepted: 04/11/2021] [Indexed: 11/17/2022] Open
Abstract
The aim of this article is to present results of theoretical study on the properties of C⋯M bonds, where C is either a carbene or carbodiphosphorane carbon atom and M is an acidic center of MX2 (M = Be, Mg, Zn). Due to the rarity of theoretical data regarding the C⋯Zn bond (i.e., the zinc bond), the main focus is placed on comparing the characteristics of this interaction with C⋯Be (beryllium bond) and C⋯Mg (magnesium bond). For this purpose, theoretical studies (ωB97X-D/6-311++G(2df,2p)) have been performed for a large group of dimers formed by MX2 (X = H, F, Cl, Br, Me) and either a carbene ((NH2)2C, imidazol-2-ylidene, imidazolidin-2-ylidene, tetrahydropyrymid-2-ylidene, cyclopropenylidene) or carbodiphosphorane ((PH3)2C, (NH3)2C) molecule. The investigated dimers are characterized by a very strong charge transfer effect from either the carbene or carbodiphosphorane molecule to the MX2 one. This may even be over six times as strong as in the water dimer. According to the QTAIM and NCI method, the zinc bond is not very different than the beryllium bond, with both featuring a significant covalent contribution. However, the zinc bond should be definitely stronger if delocalization index is considered.
Collapse
Affiliation(s)
- Mirosław Jabłoński
- Faculty of Chemistry, Nicolaus Copernicus University, 87-100 Toruń, Poland
| |
Collapse
|
8
|
Ismail TM, Mohan N, Sajith PK. Theoretical study of hydrogen bonding interactions in substituted nitroxide radicals. NEW J CHEM 2021. [DOI: 10.1039/d0nj05362g] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Interaction energy (Eint) of hydrogen bonded complexes of nitroxide radicals can be assessed in terms of the deepest minimum of molecular electrostatic potential (Vmin).
Collapse
Affiliation(s)
| | - Neetha Mohan
- Department of Chemistry
- Christopher Ingold Building
- University College London (UCL)
- London WC1H 0AJ
- UK
| | - P. K. Sajith
- Department of Chemistry
- Farook College
- Kozhikode
- India
| |
Collapse
|
9
|
Jabłoński M, Krygowski TM. Study of the influence of intermolecular interaction on classical and reverse substituent effects in para-substituted phenylboranes. NEW J CHEM 2020. [DOI: 10.1039/d0nj01334j] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The substituent effect and the reverse substituent effect in para-substituted phenylboranes and the influence of the intermolecular interaction of H⋯B type with either silane or methylsilane on the latter of these effects are extensively studied.
Collapse
Affiliation(s)
- Mirosław Jabłoński
- Faculty of Chemistry
- Nicolaus Copernicus University in Toruń
- 87-100 Toruń
- Poland
| | | |
Collapse
|
10
|
|
11
|
Yourdkhani S, Jabłoński M. Physical nature of silane⋯carbene dimers revealed by state-of-the-art ab initio calculations. J Comput Chem 2019; 40:2643-2652. [PMID: 31441520 DOI: 10.1002/jcc.26043] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2019] [Revised: 07/24/2019] [Accepted: 07/25/2019] [Indexed: 11/06/2022]
Abstract
Using the SAPT2 + 3(CCD)δMP2 method in complete basis set (CBS) limit, it is shown that the interactions in the recently studied silane⋯carbene dimers are mainly dispersive in nature. Consequently, slow convergence of dispersion energy also forces slow convergence of the interaction energy. Therefore, obtaining very accurate values requires extrapolation of the correlation part to the CBS limit. The most accurate values obtained at the CCSD(T)/CBS level of theory show that the studied silane⋯carbene dimers are rather weakly bound, with interaction energies ranging from about -1.9 to -1.3 kcal/mol. Comparing to CCSD(T)/CBS, it will be shown that SCS-MP2 and MP2C methods clearly underestimate and methods based on SAPT2+ and having some third-order corrections, as well as the MP2 method, overestimate values of interaction energies. Popular SAPT(DFT) method performs better than SCS-MP2 and MP2C; nevertheless, underestimation is still considerable. The underestimation is slightly quenched if third-order dispersion energy and its exchange counterpart is added to the SAPT(DFT). The closest value of CCSD(T)/CBS has been given by the SAPT2 + (3)(CCD)δMP2 method in quadruple-ζ basis set. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Sirous Yourdkhani
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3,CZ-12116, Prague 2, Czech Republic
| | - Mirosław Jabłoński
- Department of Quantum Chemistry, Faculty of Chemistry, Nicolaus Copernicus University in Toruń, 7-Gagarina St, 87-100, Toruń, Poland
| |
Collapse
|
12
|
Affiliation(s)
- Mirosław Jabłoński
- Department of Quantum Chemistry, Faculty of ChemistryNicolaus Copernicus University in Toruń 7‐Gagarina St. Toruń 87‐100 Poland
| |
Collapse
|