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.     

Influence of metal ion on chelate–aryl stacking interactions

CCSD(T)/CBS and DFT methods are employed to study the stacking interactions of acetylacetonate‐type (acac‐type) chelates of nickel, palladium, and platinum with benzene. The strongest chelate–aryl stacking interactions are formed by nickel and palladium chelate, with interaction energies of −5.75 kcal mol⁻¹ and −5.73 kcal mol⁻¹, while the interaction of platinum chelate is weaker, with interaction energy of −5.36 kcal mol⁻¹. These interaction energies are significantly stronger than stacking of two benzenes, −2.73 kcal mol⁻¹. The strongest nickel and palladium chelate–aryl interactions are with benzene center above the metal area, while the strongest platinum chelate–aryl interaction is with the benzene center above the C2 atom of the acac‐type chelate ring. These preferences arise from very different electrostatic potentials above the metal ions, ranging from very positive above nickel to slightly negative above platinum. While the differences in electrostatic potentials above metal atoms cause different geometries with the most stable interaction among the three metals, the dispersion (correlation energy) component is the largest contribution to the total interaction energy for all three metals.     

Enhancement of Iodine–Hydride Interaction by Substitution and Cooperative Effects in NCX–NCI–HMY Complexes

The NCX‐NCI‐HMY (X=H, Cl, Br, I, Li; M=Be, Mg; Y=H, Li, Na) trimers are investigated to find ways to enhance the iodine–hydride interaction. The interaction energy in the NCI–HMH dimer is ?2.87 and ?5.87 kcal mol^?1 for M=Be and Mg, respectively. When the free H atom in the NCI–HMH dimer is replaced with an alkali atom, the interaction energy is enhanced greatly. When NCX is added into this dimer, the interaction energy of the iodine–hydride interaction is increased by 9–45 % and its increased percentage follows the order X=Cl<Br<H<I<Li and M=Be<Mg. The combination of the alkali substitution and the cooperativity results in a more prominent enhancing effect. The largest interaction energy is found for the NCLi–NCI–HMgLi trimer (?7.03 kcal mol^?1). The influence of the I???H interaction on the X???N interaction is also studied in the trimers. Both types of interactions are analyzed with NBO, AIM, and MEP. The interaction energy in the trimer is also unveiled by a many‐body analysis.     

Ab initio calculations of the Ar–ethane intermolecular potential energy surface using bond function basis sets

The intermolecular potential energy surface (PES) of argon with ethane has been studied by ab initio calculations at the levels of second‐order Møller–Plesset perturbation (MP2) theory and coupled‐cluster theory with single, double, and noniterative triple configurations (CCSD(T)) using a series of augmented correlation‐consistent basis sets. Two sets of bond functions, bf1 (3s3p2d) and bf2 (6s6p4d2f), have been added to the basis sets to show a dramatic and systematic improvement in the convergence of the entire PES. The PES of Ar–ethane is characterized by a global minimum at a near T‐shaped configuration with a well depth of 0.611 kcal mol^?1, a second minimum at a collinear configuration with a well depth of 0.456 kcal mol^?1, and a saddle point connecting the two minima. It is shown that an augmented correlation‐consistent basis set with a set of bond functions, either bf1 or bf2, can effectively produce results equivalent to the next larger augmented correlation‐consistent basis set, that is, aug‐cc‐pVDZ‐bf1 ≈ aug‐cc‐pVTZ, aug‐cc‐pVTZ‐bf1 ≈ aug‐cc‐pVQZ. Very importantly, the use of bond functions improves the PES globally, resulting accurate potential anisotropy. Finally, MP2 method is inadequate for accurate calculations, because it gives a potentially overestimated well depth and, more seriously, a poor potential anisotropy. © 2012 Wiley Periodicals, Inc.     

Hydrogen Bonding in Phosphine Oxide/Phosphate–Phenol Complexes

