H_2O与NO,CN,OH自由基及负离子相互作用的CP研究

采用量子化学从头算方法在均衡（CounterPolse）校正和非校正势能面上研究 了珑O与NO，NO~-，CN，CNoH~-之间的相互作用比较了CP-梯度优化和非校正梯度优 化以及基组的选择对超分子结构和能量的影响．研究311++g**基组对于这些体系的 研究有很高的效率．这些自由基和负离子均能与H_2O形成强弱不同的氢强弱次序依 次为OH~-，CN，NO~_，OH，CN，NO，由相互作用能．△E~(CP)，成键临界点电荷密 E（2）分析均得到同样结果CF梯度优化和非校正梯度优化所得能量及正BSSE相差不 大。     

Theoretical Study of Hydrogen Bonding Interaction of 1:1 Dimer of HNO with HArF

The hydrogen bonding interaction of 1:1 dimer formed between HNO and HArF molecule has been completely investigated in the present study using Second-order M?ller-Plesset Perturbation (MP2) method in conjunction with 6-311+G**, 6-311++G** and 6-311++G(2d,2p) basis sets. The standard and CP-corrected calculations have been employed to determine the equilibrium structures, the vibrational frequencies and interaction energies. The interaction energies of the dimers were also calculated at G2MP2 level. Two stable structures are found as the minima. Dimer I(H···F)is a five-membered cyclic hydrogen bonded structure and is more stable than the Dimer II(H···O). The blue-shifted N-H···F hydrogen bond is confirmed with standard and CP-corrected calculations by the MP2 and DFT methods in conjunction with different basis sets. The results obtained at MP2 in conjunction with different basis sets show there is a red-shifted hydrogen bond (Ar-H···O) in the Dimer II(H···O). The topological and electronic properties, the origin of red- and blue-shifted hydrogen bonds were investigated at MP2/6-311++G(2d,2p) with CP corrected calculations. From the NBO analysis, the reasonable explanations for the red- and blue-shifted hydrogen bonds were proposed.     

Counterpoise-corrected potential energy surfaces of simple H-bonded systems

Geometries and stabilization energies of various simple H-bonded complexes (water dimer, hydrogen fluoride dimer, formamide dimer, formic acid dimer) have been determined by a gradient optimization that eliminates the basis set superposition error (BSSE) by the counterpoise (CP) method in each gradient cycle as well as by the standard gradient optimization. Both optimization methods lead to different potential energy surfaces (PES). The difference depends on the theoretical level used and is larger if correlation energy is considered. Intermolecular distances from the CP-corrected PES are consistently longer, and this difference might be significant (∼0.1 ?); also angular characteristics determined from both surfaces differ significantly. Different geometries were obtained even when passing to larger basis sets (aug-cc-pVDZ). The standard optimization procedure can result in a completely wrong structure. For example, the “quasi-linear” structure of the (HF)₂ (global minimum) does not exist at the standard MP2/ 6-31G** PES (where only cyclic structure was detected) and is found only at the CP-corrected PES. Stabilization energies obtained from the CP-corrected PES are always larger than these from the standard PES where the BSSE is added only a posteriori for the final optimized structure; both energies converge only when passing to a larger basis set (aug-cc-pVDZ). Received: 11 March 1998 / Accepted: 19 June 1998 / Published online: 4 September 1998 RID=" ID=" <E6>Acknowledgements.</E6> The project was supported by the Grant Agency of the Czech Republic (Grant No. 203/98/1166). RID=" ID=" <E5>Correspondence to</E5>: P. Hobza     

Ab initio study of the vibrational spectra of the hydrogen-bonded nitric acid dimer.

