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.     

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.     

N-methylformamide-benzene complex as a prototypical peptide N-H...pi hydrogen-bonded system: density functional theory and MP2 studies

Although the existence of peptide N-H...pi hydrogen bonds recently has been reported in protein structures, little is known about their strength and binding nature and, therefore, the relative importance of the interaction. To shed light on this binding, the N-methylformamide-benzene complex, as a model of the peptide N-H...pi hydrogen bonding, was studied by using density functional theory and M?ller-Plesset second-order perturbation (MP2) methods. The geometry of the complex was fully optimized at the B3LYP/6-31G(d,p) and MP2/6-31G(d,p) levels. The optimized interaction distances are about 3.6 and 3.2 A, respectively, at the two levels. The binding energy corrected by basis set superposition error with the MP2/cc-pVTZ method based on the MP2/6-31G geometry is -4.37 kcal/mol, which is as strong as the conventional hydrogen bonding. The calculated results suggest that the peptide N-H...pi hydrogen bonding is of sufficient strength to play an important role in the stabilization of protein structures. These results are helpful to better understand the characteristics and nature of the peptide N-H...pi interaction as well as to modify current force fields to better represent this special interaction.     

Quasi-relativistic density functional study of aurophilic interactions

Fifteen molecules containing the Au(I) species have been calculated by ab initio HF and MP2 methods and by five different density functional approaches. The aurophilic Au(d10)-Au(d10) bonding mechanism has been investigated. Both, one-electron interactions (i.e., electrostatic, polarization, charge transfer, and orbital interference) and two-electron effects (i.e., correlation, dispersion) contribute significantly to the so-called 'secondary' or metallophilic bonds representing the Au-Au interaction. Second, the applicability of density functional approaches to this type of bonding has been tested. It is well-known that present day density functionals are not yet designed to simulate the long-range London dispersion forces between nonoverlapping systems, whereas they approximately reproduce the short range dynamical electron correlations of strongly overlapping chemically bonded nondegenerate species. It is found here empirically for the investigated groups of gold(I) cluster compounds that simple local density functionals (LDF) of the Slater (or Slater plus Vosko) type yield rather reasonable estimates for the equilibrium distances, and (on the average) also for the aurophilic interaction energies, though with rather large standard deviations. Still LDF are useful for survey investigations of Au cluster compounds. Common gradient corrected DF are not recommended here, nor are the large core pseudopotentials for Au.     

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.     

Exploring the lower limit in hydrogen bonds: analysis of weak C-H...O and C-H...pi interactions in substituted coumarins from charge density analysis

Experimental electron densities in coumarin, 1-thiocoumarin, and 3-acetylcoumarin have been analyzed based on the X-ray diffraction data at 90 K. These compounds pack in the crystal lattice with weak C-H...O and C-H...pi interactions, and variations in charge density properties and derived local energy densities have been investigated in the regions of intermolecular interactions. Theoretical charge density calculations on crystals using the B3LYP/6-31G* method show remarkable agreement with the derived properties and energy densities from the experiment. The intermolecular interactions follow an exponential dependence of electron density and energy densities at the bond critical points. The Laplacian follows a "Morse-like" dependence on the length of the interaction line. Based on the set of criteria defined using the theory of "atoms in molecules", it has become possible to distinguish between a hydrogen bond (C-H...O) and a van der Waals interaction (C-H...pi). This has resulted in the identification of a "region of overlap" in terms of electron densities, energy densities, and mutual penetration of the hydrogen and acceptor atoms with respect to the interaction length. This approach suggests a possible tool to distinguish between the two types of interactions.     

Oxa-bicyclocalixarenes: a new cage for anions via C-H...anion hydrogen bonds and anion...pi interactions

Density functional theory calculations indicate that the cage molecule 4 can trap F- in the gas phase (-80.5 kcal/mol) as well as in CH2Cl2 (-14.7 kcal/mol) via strong C-H...F- hydrogen bonds and pi...F- interaction.     

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.     

Coordination complexes exhibiting anion...pi interactions: synthesis, structure, and theoretical studies

The polydentate ligand 2,4,6-tris(dipyridin-2-ylamino)-1,3,5-triazine (dpyatriz) in combination with the Cu(ClO 4) 2/CuX 2 salt mixtures (X (-) = Cl (-), Br (-), or N 3 (-)) leads to the formation of molecular coordination aggregates with formulas [Cu 3Cl 3(dpyatriz) 2](ClO 4) 3 ( 2), [Cu 3Br 3(dpyatriz) 2](ClO 4) 3 ( 3), and [Cu 4(N 3) 4(dpyatriz) 2(DMF) 4(ClO 4) 2](ClO 4) 2 ( 4). These complexes consist of two dpyatriz ligands bridged via coordination to Cu (II) and disposed either face-to-face in an eclipsed manner ( 2 and 3) or parallel and mutually shifted in one direction. The copper ions complete their coordination positions with Cl (-) ( 2), Br (-) ( 3), or N 3 (-), ClO 4 (-), and N, N-dimethylformamide (DMF) ( 4) ligands. All complexes crystallize together with noncoordinate ClO 4 (-) groups that display anion...pi interactions with the triazine rings. These interactions have been studied by means of high level ab initio calculations and the MIPp partition scheme. These calculations have proven the ClO 4 (-)...[C 3N 3] interactions to be favorable and have revealed a synergistic effect from the combined occurrence of pi-pi stacking of triazine rings and the interaction of these moieties with perchlorate ions, as observed in the experimental systems.     