To develop a new solvent‐impregnated resin (SIR) system for the removal of phenols and thiophenols from water, complex formation by hydrogen bonding of phosphine oxides and phosphates is studied using isothermal titration calorimetry (ITC) and quantum chemical modeling. Six different computational methods are used: B3LYP, M06‐2X, MP2, spin component‐scaled (SCS) MP2 [all four with 6‐311+G(d,p) basis set], a complete basis set extrapolation at the MP2 level (MP2/CBS), and the composite CBS‐Q model. This reveals a range of binding enthalpies (ΔH) for phenol–phosphine oxide and phenol–phosphate complexes and their thio analogues. Both structural (bond lengths/angles) and electronic elements (charges, bond orders) are studied. Furthermore, solvent effects are investigated theoretically by the PCM solvent model and experimentally via ITC. From our calculations, a trialkylphosphine oxide is found to be the most promising extractant for phenol in SIRs, yielding ΔH=?14.5 and ?9.8 kcal mol^?1 with phenol and thiophenol, respectively (MP2/CBS), without dimer formation that would hamper the phenol complexation. In ITC measurements, the ΔH of this complex was most negative in the noncoordinating solvent cyclohexane, and slightly less so in π–π interacting solvents such as benzene. The strongest binding is found for the dimethyl phosphate–phenol complex [?15.1 kcal mol^?1 (MP2/CBS)], due to the formation of two H‐bonds (P?O???H‐O‐ and P‐O‐H???O‐H); however, dimer formation of these phosphates competes with complexation of phenol, and would thus hamper their use in industrial extractions. CBS‐Q calculations display erroneous trends for sulfur compounds, and are found to be unsuitable. Computationally relatively cheap SCS‐MP2 and M06‐2X calculations did accurately agree with the much more elaborate MP2/CBS method, with an average deviation of less than 1 kcal mol^?1.     

Characterization of the potential energy surfaces of two small but challenging noncovalent dimers: (P2)2 and (PCCP)2

This work characterizes eight stationary points of the P₂ dimer and six stationary points of the PCCP dimer, including a newly identified minimum on both potential energy surfaces. Full geometry optimizations and corresponding harmonic vibrational frequencies were computed with the second‐order Møller–Plesset (MP2) electronic structure method and six different basis sets: aug‐cc‐pVXZ, aug‐cc‐pV(X+d)Z, and aug‐cc‐pCVXZ where X = T, Q. A new L‐shaped structure with C₂ symmetry is the only minimum for the P₂ dimer at the MP2 level of theory with these basis sets. The previously reported parallel‐slipped structure with C_2h symmetry and a newly identified cross configuration with D₂ symmetry are the only minima for the PCCP dimer. Single point energies were also computed using the canonical MP2 and CCSD(T) methods as well as the explicitly correlated MP2‐F12 and CCSD(T)‐F12 methods and the aug‐cc‐pVXZ (X = D, T, Q, 5) basis sets. The energetics obtained with the explicitly correlated methods were very similar to the canonical results for the larger basis sets. Extrapolations were performed to estimate the complete basis set (CBS) limit MP2 and CCSD(T) binding energies. MP2 and MP2‐F12 significantly overbind the P₂ and PCCP dimers relative to the CCSD(T) and CCSD(T)‐F12 binding energies by as much as 1.5 kcal mol^?1 for the former and 5.0 kcal mol^?1 for the latter at the CBS limit. The dominant attractive component of the interaction energy for each dimer configuration was dispersion according to several symmetry‐adapted perturbation theory analyses. © 2014 Wiley Periodicals, Inc.     

First Principles (DFT) Characterization of RhI/dppp‐Catalyzed CH Activation by Tandem 1,2‐Addition/1,4‐Rh Shift Reactions of Norbornene to Phenylboronic Acid

