Theoretical Study on Intermolecular Interactions and Thermodynamic Properties of Nitroamine Dimers

Ab initio self-consistent field(SCF) and Mφller-Plesset correlation correction methods employing 6-31G^** basis set have been applied to the optimizations of nitroamine dimers.The binding energies have been corrected for the basis set superposition error (BSSE) and the zero-point energy.Theree optimized dimers have been obtained.The BSSE corrected binding energy of the most stable dimer is predicted to be -31.85kJ/mol at the MP4/6-31G^**//MP2/6-31G^** level.The energy barriers of the Walden conversion for -NH2 group are 19.7kJ/mol and 18.3kJ/mol for monomer and the most stable dimer,respectively.The molecular interaction makes the internal rotation around N1-N2 even more difficult.The thermodynamic properties of nitroamine and its dimers at different temperatures have been calculated on the basis of vibrational analyses.The change of the Gibbs free energy for the aggregation from monomer to the most stable dimer at standard pressure and 298.2 K is predicted to be 14.05kJ/mol.     

A theoretical study on the intermolecular interaction of energetic system-Nitromethane dimer

Three optimized geometries of nitromethane dimer have been obtained at the HF/6-31G level.Dimer binding energies have been corrected for the basis set superposition error (BSSE) and the zero point energy.Computed results indicate that the cyclic structure of (CH3NO2)2 is the most stable of three optimized geometries,whose corrected binding energyis 17.29 kJ mol-1 at the MP4SDTQ/6-31G//HF/6-31G level.In the optimized structures of nitromethane dimer,the inter-molecular hydrogen bond has not been found; and the charge-transfer interaction between CH3NO2 subsystems is weak; and the correlation interaction energy makes a little contribution to the intermolecular interaction energy of the dimer.     

硝酸甲酯分子间相互作用的DFT和ab initio比较

用密度泛函理论(DFT)和从头算(ab initio)方法,分别在B3LYP/6 31G和HF/6 31G水平上求得硝酸甲酯三种二聚体的全优化几何构型和电子结构,并用6 311G和6 311＋＋G基组进行总能量计算.对HF/6 31G计算结果进行MP4SDTQ电子相关校正.在各基组下均进行基组叠加误差(BSSE)和零点能(ZPE)校正求得结合能.对6 31G优化构型作振动分析并基于统计热力学求得200～600 K温度下单体和二聚体的热力学性质.详细比较两种方法的相应计算结果,发现DFT求得的分子间距离较短,分子内键长较长,所得结合能均小于相应ab initio计算值.     

Pi-pi interaction in pyridine

pi-pi Interaction in pyridine dimer and trimer has been investigated in different geometries and orientations at the ab initio (HF, MP2) and DFT (B3LYP) levels of theory using various basis sets (6-31G, 6-31G, 6-311++G) and corrected for basis set superposition error (BSSE). While the HF and DFT calculations show the pyridine dimer and the trimer to be unstable with respect to the monomer, the MP2 calculations show them to be clearly stable, thus emphasizing the need to include electron correlation while determining stacking interaction in such systems. The calculated MP2/6-311++G binding energy (100% BSSE corrected) of the parallel-sandwich, antiparallel-sandwich, parallel-displaced, antiparallel-displaced, T-up and T-down geometries for pyridine dimer are 1.53, 3.05, 2.39, 3.97, 1.91, 1.47 kcal/mol, respectively. The results show the antiparallel-displaced geometry to be the most stable. The binding energies for the trimer in parallel-sandwich, antiparallel-sandwich, and antiparallel-displaced geometry are found to be 3.18, 6.14, and 8.04 kcal/mol, respectively.     

Theoretical study on the hydrogen bonding of five-membered heteroaromatics with water