Computational study of water and ammonia dimers with density functional theory methods

The structures and energies for the dimerization of water and ammonia molecules were computed with density functional theory (DFT) and ab initio methods. For all studies the same 6-311+G(2d,2p) basis set was used. Two linear hydrogen-bonded and cyclic ammonia dimer structures were computed and their relative stability is discussed. From the systematic studies, hybrid DFT methods were selected as reliable for computing the parameters of these types of van der Waals' complex.     

N-H...pi interactions in pyrroles: systematic trends from the vibrational spectroscopy of clusters

Pyrrole and some of its methylated derivatives are aggregated in a controlled way in pulsed supersonic jet expansions. The cluster N-H stretching dynamics is studied using FTIR and Raman spectroscopy. Dimers, trimers and tetramers can be differentiated. Systematic trends in the dimer N-H...pi interaction as a function of methyl substitution are identified and explored for predictions. Overtone jet absorption spectroscopy is used to extract anharmonicities for the N-H bond in different environments. The N-H anharmonicity constant increases by 10% upon dimerization. Bulk matrix shifts can be emulated by the formation of Ar-decorated clusters. The experimental results are expected to serve as benchmarks for an accurate ab initio characterization of the N-H...pi hydrogen bond.     

Interplay of phenyl-perfluorophenyl stacking,C-H...F,C-F...pi and F...F interactions in some crystalline aromatic azines

Analysis of phenyl-perfluorophenyl stacking synthon, C-H...F, C-F...pi interactions, and F...F tetramer in three closely related azine crystal structures shows the dominance of Ar-ArF synthon while other interactions are turned on/off depending on the H/F stoichiometry in the molecule.     

Reversible anion exchange and sensing in large porous materials built from 4,4'-bipyridine via pi...pi and H-bonding interactions