The C?H activation in the tandem, “merry‐go‐round”, [(dppp)Rh]‐catalyzed (dppp=1,3‐bis(diphenylphosphino)propane), four‐fold addition of norborene to PhB(OH)₂ has been postulated to occur by a C(alkyl)?H oxidative addition to square‐pyramidal Rh^III?H species, which in turn undergoes a C(aryl)?H reductive elimination. Our DFT calculations confirm the Rh^I/Rh^III mechanism. At the IEFPCM(toluene, 373.15 K)/PBE0/DGDZVP level of theory, the oxidative addition barrier was calculated to be 12.9 kcal mol^?1, and that of reductive elimination was 5.0 kcal mol^?1. The observed selectivity of the reaction correlates well with the relative energy barriers of the cycle steps. The higher barrier (20.9 kcal mol^?1) for norbornyl–Rh protonation ensures that the reaction is steered towards the 1,4‐shift (total barrier of 16.3 kcal mol^?1), acting as an equilibration shuttle. The carborhodation (13.2 kcal mol^?1) proceeds through a lower barrier than the protonation (16.7 kcal mol^?1) of the rearranged aryl–Rh species in the absence of o‐ or m‐substituents, ensuring multiple carborhodations take place. However, for 2,5‐dimethylphenyl, which was used as a model substrate, the barrier for carborhodation is increased to 19.4 kcal mol^?1, explaining the observed termination of the reaction at 1,2,3,4‐tetra(exo‐norborn‐2‐yl)benzene. Finally, calculations with (Z)‐2‐butene gave a carborhodation barrier of 20.2 kcal mol^?1, suggesting that carborhodation of non‐strained, open‐chain substrates would be disfavored relative to protonation.     

Theoretical Study of the Benzene Dimer by the Density‐Functional‐Theory Formalism Based on Electron‐Density Partitioning

The density‐functional approach based on the partition into subsystems was applied to study the benzene dimer. For several structures, the calculated interaction energy and intermolecular distance were compared with the previous theoretical results. A good agreement with high level ab initio correlated methods was found. For instance, the interaction energies obtained in this work and the CCSD(T) method agree within 0.1 – 0.6 kcal/mol depending on the structure of the dimer. The structure with the largest interaction energy is T‐shaped, in agreement with CCSD(T) results. The T‐shaped structure of benzene dimer was suggested by several experimental measurements. The calculated interaction energy of 2.09 kcal/mol agrees also well with experimental estimates based on the dissociation energy which ranges from 1.6±0.2 to 2.4±0.4 kcal/mol and the estimated zero‐point vibration energy of 0.3 – 0.5 kcal/mol.     

Dimers of N‐Heterocyclic Carbene Copper,Silver, and Gold Halides: Probing Metallophilic Interactions through Electron Density Based Concepts

Homobimetallic metallophilic interactions between copper, silver, and gold‐based [(NHC)MX]‐type complexes (NHC=N‐heterocyclic carbene, i.e, 1,3,4‐trimethyl‐4,5‐dihydro‐1H‐1,2,4‐triazol‐5‐ylidene; X=F, Cl, Br, I) were investigated by means of ab initio interaction energies, Ziegler–Rauk‐type energy‐decomposition analysis, the natural orbital for chemical valence (NOCV) framework, and the noncovalent interaction (NCI) index. It was found that the dimers of these complexes predominantly adopt a head‐to‐tail arrangement with typical M ??? M distance of 3.04–3.64 Å, in good agreement with the experimental X‐ray structure determined for [{(NHC)AuCl}₂], which has an Au ??? Au distance of 3.33 Å. The interaction energies between silver‐ and gold‐based monomers are calculated to be about ?25 kcal mol^?1, whereas that for the Cu congener is significantly lower (?19.7 kcal mol^?1). With the inclusion of thermal and solvent contributions, both of which are destabilizing, by about 15 and 8 kcal mol^?1, respectively, an equilibrium process is predicted for the formation of dimer complexes. Energy‐decomposition analysis revealed a dominant electrostatic contribution to the interaction energy, besides significantly stabilizing dispersion and orbital interactions. This electrostatic contribution is rationalized by NHC(δ⁺) ??? halogen(δ^?) interactions between monomers, as demonstrated by electrostatic potentials and derived charges. The dominant NOCV orbital indicates weakening of the π backdonation in the monomers on dimer formation, whereas the second most dominant NOCV represents an electron‐density deformation according to the formation of a very weak M ??? M bond. One of the characteristic signals found in the reduced density gradient versus electron density diagram corresponds to the noncovalent interactions between the metal centers of the monomers in the NCI plots, which is the manifestation of metallophilic interaction.     

