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

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

(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键伸长, 其振动频率发生红移. 利用电子密度拓扑方法对卤键的电子密度拓扑性质进行了分析.     

A computational study on the nature of the halogen bond between sulfides and dihalogen molecules

The characteristics and nature of the halogen bonding in a series of B···XY (B = H₂S, H₂CS, (CH₂)₂S; XY = ClF, Cl₂, BrF, BrCl, Br₂) complexes were analyzed by means of the quantum theory of “atoms in molecules” (QTAIM) and “natural bond orbital” (NBO) methodology at the second-order Møller-Plesset (MP2) level. Electrostatic potential, bond length, interaction energy, topological properties of the electron density, the dipole moment, and the charge transfer were investigated systematically. For the same electron donor, the interaction energies follows the B···BrF > B···ClF > B···BrCl > B···Br₂ > B···Cl₂ > B···ClBr order. For the same electron acceptor, the interaction energies increase in the sequence of H₂S, H₂CS, and (CH₂)₂S. Topological analyses show these halogen bonding interactions belong to weak interactions with an electrostatic nature. It was found that the strength of the halogen-bonding interaction correlates well with the electrostatic potential associated with halogen atom and the amount of charge transfer from sulfides to dihalogen molecules, indicating that electrostatic interaction plays an important role in these halogen bonds. Charge transfer is also an important factor in the halogen bonds involved with dihalogen molecules.     

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     

S···X halogen bonds and H···X hydrogen bonds in H2CS-XY (XY = FF, ClF, ClCl, BrF, BrCl, and BrBr) complexes: cooperativity and solvent effect

Using ab initio calculations, we have studied the structures, properties, and nature of halogen bonds in H(2)CS-XY (XY = FF, ClF, ClCl, BrF, BrCl, and BrBr) complexes. The results show that the ring-shaped complexes are formed by a halogen bond (S···X) and a secondary hydrogen bond (H···X). We also analyzed the H(2)CS-ClF-ClF and FCl-H(2)CS-ClF complexes to investigate the cooperative and diminutive halogen bonding. The cooperative effect of halogen bonding is found in the former, while the diminutive effect is present in the latter. We finally considered the solvent effect on the halogen bond in H(2)CS-BrCl complex and found that the solvent has a prominent enhancing effect on it. The complexes have also been analyzed with natural bond orbital, atoms in molecules, and symmetry adapted perturbation theory method.     

Regular/abnormal variation in the strength and nature of the halogen bond between H2Te and the dihalogens: Prominent effect of methyl substituents

The halogen-bonded complexes between H₂Te/Me₂Te and the dihalogen molecules XY (XY = F₂, Cl₂, Br₂, I₂, ClF, ClBr, BrF, BrCl, BrI, IF, ICl, IBr) have been studied to investigate the dependence of its strength and nature on the halogen donor X and its adjoining atom Y, as well as the methyl groups in the electron donor. The interaction energy varies between −1.7 and − 43.5 kcal/mol, indicating that the Te atom in H₂Te/Me₂Te has a strong affinity for the dihalogen molecules. For the H₂Te-XY complex, the halogen bond is stronger for the heavier halogen donor X atom and the strong electron-withdrawing group Y. However, for Me₂Te-XY, the halogen bond is stronger for the lighter halogen donor X atom. The H₂Te/Me₂Te-F₂ complex has the largest interaction energy, although the σ-hole on F₂ is the smallest in magnitude. In most of the complexes, the electrostatic and polarization contributions to the binding strength are similar in magnitude. However, for H₂Te/Me₂Te-F₂, the polarization contribution is much larger than the electrostatic contribution, with a significant contribution from charge transfer.     

Comparative study of XO···ClF and XS···ClF (X = H,CH3, and F) halogen‐bonded complexes

Quantum chemical calculations have been performed to study the single‐electron halogen bonds in HO···ClF and HS···ClF complexes. The calculation methods have a larger effect on the S···Cl halogen bond than on the O···Cl one. The interaction strength in HO···ClF complex is stronger than that in HS···ClF one, but the presence of methyl group in the halogen acceptor makes the sequence reverse. The methyl group has a greater effect on the S···Cl halogen bond than on the O···Cl one. The charge analyses indicate that the methyl group is electron‐donating and the electron‐donating role in the H₃CS? ClF complex is larger than in the H₃CO? ClF one. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Structure,properties, and nature of the BrF‐HX complexes: An Ab initio study

Structure and properties of complexes (energies and charge transfer) of complexes BrF‐HX (X = F, Cl, Br, I) have been investigated at the MP2/aug‐cc‐pVDZ (aug‐cc‐pVDZ‐pp basis sets for I) level. Two types of geometries (hydrogen‐bonded and halogen‐bonded) are observed. The calculated interaction energies show that the halogen bonded structures are more stable than the corresponding hydrogen‐bonded structures. To study the nature of the intermolecular interactions, symmetry‐adapted perturbation theory (SAPT) energy decomposition analysis reveals that the BrF‐HX complexes are dominantly electrostatic in nature. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

A new unconventional halogen bond C?X···H?M between HCCX (X = Cl and Br) and HMH (M = Be and Mg): An ab initio study

In this article, a new type of halogen‐bonded complex YCCX···HMY (X = Cl, Br; M = Be, Mg; Y = H, F, CH₃) has been predicted and characterized at the MP2/aug‐cc‐pVTZ level. We named it as halogen‐hydride halogen bonding. In each YCCX···HMY complex, a halogen bond is formed between the positively charged X atom and the negatively charged H atom. This new kind of halogen bond has similar characteristics to the conventional halogen bond, such as the elongation of the C? X bond and the red shift of the C? X stretch frequency upon complexation. The interaction strength of this type of halogen bond is in a range of 3.34–10.52 kJ/mol, which is smaller than that of dihydrogen bond and conventional halogen bond. The nature of the electrostatic interaction in this type of halogen bond has also been unveiled by means of the natural bond orbital, atoms in molecules, and energy decomposition analyses. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

乙烯、乙炔与双卤分子间π型卤键的电子密度拓扑研究

运用DFT和MP2(full)在6-311++G(d, p)和aug-cc-pvdz基组水平上, 对一系列简单的分子间π型卤键体系C2H4(C2H2)-XY(XY= F2、Cl2、Br2、ClF、BrF、BrCl) 进行构型全优化, 得到了T型卤键复合物. 结果表明MP2(full)/ 6-311++G(d, p)计算结果与实验结果较吻合. 并在MP2水平上计算了分子间的相互作用能, 用标准Counterpoise procedure (CP)方法对基函数迭加误差(BSSE)进行了校正. 利用电子密度拓扑分析方法对卤键复合物的拓扑性质进行了分析研究.     

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.     

Some measures for mediating the strengths of halogen bonds with the B–B bond in diborane(4) as an unconventional halogen acceptor

A series of halogen‐bonded complexes with diborane(4) 1 and its derivatives (Li 2 , methyl 3 , CN 4 ) as the halogen acceptors as well as with XCN, XCCH, XCF₃, XF (X = Cl, Br, I) as the halogen donors have been investigated by means of quantum chemical calculations at the MP2/aug‐cc‐pVTZ level. The result shows that the B?B bond is a good electron donor in halogen bonding, particularly for the halogen donor XF. Interestingly, for the halogen donor XF, the halogen bond becomes stronger in order of IF < BrF < ClF. It is found that the addition of electron‐donating substituents greatly strengthens the halogen bonding interaction to the point where it exceeds that of the majority of H‐bonds. When the N atom in 2 ‐BrCN is combined with another interaction, its strength has a further increase due to the cooperative effect. This study combines the boron compounds with halogen bonds and would be significant for expanding their applied fields. © 2013 Wiley Periodicals, Inc.     

The Role of Molecular Electrostatic Potentials in the Formation of a Halogen Bond in Furan⋅⋅⋅XY and Thiophene⋅⋅⋅XY Complexes

The halogen bonding of furan???XY and thiophene???XY (X=Cl, Br; Y=F, Cl, Br), involving σ‐ and π‐type interactions, was studied by using MP2 calculations and quantum theory of “atoms in molecules” (QTAIM) studies. The negative electrostatic potentials of furan and thiophene, as well as the most positive electrostatic potential (V_S,max) on the surface of the interacting X atom determined the geometries of the complexes. Linear relationships were found between interaction energy and V_S,max of the X atom, indicating that electrostatic interactions play an important role in these halogen‐bonding interactions. The halogen‐bonding interactions in furan???XY and thiophene???XY are weak, “closed‐shell” noncovalent interactions. The linear relationship of topological properties, energy properties, and the integration of interatomic surfaces versus V_S,max of atom X demonstrate the importance of the positive σ hole, as reflected by the computed V_S,max of atom X, in determining the topological properties of the halogen bonds.     

Cooperativity between the Halogen Bond and the Hydrogen Bond in H3N⋅⋅⋅XY⋅⋅⋅HF Complexes (X,Y=F,Cl, Br)

Ab initio calculations are used to provide information on H₃N???XY???HF triads (X, Y=F, Cl, Br) each having a halogen bond and a hydrogen bond. The investigated triads include H₃N???Br₂‐HF, H₃N???Cl₂???HF, H₃N???BrCI???HF, H₃N???BrF???HF, and H₃N???ClF???HF. To understand the properties of the systems better, the corresponding dyads are also investigated. Molecular geometries, binding energies, and infrared spectra of monomers, dyads, and triads are studied at the MP2 level of theory with the 6‐311++G(d,p) basis set. Because the primary aim of this study is to examine cooperative effects, particular attention is given to parameters such as cooperative energies, many‐body interaction energies, and cooperativity factors. The cooperative energy ranges from ?1.45 to ?4.64 kcal mol^?1, the three‐body interaction energy from ?2.17 to ?6.71 kcal mol^?1, and the cooperativity factor from 1.27 to 4.35. These results indicate significant cooperativity between the halogen and hydrogen bonds in these complexes. This cooperativity is much greater than that between hydrogen bonds. The effect of a halogen bond on a hydrogen bond is more pronounced than that of a hydrogen bond on a halogen bond.     

过渡金属M (M=Cu,Ag, Au)对X…Cl (X = F,Cl, Br)卤键相互作用强度的影响

采用量子化学方法,通过MCH2X…ClF(M=Cu,Ag,Au;X=F,Cl,Br)和CH3X…ClF两类复合物的对比,探讨了过渡金属对卤键相互作用强度的影响.CH3X…ClF复合物只有卤键相互作用,而优化MCH2X…ClF复合物除了得到一种只含有卤键相互作用的构型外,还得到一种含有过渡金属和Cl原子相互作用的稳定构型.含有过渡金属的复合物稳定性明显增加,Ag取代的复合物稳定性增加最为明显,Cu次之,Au最不明显.X原子最负分子表面静电势(MEP)减小是复合物稳定性增加的根本原因.利用自然键轨道(NBO)及分子中原子(AIM)分析进一步对体系的分子间相互作用进行了探讨.二阶稳定化能与键鞍点处拓扑性质的计算结果与相互作用能符合得很好.     

Electronic Structure and Properties of Trihalogen Cations X(3)(+) and XY(2)(+) (X, Y = F, Cl, Br, I)

Electronic structures, charge distributions, geometries, valence force constants, and vibrational frequencies of the homoatomic clusters F(3)(+), Cl(3)(+), Br(3)(+), and I(3)(+) and of the heteroatomic clusters ClF(2)(+), BrF(2)(+), IF(2)(+), BrCl(2)(+), ICl(2)(+), and IBr(2)(+) were determined. The self-consistent field approach extended by MP2-correlation energy or density-functional corrections was applied using various basis sets. It was found that d- and f-type polarization functions play a crucial role as in some other halogen compounds. The MP2 approach yields the most satisfactory results. The effect of the crystalline environment surrounding the Cl(3)(+), Br(3)(+), and I(3)(+) species is successfully simulated by a Madelung potential. Frequencies calculated in the crystal field are in reasonable agreement with the more reliable ones among the experimental results. Coupling force constants were determined. They are not consistent with some empirical rules. Bonding and charge distributions of the formally mixed-valence systems X(+)Y(2)(0) are discussed. X(+) behaves like a divalent chalcogen with high electronegativity. Each of the X(+)-Y bonds in XY(2)(+) is very similar to the bond in X-Y. We predict the experimentally unknown F(3)(+) to be stable in vacuum but not in the solid state. Structures and frequencies of XY(2)(+) species, which are as yet unknown, are also predicted.     

Nature of MoH···I bonds in Cp2Mo(L)H···I‐C≡C‐R Complexes (L=H,CN, PPh2, C(CH3)3; R=NO2, Cl,Br, H,OH, CH3, NH2)

The nature of the MoH···I bond in Cp₂Mo(L)H···I‐C≡C‐R (L= H, CN, PPh₂, C(CH₃)₃; R=NO₂, Cl, Br, H, OH, CH₃, NH₂) was investigated using electrostatic potential analysis, topological analysis of the electron density, energy decomposition analysis and natural bond orbital analysis. The calculated results show that MoH···I interactions in the title complexes belong to halogen‐hydride bond, which is similar to halogen bonds, not hydrogen bonds. Different to the classical halogen bonds, the directionality of MoH···I bond is low; Although electrostatic interaction is dorminant, the orbital interactions also play important roles in this kind of halogen bond, and steric interactions are weak; the strength of H···I bond can tuned by the most positive electrostatic potential of the I atom. As the electron‐withdrawing ability of the R substituent in the alkyne increases, the electrostatic potential maximum of the I atom increases, which enhances the strength of the H···I halogen bond, as well as the electron transfer.     

Computational studies of σ-type weak interactions between NCO/NCS radicals and XY(X = H,Cl;Y = F,Cl,and Br)

The triatomic radicals NCO and NCS are of interest in atmospheric chemistry,and both the ends of these radicals can potentially serve as electron donors during the formation of σ-type hydrogen/halogen bonds with electron acceptors XY(X = H,Cl;Y = F,Cl,and Br).The geometries of the weakly bonded systems NCO/NCS···XY were determined at the MP2/aug-cc-pVDZ level of calculation.The results obtained indicate that the geometries in which the hydrogen/halogen atom is bonded at the N atom are more stable than those where it is bonded at the O/S atom,and that it is the molecular electrostatic potential(MEP)-not the electronegativity-that determines the stability of the hydrogen/halogen bond.For the same electron donor(N or O/S) in the triatomic radical and the same X atom in XY,the bond strength decreases in the order Y = F > Cl > Br.In the hydrogen/halogen bond formation process for all of the complexes studied in this work,transfer of spin electron density from the electron donor to the electron acceptor is negligible,but spin density rearranges within the triatomic radicals,being transferred to the terminal atom not interacting with XY.     

Aromatization of triafulvene and its exocyclic Si,Ge, and Sn derivations by complexation with halogen atoms

The geometries of triafulvene (TF) and its exocyclic Si, Ge, and Sn analogues complexes with F, Cl, Br, and I halogen atoms (TF(X)···Y, X═C, Si, Ge, and Sn; Y═F, Cl, Br, and I) were studied. The complexes were optimized at DFT(B3LYP)/6–311+G(d,p) level of theory. To assess the aromaticity of the considered complexes the geometry-based (HOMA), magnetism-based (NICS), and recently introduced electronic-based (electric field gradient (EFG(0); Shannon aromaticity (SA)) aromaticity indices were employed. The increasing tendency of aromaticity in each complex species was noted as the series of TF(X)···F > TF(X)···Cl > TF(X)···Br > TF((X)···I. Then, the binding energies corrected by basis set super position error (BSSE) were calculated by single point energy calculations at M06-2X/6-311+G(d,p) level. Natural bond orbital (NBO) analysis confirmed that the charge transfer takes place from TF(X) to the halogen atoms. Some topological parameters, within the framework of the quantum theory of atoms in molecules (QTAIM), were also calculated to estimate the aromaticity of the complexes. It was seen that there are some important correlations between the topological parameters and aromaticity indices. In addition the most striking finding was that all the TF(X) molecules are connected with the halogen atoms through Y···C1═C2 (π) noncovalent interaction. This interaction was also investigated through noncovalent interaction (NCI) analysis.