Complexes between lithium cation and diphenylalkanes in the gas phase: the pincer effect

The gas-phase lithium cation basicities (LCB values, Gibbs free energies of binding) of alpha,omega-diphenylalkanes Ph-(CH(2))(n)-Ph (n=2, 3, or 7) and 1,1-diphenylethane Ph-CH(Me)-Ph were investigated by means of Fourier-transform ion cyclotron resonance (FTICR) mass spectrometry. Their structures, and those of the corresponding Li(+) complexes were optimized at the B3LYP/6-31G(d) level and their relative stabilities calculated at the B3LYP/6-311+G(3df,2p)//B3LYP/6-31G(d) level. Whereas the most stable conformers of the free diphenylalkanes were found to adopt a completely stretched aliphatic chain connecting the two benzene rings, the most stable Li(+) complexes correspond to conformers in which the alkali metal cation interacts simultaneously with both benzene rings through the folding of the aliphatic chain ("pincer effect"). This chelation brings about a significant enhancement of the Li(+) binding enthalpies (LBE values), which were calculated to be approximately 75 kJ mol(-1) higher than those evaluated for conventional (singly coordinated) pi complexes in which the metal cation interacts with only one of the benzene rings. The increase of the corresponding lithium cation basicities, however, (Gibbs free energies of Li(+) binding, LCB values) was calculated to be smaller by approximately 15 kJ mol(-1) as the pincer effect is entropically disfavored. The good agreement between the calculated LCB values, assuming a statistical distribution of the different conformers present in the gas phase, and the experimental LCB values measured by means of FTICR mass spectrometry are considered indirect evidence of the existence of the pincer effect.     

Halogen‐Bonding Interactions with π Systems: CCSD(T), MP2, and DFT Calculations

Halogen bonding is a noncovalent interaction between a halogen atom and a nucleophilic site. Interactions involving the π electrons of aromatic rings have received, up to now, little attention, despite the large number of systems in which they are present. We report binding energies of the interaction between either NCX or PhX (X=F, Cl, Br, I) and the aromatic benzene system as determined with the coupled cluster with perturbative triple excitations method [CCSD(T)] extrapolated at the complete basis set limit. Results are compared with those obtained by Møller–Plesset perturbation theory to second order (MP2) and density functional theory (DFT) calculations by using some of the most common functionals. Results show the important role of DFT in studying this interaction.     

Competition Between π–π and CH/π Interactions: A Comparison of the Structural and Electronic Properties of Alkoxy‐Substituted 1,8‐Bis((propyloxyphenyl)ethynyl)naphthalenes

The structural and electronic consequences of π–π and C?H/π interactions in two alkoxy‐substituted 1,8‐bis‐ ((propyloxyphenyl)ethynyl)naphthalenes are explored by using X‐ray crystallography and electronic structure computations. The crystal structure of analogue 4 , bearing an alkoxy side chain in the 4‐position of each of the phenyl rings, adopts a π‐stacked geometry, whereas analogue 8 , bearing alkoxy groups at both the 2‐ and the 5‐positions of each ring, has a geometry in which the rings are splayed away from a π‐stacked arrangement. Symmetry‐adapted perturbation theory analysis was performed on the two analogues to evaluate the interactions between the phenylethynyl arms in each molecule in terms of electrostatic, steric, polarization, and London dispersion components. The computations support the expectation that the π‐stacked geometry of the alkoxyphenyl units in 4 is simply a consequence of maximizing π–π interactions. However, the splayed geometry of 8 results from a more subtle competition between different noncovalent interactions: this geometry provides a favorable anti‐alignment of C?O bond dipoles, and two C?H/π interactions in which hydrogen atoms of the alkyl side chains interact favorably with the π electrons of the other phenyl ring. These favorable interactions overcome competing π–π interactions to give rise to a geometry in which the phenylethynyl substituents are in an offset, unstacked arrangement.     

Enhanced Aerogen–π Interaction by a Cation–π Force