The changes in the vibrational characteristics characterizing the dimerization of nitric acid have been investigated by ab initio calculations at the MP2 level, with 6-31G(d,p) and 6-31 + G(d,p) basis sets, and B3LYP/6-31G(d,p) calculations. The most consistent agreement between the computed values of the frequency shifts for the planar fully symmetric structure (2A) and those experimentally observed suggests that this structure is preferred. It was established that the most sensitive to the complexation is the stretching O-H vibration. The values of the frequency shift (-306 cm(-1)) is indicative for the formation of the relatively strong hydrogen bonds. The calculations predict an increase of the infrared intensity of the stretching O-H vibration in the nitric acid dimer more than 26 times.     

Comparison of some representative density functional theory and wave function theory methods for the studies of amino acids

Energies of different conformers of 22 amino acid molecules and their protonated and deprotonated species were calculated by some density functional theory (DFT; SVWN, B3LYP, B3PW91, MPWB1K, BHandHLYP) and wave function theory (WFT; HF, MP2) methods with the 6-311++G(d,p) basis set to obtain the relative conformer energies, vertical electron detachment energies, deprotonation energies, and proton affinities. Taking the CCSD/6-311++G(d,p) results as the references, the performances of the tested DFT and WFT methods for amino acids with various intramolecular hydrogen bonds were determined. The BHandHLYP method was the best overall performer among the tested DFT methods, and its accuracy was even better than that of the more expensive MP2 method. The computational dependencies of the five DFT methods and the HF and MP2 methods on the basis sets were further examined with the 6-31G(d,p), 6-311++G(d,p), aug-cc-pVDZ, 6-311++G(2df,p), and aug-cc-pVTZ basis sets. The differences between the small and large basis set results have decreased quickly for the hybrid generalized gradient approximation (GGA) methods. The basis set convergence of the MP2 results has been, however, very slow. Considering both the cost and the accuracy, the BHandHLYP functional with the 6-311++G(d,p) basis set is the best choice for the amino acid systems that are rich in hydrogen bonds.     

A comparison of the behavior of functional/basis set combinations for hydrogen-bonding in the water dimer with emphasis on basis set superposition error

We evaluate the performance of ten functionals (B3LYP, M05, M05-2X, M06, M06-2X, B2PLYP, B2PLYPD, X3LYP, B97D, and MPWB1K) in combination with 16 basis sets ranging in complexity from 6-31G(d) to aug-cc-pV5Z for the calculation of the H-bonded water dimer with the goal of defining which combinations of functionals and basis sets provide a combination of economy and accuracy for H-bonded systems. We have compared the results to the best non-density functional theory (non-DFT) molecular orbital (MO) calculations and to experimental results. Several of the smaller basis sets lead to qualitatively incorrect geometries when optimized on a normal potential energy surface (PES). This problem disappears when the optimization is performed on a counterpoise (CP) corrected PES. The calculated interaction energies (ΔEs) with the largest basis sets vary from -4.42 (B97D) to -5.19 (B2PLYPD) kcal/mol for the different functionals. Small basis sets generally predict stronger interactions than the large ones. We found that, because of error compensation, the smaller basis sets gave the best results (in comparison to experimental and high-level non-DFT MO calculations) when combined with a functional that predicts a weak interaction with the largest basis set. As many applications are complex systems and require economical calculations, we suggest the following functional/basis set combinations in order of increasing complexity and cost: (1) D95(d,p) with B3LYP, B97D, M06, or MPWB1k; (2) 6-311G(d,p) with B3LYP; (3) D95++(d,p) with B3LYP, B97D, or MPWB1K; (4) 6-311++G(d,p) with B3LYP or B97D; and (5) aug-cc-pVDZ with M05-2X, M06-2X, or X3LYP.     

Gas-phase acidity of D-glucose. A density functional theory study