Remarks on the analysis of torsional energy surfaces of cycloheptane and cyclooctane by molecular mechanics

A modified version (MM 2′) of the Allinger's 1977 force field is checked against cycloheptane and cyclooctane. Cycloheptane is characterized by two pseudorotating itineraries, chair/twist-chair and boat/twist-boat, separated by a barrier of 8.5 kcal mol^?1. The activation energy in the C/TC pseudorotation is estimated to be 0.96 kcal mol^?1, while B and TB transform into each other freely at an energy level 3.8 kcal mol^?1 above the global energy minimum (TC). With cyclooctane the lowest energy is calculated for the boat-chair form which participates in a pseudorotational process with TBC through a saddle point lying 3.5 kcal mol^?1 above BC. The chair/chair and boat/boat families contain only one local minimum, crown and BB, respectively, on the MM 2′ surface. The results are presented as an illustration for quick coverage of torsional energy surface by two-bond driver calculation with the block-diagonal Newton–Raphson minimization, followed by the force search of stationary points by full-matrix Newton–Raphson optimization.     

A Stable All‐Thiophene‐Based Core‐Modified [38]Octaphyrin Diradicaloid: Conformation and Aromaticity Switch at Different Oxidation States

A soluble and stable core‐modified [38]octaphyrin, MC‐T8 , containing eight thiophene rings was synthesized by Yamamoto coupling followed by oxidative dehydrogenation. X‐ray crystallographic analysis revealed a nearly planar backbone, and the molecule is globally aromatic with a dominant 38π conjugation pathway. The dication MC‐T8²⁺ is antiaromatic, and the backbone is distorted, with a different orientation of the thiophene rings. The tetracation MC‐T8⁴⁺ becomes aromatic again, with a shallow‐bowl‐shaped geometry. Both the neutral compound and the dication demonstrated open‐shell diradical character with a small singlet–triplet energy gap (?2.70 kcal mol^?1 for MC‐T8 and ?3.78 kcal mol^?1 for MC‐T8²⁺ ), and they are stable owing to effective spin delocalization.     

Structures and Dynamic Solution Behavior of Cationic,Two‐Coordinate Gold(I)–π‐Allene Complexes

A family of seven cationic gold complexes that contain both an alkyl substituted π‐allene ligand and an electron‐rich, sterically hindered supporting ligand was isolated in >90 % yield and characterized by spectroscopy and, in three cases, by X‐ray crystallography. Solution‐phase and solid‐state analysis of these complexes established preferential binding of gold to the less substituted C?C bond of the allene and to the allene π face trans to the substituent on the uncomplexed allenyl C?C bond. Kinetic analysis of intermolecular allene exchange established two‐term rate laws of the form rate=k₁[complex]+k₂[complex][allene] consistent with allene‐independent and allene‐dependent exchange pathways with energy barriers of ΔG^≠₁=17.4–18.8 and ΔG^≠₂=15.2–17.6 kcal mol^?1, respectively. Variable temperature (VT) NMR analysis revealed fluxional behavior consistent with facile (ΔG^≠=8.9–11.4 kcal mol^?1) intramolecular exchange of the allene π faces through η¹‐allene transition states and/or intermediates that retain a staggered arrangement of the allene substituents. VT NMR/spin saturation transfer analysis of [{P(tBu)₂o‐binaphthyl}Au(η²‐4,5‐nonadiene) ]⁺SbF₆^? ( 5 ), which contains elements of chirality in both the phosphine and allene ligands, revealed no epimerization of the allene ligand below the threshold for intermolecular allene exchange (ΔG^≠_298K=17.4 kcal mol^?1), which ruled out the participation of a η¹‐allylic cation species in the low‐energy π‐face exchange process for this complex.     