The interaction between a noble gas atom and an aromatic π‐electron system, which mainly originates from the London dispersion force, is very weak and has not attracted enough attention yet. Herein, we reported a type of notably enhanced aerogen–π interaction between cation–π systems and noble gas atoms. The binding strength of a divalent cation–π system with a xenon atom is comparable to a moderate hydrogen bond (up to ca. 7 kcal mol^?1), whereas krypton and argon atoms produce slightly weaker interactions. Energy‐decomposition analysis reveals that the induction interaction is responsible for the stabilization of divalent cation–π?Xe species besides the dispersion interaction. Our results might be helpful to increase the understanding of some unsolved mysteries of aerogens.     

Influence of the Li⋅⋅⋅π Interaction on the H/X⋅⋅⋅π Interactions in HOLi⋅⋅⋅C6H6⋅⋅⋅HOX/XOH (X=F,Cl, Br,I) Complexes

The influences of the Li???π interaction of C₆H₆???LiOH on the H???π interaction of C₆H₆???HOX (X=F, Cl, Br, I) and the X???π interaction of C₆H₆???XOH (X=Cl, Br, I) are investigated by means of full electronic second‐order Møller–Plesset perturbation theory calculations and “quantum theory of atoms in molecules” (QTAIM) studies. The binding energies, binding distances, infrared vibrational frequencies, and electron densities at the bond critical points (BCPs) of the hydrogen bonds and halogen bonds prove that the addition of the Li???π interaction to benzene weakens the H???π and X???π interactions. The influences of the Li???π interaction on H???π interactions are greater than those on X???π interactions; the influences of the H???π interactions on the Li???π interaction are greater than X???π interactions on Li???π interaction. The greater the influence of Li???π interaction on H/X???π interactions, the greater the influences of H/X???π interactions on Li???π interaction. QTAIM studies show that the intermolecular interactions of C₆H₆???HOX and C₆H₆???XOH are mainly of the π type. The electron densities at the BCPs of hydrogen bonds and halogen bonds decrease on going from bimolecular complexes to termolecular complexes, and the π‐electron densities at the BCPs show the same pattern. Natural bond orbital analyses show that the Li???π interaction reduces electron transfer from C₆H₆ to HOX and XOH.     

Orbital Effect‐Induced Anomalous Anion–π Interactions between Electron‐Rich Aromatic Hydrocarbons and Fluoride

Anion–π interactions generally exist between an anion and an electron‐deficient π‐ring because of the electron‐accepting character of the ring. In this paper, we report orbital effect‐induced anomalous binding between electron‐rich π systems and F^? through anion–π interactions calculated at the MP2/6‐31+G(d,p) and ωB97X‐D/6‐31+G(d,p) levels of theory. We find that anion–π interactions between F^? and electron‐rich π rings increase markedly with increasing number of π electrons and size of the π rings. This is contrary to intuition because anion–π interactions would be expected to gradually decrease because of gradually increasing Coulombic repulsion between the negative charge of the anions and gradually increasing number of π electrons of the aromatic surfaces. Energy decomposition analysis showed that the key to this anomalous effect is the more effective delocalization of negative charge to the unoccupied π* orbitals of larger π rings, which is termed an “orbital effect”.     

Rational Design of a Nonbasic Molecular Receptor for Selective NH4+/K+ Complexation in the Gas Phase

A new molecular receptor ( 1 ) for ammonium recognition has been designed and constructed by using only carbon atoms. This molecular receptor can co‐exist in two different isoenergetic conformations but, upon complexation, the conformers are no longer isoenergetic, and a basket‐shaped conformation becomes clearly more stable. The pre‐organised tetrahedral structure of this basket‐shaped molecule favours the complexation of ammonium ions by N? H???π interactions with the four phenyl groups of the host. A similar behaviour is not observed in a similar, but less pre‐organised, reference molecule. ESI‐MS competition experiments show that 1 is able to bind NH₄⁺ over K⁺ selectively. This is the first example of a neutral molecular receptor that shows a remarkable NH₄⁺/K⁺ selectivity. DFT‐calculations provide insight into the nature of host–guest interactions of both 1? NH₄⁺ and 1? K⁺ complexes as well as in the mechanism involved in multiple cation–π interactions and the influence of these interactions on the conformational stability and the selective binding of the host.     