The gas-phase acidity of D-glucopyranose was studied by means of B3LYP calculations combined with 6-31G(d,p) or 6-31+G(d,p) standard basis sets. For each anomer, deprotonation of the various primary and secondary hydroxyl groups was considered. As in solution, the anomeric hydroxyl is found to be the most acidic for both anomers, but only when the 6-31+G(d,p) basis set is used for geometry optimization. Deprotonation of the anomeric hydroxyl induces an important C(1)--O endocyclic bond elongation and subsequently promotes an energetically favored ring-opening process as attested by the very small calculated activation barriers. The results also suggest that interconversion between the various deprotonated alpha- and beta-anomers may easily occur under slightly energetic conditions. B3LYP/6-311+G(2df,2p) calculations led to the an absolute gas-phase acidity of deltaacidGo(298)(alpha-D-glucose) = 1398 kJ mol(-1). This estimate matches well the only experimental value available to date. Finally, this study again confirms that the use of diffuse functions on heavy atoms is necessary to describe anionic systems properly and to achieve good relative and absolute gas-phase acidities.     

Hydrogen peroxide clusters: the role of open book motif in cage and helical structures

Hartree-Fock (HF) calculations using 6-31G*, 6-311++G(d,p), aug-cc-pVDZ, and aug-cc-pVTZ basis sets show that hydrogen peroxide molecular clusters tend to form hydrogen-bonded cyclic and cage structures along the lines expected of a molecule which can act as a proton donor as well as an acceptor. These results are reiterated by density functional theoretic (DFT) calculations with B3LYP parametrization and also by second-order M?ller-Plesset perturbation (MP2) theory using 6-31G* and 6-311++G(d,p) basis sets. Trends in stabilization energies and geometrical parameters obtained at the HF level using 6-311++G(d,p), aug-cc-pVDZ, and aug-cc-pVTZ basis sets are similar to those obtained from HF/6-31G* calculation. In addition, the HF calculations suggest the formation of stable helical structures for larger clusters, provided the neighbors form an open book structure.     

A density functional theory study of the structure and vibrational spectra of beta-carotene, capsanthin, and capsorubin

A theoretical study of the structure and the vibrational spectra of the beta-carotene molecule and its derivatives capsanthin and capsorubin is carried out. We first investigate systematically the theoretical method which provides the best results for beta-carotene by performing ab initio calculations at the HF/6-31G(d), SVWN/6-31G(d), PBE0/6-31G(d), BLYP/6-31G(d), B3LYP/6-31G(d), B3LYP/6-31G(d,p), B3LYP/6-311G(d), and B3LYP/6-311G(d,p) levels and by using previous theoretical results available in the literature obtained at the AM1 and BPW91/6-31G(d) levels. The influence of both the level of calculation and the size of the basis set used in the geometry optimization and in the determination of the IR and Raman spectra of this molecule is thus analyzed. It is confirmed that the hybrid functional B3LYP with the basis 6-31G(d) is the method that gives the best results as a whole. By use of this level of calculation, we next optimize the molecular geometries of related molecules of capsanthin and capsorubin, which to the best of our knowledge have only been studied at the semiempirical AM1 level. In addition we calculate the IR and Raman spectra of these molecules at the B3LYP/6-31G(d) level of theory. The results obtained for capsanthin show on the one hand that the double bond of the beta-ionone ring is outside the polyene chain plane, due to the repulsion between the hydrogen atoms of the ring methyl groups and the hydrogen atoms of the polyene chain, and on the other hand that the carbonyl double bond in the other headgroup is very close to planarity with the polyene chain, since in this case such a repulsion does not exist. For the molecule of capsorubin the two carbonyl groups also take the same coplanar orientation relative to the polyene chain. The IR and Raman spectra theoretically computed for these two molecules are finally compared with their experimental spectra and the vibrational normal modes of the main signals are interpreted.     

Theoretical study of the structures and vibrational spectra of the hydrogen-bonded systems of 4-cyanophenol with N-bases

