The correlation-consistent composite approach: application to the G3/99 test set

The correlation-consistent composite approach (ccCA), an ab initio composite technique for computing atomic and molecular energies, recently has been shown to successfully reproduce experimental data for a number of systems. The ccCA is applied to the G3/99 test set, which includes 223 enthalpies of formation, 88 adiabatic ionization potentials, 58 adiabatic electron affinities, and 8 adiabatic proton affinities. Improvements on the original ccCA formalism include replacing the small basis set quadratic configuration interaction computation with a coupled cluster computation, employing a correction for scalar relativistic effects, utilizing the tight-d forms of the second-row correlation-consistent basis sets, and revisiting the basis set chosen for geometry optimization. With two types of complete basis set extrapolation of MP2 energies, ccCA results in an almost zero mean deviation for the G3/99 set (with a best value of -0.10 kcal mol(-1)), and a 0.96 kcal mol(-1) mean absolute deviation, which is equivalent to the accuracy of the G3X model chemistry. There are no optimized or empirical parameters included in the computation of ccCA energies. Except for a few systems to be discussed, ccCA performs as well as or better than Gn methods for most systems containing first-row atoms, while for systems containing second-row atoms, ccCA is an improvement over Gn model chemistries.     

B3LYP and MP2 calculations of the enthalpies of hydrogen-bonded complexes of methanol with neutral bases and anions: comparison with experimental data

To perfect a method for building a theoretical hydrogen-bond basicity scale, the enthalpy of hydrogen bonding between methanol and thirteen neutral and anionic bases (MeOH, MeNH2, Me2NH, Et2NH, Me3N, Et3N, Br-, CN-, SH-, Cl-, HCOO-, MeO-, F-) was calculated by DFT and ab initio methods. The theoretical results were compared to selected experimental ones. It appears that B3LYP/6-31+G(d,p) calculations are satisfactory for optimizing the geometry of complexes and giving a general order of basicity. However, they are deficient for reproducing the large effect of alkyl groups on the hydrogen-bond basicity of amines. This deficiency is explained by intermolecular perturbation theory calculations, which show that the alkylation of nitrogen dramatically increases the dispersion energy component not taken into account by the B3LYP functional. Of the methods considered, only MP2/aug-cc-pVTZ calculations are capable of reproducing the binding enthalpy within the experimental error for the first-row acceptor atoms N, O, and F, and of accounting for dispersion effects created by alkylation at the hydrogen-bond acceptor site.     

Ab initio study of hydrogen-bond formation between cyclic ethers and selected amino acid side chains

Binding energies for hydrogen-bonded complexes of six cyclic ethers with five hydrogen-bond donor molecules that mimic selected amino acid side chains have been calculated at the MP2/6-31G*, MP2/6-31+G*, MP2/6-311++G**(single point), and MP2/aug-cc-pvtz levels, using geometries obtained with or without counterpoise corrections throughout the geometry optimization. The calculated basis set superposition error (BSSE) amounts to 10-20% and 5-10% of the uncorrected binding energies for the neutral and ionic species, respectively, at the MP2/aug-cc-pvtz level. The authors conclude that the O...H distances in the hydrogen bonds and binding energies for the studied systems may be determined with uncertainties of up to 0.08 A and 1-2 kcal/mol, respectively, in comparison with the MP2/aug-cc-pvtz values at a reasonable computational cost by performing standard geometry optimization at the MP2/6-31+G* level. Hydrogen-bond formation energies are more negative for cyclic ethers compared to their counterparts with a C=C double bond in the ring next to the oxygen atom. The less negative hydrogen-bonding energy and the increased O...H separation have been attributed to the reduced basicity of the ether oxygen when the lone pairs can enter conjugation with the pi-electrons of the Calpha=Cbeta double bond. The present study is the first step toward the development of an affordable computational level for estimating the binding energies of natural product, fused ring ether systems to the human estrogen receptor.     

Proton-bound homodimers: how are the binding energies related to proton affinities?