The Nature of the [TTF].+⋅⋅⋅[TTF].+ Interactions in the [TTF]22+ Dimers Embedded in Charged [3]Catenanes: Room‐Temperature Multicenter Long Bonds

The properties of tetrathiafulvalene dimers ([TTF]₂²⁺) and the functionalized ring‐shaped bispropargyl (BPP)‐functionalized TTF dimers, [BPP–TTF]₂²⁺, found at room temperature in charged [3]catenanes, were evaluated by M06L calculations. The results showed that their isolated [TTF]₂²⁺ and [BPP–TTF]₂²⁺ dimers are energetically unstable towards dissociation. When enclosed in the 4⁺‐charged central cyclophane ring of charged [3]catenanes (CBPQT⁴⁺), [TTF]₂²⁺ and [BPP–TTF]₂²⁺ dimers are also energetically unstable with respect to leaving the CBPQT⁴⁺ ring; since the barrier for the exiting process is only about 3 kcal mol^?1, that is, within the reach of thermal energies at room temperature (neutral [TTF]₂⁰ dimers are stable within the CBPQT⁴⁺ ring). However, the [BPP–TTF]₂²⁺ dimers in charged [3]catenanes cannot exit, because this would imply breaking the covalent bonds of the BPP–TTF⁺ macrocycle. Finally, it was shown that the [TTF]₂²⁺, [BPP–TTF]₂²⁺ dimers, and charged [3]catenanes are energetically stable in solution and in crystals of their salts, in the first case due to the interactions with the solvent, and in the second case mostly due to cation–anion interactions. In these environmental conditions at room temperature the TTF units of the [BPP–TTF]₂²⁺ dimers make short contacts, thus allowing their SOMO orbitals to overlap: a room‐temperature multicenter long bond is formed, similar to those previously found in other [TTF]₂²⁺ salts and their solutions.     

Understanding the Fundamental Role of π/π, σ/σ, and σ/π Dispersion Interactions in Shaping Carbon‐Based Materials

Noncovalent interactions involving aromatic rings, such as π‐stacking and CH/π interactions, are central to many areas of modern chemistry. However, recent studies proved that aromaticity is not required for stacking interactions, since similar interaction energies were computed for several aromatic and aliphatic dimers. Herein, the nature and origin of π/π, σ/σ, and σ/π dispersion interactions has been investigated by using dispersion‐corrected density functional theory, energy decomposition analysis, and the recently developed noncovalent interaction (NCI) method. Our analysis shows that π/π and σ/σ stacking interactions are equally important for the benzene and cyclohexane dimers, explaining why both compounds have similar boiling points. Also, similar dispersion forces are found in the benzene???methane and cyclohexane???methane complexes. However, for systems larger than naphthalene, there are enhanced stacking interactions in the aromatic dimers adopting a parallel‐displaced configuration compared to the analogous saturated systems. Although dispersion plays a decisive role in stabilizing all the complexes, the origin of the π/π, σ/σ, and σ/π interactions is different. The NCI method reveals that the dispersion interactions between the hydrogen atoms are responsible for the surprisingly strong aliphatic interactions. Moreover, whereas σ/σ and σ/π interactions are local, the π/π stacking are inherently delocalized, which give rise to a non‐additive effect. These new types of dispersion interactions between saturated groups can be exploited in the rational design of novel carbon materials.     

Accurate binding energies of hydrogen,halogen, and dihydrogen bonded complexes and cation enhanced binding strengths

