Reliable structures and energetics for two new delocalized pi...pi prototypes: cyanogen dimer and diacetylene dimer

Two new prototype delocalized pi[dot dot dot]pi complexes are introduced: the dimers of cyanogen, (N[triple bond]C-C[triple bond]N)(2), and diacetylene, (HC[triple bond]C-C[triple bond]CH)(2). These dimers have properties similar to larger delocalized pi...pi systems such as benzene dimer but are small enough that they can be probed in far greater detail with high accuracy electronic structure methods. Parallel-slipped and T-shaped structures of both cyanogen dimer and diacetylene dimer have been optimized with 15 different procedures. The effects of basis set size, theoretical method, counterpoise correction, and the rigid monomer approximation on the structure and energetics of each dimer have been examined. MP2 and CCSD(T) optimized geometries for all four dimer structures are reported, as well as estimates of the CCSD(T) complete basis set (CBS) interaction energy for every optimized geometry. The data reported here suggest that future optimizations of delocalized pi[dot dot dot]pi clusters should be carried out with basis sets of triple-zeta quality. Larger basis sets and the expensive counterpoise correction to the molecular geometry are not necessary. The rigid monomer approximation has very little effect on structure and energetics of these dimers and may be used without consequence. Due to a consistent cancellation of errors, optimization with the MP2 method leads to CCSD(T)/CBS interaction energies that are within 0.2 kcal mol(-1) of those for structures optimized with the CCSD(T) method. Future studies that aim to resolve structures separated by a few tenths of a kcal mol(-1) should consider the effects of optimization with the CCSD(T) method.     

Noncovalent interactions between cytosine and SWCNT: curvature dependence of complexes via pi...pi stacking and cooperative CH...pi/NH...pi

Fragments of C24H12, adapted from a variety of armchair [(n,n), (n = 5, 7, and 8)] and zigzag [(m,0) (m = 8, 10, and 12)] single-walled carbon nanotube (SWCNT), are used to model corresponding SWCNTs with different diameters and electronic structures. The parallel binding mainly through pi...pi stacking interaction, as well as the perpendicular binding via cooperative NH...pi and CH...pi between cytosine and the fragments of SWCNT have been extensively investigated with a GGA type of DFT, PW91LYP/6-311++G(d,p). The eclipsed tangential (ET) conformation with respect to the six-membered ring of cytosine and the central ring of SWCNT fragments is less stable than the slipped tangential (ST) conformation for the given fragment; perpendicular conformations with NH2 and CH ends have higher negative binding energy than those with NH and CH ends. At PW91LYP/6-311++G(d,p) level, two tangential complexes are less bound than perpendicular complexes. However, as electron correlation is treated with MP2/6-311G(d,p) for PW91LYP/6-311++G(d,p) optimized complexes, it turns out there is an opposite trend that two tangential complexes become more stable than three perpendicular complexes. This result implies that electron correlation, a primary source to dispersion energy, has more significant contributions to the pi...pi stacking complexes than to the complexes via cooperative NH...pi and CH...pi interactions. In addition, it was found for the first time that binding energies for two tangential complexes become more negative with increasing nanotube diameter, while those for three perpendicular complexes have a weaker dependence on the curvature; i.e., binding energies are slightly less and less negative. The performance of a novel hybrid DFT, MPWB1K, was also discussed.     

I...pi halogen-bonding in the structures of platinum-alkynyl compounds: an overlooked supramolecular interaction

The crystal structures of both a symmetric and an unsymmetric platinum--alkynyl compound with terminal iodophenyl functionalities have been determined, and found to exhibit an attractive interaction between the iodine atom and the carbon-carbon triple bond of the alkynyl ligand.     

Separation of sigma and pi energies

This work presents an all-electron density functional theory implementation of the separation of sigma and pi energies. On the basis of the separation of the electronic density, rho, into sigma and pi parts, an ansatz for the separation of the exchange-correlation energy is proposed. The behavior of the sigma and pi energy parts in benzene is investigated under different distortions. The effect of local and nonlocal functionals on the separation of the exchange-correlation energy is studied, too.     

