(CH2)2O, (CH2)2S与双卤分子间卤键的理论研究

运用量子化学密度泛函B3LYP方法, 采用6-311＋＋G(d,p)及aug-cc-pVDZ基组, 通过CP校正的几何梯度优化对(CH2)2O和(CH2)2S与双卤分子XY (XY＝Cl2, Br2, ClF, BrF, BrCl)形成的卤键复合物的几何构型、振动频率和相互作用能等进行了研究. 利用电子密度拓扑分析理论方法对卤键复合物的拓扑性质进行了分析研究, 探讨了该类分子间卤键的作用本质. 结果表明, (CH2)2O和(CH2)2S与双卤分子间的卤键介于共价键与离子键之间, 偏于静电作用成分为主. 形成卤键后, 双卤分子的键长增加, 振动频率减小, 原子积分性质发生改变. 卤键键长的变化、键能的强弱、键鞍点处的电子密度值与双卤分子的电负性有关.     

硫代双烯酮与异氰酸环加成反应的理论研究

采用从头算ＲＨＦ／６－３１Ｇ方法研究了硫代双烯酮与异氰酸之间两种可能的环加成反应的机理,并对反应各驻点进行了电子密度拓扑分析研究．结果表明,这两个生成不同四元杂环产物的平行反应均为非同步的协同反应,但两个反应进行的难易程度不同,形成硫氮杂环的反应更容易一些,而形成氮杂环反应的产物在热力学上更稳定一些     

呋喃与双卤分子间卤键的电子密度拓扑研究

用DFT和MP2(full)方法研究了呋喃-XY (XY＝ClF, BrF, BrCl)体系分子间的相互作用, 讨论了两种卤键类型 Y—X…O和Y—X…π的作用. 研究表明: Y—X…π的相互作用比Y—X…O的作用强, 计算得到的分子间相互作用能按照呋喃-BrF＞呋喃-ClF＞呋喃-BrCl的顺序依次降低. 形成卤键后, 两种类型卤键复合物中的电子受体X—Y键伸长, 其振动频率发生红移. 利用电子密度拓扑方法对卤键的电子密度拓扑性质进行了分析.     

蒽与苯甲酸及其衍生物之间的作用机理研究

用紫外差光谱、红外光谱法和荧光猝灭法研究了典型多环芳烃蒽与苯甲酸及其羟基取代衍生物邻羟基苯甲酸和对羟基苯甲酸之间的作用机理。实验结果表明,在此芳香羧酸与蒽之间存在定向的、特异性的作用,其作用方式受到反应物结构和环境酸度的影响。对苯甲酸和邻羟基苯甲酸而言,当pH<pKa时,二者之间以蒽离域大π电子与羧基质子之间的π-H氢键作用为主,pH>pKa时,π-π电子给体受体作用逐渐成为主要结合方式。邻位羟基的存在使得苯甲酸与蒽的作用强度明显降低。对羟基苯甲酸特殊的D-π-A型分子结构使得它在溶液中形成平面多分子聚集体,这种多分子聚集体的生成使得对羟基苯甲酸与蒽的结合方式不随酸度变化,在pH2.0~10.0的范围内均以π-π电子给体受体作用相结合,且结合强度大于苯甲酸和邻羟基苯甲酸。     

Ab initio and AIM studies on typical π‐type and pseudo‐π‐type halogen bonds: Comparison with hydrogen bonds

Series of typical π‐type and pseudo‐π‐type halogen‐bonded complexes B ··· ClY and B ··· BrY and hydrogen‐bonded complex B ··· HY (B = C₂H₄, C₂H₂, and C₃H₆; Y = F, Cl, and Br) have been investigated using the MP2/aug‐cc‐pVDZ method. A striking parallelism was found in the geometries, vibrational frequencies, binding energies, and topological properties between B ··· XY and B ··· HY (X = Cl and Br). It has been found that the lengths of the weak bond d(X ··· π)/d(H ··· π), the frequencies of the weak bond ν(X ··· π)/ν(H ··· π), the frequency shifts Δν(X? Y)/Δν(H? Y), the electron densities at the bond critical point of the weak bonds ρ_c(X ··· π)/ρ_c(H ··· π), and the electron density changes Δρ_c(X? Y)/Δρ_c(H? Y) could be used as measures of the strengths of typical π‐type and pseudo‐π‐type halogen/hydrogen bonds. The typical π‐type and pseudo‐π‐type halogen bond and hydrogen bond are noncovalent interactions. For the same Y, the halogen bond strengths are in the order B ··· ClY < B ··· BrY. For the same X, the halogen bond strength decreases according to the sequence F > Cl > Br that is in agreement with the hydrogen bond strengths B ··· HF > B ··· HCl > B ··· HBr. All of these typical π‐type and pseudo‐π‐type hydrogen‐bonded and halogen‐bonded complexes have the “conflict‐type” structure. Contour maps of the Laplacian of π electron density indicate that the formation of B ··· XY halogen‐bonded complex and B ··· HY hydrogen‐bonded complex is very similar. Charge transfer is observed from B to XY/HY and both the dipolar polarization and the volume of the halogen atom or hydrogen atom decrease on B ··· XY/B ··· HY complex formation. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

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