The hydrogen-bonding ability of five-membered heteroaromatic molecules containing one chalcogen and two heteroatoms with nitrogen in addition to chalcogen, respectively, have been analyzed using density functional and molecular orbital methods through adduct formation with water. The stabilization energies for all the adducts are established at B3LYP/6-31+G* and MP2/6-31+G* levels after correcting for the basis set superposition error by using the counterpoise method and also corrected for zero-point vibrational energies. A natural bond orbital analysis at B3LYP/6-31+G* level and natural energy decomposition analysis at HF/6-31+G* using MP2/6-31+G* geometries have been carried out to understand the nature of hydrogen-bonding interaction in monohydrated heterocyclic adducts. Nucleus-independent chemical shift have been evaluated to understand the correlation between hydrogen bond formation and aromaticity.     

Accuracy of the three-body fragment molecular orbital method applied to Møller-Plesset perturbation theory

The three‐body energy expansion in the fragment molecular orbital method (FMO) was applied to the 2nd order Møller–Plesset theory (MP2). The accuracy of both the two and three‐body expansions was determined for water clusters, alanine n‐mers (α‐helices and β‐strands) and one synthetic protein, using the 6‐31G* and 6‐311G* basis sets. At the best level of theory (three‐body, two molecules/residues per fragment), the absolute errors in energy relative to ab initio MP2 were at most 1.2 and 5.0 mhartree, for the 6‐31G* and 6‐311G* basis sets, respectively. The relative accuracy was at worst 99.996% and 99.96%, for 6‐31G* and 6‐311G*, respectively. A three‐body approximation was introduced and the optimum threshold value was determined. The protein calculation (6‐31G*) at the production level (FMO2/2) took 3 h on 36 3.2‐GHz Pentium 4 nodes and had the absolute error in the MP2 correlation energy of only 2 kcal/mol. © 2007 Wiley Periodicals, Inc. J Comput Chem, 2007     

硝酸乙酯分子间相互作用的ab initio研究

在abinitio-HF/6-31G水平上求得硝酸乙酯二聚体势能面上的四种优化构型和电子结构。经MP2电子相关校正和基组叠加误差(BSSE)以及零点能(ZPE)校正，求得二聚体的最大结合能为11.46kJ.mol^-^1，还进行HF/6-311G和HF/6-311++G水平的总能量比较计算，发现6-31G基组对计算结合能比较适合，二子体系间的电荷转移很少，对优化构型进行振动分析，并基于统计热力学求得从单体形成二聚体的热力学性质变化。     

Theoretical Study on Intermolecular Interactions and Thermodynamic Properties of Difluoroamine Complex

Introduction　　ApplicationsofabinitiocalculationstointermolecularinteractionsincludingeitherweakvanderWaalsorstrongerhydrogenbondinghavedrawnmuchattentioninthepastdecadesbecausetheyareimportantinawiderangeofphysical,chemicalandbiologicalfields .1 5Inre centyears ,wehaveappliedtheintermolecularinteractionstoenergeticsystemsandobtainedsomemeaningfulinfor mationthatisvaluableforthestudyofenergeticmateri als .6 14 Thebehaviorofmolecularcomplexesisusuallybe tweentwoextremes :thegasphaseandthecrys…     

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.     

Density functional theory study of formaldehyde oligomers

Hydrogen‐bonded formaldehyde oligomers (dimer to pentamer) are studied using density functional theory (DFT), the B3LYP method, and the 6‐311+G* basis set. Many‐body interaction energies are obtained to study the contribution of many‐body terms to binding energy. The basis set superposition error (BSSE)‐corrected total energies are ?229.08170, ?343.61410, ?458.16660, and ?572.70901 hartrees for dimer, trimer, tetramer, and pentamer, respectively, with corresponding binding energies ?2.55, ?4.86, ?6.99, and ?9.49 kcal/mol. Two‐body energies have been found to contribute significantly to the total binding energy in dimer to pentamer, whereas higher‐order interaction energies are negligible. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Comparison of levels of electronic structure theory in direct dynamics simulations of C2H5F --> HF + C2H4 product energy partitioning

