Density Functional Theory Studies on the Intermolecular Interactions of Five Aza-calix [6] arene Host with HMX

Five fully optimized structures of complexes between aza-calix[6]arene host monomers(Ma～Me) and complexes(a～e) have been obtained at the B3LYP/6-31G(d) level.Natural bond orbital(NBO) analysis was performed to reveal the origin of the interaction.The intermolecular interaction energy was evaluated with basis set superposition error correction(BSSE) and zero point energy correction(ZPEC).The B3LYP/6-31G(d) calculations on the five complexes have shown that the greatest interaction(-13.98 kJ/mol) is found in the complex between HMX and hexa-aza-calix[3]-p-tri-arene[3]-2-amido-1,3,5-tri-azine.The results have indicated that intermolecular interaction energies of aza-calix[6]arenes with substituted group are stronger than those without substituted group,and those with amido are greater than with nitryl.Thus,hexa-azacalix[3]-p-tri-arene[3]-2-amido-1,3,5-tri-azine is rather equal to eliminate HMX from explosive waste water.     

Theory Studies on the Intermolecular Interactions of Urea Nitrate with RDX

Six fully optimized geometries of urea nitrate cation and RDX complexes have been obtained with DFT-B3LYP and MP2 methods at the 6-311++G** level. The intermolecular interaction energies have been calculated with basis set superposition error (BSSE) and zero point energy (ZPE) correction. The nature of intermolecular interaction has been revealed by the analysis of AIM and NBO. The results indicate that the greatest binding energy of urea nitrate with RDX is –82.47kJ/mol. The O–H…O and N–H…O hydrogen bonds are important intermolecular interactions of urea nitrate cation with RDX, and the origin of hydrogen bonds is the oxygen atom offering its lone-pair electrons to the σ(O-H)* or σ(O-H)* antibonding orbital. The intermolecular interactions strengthen the N–NO2 bond, leading to the reduced sensitivity of urea nitrate and RDX mixture explosive.     

An ab initio study of intermolecular interactions of nitromethane dimer and nitromethane trimer

Different geometries of nitromethane dimer and nitromethane trimer have been fully optimized employing the density functional theory B3LYP method and the 6-31++G** basis set. Three-body interaction energy has been obtained with the ab initio supermolecular approach at the levels of MP2/6-31++G**//B3LYP/6-31++G** and MP2/aug-cc-pVDZ//B3LYP/6-31++G**. The internal rotation of methyl group induced by intermolecular interaction has been observed theoretically. For the optimized structures of nitromethane dimer, the strength of C--H...O--N H-bond ranges from -9.0 to -12.4 kJ mol(-1) at the MP2/aug-cc-pVDZ//B3LYP/6-31++G** level, and the B3LYP method underestimates the interaction strength compared with the MP2 method, while MP2/6-31++G**//B3LYP/6-31++G** calculated DeltaE(C) is within 2.5 kJ mol(-1) of the corresponding value at the MP4(SDTQ)/6-31G**//B3LYP/6-31++G** level. The analytic atom-atom intermolecular potential has been successfully regressed by using the MP2/6-31++G**//B3LYP/6-31++G** calculated interaction energies of nitromethane dimer. For the optimized structures of nitromethane trimer the three-body interaction energies occupy small percentage of corresponding total binding energies, but become important for the compressed nitromethane explosive. In addition, it has been discovered that the three-body interaction energy in the cyclic nitromethane trimer is more and more negative as intermolecular distances decrease from 2.2 to 1.7 A.     

Topomerization of 2-Halo-[9] Annulen Anion at the DFT Level

Topomerization of [9] annulen anion(1)and its 2-fluoro-,2-chloro-and 2-bromo-derivatives(2,3 and 4,respectively)are studied at the HF/6-31G* and B3LYP/6-311++G** levels of theory.The relative ease of topomerization is dependent on the charge distribution and planarity of the ground state and the transition state of 9-membered rings as well as the size and electronegativity of halogen substituent.Consequently,the endo-2-halo-[9] annulen anion topomers become unstable and easily convert to related exo-topomers.Hence,according to the DFT calcu-lations,the order of topomerization energy barrier for endo exo topomerization is 1endo > 2 endo > 3endo.     