Different catalysis role of in‐loop and out‐of‐loop waters in assisting HNS/HSN proton transfer isomerizations: Bridging vs. surrounding effect

In this work, a density function theory (DFT) study is presented for the HNS/HSN isomerization assisted by 1–4 water molecules on the singlet state potential energy surface (PES). Two modes are considered to model the catalytic effect of these water molecules: (i) water molecule(s) participate directly in forming a proton transfer loop with HNS/HSN species, and (ii) water molecules are out of loop (referred to as out‐of‐loop waters) to assist the proton transfer. In the first mode, for the monohydration mechanism, the heat of reaction is 21.55 kcal · mol^?1 at the B3LYP/6‐311++G** level. The corresponding forward/backward barrier lowerings are obtained as 24.41/24.32 kcal · mol^?1 compared with the no‐water‐assisting isomerization barrier T (65.52/43.87 kcal · mol^?1). But when adding one water molecule on the HNS, there is another special proton‐transfer isomerization pathway with a transition state 10T′ in which the water is out of the proton transfer loop. The corresponding forward/backward barriers are 65.89/65.89 kcal · mol^?1. Clearly, this process is more difficult to follow than the R–T–P process. For the two‐water‐assisting mechanism, the heat of reaction is 19.61 kcal · mol^?1, and the forward/backward barriers are 32.27/12.66 kcal · mol^?1, decreased by 33.25/31.21 kcal · mol^?1 compared with T. For trihydration and tetrahydration, the forward/backward barriers decrease as 32.00/12.60 (30T) and 37.38/17.26 (40T) kcal · mol^?1, and the heat of reaction decreases by 19.39 and 19.23 kcal · mol^?1, compared with T, respectively. But, when four water molecules are involved in the reactant loop, the corresponding energy aspects increase compared with those of the trihydration. The forward/backward barriers are increased by 5.38 and 4.66 kcal · mol^?1 than the trihydration situation. In the second mode, the outer‐sphere water effect from the other water molecules directly H‐bonded to the loop is considered. When one to three water molecules attach to the looped water in one‐water in‐loop‐assisting proton transfer isomerization, their effects on the three energies are small, and the deviations are not more than 3 kcal · mol^?1 compared with the original monohydration‐assisting case. When adding one or two water molecules on the dihydration‐assisting mechanism, and increasing one water molecule on the trihydration, the corresponding energies also are not obviously changed. The results indicate that the forward/backward barriers for the three in‐loop water‐assisting case are the lowest, and the surrounding water molecules (out‐of‐loop) yield only a small effect. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

The Beryllium Pentamer: Trailing an Uneven Sequence of Dissociation Energies

Recent high‐resolution spectroscopic studies by Merritt, Bondybey, and Heaven (Science 2009 , 324, 1548) have heightened the anticipation that small beryllium clusters will soon be observed in the laboratory. Beryllium clusters are important discrete models for the theoretical study of metals. The trigonal bipyramidal Be₅ molecule is studied using high‐level coupled cluster methods. We obtain the optimized geometry, atomization and dissociation energies, and vibrational frequencies. The c～CCSDT(Q) method is employed to compute the atomization and dissociation energies. In this approach, complete basis set (CBS) extrapolations at the CCSD(T) level of theory are combined with an additive correction for the effect of iterative triple and perturbative quadruple excitations. Harmonic vibrational frequencies are obtained using analytic gradients computed at the CCSD(T) level of theory. We report an atomization energy of 129.6 kcal mol^?1 at the trigonal bipyramid global minimum geometry. The Be₅→Be₄+Be dissociation energy is predicted to be 39.5 kcal mol^?1. The analogous dissociation energies for the smaller beryllium clusters are 64.0 kcal mol^?1 (Be₄→Be₃+Be), 24.2 kcal mol^?1 (Be₃→Be₂+Be), and 2.7 kcal mol^?1 (Be₂→Be+Be). The trigonal bipyramidal Be₅ structure has an equatorial–equatorial bond length of 2.000 Å and an axial–equatorial distance of 2.060 Å. Harmonic frequencies of 730, 611, 456, 583, 488, and 338 cm^?1 are obtained at the CCSD(T)/cc‐pCVQZ level of theory. Quadruple excitations are found to make noticeable contributions to the energetics of the pentamer, which exhibits a significant level of static correlation.     