The structural and vibrational features of the hydrogen bonded complexes of 1,5,7-triazabicyclo [4.4.0] dec-5-ene (TBD) with one and two 4-CNPhOH molecules have been studied extensively by ab initio SCF/6-31G(d,p) and BLYP calculations with various basis sets: 6-31G(d,p), 6-31+G(d,p) and 6-31++G(d,p). Full geometry optimization was made for the complexes studied. The nature of the hydrogen bonding and the influence of the hydrogen bonding on the structural and vibrational characteristics of the monomers have been investigated. The corrected values of the dissociation energy for the hydrogen-bonded complexes have been calculated in order to estimate their stability. The calculated values of the dissociation energy per phenol molecule indicate that the complex: TBD: 4-CNPhOH (1:1) is more stable than the complex: TBD: 4-CNPhOH (1:2). The changes in the structural and vibrational characteristics upon hydrogen bonding depend on the strength of the hydrogen bonds. In agreement with the experiment, the calculations show that the complexation between TBD and 4-CNPhOH leads to considerably changes in the vibrational characteristics of the stretching O-H vibration. The vibrational frequency of the O-H stretching vibration is shifted to lower wave numbers upon hydrogen bonding. The predicted frequency shifts Deltanu(O-H) for the complexes--TBD: 4-CNPhOH (1:1) and TBD: 4-CNPhOH (1:2) are in the range from -190 cm(-1) to -586 cm(-1). In the same time the IR intensity of the O-H stretching vibration increases dramatically in the hydrogen-bonded complexes.     

Methyl hydrogen peroxide dimer: A structural study

Methyl hydrogen peroxide (MHP) exhibits a tendency to form a stable dimer by hydrogen-bonding. Ab initio theoretical investigations on methyl hydrogen peroxide dimer (MHPD) carried out herein lead to several energetically stable structures that have a direct bearing on the reactivity of the monomer in terms of its molecular electrostatic potential (MESP). To gauge the role played by the electron-correlation in lending stability to MHP and its dimer, we employ the density functional theory (DFT) (as implemented by B3LYP-functional), and subsequently second order Møller-Plesset (MP2) perturbation theory, using the basis sets 6-31G(d, p) and 6-311++G(2d, 2p). Simulated infra-red vibrational spectra lead to spectral intensity redistribution upon dimerization. Energetically the lowest MHPD is endowed with inversion symmetry and has two hydrogen bonds, while three other structures emerge: one energetically very close with two H-bonds, and the two others, with three H-bonds each, yet higher by about 2 kcal mol⁻¹.     

分子几何构型优化方法的系统性比较

对《ＣＲＣ物理化学手册》第７７版中收集的第三周期以前的所有已知实验构型的无机分子的构型,以ＭＰ２、Ｂ３ＬＹＰ、Ｂ３ＰＷ９１级别上进行了构型优化的系统性比较,优化采用基组为６－３１Ｇ（ｄ,ｐ）,６－３１１Ｇ（ｄ,ｐ）,６－３１１Ｇ（２ｄ,ｐ）。对大多数分子另比较了ＱＣＩＳＤ（Ｔ）方法,对多原子分子比较了ＢＰＷ９１方法,对含氢双原子分子也使用ＱＣＩＳＤ方法。结果表明,键长的平均绝对偏差（单位：ｐｐｍ）     

Estimation on the intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies in glycine and alanine peptides

Computation of accurate intramolecular hydrogen-bonding energies for peptides is of great importance in understanding the conformational stabilities of peptides and developing a more accurate force field for proteins. We have proposed a method to determine the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies in glycine and alanine peptides. In this article, the method is further applied to evaluate the intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies in peptides. The optimal structures of the intramolecular 10-membered ring N-H...O=C hydrogen bonds in glycine and alanine tripetide molecules are obtained at the MP2 level with 6-31G(d), 6-31G(d,p), and 6-31+G(d,p) basis sets. The intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies are then evaluated based on our method at the MP2/6-311++G(3df,2p) level with basis set superposition error correction. The intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies are calculated to be in the range of -6.84 to -7.66, -4.44 to -4.98, and -6.95 to -7.88 kcal/mol. The method is also applied to estimate the individual intermolecular hydrogen-bonding energies in the dimers of amino-acetaldehyde, 2-amino-acetamide, formamide, and oxalamide, each dimer having two identical intermolecular hydrogen bonds. According to our method, the individual intermolecular hydrogen-bonding energies in the four dimers are calculated to be -1.77, -1.67, -6.35, and -4.82 kcal/mol at the MP2/6-311++G(d,p) level, which are in good agreement with the values of -1.84, -1.72, -6.23, and -4.93 kcal/mol predicted by the supermolecular method.     