Interaction energy (E_int) values of a variety of hydrogen, halogen, and dihydrogen bonded complexes in the weak, medium, and strong regimes have been computed using W1BD, MP2, M06L density functional theory, and hybrid methods MP4//MP2, MP4//M06L, and CCSD(T)//MP2. W1BD level E_int and CCSD(T) results reported in the literature show very good agreement (mean absolute deviation = 0.19 kcal/mol). MP2 underestimates E_int while M06L shows accurate behavior for all except halogen and charge‐assisted hydrogen bonds. MP4//MP2, MP4//M06L, and CCSD(T)//MP2 yield E_int very close to those obtained from W1BD. The high accuracy energy data at MP4/MP2 is used to study the effect of a cation (Li⁺, NH₄⁺) on the E_int. The cation enhances electron donation from the donor to noncovalent bonding region leading to substantial enhancement in E_int (～141–566% for Li⁺ and ～105–539% for NH₄⁺) and promotes a noncovalent bond in the weak regime to medium regime and that in the medium regime to strong regime. © 2014 Wiley Periodicals, Inc.     

Intermolecular Interactions in Weak Donor–Acceptor Complexes from Symmetry‐Adapted Perturbation and Coupled‐Cluster Theory: Tetracyanoethylene–Benzene and Tetracyanoethylene–p‐Xylene

The interactions in the complexes of tetracyanothylene (TCNE) with benzene and p‐xylene, often classified as weak electron donor–acceptor (EDA) complexes, are investigated by a range of quantum chemical methods including intermolecular perturbation theory at the DFT‐SAPT (symmetry‐adapted perturbation theory combined with density functional theory) level and explicitly correlated coupled‐cluster theory at the CCSD(T)‐F12 level. The DFT‐SAPT interaction energies for TCNE–benzene and TCNE–p‐xylene are estimated to be ?35.7 and ?44.9 kJ mol^?1, respectively, at the complete basis set limit. The best estimates for the CCSD(T) interaction energy are ?37.5 and ?46.0 kJ mol^?1, respectively. It is shown that the second‐order dispersion term provides the most important attractive contribution to the interaction energy, followed by the first‐order electrostatic term. The sum of second‐ and higher‐order induction and exchange–induction energies is found to provide nearly 40 % of the total interaction energy. After addition of vibrational, rigid‐rotor, and translational contributions, the computed internal energy changes on complex formation approach results from gas‐phase spectrophotometry at elevated temperatures within experimental uncertainties, while the corresponding entropy changes differ substantially.     

(Isocyano Group π-Hole)⋅⋅⋅[d -MII] Interactions of (Isocyanide)[MII] Complexes,in which Positively Charged Metal Centers (d8-M=Pt,Pd) Act as Nucleophiles

Inspection of the X-ray structures of the newly prepared trans-[M^II(CNXyl)₂(DAPT)₂]Cl(BF₄) (M=Pd, Pt; Xyl=2,6-Me₂C₆H₃; DAPT=4,6-diaminopyrimidine-2(1H)-thione) complexes and the appropriate Hirshfeld molecular surface analysis allowed the recognition of the previously unknown π-hole⋅⋅⋅metal interactions between a ligated isocyano group (acting as a π-hole donor) and the positively charged d⁸-Pt^II and d⁸-Pd^II metal centers (acting as nucleophiles); this is the first identification of π-hole⋅⋅⋅metal interactions with triple-bond species. Results of DFT calculations followed by the topological analysis of the electron density distribution within the framework of Bader's theory (quantum theory of atoms in molecules, QTAIM) confirmed the presence of these contacts. The electrostatic surface potential calculations indicated that π-hole⋅⋅⋅metal contacts are formed upon interaction between the electrophilic isocyano C atom (π-hole donor) and the nucleophilic d orbital of the metal centers, which act as π-hole acceptors. Available CCDC data were processed from the perspective of the π-hole⋅⋅⋅metal interactions with isocyanide ligands, and their analysis disclosed the role of metal nucleophilicity in the corresponding π-hole acceptor ability.     

DFT Mechanistic Study on Diene Metathesis Catalyzed by Ru‐Based Grubbs–Hoveyda‐Type Carbenes: The Key Role of π‐Electron Density Delocalization in the Hoveyda Ligand