High-level quantum chemical calculations [G3(MP2)-RAD//MP2/6-31+G(d,p)] have been employed to investigate the relationship between the binding energy (BE) of a substrate (X) and its protonated form [H-X]+ with the proton affinity (PA) of the substrate (X) in several series of protonated homodimers ([X...H-X]+). We find that for each series of closely related substrates, the binding energy (BE) is correlated with the proton affinity (PA) in an approximately quadratic manner. Thus, for a given series, the BE initially increases in magnitude with increasing PA, reaches a point of maximum binding, and then becomes smaller as the PA increases further. This behavior can be attributed to the competing effects of the exothermic partial protonation of the substrate and the endothermic partial deprotonation of the protonated substrate. As the PA increases, protonation of X contributes to increased binding but the penalty for partial deprotonation of [H-X]+ also increases. Once the PA becomes sufficiently high, the penalty for the partial deprotonation of [H-X]+ dominates, leading to maximum binding occurring at intermediate PA.     

Accurate atomic Gaussian basis functions for second-row atoms: Small split-valence 3-21SP and 4-22SP basis sets

Small split-valence Gaussian 3-21SP and 4-22SP basis sets, previously reported for the first-row atoms [Chem. Phys. Lett., 229 , 151 (1996)], have been extended for the second-row elements of the Periodic Table. The total energies of the ground states of the second-row atoms calculated with the new basis sets are significantly lower than those obtained with the well-known 3-21G (J. Am. Chem. Soc., 104 , 2797 (1982)] and 4-31G [J. Chem. Phys., 56 , 5255 (1972)] basis sets. This is because, as first noted in our previous work for first-row atoms, that the 3-21G and 4-31G basis sets only correspond to a local minimum of the Hartree–Fock energy functional, which is relatively far from its global minimum. The proposed basis sets have been tested by performing geometry optimizations and calculations of normal frequencies in the harmonic approximation of some diatomic and polyatomic molecules at the Hartree–Fock level. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1200–1210     

Toward true DNA base-stacking energies: MP2, CCSD(T), and complete basis set calculations

Stacking energies in low-energy geometries of pyrimidine, uracil, cytosine, and guanine homodimers were determined by the MP2 and CCSD(T) calculations utilizing a wide range of split-valence, correlation-consistent, and bond-functions basis sets. Complete basis set MP2 (CBS MP2) stacking energies extrapolated using aug-cc-pVXZ (X = D, T, and for pyrimidine dimer Q) basis sets equal to -5.3, -12.3, and -11.2 kcal/mol for the first three dimers, respectively. Higher-order correlation corrections estimated as the difference between MP2 and CCSD(T) stacking energies amount to 2.0, 0.7, and 0.9 kcal/mol and lead to final estimates of the genuine stacking energies for the three dimers of -3.4, -11.6, and -10.4 kcal/mol. The CBS MP2 stacking-energy estimate for guanine dimer (-14.8 kcal/mol) was based on the 6-31G(0.25) and aug-cc-pVDZ calculations. This simplified extrapolation can be routinely used with a meaningful accuracy around 1 kcal/mol for large aromatic stacking clusters. The final estimate of the guanine stacking energy after the CCSD(T) correction amounts to -12.9 kcal/mol. The MP2/6-31G(0.25) method previously used as the standard level to calculate aromatic stacking in hundreds of geometries of nucleobase dimers systematically underestimates the base stacking by ca. 1.0-2.5 kcal/mol per stacked dimer, covering 75-90% of the intermolecular correlation stabilization. We suggest that this correction is to be considered in calibration of force fields and other cheaper computational methods. The quality of the MP2/6-31G(0.25) predictions is nevertheless considerably better than suggested on the basis of monomer polarizability calculations. Fast and very accurate estimates of the MP2 aromatic stacking energies can be achieved using the RI-MP2 method. The CBS MP2 calculations and the CCSD(T) correction, when taken together, bring only marginal changes to the relative stability of H-bonded and stacked base pairs, with a slight shift of ca. 1 kcal/mol in favor of H-bonding. We suggest that the present values are very close to ultimate predictions of the strength of aromatic base stacking of DNA and RNA bases.     

An ab initio molecular orbital study of the structures and energies of neutral and charged bimolecular complexes of NH3 with the hydrides AHn (A = N,O, F,P, S,and Cl)

Hartree-Fock 6-31G(d) structures for the neutral, positive ion, and negative ion bimolecular complexes of NH₃ with the first- and second-row hydrides AH_n (AH_n = NH₃, OH₂, FH, PH₃, SH₂, and ClH) have been determined. All of the stable neutral complexes except (NH₃)₂, the positive ion complexes with NH₃ as the proton acceptor, and the negative ion complexes containing first-row anions exhibit conventional hydrogen bonded structures with essentially linear hydrogen bonds and directed lone pairs of electrons. The positive ion complex NH₄⁺ …? OH₂ has the dipole moment vector of H₂O instead of a lone pair directed along the intermolecular line, while the complexes of NH₄⁺ with SH₂, FH, and ClH have structures intermediate between the lone-pair directed and dipole directed forms. The negative ion complexes containing second-row anions have nonlinear hydrogen bonds. The addition of diffuse functions on nonhydrogen atoms to the valence double-split plus polarization 6-31G(d,p) basis set usually decreases the computed stabilization energies of these complexes. Splitting d polarization functions usually destabilizes these complexes, whereas splitting p polarization functions either has no effect or leads to stabilization. The overall effect of augmenting the 6-31G(d,p) basis set with diffuse functions on nonhydrogen atoms and two sets of polarization functions is to lower computed stabilization energies. Electron correlation stabilizes all of these complexes. The second-order Møller–Plesset correlation term is the largest term and always has a stabilizing effect, whereas the third and fourth-order terms are smaller and often of opposite sign. The recommended level of theory for computing the stabilization energies of these complexes is MP2/6-31+G(2d,2p), although MP2/6-31+G(d,p) is appropriate for the negative ion complexes.     

Rapid evaluation of the binding energies between peptide amide and DNA base

The binding energies and the equilibrium hydrogen bond distances as well as the potential energy curves of 20 hydrogen‐bonded amide–base dimers are evaluated from the analytic potential energy function established in our laboratory recently. The analytic potential energy function is used to calculate the N? H···N, N? H···O?C, C? H···N, and C? H···O?C dipole–dipole attractive interaction energies and C?O···O?C, N? H···H? N, and N? H···H? C dipole–dipole repulsive interaction energies in the 20 dimers composed of DNA bases adenine, guanine, cytosine, or thymine and peptide amide. The calculation results show that the potential energy curves obtained from the analytic potential energy function are in good agreement with those obtained from MP2/6‐311+G** calculations by including the basis set superposition error (BSSE) correction. For all the 20 dimers, the analytic potential energy function yields the binding energies of the MP2/6‐311+G** with BSSE correction within the error limits of 0.50 kcal/mol for 19 dimers, only one difference is larger than 0.50 kcal/mol and the difference is only 0.61 kcal/mol. The analytic potential energy function produces the equilibrium hydrogen bond distances of the MP2/6‐311+G** with BSSE correction within the error limits of 0.030 Å for all the 20 dimers. The analytic potential energy function is further applied to four more complicated DNA base‐peptide amide systems involving amino acid side chain and β‐sheet. The values of the binding energies and equilibrium hydrogen bond distances obtained from the analytic potential energy function are also in good agreement with those obtained from MP2 calculations with the BSSE correction. These results demonstrate that the analytic potential energy function can be used to evaluate the binding energies in hydrogen‐bonded peptide amide–DNA base dimers quickly and accurately. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Gaussian-4 theory

The Gaussian-4 theory (G4 theory) for the calculation of energies of compounds containing first- (Li-F), second- (Na-Cl), and third-row main group (K, Ca, and Ga-Kr) atoms is presented. This theoretical procedure is the fourth in the Gaussian-n series of quantum chemical methods based on a sequence of single point energy calculations. The G4 theory modifies the Gaussian-3 (G3) theory in five ways. First, an extrapolation procedure is used to obtain the Hartree-Fock limit for inclusion in the total energy calculation. Second, the d-polarization sets are increased to 3d on the first-row atoms and to 4d on the second-row atoms, with reoptimization of the exponents for the latter. Third, the QCISD(T) method is replaced by the CCSD(T) method for the highest level of correlation treatment. Fourth, optimized geometries and zero-point energies are obtained with the B3LYP density functional. Fifth, two new higher level corrections are added to account for deficiencies in the energy calculations. The new method is assessed on the 454 experimental energies in the G305 test set [L. A. Curtiss, P. C. Redfern, and K. Raghavachari, J. Chem. Phys. 123, 124107 (2005)], and the average absolute deviation from experiment shows significant improvement from 1.13 kcal/mol (G3 theory) to 0.83 kcal/mol (G4 theory). The largest improvement is found for 79 nonhydrogen systems (2.10 kcal/mol for G3 versus 1.13 kcal/mol for G4). The contributions of the new features to this improvement are analyzed and the performance on different types of energies is discussed.     

Rapid evaluation of the binding energies in hydrogen-bonded amide-thymine and amide-uracil dimers in gas phase

The binding energies and the equilibrium hydrogen bond distances as well as the potential energy curves of 48 hydrogen‐bonded amide–thymine and amide–uracil dimers are evaluated from the analytic potential energy function established in our lab recently. The calculation results show that the potential energy curves obtained from the analytic potential energy function are in good agreement with those obtained from MP2/6‐311+G** calculations by including the BSSE correction. For all the 48 dimers, the analytic potential energy function yields the binding energies of the MP2/6‐311+G** with BSSE correction within the error limits of 0.50 kcal/mol for 46 dimers, only two differences are larger than 0.50 kcal/mol and the largest one is only 0.60 kcal/mol. The analytic potential energy function produces the equilibrium hydrogen bond distances of the MP2/6‐311+G** with BSSE correction within the error limits of 0.050 Å for all the 48 dimers. The analytic potential energy function is further applied to four more complicated hydrogen‐bonded amide–base systems involving amino acid side chain and β‐sheet. The values of the binding energies and equilibrium hydrogen bond distances obtained from the analytic potential energy function are also in good agreement with those obtained from MP2 calculations with the BSSE correction. These results demonstrate that the analytic potential energy function can be used to evaluate the binding energies in hydrogen‐bonded amide–base dimers quickly and accurately. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

A computational study of intramolecular proton transfer in gaseous protonated glycine

The minimum energy paths for intramolecular proton transfer between the amino nitrogen and carbonyl oxygen atoms in gaseous protonated glycine were estimated at the Hartree-Fock (HF) and second-order M?ller-Plesset Perturbation (MP2) levels of theory. Potential energy profiles and their associated reactant, transition state, and product species calculated at the MP2/6-31G* level were shown to differ significantly from those obtained at the HF/6-31G* level. Effects of electron correlation and basis functions on the calculated geometries and energies of relevant species were examined at the HF, MP2, MP4, CCSD, and B3LYP levels using the 6-31G*, 6-31G**, 6-31+G**, 6-311+G**, 6-31+G(2d,2p), 6-311+G(3df,2p), cc-pVDZ, aug-cc-pVDZ, and cc-pVTZ basis sets. The HF and MP2 optimized levels with the 6-31G*, 6-31G**, 6-31+G**, and 6-311+G** bases were used to calculate the thermodynamic and kinetic properties of the proton transfer reaction at 298.15 K and 1 atm, which include enthalpy, entropy, Gibbs free energy, equilibrium constant, potential energy barriers, tunneling transmission coefficients, and rate constants. Results indicate that the proton in a carbonyl O-protonated glycine undergoes a rapid migration to the amino nitrogen atom, while the reverse process is extremely unfavorable. The objective of this work is to develop practical theoretical procedures for studying proton transfer reactions in amino acids and peptides and to assemble physical data from these model calculations for future references.     

羟基二酰亚胺二聚体分子间相互作用的量子化学研究

采用从头算MP2方法在6-311++G~**基组水平上讨论了CP梯度校正对两种羟基 二酰亚胺异构体所的相互作用能和几何结构的影响，并利用分子中的原子理论（ Atoms in molecules，AIM）计算了五个拓扑参数：键临荷密度、电荷密度的 Iaplacian值、氢键中氢原子的体积、氢原子集居数、氢原子能量来表征氢键的形 成．种构型氢键体系中还讨论了二聚体的相互作用能与氢键临界点的电荷密度、质 子供体X-H键长的线性相关性问题．表明这种线性相关性的存在有范围限制，复合 物和其中单体的构型能够影响这种关系的存在.     

Binding energies of the proton-bound amino Acid dimers gly.gly, ala.ala, gly.ala, and lys.lys measured by blackbody infrared radiative dissociation

Arrhenius activation energies in the zero-pressure limit for dissociation of gas-phase proton-bound homodimers of N,N-dimethylacetamide (N,N-DMA), glycine, alanine, and lysine and the heterodimer alanine.glycine were measured using blackbody infrared radiative dissociation (BIRD). In combination with master equation modeling of the kinetic data, binding energies of these dimers were determined. A value of 1.25 +/- 0.05 eV is obtained for N,N-DMA and is in excellent agreement with that reported in the literature. The value obtained from the truncated Boltzmann model is significantly higher, indicating that the assumptions of this model do not apply to these ions. This is due to the competitive rates of photon emission and dissociation for these relatively large ions. The binding energies of the amino acid dimers are ~1.15 +/- 0.05 eV and are indistinguishable despite the difference in their gas-phase basicity and structure. The threshold dissociation energies can be accurately modeled using a range of dissociation parameters and absorption/emission rates. However, the absolute values of the dissociation rates depend more strongly on the absorption/emission rates. For N,N-DMA and glycine, an accurate fit was obtained using frequencies and transition dipole moments calculated at the ab initio RHF/2-31G* and MP2/2-31G* level, respectively. In order to obtain a similar accuracy using values obtained from AM1 semiempirical calculations, it was necessary to multiply the transition dipole moments by a factor of 3. These results demonstrate that in combination with master equation modeling, BIRD can be used to obtain accurate threshold dissociation energies of relatively small ions of biological interest.     

Revised model core potentials for second-row transition metal atoms from Y to Cd

We have produced new relativistic model core potentials (spdsMCPs) for the second-row transition-metal atoms from Y to Cd treating explicitly 4s and 4p electrons in addition to 4d and 5s electrons in the same manner as for the first-row transition-metal atoms given in [Y. Osanai, M.S. Mon, T. Noro, H. Mori, H. Nakashima, M. Klobukowski, E. Miyoshi, Chem. Phys. Lett. 452 (2008) 210]. Using suitable correlating functions together with the split valence MCP functions, we demonstrate that the present MCP basis sets show reasonable performance in predicting the electronic structures of atoms and molecules, bringing about accurate excitation energies for atoms and reasonable spectroscopic constants for AgH.     

Theoretical study of hydrogen bonding in homodimers and heterodimers of amide, boronic acid, and carboxylic acid, free and in encapsulation complexes

The homodimers and the heterodimers of two amides, two boronic acids, and two carboxylic acids have been calculated in the gas phase and in N,N-dimethylformamide (DMF) and CCl(4) solvents using the DFT (M06-2X and M06-L) and the MP2 methods in conjunction with the 6-31G(d,p) and 6-311+G(d,p) basis sets. Furthermore, their pairwise coencapsulation was studied to examine its effect on the calculated properties of the hydrogen bonds at the ONIOM[M06-2X/6-31G(d,p);PM6], ONIOM[MP2/6-31G(d,p); PM6], and M06-2X/6-31G(d,p) levels of theory. The present work is directed toward the theoretical rationalization and interpretation of recent experimental results on hydrogen bonding in encaptulation complexes [D. Ajami et al. J. Am. Chem. Soc. 2011, 133, 9689-9691]. The calculated dimerization energy (ΔE) values range from 0.74 to 0.35 eV for the different dimers in the gas phase, with the ordering carboxylic homodimers > amide-carboxylic dimers > amide homodimers > boronic-carboxylic dimers > amide-boronic dimers > boronic homodimers. In solvents, generally smaller ΔE values are calculated with only small variations in the ordering. In the capsule, the ΔE values range between 0.67 and 0.33 eV with practically the same ordering as in the gas phase. The calculated % distributions of the encapsulated dimers, taking into account statistical factors, are in agreement with the experimental distribution, where the occurrence of boronic homodimer dominates, even though it is calculated to have the smallest ΔE.     

The simulated ab initio molecular orbital (samo) method: Application To Second-Row Elements

The simulated ab initio molecular orbital (SAMO) method is extended to molecules containing second-row atoms with d-orbitals included in the basis. Results for 1-butanethiol are as successful as previous studies for molecules containing only first-row atoms. The role of hybrid orbitais is critically discussed.     

Ab initio study of hydrogen-bond formation between aliphatic and phenolic hydroxy groups and selected amino acid side chains

Hydrogen bonding was studied in 24 pairs of isopropyl alcohol and phenol as one partner, and water and amino-acid mimics (methanol, acetamide, neutral and protonated imidazole, protonated methylalamine, methyl-guanidium cation, and acetate anion) as the other partner. MP2/6-31+G* and MP2/aug-cc-pvtz calculations were conducted in the gas phase and in a model continuum dielectric environment with dielectric constant of 15.0. Structures were optimized in the gas phase with both basis sets, and zero-point energies were calculated at the MP2/6-31+G* level. At the MP2/aug-cc-pvtz level, the BSSE values from the Boys-Bernardi counterpoise calculations amount to 10-20 and 5-10% of the uncorrected binding energies of the neutral and ionic complexes, respectively. The geometry distortion energy upon hydrogen-bond formation is up to 2 kcal/mol, with the exception of the most strongly bound complexes. The BSSE-corrected MP2/aug-cc-pvtz binding energy of -27.56 kcal/mol for the gas-phase acetate...phenol system has been classified as a short and strong hydrogen bond (SSHB). The CH3NH3+...isopropyl alcohol complex with binding energy of -22.54 kcal/mol approaches this classification. The complete basis set limit (CBS) for the binding energy was calculated for twelve and six complexes on the basis of standard and counterpoise-corrected geometry optimizations, respectively. The X...Y distances of the X-H...Y bridges differ by up to 0.03 A as calculated by the two methods, whereas the corresponding CBS energy values differ by up to 0.03 kcal/mol. Uncorrected MP2/aug-cc-pvtz hydrogen-bonding energies are more negative by up to 0.35 kcal/mol than the MP2/CBS values, and overestimate the CCSD(T)/CBS binding energies generally by up to 5% for the eight studied complexes in the gas phase. The uncorrected MP2/aug-cc-pvtz binding energies decreased (in absolute value) by 11-18 kcal/mol for the ionic species and by up to 5 kcal/mol for the neutral complexes when the electrostatic effect of a polarizable model environment was considered. The DeltaECCSD(T) - DeltaEMP2 corrections still remained close to their gas-phase values for four complexes with 0, +/-1 net charges. Good correlations (R2 = 0.918-0.958) for the in-environment MP2/aug-cc-pvtz and MP2/6-31+G* hydrogen-bonding energies facilitate the high-level prediction of these energies on the basis of relatively simple MP2/6-31+G* calculations.     

Basis sets for geometry optimizations of second-row transition metal inorganic and organometallic complexes

Modest-sized basis sets for the second-row transition metal atoms are developed for use in geometry optimization calculations. Our method is patterned after previous work on basis sets for first-row transition metal atoms. The basis sets are constructed from the minimal basis sets of Huzinaga and are augmented with a set of diffuse p and d functions. The exponents of these diffuse functions are chosen to minimize both the difference between the calculated and experimental equilibrium geometries and the total molecular energies for several second-row transition metal inorganic and organon etallic complexes. Slightly smaller basis sets, based on the same Huzinaga minimal sets but augmented with a set of diffuse s and p functions rather than diffuse p and d functions, are also presented. The performance of these basis sets is tested on a wide variety of second-row transition metal inorganic and organometallic complexes and is compared to pseudopotential basis sets incorporating effective core potentials.     

Improvement of polarized double-zeta basis sets for molecular interactions. Complexes of NH3, OH2, and FH with H+ and Li+

Interaction energies of a proton and lithium cation with NH₃, OH₂, and HF axe computed at the SCF and Møller-Plesset levels. It is found that the basis set superposition errors at all levels with the 6-3 1G** basis set may be greatly reduced by including an additional diffuse sp shell on first-row atoms; the exponent of this shell has been optimized to 0.1.     

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.