Two unusual d10 compounds, [Zn2(bipy)3(H2O)8(ClO4)2(paba)2].2(bipy).4H2O (1) and [Cd2(bipy)4(H2O)6(ClO4)2(paba)2].(bipy).5H2O (2) (bipy = 4,4'-bipyridine, paba = p-aminobenzoate), were obtained from reaction of the metal salts, bipy and paba in an EtOH/H2O mixture. The bipy ligands in the two compounds exhibit two new modes of coordinating to transition metal ions, resulting in the formation of large porous frameworks. Immersion of single crystals of 1 in an aqueous solution of NH4PF6 results in the formation of its hexafluorophosphate derivative 3 as shown by single crystal diffraction. Immersion of crystals of 3 in NaClO4 regenerates 1. Furthermore, compound 1 also shows interesting anion sensing properties in an EtOH/H2O mixture.     

Encapsulation of small base molecules and tetrahedral/cubane-like clusters of group V atoms in the boron buckyball: a density functional theory study

A density functional theory study of small base molecules and tetrahedral and cubane-like group V clusters encapsulated in B(80) shows that the boron buckyball is a hard acid and prefers hard bases like NH(3) or N(2)H(4) to form stable off-centered complexes. In contrast, tetrahedral and cubane-like clusters of this family are metastable in the cage. The most favorable clusters are the mixed tetrahedral and cubane clusters formed by nitrogen and phosphorus atoms such as P(2)N(2)@B(80), P(3)N@B(80), and P(4)N(4)@B(80). The boron cap atoms are electrophilic centers, and prefer mainly to react with electron rich nucleophilic sites. The stability of the complexes will be governed by the size and electron donating character of the encapsulated clusters. B(80) forms stable complexes with hard materials where a bidentate interaction of the encapsulated molecule with two boron cap atoms is preferred over a single direct complex toward a single endohedral boron.     

Comparison of C-H...pi and hydrophobic interactions in a beta-hairpin peptide: impact on stability and specificity

We have examined the impact of C-H...pi and hydrophobic interactions in the diagonal position of a beta-hairpin peptide through comparison of the interaction of Phe, Trp, or Cha (cyclohexylalanine) with Lys or Nle (norleucine). NMR studies, including NOESY and chemical shift perturbation studies, of the Lys side chain indicates that Lys interacts in a specific geometry with Phe or Trp through the polarized C epsilon. In contrast, Nle does not interact in a specific manner with the diagonal aromatic residue. Thermal denaturation provides additional support that Lys and Nle interact in fundamentally different manners. Folding of the peptide with a diagonal Trp...Lys interaction was found to be enthalpically driven, whereas the peptide with a diagonal Trp...Nle interaction displayed cold denaturation, as did the control peptide with a diagonal Cha...Nle interaction, indicating different driving forces for interaction of Lys and Nle with Trp. These findings have significant implications for specificity in protein folding and de novo protein design.     

Design of a magnetic bistability molecular system constructed by H-bonding and pi...pi-stacking interactions

Crystal structures and magnetic properties were determined for two novel polymorphs of the complex [H2DABCO][Ni(mnt)2] [(H2DABCO)2+ = diprotonated 1,4-diazabicyclo[2.2.2]octane; mnt2- = maleonitriledithiolate]. For each polymorph, anions form a layered structure in which two kinds of dimers were observed. The adjacent anionic sheets are held together by cations via H-bonding interactions between protons of cations and CN groups of anions. Two polymorphs possess spin bistability; namely, upon cooling, a magnetic transition happens at around 120 K with about 1 K hysteresis on heating for the alpha phase and at 112 K with about 10 K hysteresis for the beta phase. Above the transition, the magnetic behaviors of two polymorphs can be approximately interpreted by a singlet-triplet model of an antiferromagnetically coupled S = 1/2 dimer, which is supported by the crystal structures and spin dimer analyses based on extended Hückel molecular orbital calculations.     

The response of selected isomers of B80 buckyball toward NH3 adsorption: a density functional theory investigation

Density functional theory calculations at the B3LYP/6-31G(d) and B3LYP/6-31+G(d) levels were carried out for the adsorption of NH₃ on three symmetric isomers of B₈₀ {C ₁, T _h, I _h}. To investigate the binding features of B₈₀ isomers with NH₃, different studies including the structural and electronic parameters, the ¹⁴N electric field gradient tensors and the atoms in molecules (AIM) properties were considered. The calculated parameters by these investigations can be used as powerful tools to find out some of the unknown aspects of electronic structures of the boron buckyball and its isomers. According to previous studies, boron buckyball as an amphoteric and a hard molecule has two distinct reactive sites defined as cap and frame which act as an acid and a base, respectively. Regarding the obtained results in this study, all the isomers had the same exposure when NH₃ molecule reacted with the external wall of B₈₀. For instance, the stability of N–B bond in the cap site was significantly more than the stability of N–B bond in the frame. Moreover, the adsorption of NH₃ on frame site showed a considerable reduction in HOMO–LUMO energy gap. According to AIM theory, an electrostatic nature was observed for N–B interaction. Concerning the selected isomers of buckyball, the capability of the NH₃–B₈₀ complexes to localize electron at the N–B bond critical points depend on the reaction sites significantly. In general, ¹⁴N nuclear quadruple coupling constants and asymmetry parameter reveal a remarkable effect of NH₃ adsorption on electronic structure of the B₈₀.     

NMR and quantum chemical study on the OH...pi and CH...O interactions between trehalose and unsaturated fatty acids: implication for the mechanism of antioxidant function of trehalose

Trehalose is a disaccharide that attracts much attention as a stress protectant. In this study, we investigated the mechanism of the antioxidant function of trehalose. The spin-lattice relaxation times (T(1)) of (1)H and (13)C NMR spectra were measured to investigate the interaction between trehalose and unsaturated fatty acid (UFA). We selected several kinds of UFA that differ in the number of double bonds and in their configurations (cis or trans). Several other disaccharides (sucrose, maltose, neotrehalose, maltitol, and sorbitol) were also analyzed by NMR. The T(1) values for the (1)H and (13)C signals assigned to the olefin double bonds in UFA decrease with increasing concentration of trehalose and the changes reaches plateaus at integer ratios of trehalose to UFA. The characteristic T(1) change is observed only for the combination of trehalose and UFA with cis double bond(s). On the other hand, from the (13)C-T(1) measurements for trehalose, the T(1) values of the C-3 (C-3') and C-6' (C-6) are found to change remarkably by addition of UFA. (1)H[bond](1)H NOESY measurements provide direct evidence for complexation of trehalose with linoleic acid. These results indicate that one trehalose molecule stoichiometrically interacts with one cis-olefin double bond of UFA. Computer modeling study indicates that trehalose forms a stable complex with an olefin double bond through OH...pi and CH...O types of hydrogen bonding. Furthermore, a significant increase in the activation energy is found for hydrogen abstraction reaction from the methylene group located between the double bonds that are both interacting with the trehalose molecules. Therefore, trehalose has a significant depression effect on the oxidation of UFA through the weak interaction with the double bond(s). This is the first study to elucidate the antioxidant function of trehalose.     

C-H...pi interactions in the low-temperature crystal structures of alpha,omega-unsaturated linear hydrocarbons

The low-temperature crystal structures of 1,7-octadiene (C8H14), 1,9-decadiene (C10H18), 1,7-octadiyne (C8H10) and 1,9-decadiyne (C10H14), determined following in situ crystal growth from the liquid, reveal the role of C-H...pi interactions in alpha,omega-unsaturated linear hydrocarbons in the solid state.