NMR spectra, GIAO and charge density calculations of five-membered aromatic heterocycles

The B3LYP/6-31+G(d) molecular geometry optimized structures of 17 five-membered heterocycles were employed together with the gauge including atomic orbitals (GIAO) density functional theory (DFT) method at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p) and B3LYP/6-311+G(2d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants. The method of geometry optimization for pyrrole (1), N-methylpyrrole (2) and thiophene (7) using the larger 6-311++G(d,p) basis sets at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,p) and B3LYP/cc-pVTZ levels of theory gave little difference between calculated and experimental values of coupling constants. In general, the (1)H and 13C chemical shifts for all compounds are in good agreement with theoretical calculations using the smaller 6-31 basis set. The values of nJHH(n=3, 4, 5) and rmnJ(CH)(n=1, 2, 3, 4) were predicted well using the larger 6-31+G(d,p) and 6-311++G(d,p) basis sets and at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,2p) levels of theory. The computed atomic charges [Mülliken; Natural Bond Orbital Analysis (NBO); Merz-Kollman (MK); CHELP and CHELPG] for the B3LYP/6-311++G(d,p) geometry optimized structures of 1-17 were used to explore correlations with the experimental proton and carbon chemical shifts.     

Theoretical study of the structures, stability and vibrational spectra of the nitrous acid complexes with CH4

The structures, stability and vibrational spectra of the binary complexes CH4...HONO-trans and CH4...HONO-cis have been investigated using ab initio calculations at the SCF and MP2 levels with 6-311++G(d,p) basis set and B3LYP calculations with 6-31G(d,p) and 6-31+G(d,p) basis sets. Full geometry optimization was made for the complexes studied. It was established that the complex CH4...HONO-trans is more stable by 0.41 kcal mol(-1) than the complex CH4...HONO-cis. The accuracy of the ab initio calculations have been estimated by comparison between the predicted values of the vibrational characteristics (vibrational frequencies and infrared intensities) and the available experimental data. It was established, that the methods, used in this study are well adapted to the problem under examination. The predicted values with the B3LYP calculations are very near to the results, obtained with 6-311++G(d,p)/MP2. The changes in the vibrational characteristics of methane and trans-, cis-nitrous acid upon formation of the hydrogen bond show that the complexes CH4...HONO-trans and CH4...HONO-cis have geometry in which the OH group interacts with a methane molecule forming a single hydrogen bond. This fact is confirmed by relatively strong perturbation of the OH stretching vibration to lower frequencies and an increase of the infrared intensity of this vibration up to three times upon hydrogen bonding.     

态定(线性响应)-极化连续模型/含时密度泛函方法研究一种有机发光材料的吸收和发射光谱

3-(二氰亚甲基)-5,5-二甲基-1-(3-[9-(2-乙基-己基)-咔唑基]-乙烯基)环己烷(DCDHCC)是一种用于光电器件中的有机染料,它具有良好的发光特性.我们使用含时密度泛函方法(TD-PBE0,TD-BMK和TD-M06)以及极化连续模型(PCM)计算了该材料在溶剂中的吸收和发射特性.计算中使用了线性响应(LR)、态定(SS)两种溶剂模型和6-31G(d)、6-31+G(d,p)两种基组.计算了DCDHCC在苯、四氢呋喃和丙酮溶剂中的吸收和发射光谱,并与实验观测进行了比较.结果表明:对于吸收光谱的计算,杂化函数的影响大于基组和溶剂模型,在三种函数中BMK更适于研究DCDHCC的吸收光谱;而对于发射光谱,基组的影响最大,基组通过影响激发态构型从而影响发射光谱,对于激发态构型的优化需要使用6-31+G(d,p)基组.我们希望这些研究能对今后设计类似的发光分子有帮助.     