Theoretical Studies on Intermolecular Hydrogen-bond Interactions between Hexamethylenetetramine and Nitric Acid

The structures of the complexes generated by hexamethylenetetramine and nitric acid have been fully optimized by B3LYP method at the 6-311++G** and aug-cc-pVTZ levels. The intermolecular hydrogen-bonding interactions have been calculated by the B3LYP/6-311++G**, B3LYP/aug-cc-pVTZ, MP2(full)/6-311++G** and CCSD(T)/6-311++G** methods, respectively. The NBO (nature bond orbital), AIM (atom in molecule), temperature effect and solvation effect have been analyzed to reveal the origin of the interactions. The results indicate that the stable hydrogen-bonded complexes could be generated by hexamethylenetetramine and nitric acid. The interactions follow the order of (a)>>(e)>(b)>(c)>(d)>(f)>(g). The C-N bonds which are adjacent to the methylene involving the hydrogen bonds tend to break in the chemical reaction. Due to the exothermic process, low temperature is conducive to the formation of the composition, which tallies with the experimental result.     

气相中O3与HSO自由基间的氢键复合物

气相中O3与HSO自由基之间的相互作用及其反应在大气化学中非常重要.在DFT-B3LYP/6-311++G**和MP2/6-311++G**水平上求得O3+HSO复合物势能面上的稳定构型,B3LYP方法得到了三种构型(复合物Ⅰ,Ⅱ和Ⅲ),而MP2方法只能得到一种构犁(复合物Ⅱ).在复合物Ⅰ和Ⅲ中,HSO单元中的1H原子作为质子供体.与O3分子中的端基O原子作为质子受体相互作用,形成红移氢键复合物;而在复合物Ⅱ中,虽与复合物Ⅰ和Ⅲ中具有相间的质子供体和质子受体,却形成了蓝移氢键复合物.B3LYP/6-311++G**水平上计算的单体间相互作用能的计算考虑了基组重甍误差(BSSE)和零点振动能(ZPVE)校正,其值在-3.37到-4.55 kJ·mol-1之间.采用自然键轨道理论(NBO)对单体间相互作用的本质进行了考查,并通过分子中原子理论(AIM)分析了三种复合物中氢键的电子密度拓扑性质.     

Interaction of anions with perfluoro aromatic compounds

The complexes formed by a variety of anions with perfluoro derivatives of benzene, naphthalene, pyridine, thiophene, and furan have been calculated using DFT (B3LYP/6-31++G**) and MP2 (MP2/6-31++G** and MP2/6-311++G**) ab initio methods. The minimum structures show the anion interacting with the pi-cloud of the aromatic compounds. The interaction energies obtained range between -8 and -19 kcal mol(-1). The results obtained at the MP2/6-31++G** and MP2/6-311++G** levels are similar. However, the B3LYP/6-31++G** results provide longer interaction distances and smaller interaction energies than do the MP2 results. The interaction energies have been partitioned using an electrostatic, polarization, and van der Waals scheme. The AIM analysis of the electron density shows a variety of topologies depending on the aromatic system considered.     

气相中O3与HSO自由基间的氢键复合物

气相中O3与HSO自由基之间的相互作用及其反应在大气化学中非常重要. 在DFT-B3LYP/6-311++G**和MP2/6-311++G**水平上求得O3+HSO复合物势能面上的稳定构型, B3LYP方法得到了三种构型(复合物I, II和III), 而MP2方法只能得到一种构型(复合物II). 在复合物I和III中, HSO单元中的1H原子作为质子供体, 与O3分子中的端基O原子作为质子受体相互作用, 形成红移氢键复合物; 而在复合物II中, 虽与复合物I和III中具有相同的质子供体和质子受体, 却形成了蓝移氢键复合物. B3LYP/6-311++G**水平上计算的单体间相互作用能的计算考虑了基组重叠误差(BSSE)和零点振动能(ZPVE)校正, 其值在-3.37到-4.55 kJ·mol-1之间. 采用自然键轨道理论(NBO)对单体间相互作用的本质进行了考查, 并通过分子中原子理论(AIM)分析了三种复合物中氢键的电子密度拓扑性质.     

