skip to main content


Title: Weak and strong π interactions between two monomers—assessed with local vibrational mode theory
We introduce in this work a unique parameter for the quantitative assessment of the intrinsic strength of the π interaction between two monomers forming a complex. The new parameter is a local intermonomer stretching force constant, based on the local mode theory, originally developed by Konkoli and Cremer, and derived from the set of nine possible intermonomer normal vibrational modes. The new local force constant was applied to a diverse set of more than 70 molecular complexes, which was divided into four groups. Group 1 includes atoms, ions, and small molecules interacting with benzene and substituted benzenes. Group 2 includes transition metal hydrides and oxides interacting with benzene while Group 3 involves ferrocenes, chromocenes, and titanium sandwich compounds. Group 4 presents an extension to oxygen π–hole interactions in comparison with in-plane hydrogen bonding. We found that the strength of the π interactions in these diverse molecular complexes can vary from weak interactions with predominantly electrostatic character, found, e.g., for argon–benzene complexes, to strong interactions with a substantial covalent nature, found, e.g., for ferrocenes; all being seamlessly described and compared with the new intermonomer local mode force constant, which also outperforms other descriptors such as an averaged force constant or a force constant guided by the electron density bond paths. We hope that our findings will inspire the community to apply the new parameter also to other intermonomer π interactions, enriching in this way the broad field of organometallic chemistry with a new efficient assessment tool.  more » « less
Award ID(s):
2102461
NSF-PAR ID:
10430209
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Canadian Journal of Chemistry
ISSN:
0008-4042
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Density functional theory (DFT) calculations were used to examine the binding strength of one and two methane molecule(s) with graphene (62 and 186 carbon atoms) and model systems of aromatic hydrocarbons (benzene, pyrene, and coronene). We explored different possibilities of binding modes of methane such as one, two, and three C-H interacting with small π-systems. Two methane molecules were considered to bind from the same as well as opposite sides of the plane of benzene and other πsystems including graphene models. Our results show that methane molecule prefers to bind with three C-H…π interactions with all the π-systems except benzene. The preference of tripod configuration of methane on the surface of graphene systems strongly agrees with the neutron diffraction experiment of methane on graphitized carbon black. The binding strength is almost doubled by increasing the number of methane molecules from one to two. Importantly, two methane molecules prefer to bind on the same side rather than opposite sides of the plane of graphene due to stabilizing CH…HC interactions between them in addition to six CH…π interactions. Interestingly, binding strength contributions from CH…HC interactions (approx. 0.4–0.5 kcal/mol) of two methane molecules on the surface are analogous to methane dimer complex free from the surface of graphene. C-H stretching frequency shifts, bond lengths, and binding distances support the presence of CH…HC interactions between two methane molecules. Structures of complexes, binding energies, and C-H stretching frequency shifts agree with available experimental data 
    more » « less
  2. We present a comprehensive investigation on the different role of CO in carboxy- neuroglobin i) as ligand of the heme group in the active site forming a bond with the heme iron and ii) dissociated from the heme group but still trapped inside the active site, focusing on two specific orientations, one with CO perpendicular to the plane defined by the distal histidine of the enzyme (form A) and one with CO located parallel to that plane (form B). Our study includes wild type carboxy-neuroglobin and nine known protein mutations. Considering that the distal histidine interacting with the heme group can adapt two different tautomeric forms and the two possible orientations of the dissociated CO, a total of 36 protein systems were analyzed in this study. Fully optimized geometries and vibrational frequencies were calculated at the QM/MM level, followed by the local mode analysis, to decode CO bond properties. The intrinsic bond strengths derived from the local mode analysis, complemented with NBO and QTAIM data, reveal that the strength of the CO bond, in the hexacoordinate (where CO is a ligand of the heme group) and pentacoordinate (where CO is dissociated from the heme group) scenarios, is dominated by through bond and through space charge transfer between CO and Fe, fine-tuned by electrostatic and dispersion interactions with the side chain amino acids in the distal heme pocket. Suggestions are made as to advise on how protein modifications can influence the molecular properties of the coordinated or dissociated CO, which could serve the fine-tuning of existing and the design of new neuroglobin models with specific FeC and CO bond strengths. 
    more » « less
  3. null (Ed.)
    The exceptionally π-basic metal fragments {MoTp(NO)(DMAP)} and {WTp(NO)(PMe3)} (Tp = tris(pyrazolyl)borate; DMAP = 4-(N,N-dimethylamino)pyridine) form thermally stable η2-coordinated complexes with a variety of electron-deficient arenes. The tolerance of substituted arenes with fluorine-containing electron withdrawing groups (EWG; −F, −CF3, −SF5) is examined for both the molybdenum and tungsten systems. When the EWG contains a π bond (nitriles, aldehydes, ketones, ester), η2 coordination occurs predominantly on the nonaromatic functional group. However, complexation of the tungsten complex with trimethyl orthobenzoate (PhC(OMe)3) followed by hydrolysis allows access to an η2-coordinated arene with an ester substituent. In general, the tungsten system tolerates sulfur-based withdrawing groups well (e.g., PhSO2Ph, MeSO2Ph), and the integration of multiple electron-withdrawing groups on a benzene ring further enhances the π-back-bonding interaction between the metal and aromatic ligand. While the molybdenum system did not form stable η2-arene complexes with the sulfones or ortho esters, it was capable of forming rare examples of stable η2-coordinated arene complexes with a range of fluorinated benzenes (e.g., fluorobenzene, difluorobenzenes). In contrast to what has been observed for the tungsten system, these complexes formed without interference of C–H or C–F insertion. 
    more » « less
  4. In this work, we analyzed five groups of different dihydrogen bonding interactions and hydrogen clusters with an H3+ kernel utilizing the local vibrational mode theory, developed by our group, complemented with the Quantum Theory of Atoms–in–Molecules analysis to assess the strength and nature of the dihydrogen bonds in these systems. We could show that the intrinsic strength of the dihydrogen bonds investigated is primarily related to the protonic bond as opposed to the hydridic bond; thus, this should be the region of focus when designing dihydrogen bonded complexes with a particular strength. We could also show that the popular discussion of the blue/red shifts of dihydrogen bonding based on the normal mode frequencies is hampered from mode–mode coupling and that a blue/red shift discussion based on local mode frequencies is more meaningful. Based on the bond analysis of the H3+(H2)n systems, we conclude that the bond strength in these crystal–like structures makes them interesting for potential hydrogen storage applications. 
    more » « less
  5. High-level quantum chemical computations have provided significant insight into the fundamental physical nature of non-covalent interactions. These studies have focused primarily on gas-phase computations of small van der Waals dimers; however, these interactions frequently take place in complex chemical environments, such as proteins, solutions, or solids. To better understand how the chemical environment affects non-covalent interactions, we have undertaken a quantum chemical study of π– π interactions in an aqueous solution, as exemplified by T-shaped benzene dimers surrounded by 28 or 50 explicit water molecules. We report interaction energies (IEs) using second-order Møller–Plesset perturbation theory, and we apply the intramolecular and functional-group partitioning extensions of symmetry-adapted perturbation theory (ISAPT and F-SAPT, respectively) to analyze how the solvent molecules tune the π– π interactions of the solute. For complexes containing neutral monomers, even 50 explicit waters (constituting a first and partial second solvation shell) change total SAPT IEs between the two solute molecules by only tenths of a kcal mol −1 , while significant changes of up to 3 kcal mol −1 of the electrostatic component are seen for the cationic pyridinium–benzene dimer. This difference between charged and neutral solutes is attributed to large non-additive three-body interactions within solvated ion-containing complexes. Overall, except for charged solutes, our quantum computations indicate that nearby solvent molecules cause very little “tuning” of the direct solute–solute interactions. This indicates that differences in binding energies between the gas phase and solution phase are primarily indirect effects of the competition between solute–solute and solute–solvent interactions. 
    more » « less