分别在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型锂键.     

乙炔基自由基C2H与氧气反应的密度泛函理论研究

应用量子化学从头算和密度泛函理论(DFT)对C₂H自由基和O₂的反应进行了研究.在B3LYP/6-311G^**水平上优化了反应通道上各驻点(反应物、中间体、过渡态和产物)的几何构型,并计算出它们的振动频率和零点振动能(ZPVE).各物种的总能量由CCSD(T)/6-311G^**//B3LYP/6-311G^**给出,并对能量进行了零点能校正.计算结果表明,反应物中自由基C₂H中的边端C进攻O₂形成了中间体1 (HCCOO),中间体1是一个加合产物.由中间体1经过不同的反应通道可以生成不同的产物P₁ (HCO+CO), P₂ (HCCO+O), P₃(CO₂+CH), P₄ (C₂O+OH)和P₅ (2CO+H).反应通道之间存在着竞争机制.其中P₁, P₂是主要产物,其次还有一定比例的P5生成,而产物P₃, P₄的生成几率较低.各条反应通道化学反应热的计算与实验吻合较好.     

H2O, H2S与双卤分子间相互作用的电子密度拓扑研究

运用量子化学微扰理论MP2和密度泛函B3LYP方法, 采用6-311＋＋G(d,p)基组, 对H₂O, H₂S与双卤分子XY (XY＝F₂, Cl₂, Br₂, ClF, BrF, BrCl)形成的卤键复合物进行构型全优化, 并计算得到了这些体系的分子间相互作用能. 利用电子密度拓扑分析方法对卤键复合物的拓扑性质进行了分析研究, 探讨了该类分子间卤键的作用本质. 结果表明, 形成卤键后, 作为电子受体的双卤分子X—Y键长增长, 振动频率减小. 复合物体系中的卤键介于共价键与离子键之间, 偏于静电作用成分为主.     

Theoretical study on the interactions of halogen‐bonds and pnicogen‐bonds in phosphine derivatives with Br2, BrCl,and BrF

The MP2 ab initio quantum chemistry methods were utilized to study the halogen‐bond and pnicogen‐bond system formed between PH₂X (X = Br, CH₃, OH, CN, NO₂, CF₃) and BrY (Y = Br, Cl, F). Calculated results show that all substituent can form halogen‐bond complexes while part substituent can form pnicogen‐bond complexes. Traditional, chlorine‐shared and ion‐pair halogen‐bonds complexes have been found with the different substituent X and Y. The halogen‐bonds are stronger than the related pnicogen‐bonds. For halogen‐bonds, strongly electronegative substituents which are connected to the Lewis acid can strengthen the bonds and significantly influenced the structures and properties of the compounds. In contrast, the substituents which connected to the Lewis bases can produce opposite effects. The interaction energies of halogen‐bonds are 2.56 to 32.06 kcal·mol^?1; The strongest halogen‐bond was found in the complex of PH₂OH???BrF. The interaction energies of pnicogen‐bonds are in the range 1.20 to 2.28 kcal·mol^?1; the strongest pnicogen‐bond was found in PH₂Br???Br₂ complex. The charge transfer of lp(P) ? σ*(Br? Y), lp(F) ? σ*(Br? P), and lp(Br) ? σ*(X? P) play important roles in the formation of the halogen‐bonds and pnicogen‐bonds, which lead to polarization of the monomers. The polarization caused by the halogen‐bond is more obvious than that by the pnicogen‐bond, resulting in that some halogen‐bonds having little covalent character. The symmetry adapted perturbation theory (SAPT) energy decomposition analysis showes that the halogen‐bond and pnicogen‐bond interactions are predominantly electrostatic and dispersion, respectively.     