Stabilization of HHeF by Complexation: Is it a Really Viable Strategy?

Ab initio calculations at the MP2 and CCSD(T) levels of theory have disclosed the conceivable existence of fluorine‐coordinated complexes of HHeF with alkali‐metal ions and molecules M⁺ (M⁺=Li⁺–Cs⁺), M⁺–OH₂, M⁺–NH₃ (M⁺=Li⁺, Na⁺), and MX (M=Li, Na; X=F, Cl, Br). All these ligands L induce a shortening of the H? He distance and a lengthening of the He? F distance accompanied by consistent blue‐ and redshifts, respectively, of the H? He and He? F stretching modes. These structural effects are qualitatively similar to those predicted for other investigated complexes of the noble gas hydrides HNgY, but are quantitatively more pronounced. For example, the blueshifts of the H? He stretching mode are exceptionally large, ranging between around 750 and 1000 cm^?1. The interactions of HHeF with the ligands investigated herein also enhance the (HHe)⁺F^? dipole character and produce large complexation energies of around 20–60 kcal mol^?1. Most of the HHeF–L complexes are indeed so stable that the three‐body dissociation of HHeF into H+He+F, exothermic by around 25–30 kcal mol^?1, becomes endothermic. This effect is, however, accompanied by a strong decrease in the H? He? F bending barrier. The complexation energies, ΔE, and the bending barriers, E*, are, in particular, related by the inverse relationship E*(kcal mol^?1)=6.9exp[?0.041ΔE(kcal mol^?1)]. Therefore the HHeF? L complexes, which are definitely stable with respect to H+He+F+L (ΔE≈25–30 kcal mol^?1), are predicted to have bending barriers of only 0.5–2 kcal mol^?1. Overall, our calculations cast doubt on the conceivable stabilization of HHeF by complexation.     

Generating and Trapping Metalla‐N‐Heterocyclic Carbenes

By means of a combined experimental and theoretical approach, the electronic features and chemical behavior of metalla‐N‐heterocyclic carbenes (MNHCs, N‐heterocyclic carbenes containing a metal atom within the heterocyclic skeleton) have been established and compared with those of classical NHCs. MNHCs are strongly basic (proton affinity and pK_a values around 290 kcal mol^?1 and 36, respectively) with a narrow singlet–triplet gap (around 23 kcal mol^?1). MNHCs can be generated from the corresponding metalla‐imidazolium salts and trapped by addition of transition‐metal complexes affording the corresponding heterodimetallic dicarbene derivatives, which can serve as carbene transfer agents.     

Triel–hydride triel bond between ZX3 (Z = B and Al; X = H and Me) and THMe3 (T = Si,Ge and Sn)

Complexes between THMe₃ (T = Si, Ge and Sn) and ZX₃ (Z = B and Al; X = H and Me) have been characterized using MP2/aug‐cc‐pVTZ calculations. These complexes are chiefly stabilized by a triel–hydride triel bond with the T–H bond pointing to the π‐hole on the triel atom. The triel–hydride interaction is mainly attributed to the charge transfer from the T–H bond orbital to the empty p orbital of the triel atom. These complexes are very stable with a large interaction energy (>10 kcal mol^?1) excluding THMe₃···BMe₃ (T = Si and Ge), indicating that the sp²‐hydridized triel atom has a strong affinity for the T–H bond. The formation of THMe₃···BH₃ results in proton transfer, characterized by conversion of orbital interaction and large charge transfer (ca 0.5e). The large deformation is primarily responsible for the abnormally greater interaction energy in THMe₃···BH₃ (>30 kcal mol^?1) than in the AlH₃ analogue. Methyl substitution on the triel atom weakens the triel–hydride interaction and causes a larger interaction energy in THMe₃···AlMe₃ with respect to its BMe₃ counterpart. Most of these interactions possess characteristics of covalent bonds. Polarization makes a contribution to the stability of most complexes nearly equivalent to the electrostatic term.     