Direct dynamics simulations at the MP2/6-311++G** level of theory were performed to study C(2)H(5)F --> HF + C(2)H(4) product energy partitioning. The simulation results are compared with experiment and a previous MP2/6-31G* simulation. The current simulation with the larger basis set releases more energy to HF vibration and less to HF + C(2)H(4) relative translation as compared to the previous simulation with the 6-31G* basis set. The HF rotation and vibration energy distributions determined from the current simulation are in overall very good agreement with previous experimental studies of C(2)H(5)F dissociation by chemical activation and IRMPA. A comparison of the simulations with experiments suggests there may be important mass effects for energy partitioning in HX elimination from haloalkanes. The transition state (TS) structures and energies calculated with MP2 and the 6-31G* and 6-311++G** basis sets are compared with those calculated using CCD, CCSD, CCSD(T), and the 6-311++G** basis set.     

Theoretical studies on the conformations and properties of phenylazonaphthalenes

The conformations of phenylazo‐2‐naphthalene (I) and phenylazo‐1‐naphthalene (II) have been studied using ab initio methods and density functional theory. The rotational potential energy surfaces about the C N bonds were calculated for both the trans and the cis forms at the PM3 and HF/STO‐3G levels. The PM3 method was found to be incapable for locating the energy minima correctly. The geometries of rotamers obtained from the HF/STO‐3G surface were fully optimized at the HF/6‐31G* and B3LYP/6‐31G* levels. Single‐point MP2/6‐31G* calculations have also been carried out at the HF/6‐31G* geometries. The vibrational spectra, standard thermodynamic functions, heats of formation, and equilibrium molar fractions have been obtained. According to the calculated results, I is comparatively more stable than II. The trans forms have much lower energies than the cis forms, implying that I and II mainly exist in the trans planar forms. The energy gap (ΔE) between the frontier molecular orbitals (MOs) of the cis forms are about 0.7 eV higher than those of the trans ones, indicating that change in geometry from trans to cis form will cause a hypsochromic effect. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 79: 25–33, 2000     

A quantum mechanical approach to the kinetics of the hydrogen abstraction reaction H2O2 + •OH → HO2 + H2O

The kinetics of the hydrogen abstraction from H₂O₂ by ^?OH has been modeled with MP2/6‐31G*//MP2/6‐31G*, MP2‐SAC//MP2/6‐31G*, MP2/6‐31+G**//MP2/6‐31+G**, MP2‐SAC// MP2/6‐31+G**, MP4(SDTQ)/6‐311G**//MP2/6‐31G*, CCSD(T)/6‐31G*//CCSD(T)/6‐31G*, CCSD(T)/6‐31G**//CCSD(T)/6‐31G**, CCSD(T)/6‐311++G**//MP2/6‐31G* in the gas phase. MD simulations have been used to generate initial geometries for the stationary points along the potential energy surface for hydrogen abstraction from H₂O₂. The effective fragment potential (EFP) has been used to optimize the relevant structures in solution. Furthermore, the IEFPCM model has been used for the supermolecules generated via MD calculations. IEFPCM/MP2/6‐31G* and IEFPCM/CCSD(T)/6‐31G* calculations have also been performed for structures without explicit water molecules. Experimentally, the rate constant for hydrogen abstraction by ^?OH drops from 1.75 × 10^?12 cm³ molecule^?1 s^?1 in the gas phase to 4.48 × 10^?14 cm³ molecule^?1 s^?1 in solution. The same trend has been reproduced best with MP4 (SDTQ)/6‐311G**//MP2/6‐31G* in the gas phase (0.415 × 10^?12 cm³ molecule^?1 s^?1) and with EFP (UHF/6‐31G*) in solution (3.23 × 10^?14 cm³ molecule^?1 s^?1). © 2005 Wiley Periodicals, Inc. Int J Chem Kinet 37: 502–514, 2005     

Conformational analysis of carbonic anhydrase inhibitors using ab initio molecular orbital methods. 1. Rotational isomerism in methane sulfonamide anion, CH3-SO2-NH−