Theoretical study on the mechanisms of the nucleophilic substitution reactions of hydrosulfide ion and halomethanes

Three species involved in the nucleophilic substitution reaction of hydrosulfide ion and halomethanes are investigated by ab initio calculations. Geometries for stationary structures along the reaction paths are fully performed with the second‐order Møller–Plesset perturbation approximation with the cc‐pVDZ basis set. The monomer geometries determined by the MP2 method match the experimental results very well. Single point energy calculations are carried out at the coupled cluster with perturbative triple excitations CCSD (T) theory with aug‐cc‐pVDZ basis set. Halomethanes have three conformers here, which lead to the three product channels, HSCH₃ + F^?1, HSCH₃ + Cl^?1, and HSCH₃ + Br^?1. The investigation encompasses the six complexes formed among three channels, respectively. By selecting the six complexes as the model, we investigate the binding energy, topological property of the electron charge density and their Laplacian in detail theoretically. Electrostatic density potential maps of halomethanes are generated for the determination of attractive interaction sites. It is proved that the similar misshaped electron clouds of the three halogen atoms result in the similar properties of the carbon‐halogen bonds, and reveals that the product ion‐dipole complexes interactions are predominantly electrostatic in nature. The calculated results predict the binding energy of the most stable complex in six complexes is ?47.06 kcal/mol at the MP2 level of theory. The second channel has the lowest energy barrier, which is ?3.63 kcal/mol at the CCSD (T) levels of theory, is expected to be the most important pathway. It occurs via C? Cl cleavage accompanied by C? S bond formation. The other two channels have higher energy barriers and are expected to have smaller rates. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Study on the Cation-π Interactions between Ammonium Ion and Aromatic π Systems

The nature and strength of the cation-π interactions between NH4^＋ and toluene, p-cresol, or Me-indole were studied in terms of the topological properties of molecular charge density and binding energy decomposition. The results display that the diversity in the distribution pattern of bond and cage critical points reflects the profound influence of the number and nature of substituent on the electron density of the aromatic rings. On the other hand, the energy decomposition shows that dispersion and repulsive exchange forces play an important role in the organic cation （NH4^＋）-π interaction, although the electrostatic and induction forces dominate the interaction. In addition, it is intriguing that there is an excellent correlation between the electrostatic energy and ellipticity at the bond critical point of the aromatic π systems, which would be helpful to further understand the electrostatic interaction in the cation-π complexes.     

Interplay between anion-pi and hydrogen bonding interactions

The interplay between two important noncovalent interactions involving aromatic rings is studied by means of high level ab initio calculations. They demonstrate that synergistic effects are present in complexes where anion-pi and hydrogen bonding interactions coexist. These synergistic effects have been studied using the "atoms-in-molecules" theory and the Molecular Interaction Potential with polarization partition scheme. The present study examines how these two interactions mutually influence each other.     

杂氮硅三环类化合物的紫外光电子能谱及化学键的理论研究

报道了10种杂氮硅三环类化合物的紫外光电子能谱(UPS).采用RHF/3-21G^*对各化合物进行了构型优化.根据各化合物UPS的谱带特征结合RHF/6-31G^*的计算结果对谱图进行了合理的解析和指认.利用6-31G^*基组进行的电子密度拓扑分析表明,各化合物的N和Si原子间均存在键鞍点,从理论上证实了N和Si原子间存在成键作用.各化合物的UPS和SCFMO计算和电子密度拓扑分析都表明,在该类体系中σ_N→Si配键是比杂氮硼三环类化合物中σ_N→B配键弱得多的成键作用;文中同时讨论了不同取代基对Si和N间成键作用的影响.     

Long‐range π‐type hydrogen bond in dimers CH2?HF,CH2O?H2O,and CH2O?NH3