Probing the Chemical Functionalization of Single‐Walled Carbon Nanotubes with Multiple Carbon Ad‐Dimer Defects

Drying‐tube‐shaped single‐walled carbon nanotubes (SWCNTs) with multiple carbon ad‐dimer (CD) defects are obtained from armchair (n,n,m) SWCNTs (n=4, 5, 6, 7, 8; m=7, 13). According to the isolated‐pentagon rule (IPR) the drying‐tube‐shaped SWCNTs are unstable non‐IPR species, and their hydrogenated, fluorinated, and chlorinated derivatives are investigated. Interestingly, chemisorptions of hydrogen, fluorine, and chlorine atoms on the drying tube‐shaped SWCNTs are exothermic processes. Compared to the reaction energies for binding of H, F, and Cl atoms to perfect and Stone–Wales‐defective armchair (5,5) nanotubes, binding of F with the multiply CD defective SWCNTs is stronger than with perfect and Stone–Wales‐defective nanotubes. The reaction energy for per F₂ addition is between 85 and 88 kcal mol^?1 more negative than that per H₂ addition. Electronic structure analysis of their energy gaps shows that the CD defects have a tendency to decrease the energy gap from 1.98–2.52 to 0.80–1.17 eV. After hydrogenation, fluorination, and chlorination, the energy gaps of the drying‐tube‐shaped SWCNTs with multiple CD defects are substantially increased to 1.65–3.85 eV. Furthermore, analyses of thermodynamic stability and nucleus‐independent chemical shifts (NICS) are performed to analyze the stability of these molecules.     

Binding Energies of the π‐Stacked Anisole Dimer: New Molecular Beam—Laser Spectroscopy Experiments and CCSD(T) Calculations

Among noncovalent interactions, π–π stacking is a very important binding motif governed mainly by London dispersion. Despite its importance, for instance, for the structure of bio‐macromolecules, the direct experimental measurement of binding energies in π–π stacked complexes has been elusive for a long time. Only recently, an experimental value for the binding energy of the anisole dimer was presented, determined by velocity mapping ion imaging in a two‐photon resonant ionisation molecular beam experiment. However, in that paper, a discrepancy was already noted between the obtained experimental value and a theoretical estimate. Here, we present an accurate recalculation of the binding energy based on the combination of the CCSD(T)/CBS interaction energy and a DFT‐D3 vibrational analysis. This proves unambiguously that the previously reported experimental value is too high and a new series of measurements with a different, more sensitive apparatus was performed. The new experimental value of 1800±100 cm^?1 (5.15±0.29 kcal mol^?1) is close to the present theoretical prediction of 5.04±0.40 kcal mol^?1. Additional calculations of the properties of the cationic and excited states involved in the photodissociation of the dimer were used to identify and rationalise the difficulties encountered in the experimental work.     

Modulation of Stacking Interactions by Transition‐Metal Coordination: Ab Initio Benchmark Studies

A series of ab initio calculations are used to determine the C? H???π and π???π‐stacking interactions of aromatic rings coordinated to transition‐metal centres. Two model complexes have been employed, namely, ferrocene and chromium benzene tricarbonyl. Benchmark data obtained from extrapolation of MP2 energies to the basis set limit, coupled with CCSD(T) correction, indicate that coordinated aromatic rings are slightly weaker hydrogen‐bond acceptors but are significantly stronger hydrogen‐bond donors than uncomplexed rings. It is found that π???π stacking to a second benzene is stronger than in the free benzene dimer, especially in the chromium case. This is assigned, by use of energy partitioning in the local correlation method, to dispersion interactions between metal d and benzene π orbitals. The benchmark data is also used to test the performance of more efficient theoretical methods, indicating that spin‐component scaling of MP2 energies performs well in all cases, whereas various density functionals describe some complexes well, but others with errors of more than 1 kcal mol^?1.