Ab initio calculations have been performed on methane sulfonamide anion. Geometries have been optimized using Hartree-Fock basis sets up to 6-31+G*, and single-point calculations employing those Hartree-Fock geometries have been performed at levels up to MP2/6-311++G**. In addition, geometry optimizations for the 0°, 90°, 150°, and 180° conformers have been carried out at the MP2/6-31G*, MP2/6-31+G*, and MP2/6-311++G** levels. Vibrational frequencies have been calculated using the HF/4-31G*, MP2/6-31G*, and MP2/6-31+G* geometries. All calculations at or above the 4-31G* level agree that H—N—S—C ˜90° is the global minimum. The H—N—S—C = 180° conformer is clearly higher in energy although the relative energy of this conformer varies from 0.36 to 1.03 kcal/mol for the post-HF calculations depending on basis set. The H—N—S—C = 180° conformer appears to be a very shallow local minimum. However, the potential energy surface is quite flat in this region, and the highest-level calculations, including MP2 optimizations and vibrational frequency analysis, are ambiguous on this point. The conformer with an H—N—S—C torsion of 0° is a transition state with a relative energy ˜8 kcal/mol. Received: 3 December 1996 / Accepted: 2 January 1997     

Ab initio Study on the Intermolecular Interaction and Thermodynamic Properties of Methyl Nitrate Dimer

Three stable dimers of methyl nitrate have been obtained and their geometries have been fully optimized at the HF/6‐31G,. level. Binding energies have been calculated with correction for the basis set superposition error (BSSE) and zero point energy (ZPE). The cyclic overlap‐type structure, the binding energy of which is 11.97 kJ/mol at the MP4SDTQ/6‐31G. / HF/6‐31G. level, is the most stable. No intermolecular hydrogen bond was found, and the charge transfer between two subsystems is minute. The thermodynamic properties of methyl nitrate and its dimers have been calculated based on the vibrational analysis and statistical thermodynamics.     

Theoretical Study on Intermolecular Interactions and Thermodynamic Properties of Dimethylnitroamine Clusters

Ab initio SCF and Mφller-Plesset correlation correction methods in combination with counterpose procedure for BSSE correction have been applied to the theroetical studying of dimethylnitroamine and its dimers and trimers.Three optimized stable dimers and two trimers have been obtained.The corrected binding energies of the most stable dimer and trimer were predicted to be -24.68kJ/mol and -47.27kJ/mol,respectively at the MP2/6-31G^*//HF/6-31G^* level.The proportion of correlated interation energies to their total interaction energies for all clusters was at least 29.3 percent,and the BSSE of ΔE(MP2) was at least 10.0kJ/mol.Dispersion and/or electrostatic force were dominant in all clusters.There exist cooperative effects in both the chain and the cyclic trimers.The vibrational frequencies associated with N-O stretches or wags exhibit slight red shifts,but the modes associated with the motion of hydrogen atoms of the methyl group show somewhat blue shifts with respect to those of monomer.Thermodynamic properties of dimethylnitroamine and its clusters at different temperatures have been calculated on the basis of vibrational analyses.The changes of the Gibbs free energies for the aggregation from monomer to the most stable dimer and trimer were predicted to be 14.37kJ/mol and 30.40kJ/mol,respectively,at 1 atm and 298.15K.     

Ab initio study of reactions of hydroxyl radicals with chloro‐ and fluoro‐substituted methanes