Ab initio quantum chemical investigation of the ground and excited states of salicylic acid dimer

The structure and stability of different forms of salicylic acid dimer have been examined by Hartree-Fock and density functional theoretic calculations using 6-31G(d,p) and 6-311++g(d,p) basis sets. Vertical excitation energies for the monomer as well as the dimer have been computed using the time-dependent density functional theory using 6-311++G(d,p) basis set. The predicted absorption maxima for the first excited singlet state of salicylic acid monomer and the dimer of the primary form are in reasonable agreement with the experimental result. There is a slight red shift (approximately 6 nm) in the absorption maximum in going from the monomer to the dimer, in accord with the experimental observation. Configuration-interaction calculations including single excitation have been carried out to map the potential-energy profile for the intra- as well as the intermolecular proton transfer in different forms of the dimer. The barrier for proton transfer in the ground state as well as the excited states makes it clear that most of the processes take place in the primary form and largely by intramolecular proton transfer.     

Coordination of methanol clusters to benzene: a computational study

Benzene-methanol cluster structures were investigated with theoretical chemistry methods to describe the microsolvation of benzene and the benzene-methanol azeotrope. Benzene-methanol (MeOH) clusters containing up to six methanol molecules have been calculated by ab initio [MP2/6-311++G(d,p)//MP2/6-31+G(d,p) + BSSE correction] method. The BSSE was found quite large with this basis set, hence, different extrapolation schemes in combination with the aug-cc-pVxZ basis sets have been used to estimate the complete basis set limit of the MP2 interaction energy [ΔE(MP2/CBS)]. For smaller clusters, n ≤ 3, DFT procedures (DFTB+, MPWB1K, M06-2X) have also been applied. Geometries obtained for these clusters by M06-2X and MP2 calculations are quite similar. Based on the MP2/CBS results, the most stable C(6)H(6)(MeOH)(3) cluster is characterized by a hydrogen bonded MeOH trimer chain interacting with benzene via π···H-O and O···H-C(benzene) hydrogen bonds. Larger benzene-MeOH clusters with n ≥ 4 consist of cyclic (MeOH)(n) subclusters interacting with benzene by dispersive forces, to be denoted by C(6)H(6) + (MeOH)(n). Interaction energies and cooperativity effects are discussed in comparison with methanol clusters. Besides MP2/CBS calculations, for selected larger clusters the M06-2X/6-311++G(d,p)//M06-2X/6-31+G(d,p) procedure including the BSSE correction was also used. Interaction energies obtained thereby are usually close to the MP2/CBS limit. To model the benzene-MeOH azeotrope, several structures for (C(6)H(6))(2)(MeOH)(3) clusters have been calculated. The most stable structures contain a tilted T-shaped benzene dimer interacting by π···H-O and O···H-C (benzene) hydrogen bonds with a (MeOH)(3) chain. A slightly less negative interaction energy results for a parallel displaced benzene sandwich dimer with a (MeOH)(3) chain atop of one of the benzene molecules.     

氮原子与卤代甲烷反应的直接氢抽提过程研究

利用B3LYP理论研究了N(~4S)+CH_3X(X = H, F, Cl)反应体系的直接氢抽提过 程,分别得到了各反应物、产物和过渡态的优化构型和谐振频率。同时应用了6- 31G(d), 6-311+G(d,p)和6-311+ + G(2d,2p)基组,考察其大小对反应体系中各物 种构型及能量的影响。理论计算表明,随着基组的增加,反应势垒逐渐降低,反应 吸热减少。对比取代甲烷的情形,结果表明反应过程中卤素原子具有典型的诱导效 应,降低了抽提势垒。