CH...pi interaction for rhenium-based rectangles: an interaction that is rarely designed into a host-guest pair

Alkoxy- and thiolato-bridged Re(I) molecular rectangles [{(CO)3Re(mu-ER)2Re(CO)3}2(mu-bpy)2] (ER = SC4H9, 1a; SC8H17, 1b; OC4H9, 2a; OC12H25, 2b; bpy = 4,4'-bipyridine) exhibit strong interactions with several planar aromatic molecules. The nature of their binding was studied by spectral techniques and verified by X-ray diffraction analysis. Standard absorption and fluorescence titrations showed that a relatively strong 1:1 interaction occurs between aromatic guests such as pyrene and these rectangles. The results of a single-crystal X-ray diffraction analysis show that the recognition of 1 with a pyrene molecule is mainly due to CH...pi interactions and the face of the guest pyrene is located over the edges of the bpy linkers of 1. This is a fairly novel example of an interaction that is rarely designed into a host-guest pair. Furthermore, the interaction of 1 with Ag+ results in the self-organization of supramolecular arrays, as revealed by solid-state data.     

Halogen bonding and pi...pi stacking control reactivity in the solid state

The halogen bonding and the pi...pi stacking interactions induce the noncovalent self-assembly of modules into photoreactive supramolecular architecture. The pi...pi interaction pre-organizes the template, and the halogen bonding aligns the olefins to conform to the topochemical principle for photoreaction. The UV irradiation of the crystal resulted in a cyclization product with quantitative yield and stereospecificity.     

Possible artifacts occurring in the calculation of intermolecular energies from delocalized pictures

The intermolecular energy between two identical subsystems may be calculated from symmetrydelocalized MO's resulting for instance from a preliminary SCF calculation of the supersystem. Then each second-order energy correction mixes intramolecular correlation,R ^–6 intermolecular dispersion energy, andR ^–3 components. TheR ^–3 components disappear through subtle cancellations. The shifted Epstein-Nesbet energy denominators introduce an artificial second-order intermolecularR ^–1 component, which would be cancelled by off-diagonal third-order terms, as well as a bad asymptotic limit at infinite distances. TheR ^–1 artifact will also occur in strong symmetrical chemical bonds calculated in the Epstein-Nesbet perturbation scheme from delocalized MO's. These defects will occur in all variational approximate CI techniques which neglect off-diagonal elements between delocalized doubly excited determinants. These artifacts are avoided when using the Moller-Plesset definition of the zeroth order Hamiltonian or when starting from (SCF)localized MO's (even in the Epstein-Nesbet perturbation). The discussion is exemplified on an accurateab initio calculation of the Ar₂ molecule.     

Multicoefficient extrapolated density functional theory studies of pi...pi interactions: the benzene dimer

We report tests of new (2005) and established (1999-2003) multilevel methods against essentially converged benchmark results for nonbonded interactions in benzene dimers. We found that the newly developed multicoefficient extrapolated density functional theory (DFT) methods (which combine DFT with correlated wave function methods) give better performance than multilevel methods such as G3SX, G3SX(MP3), and CBS-QB3 that are based purely on wave function theory (WFT); furthermore, they have a lower computational cost. We conclude that our empirical approach for combining WFT methods with DFT methods is a very efficient and effective way for describing not only covalent interactions (as shown previously) but also nonbonded interactions.     

Experimental and theoretical determination of the accurate interaction energies in benzene-halomethane: the unique nature of the activated CH/pi interaction of haloalkanes

The CH/pi interaction energies between benzene and halomethanes (CH(2)Cl(2) and CHCl(3)) were accurately determined. Two-color ionization spectroscopy was applied to the benzene-CH(2)Cl(2) and -CHCl(3) clusters, and the binding energies in the neutral ground state, i.e. the CH/pi interaction energies in these model cluster systems, were precisely evaluated on the basis of the dissociation threshold measurements of the clusters in the cationic state and the ionization potential value of the bare molecule. The experimentally determined interaction energies were 3.8 +/- 0.2 and 5.2 +/- 0.2 kcal mol(-1) for benzene-CH(2)Cl(2) and -CHCl(3) respectively, and the remarkable enhancement of the CH/pi interaction energy with chlorine-substitution was quantitatively confirmed. The experimental interaction energies were well reproduced by the high-level ab initio calculations. The theoretical calculations clarified the unique nature of the activation of the CH/pi interaction by the chlorine-substitution.     

Perturbing the structure of the 2-norbornyl cation through C-H...N and C-H...pi interactions

Quantum chemical studies on complexes of the 2-norbornyl cation with ammonia and benzene are described. These calculations demonstrate that the delocalized, bridged, nonclassical geometry that is usually favored for this carbocation can be perturbed significantly toward the normally less favorable classical geometry through appropriately oriented noncovalent interactions. Since such cations have been proposed as intermediates in enzyme and antibody catalyzed reactions, these results have implications for the nature of the cationic species that may be generated in the presence of electron-rich amino acid side chains of the sort that may be present in the active sites of biocatalysts.     

Solvent-induced aggregation through metal...metal/pi...pi interactions: large solvatochromism of luminescent organoplatinum(II) terpyridyl complexes

A dramatic color change and tremendous emission enhancement have been "switched on" upon increasing diethyl ether ratio in acetonitrile or acetone solution of [Pt(tpy)(CC-CCH)]OTf, attributed to the formation of Pt...Pt and pi-pi interactions. Two crystal forms (dark-green and red) of [Pt(tpy)(CC-CCH)]OTf, together with [Pt(tBu3-tpy)(CC-CCH)]OTf, show different crystal-packing modes as revealed by X-ray crystallography.     

Length dependence of excitation energies in linear polyenes: Localized and delocalized descriptions

The length dependence of the lowest allowed transition energy of linear polyenes is studied using delocalized SCF and localized excitonic approaches. Within the PPP SCF approximations the calculated transition energies converge to a finite values as N^?1 as the number of double bonds (N) becomes large, when the excited state contains all singly excited configurations. On the other hand, the fully localized excitonic method at the level of single excitations, although it predicts a gap in the excitation spectrum of an infinite polyene, gives results which converge to this value as N^?2. The inclusion of double and triple excitations into the excitonic method by means of perturbation theory does not appear to change this behavior. The reasons for the discrepancy between the two approaches is analyzed. Experimentally, the transition energies in solution converge to a finite value as N^?1 to a good degree of approximation. If it is assumed that the solvent shift is constant for long polyenes, the available experimental results favour the delocalized approach as the starting point in describing the length dependence of the excitation energy of long polyenes.     

A bottom-up valence bond derivation of excitation energies in 1D-like delocalized systems

Using the chemically relevant parameters hopping integral t(0) and on-site repulsion energy U, the charge gap (lowest dipolarly allowed transition energy) in 1D systems is examined through a bottom-up strategy. The method is based on the locally ionized states, the energies of which are corrected using short-range delocalization effects. In a valence bond framework, these states interact to produce an excitonic matrix which accounts for the delocalized character of excited states. The treatment, which gives access to the correlated spectrum of ionization potentials, is entirely analytical and valid whatever the U/|t(0)| ratio for such systems ruled by Peierls-Hubbard Hamiltonians. This second-order analytical derivation is finally confronted to numerical results of a renormalized excitonic treatment using larger blocks as functions of the U/|t(0)| ratio. The method is applied to dimerized chains and to fused polybenzenic 1D lattices. Such approaches complement the traditional Bloch-function based picture and deliver a conceptual understanding of the charge gap opening process based on a chemical intuitive picture.     

Electrophilic behavior of the pi delocalized azepinium ion: Friedel-Crafts reactions with benzenes and five-membered aromatic heterocycles

[reaction: see text] Although the reactivity of tropylium ion with aromatic substrates is low, the reaction of azepinium ion with aromatic substrates such as benzene, phenol, furan, and thiophene resulted in the formation of 2-aryl-2H-azepine as a major product. An exceptional result in the formation of ring-contracted product was observed in the reaction with pyrrole.     

Computational characterization and modeling of buckyball tweezers: density functional study of concave-convex pi...pi interactions

The geometries and binding energies of a recent buckyball tweezers (C(60)H(28)) and its supramolecular complexes are investigated using recently developed density functionals (M06-L and M06-2X) that include an accurate treatment of medium-range correlation energy. The pincer part of the tweezers, corannulene, has a strong attractive interaction with C(60). However, due to the entropy penalty, the calculated gas-phase free energy of association of the C(60)@corannulene supramolecule is positive 3.5 kcal mol(-1); and this entropy penalty explains why it is difficult to observe C(60)@corannulene supramolecule experimentally. By using a pi-extended tetrathiafulvalene (TTF), in particular 9,10-bis(1,3-dithiol-2-ylidene)-9,10-dihydroanthracene (TTFAQ or C(20)H(10)S(4)), as the pincer part, we modeled a new buckyball tweezers. The geometries and binding energies of the new buckyball tweezers and its supramolecular complexes are also calculated. Due to fact that the attractive interaction between TTFAQ and C(60) is weaker than that between corannulene and C(60), the gas-phase binding free energy in the C(60)@C(60)H (32)S(8) supramolecular complex is smaller than that in the C(60)@C(60)H(28) supramolecule. We also discuss solvent effects.     

Origin of attraction and directionality of the pi/pi interaction: model chemistry calculations of benzene dimer interaction.

A model chemistry for the evaluation of intermolecular interaction between aromatic molecules (AIMI Model) has been developed. The CCSD(T) interaction energy at the basis set limit has been estimated from the MP2 interaction energy near the basis set limit and the CCSD(T) correction term obtained by using a medium size basis set. The calculated interaction energies of the parallel, T-shaped,and slipped-parallel benzene dimers are -1.48, -2.46, and -2.48 kcal/mol, respectively. The substantial attractive interaction in benzene dimer, even where the molecules are well separated, shows that the major source of attraction is not short-range interactions such as charge-transfer but long-range interactions such as electrostatic and dispersion. The inclusion of electron correlation increases attraction significantly. The dispersion interaction is found to be the major source of attraction in the benzene dimer. The orientation dependence of the dimer interaction is mainly controlled by long-range interactions. Although electrostatic interaction is considerably weaker than dispersion interaction, it is highly orientation dependent. Dispersion and electrostatic interactions are both important for the directionality of the benzene dimer interaction.     

Probing the interaction mode in hydrophilic interaction chromatography

This work aims at characterizing interactions between a select set of probes and 22 hydrophilic and polar commercial stationary phases, to develop an understanding of the relationship between the chemical properties of those phases and their interplay with the eluent and solutes in hydrophilic interaction chromatography. "Hydrophilic interaction" is a somewhat inexact term, and an attempt was therefore made to characterize the interactions involved in HILIC as hydrophilic, hydrophobic, electrostatic, hydrogen bonding, dipole-dipole, π-π interaction, and shape-selectivity. Each specific interaction was quantified from the separation factors of a pair of similar substances of which one had properties promoting the interaction mode being probed while the other did not. The effects of particle size and pore size of the phases on retention and selectivity were also studied. The phases investigated covered a wide range of surface functional groups including zwitterionic (sulfobetaine and phosphocholine), neutral (amide and hydroxyl), cationic (amine), and anionic (sulfonic acid and silanol). Principal component analysis of the data showed that partitioning was a dominating mechanism for uncharged solutes in HILIC. However, correlations between functional groups and interactions were also observed, which confirms that the HILIC retention mechanism is partly contributed by adsorption mechanisms involving electrostatic interaction and multipoint hydrogen bonding. Phases with smaller pore diameters yielded longer retention of solutes, but did not significantly change the column selectivities. The particle diameter had no significant effect, neither on retention, nor on the selectivities. An increased water content in the eluent reduced the multipoint hydrogen bonding interactions, while an increased electrolyte concentration lowered the selectivities of the tested columns and made their interaction patterns more similar.