Ab initio calculations at the unrestricted Hartree–Fock (UHF) level have been performed to investigate the hydrogen abstraction reactions of ^? OH radicals with methane and nine halogen‐substituted methanes (F, Cl). Geometry optimization and vibrational frequency calculations have been performed on all reactants, adducts, products, and transition states at the UHF/6‐31G* level. Single‐point energy calculations at the MP2/6‐31++G* level using the UHF/6‐31G* optimized geometries have also been carried out on all species. Pre‐ and postreaction adducts have been detected on the UHF/6‐31G* potential energy surfaces of the studied reactions. Energy barriers, ΔE^?, reaction energies, ΔE_r, reaction enthalpies, ΔH_r, and activation energies, E_a, have been determined for all reactions and corrected for zero‐point energy effects. Both E_a and ΔH_r come into reasonable agreement with the experiment when correlation energy is taken into account and when more polarized and diffuse basis sets are used. The E_a values, estimated at the PMP2/6‐31++G* level, are found to be in good agreement with the experimental ones and correctly reproduce the experimentally observed trends in fluorine and chlorine substitution effects. A linear correlation between E_a and ΔH_r is obtained, suggesting the presence of an Evans–Polanyi type of relationship. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem 84: 426–440, 2001     

Calculated structures and relative stabilities of furoxan,some 1,2-dinitrosoethylenes and other isomers

The structures and relative stabilities of furoxan and some of its isomers, e.g., the 1,2-dinitrosoethylenes, have been determined by means of ab initio Hartee–Fock and Møller–Plesset calculations. Geometries were optimized at the HF/3-21G, HF/6-31G* and MP2/6-31G* levels, and subsequently used for computing MP2/6-31G*, MP3/6-31G*, and MP4/6-31G* energies. The results are markedly affected by the inclusion of electronic correlation, which renders three of the isomers unstable. It also emphasizes the importance of a zwitterionic contribution to the structure of furoxan, which promotes ring-opening through a cis 1,2-dinitrosoethylene intermediate/transition state that has an MP4/6-31G*//MP2/6-31G* energy that is 31.6 kcal/mol above furoxan.     

Molecular dynamics simulations of fluid methane properties using ab initio intermolecular interaction potentials

Intermolecular interaction energy data for the methane dimer have been calculated at a spectroscopic accuracy and employed to construct an ab initio potential energy surface (PES) for molecular dynamics (MD) simulations of fluid methane properties. The full potential curves of the methane dimer at 12 symmetric conformations were calculated by the supermolecule counterpoise‐corrected second‐order Møller‐Plesset (MP2) perturbation theory. Single‐point coupled cluster with single and double and perturbative triple excitations [CCSD(T)] calculations were also carried out to calibrate the MP2 potentials. We employed Pople's medium size basis sets [up to 6‐311++G(3df, 3pd)] and Dunning's correlation consistent basis sets (cc‐pVXZ and aug‐cc‐pVXZ, X = D, T, Q). For each conformer, the intermolecular carbon–carbon separation was sampled in a step 0.1 Å for a range of 3–9 Å, resulting in a total of 732 configuration points calculated. The MP2 binding curves display significant anisotropy with respect to the relative orientations of the dimer. The potential curves at the complete basis set (CBS) limit were estimated using well‐established analytical extrapolation schemes. A 4‐site potential model with sites located at the hydrogen atoms was used to fit the ab initio potential data. This model stems from a hydrogen–hydrogen repulsion mechanism to explain the stability of the dimer structure. MD simulations using the ab initio PES show quantitative agreements on both the atom‐wise radial distribution functions and the self‐diffusion coefficients over a wide range of experimental conditions. © 2008 Wiley Periodicals, Inc. J Comput Chem 2009     

乙烯、乙炔和苯与氨基锂间的π型锂键相互作用

分别在DFT-B3LYP和MP2/6-311++G**水平上求得乙烯、乙炔和苯与氨基锂锂键复合物势能面上的3个几何构型. 频率分析表明,3个构型均为稳定构型. 计算结果表明,形成锂键复合物后,质子供体N(2)- Li(4)的键长明显增大,且其伸缩振动的频率发生了不同程度的红移. 分别在乙烯…氨基锂、乙炔…氨基锂和苯…氨基锂三种复合物中,经MP2/6-311++G**水平计算的同时含基组重叠误差(BSSE)和零点振动能校正的单体间锂键相互作用能分别为-26.04、-24.86 和 -30.02 kJ·mol-1. 自然键轨道理论（NBO）分析表明单体间的弱相互作用属于π-s型锂键.