Using four basis bets, (6‐311G(d,p), 6‐31+G(d,p), 6‐31++G(2d,2p), and 6‐311++G(3df,3pd), the optimized structures with all real frequencies were obtained at the MP2 level for the dimers CH₂O? HF, CH₂O? H₂O, CH₂O? NH₃, and CH₂O? CH₄. The structures of CH₂O? HF, CH₂O? H₂O, and CH₂O? NH₃ are cycle‐shaped, which result from the larger bend of σ‐type hydrogen bonds. The bend of σ‐type H‐bond O…H? Y (Y?F, O, N) is illustrated and interpreted by an attractive interaction of a chemically intuitive π‐type hydrogen bond. The π‐type hydrogen bond is the interaction between one of the H atoms of CH₂O and lone pair(s) on the F atom in HF, the O atom in H₂O, or the N atom in NH₃. In contrast with the above three dimers, for CH₂O? CH₄, because there is not a π‐type hydrogen bond to bend its linear hydrogen bond, the structure of CH₂O? CH₄ is noncyclic shaped. The interaction energy of hydrogen bonds and the π‐type H‐bond are calculated and discussed at the CCSD (T)/6‐311++G(3df,3pd) level. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

硫代双烯酮与硫甲醛环加成反应机理的理论研究

采用从头算HF/6-31G^*和密度泛函B3LYP/6-31++G^*^*方法研究了硫代双烯酮与硫甲醛两种可能的环加成反应的机理,并对反应各驻点进行了电子密度拓扑分析研究.结果表明,这两个生成不同四元杂环产物的平行反应均为有一两性离子中间体的分步反应.两个反应均较易进行,但反应(1)更容易一些,结果与实验一致.     

Theoretical Study of the Isomerization Reaction of 1-3 H Transfer on Formamidine Substituted by Halogen

The ab initio method has been used to study the 1-3 H transfer reaction on formamidine substituted by halogen. The calculation results show that the substituted halogen has two effects on the 1-3 H transfer reaction: decreasing the activation energy and stabilizing the C=N double bond owing to the conjugative effect of p-π-p of products and transition states.     

Competition between Halogen Bonding and π‐Hole Interactions Involving Various Donors: The Role of Dispersion Effects

In this study several σ‐ and π‐hole complexes between IF and pnicogen ZO₂F (Z=P, As), chalcogen ChO₃ (Ch=S, Se) and tetrel TrOF₂ (Tr=Si, Ge) ‐bearing compounds were optimized at the RI‐MP2/def2‐TZVPD level of theory. All complexes were characterized as minima by frequency analysis calculations. In addition, a comparative CCSD(T) and DFT (with and without dispersion correction) study using the BP86, B3LYP and M06‐2X method was done in order to analyze the role of dispersion effects in the σ‐/π‐hole binding. Finally the Bader’s AIM analysis of several complexes was performed to further characterize the interactions discussed herein.     

Effects of π bond type,backbone size,and halogen on structural and spectroscopic properties of hydrogen‐bonded complexes of the X?H … π type between alkenes or alkynes and haloacids (HF and HCl)

The effects of CC bond type (double or triple), substituent (H or methyl), and halogen (F and Cl) on three properties of hydrogen‐bonded complexes formed between unsaturated hydrocarbons and HX (X?F, Cl) are studied. The properties comprise hydrogen bond distances (R_H), stabilization energies (SE), and frequency shifts (Δν). A 2³ factorial design technique, along with ab initio (HF and MP2) and DFT (B3LYP and PBE1PBE) calculations, has been employed. All three responses are mainly affected by the halogen, and when it is changed from F to Cl, R_H tends to increase, while SE tends to decrease. Surprisingly, the type of substituent is more important than the type of CC bond, for all three responses. Both effects tend to decrease R_H. Significant interaction effects are obtained for the type of CC bond along with the type of substituent, and for the type of substituent along with the type of halogen. Both interaction effects are smaller than the main effects and also tend to decrease R_H. The greatest SE values are obtained with PBE1 functional (BSSE + ZPE corrected values). Again, the next more important effect is due to the type of substituent, and the replacement of H by CH₃ group tends to increase SE. The effect due to the CC bond type is not significant, at all computational levels. The only interaction effect that is significant for SE (corrected) and Δν is between factors 1 (CC bond type) and 2 (substituent), but only at HF and B3LYP levels, and it tends to increase both properties. As the halogen changes from F to Cl, Δν tends to decrease. In contrast, changing the substituent from H to CH₃ leads to greater values of Δν. The effect of CC bond type is not significant at HF level, and when it is changed from double to triple Δν is decreased, at B3LYP and PBE1 levels. A suggestion as to how the results may point toward a better experimental detection of similar (π‐type) complexes is also given. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006