The catalytic activity and catalyst recovery of two heterogenized ruthenium‐based precatalysts ( H and NO₂(4) ) in diene ring‐closing metathesis have been studied by means of density functional calculations at the B3LYP level of theory. For comparison and rationalization of the key factors that lead to higher activities and higher catalyst recoveries, four other Grubbs–Hoveyda complexes have also been investigated. The full catalytic cycle (catalyst formation, propagation, and precatalyst regeneration) has been considered. DFT calculations suggest that either for the homogeneous and heterogenized systems the activity of the catalysts mainly depends on the ability of the precursor to generate the propagating carbene. This ability does not correlate with the traditionally identified key factor, the Ru???O interaction strength. In contrast, precatalysts with lower alkoxy‐dissociation energy barriers and lower stabilities compared with the propagating carbene also present larger C1? C2 bond length (i.e., lower π character of the C? C bond that exists between the metal–carbene (Ru?C) and the phenyl ring of the Hoveyda ligand). Catalyst recovery, regardless of whether a release–return mechanism occurs or not, is also mainly determined by the π delocalization. Therefore, future Grubbs–Hoveyda‐type catalyst development should be based on fine‐tuning the π‐electron density of the phenyl moiety, with the subsequent effect on the metalloaromaticity of the ruthenafurane ring, rather than considering the modification of the Ru???O interaction.     

Dual Role of the 1,2,3‐Triazolium Ring as a Hydrogen‐Bond Donor and Anion–π Receptor in Anion‐Recognition Processes

Several bis(triazolium)‐based receptors have been synthesized as chemosensors for anion recognition. The central naphthalene core features two aryltriazolium side‐arms. NMR experiments revealed differences between the binding modes of the two triazolium rings: one triazolium ring acts as a hydrogen‐bond donor, the other as an anion–π receptor. Receptors 9²⁺?2BF₄ ^? (C₆H₅), 11²⁺?2BF₄ ^? (4‐NO₂?C₆H₄), and 13²⁺?2BF^4? (ferrocenyl) bind HP₂O₇^3? anions in a mixed‐binding mode that features a combination of hydrogen‐bonding and anion–π interactions and results in strong binding. On the other hand, receptor 10²⁺?2 BF₄ ^? (4‐CH₃O?C₆H₄) only displays combined Csp₂?H/anion–π interactions between the two arms of the receptors and the bound anion rather than triazolium (CH)⁺???anion hydrogen bonding. All receptors undergo a downfield shift of the triazolium protons, as well as the inner naphthalene protons, in the presence of H₂PO₄^? anions. That suggests that only hydrogen‐bonding interactions exist between the binding site and the bound anion, and involve a combination of cationic (triazolium) and neutral (naphthalene) C?H donor interactions. Theoretical calculations relate the electronic structure of the substituent on the aromatic group with the interaction energies and provide a minimum‐energy conformation for all the complexes that explains their measured properties.     

Theoretical Study on the Complexes of Benzene with Isoelectronic Nitrogen‐Containing Heterocycles

The π–π interactions between benzene and the aromatic nitrogen heterocycles pyridine, pyrimidine, 1,3,5‐triazine, 1,2,3‐triazine, 1,2,4,5‐tetrazine, and 1,2,3,4,5‐pentazine are systematically investigated. The T‐shaped structures of all complexes studied exhibit a contraction of the C? H bond accompanied by a rather large blue shift (40–52 cm^?1) of its stretching frequency, and they are almost isoenergetic with the corresponding displaced‐parallel structures at reliable levels of theory. With increasing number of nitrogen atoms in the heterocycle, the geometries, frequencies, energies, percentage of s character at C, and the electron density in the C? H σ antibonding orbital of the complexes all increase or decrease systematically. Decomposition analysis of the total binding energy showed that for all the complexes, the dispersion energy is the dominant attractive contribution, and a rather large attraction originating from electrostatic contribution is compensated by its exchange counterpart.