Theoretical study of light-metal tetrahydroaluminates

Ab initio calculations of the structural, energetic, vibrational, and magnetic characteristics of the lowest-lying structures for isolated molecules and ions of light-metal tetrahydroaluminates (Li, Na, Be, Mg, and Al) have been performed by the perturbation theory (MP2), quadratic configuration interaction (QCISD(T)), coupled cluster (CCSD(T)), and density functional theory (B3LYP) methods using the 6-31G*, 6-31G**, 6-311+G**, and 6-311++G** basis sets. The trends in the behavior of the molecular characteristics have been analyzed in various related series of these compounds. The results are compared with the data on analogous light-metal tetrahydroborates calculated at the same levels of theory. The differences in structure and stability between analogous hydroborate and alanate complexes are examined. The economical approximation B3LYP/6-311++G**//B3LYP/6-311+G** has been shown to adequately reproduce the results obtained at the higher level of theory CCSD(T)/6-311++G**//MP2/6-31G* even though it requires considerably shorter CPU times and smaller amounts of memory.     

Correlated ab initio quantum chemical calculations of di- and trisaccharide conformations

High level correlated quantum chemical calculations, using MP2 and local MP2 theory, have been performed for conformations of the disaccharide, beta-maltose, and the trisaccharide, 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranose. For beta-maltose, MP2 and local MP2 calculations using the 6-311++G** basis set are in good agreement, predicting a global minimum gas-phase conformation with a counterclockwise hydrogen bond network and the experimentally-observed intersaccharide hydrogen bonding arrangement. For conformations of 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranose, MP2/6-311++G**, and local MP2/6-311++G** calculations do not provide a consensus prediction of relative energetics, with the MP2 method finding large differences in stability between extended and folded trisaccharide conformations. Local MP2 calculations, less susceptible to intramolecular basis set superposition errors, predict a narrower range of trisaccharide energetics, in line with estimates from Hartree-Fock theory and B3LYP and BP86 density functionals. All levels of theory predict compact, highly hydrogen-bonded conformations as lowest in energy on the in vacuo potential energy surface of the trisaccharide. These high level, correlated local MP2/6-311++G** calculations of di- and trisaccharide energetics constitute potential reference data in the development and testing of improved empirical and semiempirical potentials for modeling of carbohydrates in the condensed phase.     

HOCl…HCOCl复合物的结构和电子性质

在DFT-B3LYP/6-311++G**水平上求得HOCl+HCOCl复合物势能面上的四种稳定构型(S1, S2, S3和S4). 其中, 在复合物S1和S3中, HOCl单体的5H原子作为质子供体, 与HCOCl单体中作为质子受体的1O原子相互作用, 形成红移氢键复合物; 在复合物S4中, HOCl单体的7Cl原子作为质子供体, 与HCOCl单体中作为质子受体的1O原子相互作用, 形成红移卤键复合物; 而在复合物S2中, 同时存在2C—3H…6O蓝移氢键和4Cl…5O相互作用. 在MP2/6-311++G**水平上计算的单体间的相互作用能考虑了基组重叠误差(BSSE)和零点振动能(ZPVE)校正, 其值在-5.05与-14.76 kJ·mol-1之间. 采用自然键轨道理论(NBO)对两种单体间相互作用的本质进行了考查, 并通过分子中原子理论(AIM)分析了复合物中氢键和卤键键鞍点处的电子密度拓扑性质.     

Comparative study of calix[4]arene derivatives: implications for ligand design

The first comparative theoretical study of three parent calix[4]arene analogues (calix[4]arene, thiacalix[4]arene, and homooxacalix[4]arene) has been performed using molecular dynamic simulations and density functional theory (MPWB1K/6-311G^∗∗//B3LYP/6-311G^∗∗) methods. The theoretical observations herein including optimized geometry, polarity, and atomic charge data provide that homooxacalix[4]arene would offer more efficient platform for metal ion recognition compared to thiacalix[4]arene or calix[4]arene.     

The magnitude of [C-H...O] hydrogen bonding in molecular and supramolecular assemblies

Ab initio calculations at the MP2/6-311++G** level on model systems (N-methylpyridinium complexes of dimethyl ether and dimethyl phosphate anion) provide quantitative measures of the large stabilization energies that arise from [C-H...O] contacts in charged systems. These attractive interactions control (i) the self-assembly of bipyridinium-based catenanes and rotaxanes in solution, (ii) the self-organization of left-handed Z-DNA with alternating [dC-dG] sequences in the solid state, and (iii) the binding of pyridinium derivatives with single- and double-stranded DNA. Slightly attractive interactions occur between the donor ether and phosphate moieties and a neutral pyridine molecule in the gas phase. Electrostatic potential and solvation calculations demonstrate that [C-H...O] interactions which involve a cationic [C-H] donor are dominated by electrostatic terms.     

Structural studies on ethyl isovalerate by microwave spectroscopy and quantum chemical calculations

We observed the microwave spectrum of ethyl isovalerate by molecular beam Fourier transform microwave spectroscopy. The rotational and centrifugal distortion constants of the most abundant conformer were determined. Its structure was investigated by comparison of the experimental rotational constants with those obtained by ab initio methods. In a first step, the rotational constants of various conformers were calculated at the MP2/6-311++G** level of theory. Surprisingly, no agreement with the experimental results was found. Therefore, we concluded that in the case of ethyl isovalerate more advanced quantum chemical methods are required to obtain a reliable molecular geometry. Ab initio calculations carried out at MP3/6-311++G**, MP4/6-311++G**, and CCSD/6-311++G** levels and also density functional theory calculations using the B3LYP/6-311++G** method gave similar results for the rotational constants, but they were clearly distinct from those obtained at the MP2/6-311++G** level. With use of these more advanced methods, the rotational constants of the lowest energy conformer were in good agreement with those obtained from the microwave spectrum.     

Spectrophotometric study of the supramolecular complexes of [60]- and [70]Fullerenes with biscalix[6]arene and crown[4]calix[6]arene

[60]- and [70]Fullerenes have been shown to form 1:1 supramolecular complexes with bis[2-(5,11,17,23,29,35-hexa-tert-butyl-37,38,39,40,41-pentahydroxycalix[6]arenyl-oxy ethyl ether) (1) and 5,11,17,23,29,35-hexa-tert-butyl-37,38,40,41-tetra hydroxyl-39,42-(crown-4)calix[6]arene (2) in CHCl3 medium by electronic absorption spectroscopy. Formation constants (K) of the complexes of [60]- and [70]fullerenes with 1 and 2 have been determined at room temperature from which free energy of formation values of the complexes have been estimated. The very high formation constant value of [60]fullerene/1 complex (5900 dm3 mol-1) in indicative of formation of inclusion complex. Moreover, PM3 calculations reveal that intermolecular interaction between [60]fullerene and 1 proceeds through quite deep energy molecular orbital.     

Polar substituent effects in the bicyclo[1.1.1]pentane ring system: acidities of 3-substituted bicyclo[1.1.1]pentane-1-carboxylic acids

Experimental gas-phase acidities are reported for a series of 3-substituted (X) bicyclo [1.1.1]pent-1-yl carboxylic acids (1, Y = COOH). A comparison with available calculated data (MP2/6-311++G**// B3LYP/6-311+G**) reveals good agreement. The relative substituent effects are shown to be adequately described by a much lower level of theory (B3LYP/6-31+G*). Various correlations are presented which clearly point to polar field effects as being the origin of the relative acidities.     

Basis set superposition error in MP2 and density-functional theory: a case of methane-nitric oxide association

A systematic study of basis set superposition error (BSSE) behavior in H(3)C-H[ellipsis (horizontal)][NO] complexes for both -H...N- and -H...O- orientations were carried out using MP2 and density-functional theory with Pople's [6-31G(d,p),6-311++G(nd,nd), where n=1,2,3, and 6-311++G(3df,3pd)] and Dunning's augmented correlation consistent basis sets [aug-cc-pVXZ (X=D and T)]. Corrected and uncorrected counterpoise potential-energy surfaces (PESs) were explored and differences obtained between them indicate that reliable optimizations of these molecular interactions must be carried out in a PES free of BSSE, even in the case of large basis sets and popularly used functionals such as B3LYP. Although all basis used could be always considered within a margin of approximation for representing molecular orbitals and show important values of BSSE, 6-311++G(2d,2p) basis set shows the best results in uncorrected PES with respect to the corrected ones. B3LYP functional produces erratic results: complexes appear repulsive and the intermolecular distances are always large, evidencing the lack of a correct dispersive forces treatment in the original parameterization. According to the MP2 results, the -H...N- interactions appear as slightly more stable than those of the -H...O- orientation.     

AuClx (x=1, 3)催化2-炔丙基苯酮与苯炔环化反应机理

对金属元素Au采用LanL2DZ基组,对非金属元素C,H,O,Cl采用6-31G*和6-311++G**基组,用密度泛函理论的B3LYP、B3PW91、UB3LYP方法和二阶微扰理论MP2方法研究了不同氧化态的金催化剂催化2-炔丙基苯酮与苯炔环化反应的机理.结果表明:在AuCl和AuCl3的作用下,反应均能通过[4+2]和[3+2]途径生成产物.但在AuCl催化下,[4+2]反应途径比[3+2]反应途径具有更低的活化自由能,反应主要通过[4+2]途径进行;而在AuCl3催化下,[4+2]反应途径与[3+2]反应途径的活化自由能相近,反应通过两条途径竞争生成产物.比较AuCl和AuCl3的催化效果发现,不同氧化态的金催化剂改变了反应的机理,该反应的活化自由能在AuCl作用下比在AuCl3作用下低11.18 kJ·mol-1.对于该反应,AuCl表现出了更好的催化活性.这些计算结果和实验现象相吻合.     

Structures and Electronic Properties of HOCl···HCOCl Complexes

B3LYP/6-311++G** and MP2/6-311++G** calculations were used to analyze the interaction between hypochlorous acid (HOCl) and formyl chloride (HCOCl). The results showed that there were four equilibrium geometries (S1, S2, S3, and S4) optimized at B3LYP/6-311++G** level, and all the equilibrium geometries were confirmed to be in stable states by analytical frequency calculations. Complexes S1 and S3 use the 5H atom of HOCl as proton donor and the terminal 1O atom of HCOCl as acceptor to form red shift hydrogen bond systems. However, the blue-shifted hydrogen bond (2C-3H···6O) coexists with 4Cl···5O interaction in structures S2. As for S4, it uses the 7Cl atom of HOCl as proton donor and the terminal 1O atom of HCOCl as acceptor to form red shift halogen bond system. Interaction energies between monomers in the four complexes corrected with basis set superposition error (BSSE) and zero-point vibrational energy (ZPVE) lie in the range from −5.05 to −14.76 kJ·mol⁻¹ at MP2/6-311++G** level. The natural bond orbital (NBO) and atoms in molecules (AIM) theories have also been applied to explain the structures and the properties of the complexes.     

Theoretical study of light-metal tetrahydroborates

Ab initio calculations of the structural, energetic, vibrational, and magnetic characteristics of the lowest-lying structures for isolated molecules and ions of light-metal tetrahydroborates (Li, Na, Be, Mg, and Al) have been performed by the perturbation theory (MP2), quadratic configuration interaction (QCISD(T)), coupled cluster (CCSD(T)), and density functional theory (B3LYP) methods using the 6-31G*, 6-31G**, 6-311+G**, and 6-311++G** basis sets. The trends in the behavior of the structural parameters, the energies of different decomposition pathways, barriers to internal rotation of BH₄ groups, normal mode frequencies, magnetic shielding constants, and spin density distribution (in radicals) have been analyzed in various related series of these compounds. The results obtained by ab initio methods and at the DFT level are compared. The economical approximation B3LYP/6-311++G**//B3LYP/6-311+G** adequately reproduces the results obtained at the higher level of theory CCSD(T)/6-311++G*s*//MP2/6-31G* even though it requires considerably shorter CPU times and smaller amounts of memory.