Comparison of the Ability of N-Bases to Engage in Noncovalent Bonds

The lone pair of the N atom is a common electron donor in noncovalent bonds. Quantum calculations examine how various aspects of the base on which the N is located affect the strength and other properties of complexes formed with Lewis acids FH, FBr, F₂Se, and F₃As that respectively encompass hydrogen, halogen, chalcogen, and pnicogen bonds. In most cases the halogen bond is the strongest, followed in order by chalcogen, hydrogen, and pnicogen. The noncovalent bond strength increases in the sp<sp²<sp³ order of hybridization of N. Replacement of H substituents on the base by a methyl group or substituting N by C atom to which the base N is attached, strengthens the bond. The strongest bonds occur for trimethylamine and the weakest for N₂.     

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.     

H2XP…SHY复合物中磷键与硫键的理论研究

用从头算量子化学方法MP2 与CCSD(T)研究了H₂XP和SHY (X, Y=H, F, Cl, Br)分子的P与S之间形成的磷键X―P…S与硫键Y―S…P的本质与规律以及取代基X与Y对成键的影响. 计算结果表明, 硫键比磷键强, 连接在Lewis 酸上的取代基的电负性增大导致形成的磷键或硫键增强, 键能增大, 对单体的结构和性质的影响也增大; 而连接在Lewis 碱上的取代基效应则相反. 硫键键能为8.37-23.45 kJ·mol^-1, 最强的硫键结构是Y 电负性最大而X 电负性最小的HFS…PH₃, CCSD(T)计算的键能是16.04 kJ·mol^-1; 磷键键能为7.54-14.65 kJ·mol^-1, 最强的磷键结构是X 电负性最大而Y 电负性最小的H₂FP…SH₂, CCSD(T)计算的键能是12.52 kJ·mol^-1. 对磷键与硫键能量贡献较大的是交换与静电作用. 分子间超共轭lp(S)-σ^*(PX)与lp(P)-σ^*(SY)对磷键与硫键的形成起着重要作用, 它导致单体的极化, 其中硫键的极化效应较大, 从而有一定的共价特征.     

Theoretical Density Functional Theory insights into the nature of chalcogen bonding between CX2 (X = S,Se, Te) and diazine from monomer to supramolecular complexes

Chalcogen bonding is a noncovalent interaction, highly similar to halogen and hydrogen bonding, occurring between a chalcogen atom and a nucleophilic region. Two density functional theory (DFT) approaches B3LY-D3 and B97-D3 were performed on a series of complexes formed between CX₂ (X = S, Se, Te) and diazine (pyridazine, pyrimidine and pyrazine). Chalcogen atoms prefer interacting with the lone pair of a nitrogen atom rather than with the π-cloud of an aromatic ring. CTe₂ and CSe₂ form a stronger chalcogen bond than CS₂. The electrostatic potential of CX₂ (X = S, Se and Te) reveals the presence of two equivalent σ-holes, one on each chalcogen atom. These CX₂ molecules interact with diazine giving rise to supramolecular interactions. Wiberg bond index and second-order perturbation theory analysis in NBO were performed to better understand the nature of the chalcogen bond interaction.     

Activating Chalcogen Bonding (ChB) in Alkylseleno/Alkyltelluroacetylenes toward Chalcogen Bonding Directionality Control

Activation of a deep electron-deficient area on chalcogen atoms (Ch=Se, Te) is demonstrated in alkynyl chalcogen derivatives, in the prolongation of the (C≡)C−Ch bond. The solid-state structures of 1,4-bis(methylselenoethynyl)perfluorobenzene ( 1Se ) show the formation of recurrent chalcogen-bonded (ChB) motifs. Association of 1Se and the tellurium analogue 1Te with 4,4′-bipyridine and with the stronger Lewis base 1,4-di(4-pyridyl)piperazine gives 1:1 co-crystals with 1D extended structures linked by short and directional ChB interactions, comparable to those observed with the corresponding halogen bond (XB) donor, 1,4-bis(iodoethynyl)-perfluorobenzene. This “alkynyl” approach for chalcogen activation provides the crystal-engineering community with efficient, and neutral ChB donors for the elaboration of supramolecular 1D (and potentially 2D or 3D) architectures, with a degree of strength and predictability comparable to that of halogen bonding in iodoacetylene derivatives.     

Revealing substituent effects on the concerted interaction of pnicogen, chalcogen, and halogen bonds in substituted s-triazine ring

Quantum chemical calculations have been performed to gauge the effect of substituents on concerted interactions of pnicogen, chalcogen, and halogen bonds in the X–TAZ···Y complexes (X = CN, F, Cl, Br, H, CH₃, OH, and NH₂, where TAZ and Y denote s-triazine ring and P, S, and Cl atoms, respectively) at the M06-2X/aug-cc-pVDZ level. The mutual interplay of these interactions is also investigated. The results indicate that diminutive effects are observed when the three kinds of noncovalent interactions pnicogen, chalcogen, and halogen bonds are coexisted in the complexes. These effects are studied in terms of energetic and geometric features of the complexes. In addition, Bader’s theory of “atoms in molecules” is used to analyze their strength of varying electron density at bond critical points. Natural bond orbital (NBO) theory is used to characterize the orbital interactions. The results indicate that the electron-withdrawing/donating substituents decrease/increase the magnitude of the binding energies compared to the unsubstituted X–TAZ···Y (X = H) complex. Good correlations among binding energies, Hammett constants, geometrical, atoms in molecular and NBO parameters are established in X–TAZ···Y complexes. By taking advantage of all the aforementioned computational methods, this study examines how these interactions mutually influence each other.     

二元和三元复合物中阳离子硫键与磷键的理论研究

采用MP2/cc-pVDZ和cc-pVTZ基组分别对复合物XH₂S⁺…NCH₂P和NCH₂P…PyX(X=NH₂, CH₃, H, CN, F, Cl, Br)中的硫键和磷键进行了研究, 讨论了键长、 键临界点的电荷密度(ρ)、 拉普拉斯密度(▽²ρ)、 范德华表面穿透距离、 二阶稳定化能和电荷转移量对硫键和磷键相互作用能的影响. 结果表明, 当取代基X为吸电子基团时, 形成的硫键较强. 当X为给电子基团时, 形成的磷键较强. 利用能量分解方法分析了取代基—CN导致硫键稳定性反常的可能原因. 还进一步讨论了三元复合物H₃S⁺…NCH₂P…PyX(X=NH₂, CH₃, H, CN, F, Cl, Br)中硫键和磷键的协同相互作用以及取代基对复合物稳定性的影响. 并通过对比相同的2种单体在三元复合物和二元复合物中的二阶稳定化能和相互作用能的差值, 说明了硫键与磷键起到相互促进的正协同作用, 增强了三元复合物的稳定性.     

The competition of Y⋯o and X⋯n halogen bonds to enhance the group V σ‐hole interaction in the NCY⋯oPH3⋯NCX and OPH3⋯NCX⋯NCY (X,YF,Cl, and Br) complexes

The positive electrostatic potentials (ESP) outside the σ‐hole along the extension of O? P bond in O?PH₃ and the negative ESP outside the nitrogen atom along the extension of the C? N bond in NCX could form the Group V σ‐hole interaction O?PH₃?NCX. In this work, the complexes NCY?O?PH₃?NCX and O?PH₃?NCX?NCY (X, Y?F, Cl, Br) were designed to investigate the enhancing effects of Y?O and X?N halogen bonds on the P?N Group V σ‐hole interaction. With the addition of Y?O halogen bond, the V _{S, max} values outside the σ‐hole region of O?PH₃ becomes increasingly positive resulting in a stronger and more polarizable P?N interaction. With the addition of X?N halogen bond, the V _{S, min} values outside the nitrogen atom of NCX becomes increasingly negative, also resulting in a stronger and more polarizable P?N interaction. The Y?O halogen bonds affect the σ‐hole region (decreased density region) outside the phosphorus atom more than the P?N internuclear region (increased density region outside the nitrogen atom), while it is contrary for the X?N halogen bonds. © 2015 Wiley Periodicals, Inc.     

Competition Between the Two σ-Holes in the Formation of a Chalcogen Bond

A chalcogen atom Y contains two separate σ-holes when in a R₁YR₂ molecular bonding pattern. Quantum chemical calculations consider competition between these two σ-holes to engage in a chalcogen bond (ChB) with a NH₃ base. R groups considered include F, Br, I, and tert-butyl (tBu). Also examined is the situation where the Y lies within a chalcogenazole ring, where its neighbors are C and N. Both electron-withdrawing substituents R₁ and R₂ act cooperatively to deepen the two σ-holes, but the deeper of the two holes consistently lies opposite to the more electron-withdrawing group, and is also favored to form a stronger ChB. The formation of two simultaneous ChBs in a triad requires the Y atom to act as double electron acceptor, and so anti-cooperativity weakens each bond relative to the simple dyad. This effect is such that some of the shallower σ-holes are unable to form a ChB at all when a base occupies the other site.     

A Quantitative Molecular Orbital Perspective of the Chalcogen Bond

We have quantum chemically analyzed the structure and stability of archetypal chalcogen-bonded model complexes D₂Ch⋅⋅⋅A⁻ (Ch = O, S, Se, Te; D, A = F, Cl, Br) using relativistic density functional theory at ZORA-M06/QZ4P. Our purpose is twofold: (i) to compute accurate trends in chalcogen-bond strength based on a set of consistent data; and (ii) to rationalize these trends in terms of detailed analyses of the bonding mechanism based on quantitative Kohn-Sham molecular orbital (KS-MO) theory in combination with a canonical energy decomposition analysis (EDA). At odds with the commonly accepted view of chalcogen bonding as a predominantly electrostatic phenomenon, we find that chalcogen bonds, just as hydrogen and halogen bonds, have a significant covalent character stemming from strong HOMO−LUMO interactions. Besides providing significantly to the bond strength, these orbital interactions are also manifested by the structural distortions they induce as well as the associated charge transfer from A⁻ to D₂Ch.     

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.     

NCO/NCS自由基与XY（X=H,Cl;Y=F,Cl,Br）间σ-型弱作用的理论研究

NCO和NCS是大气化学中非常引人关注的自由基,它们均有三个原子并且两个端基原子均可作为电子给体形成σ-型氢/卤键.本文在MP2/aug-cc-pVDZ水平上研究了NCO/NCS...XY（X=H,Cl;Y=F,Cl,Br）体系中的弱化学键.计算结果表明,氢/卤原子与N原子相连形成的复合物比与O/S原子相连形成的复合物稳定;氢/卤键的稳定性由分子静电势决定,而非原子电负性;对相同的电子给体B（B=N,O/S）和相同的卤原子来说,化学键的强度按Y=F,Cl,Br的顺序逐渐减弱.在氢/卤键形成过程中,自旋电子密度在电子给体和电子受体间的转移较少,但它在自由基内部发生重排,就本文研究的所有复合物而言,自旋电子密度均转移向XY分子的相反位置.     

The Lewis basicity of nitrido complexes. Theoretical investigation of the structure and bonding of Cl2 (PH3)3ReN–X (X = BH3, BCl3, BBr3, AlH3, AlCl3, AlBr3, GaH3, GaCl3, GaBr3, O, S, Se, Te)

Quantum chemical calculations using gradient-corrected density functional theory (B3LYP) and ab initio methods at the MP2 level are reported for the geometries and bond energies of the nitrido complexes Cl₂ (PH₃)₃ReN–X (X = BH₃, BCl₃, BBr₃, AlH₃, AlCl₃, AlBr₃, GaH₃, GaCl₃, GaBr₃, O, S, Se, Te). The theoretical geometries are in excellent agreement with experimental values of related complexes which have larger phosphine ligands. The parent nitrido complex Cl₂(PH₃)₃ReN is a very strong Lewis base. The calculated bond dissociation energy of Cl₂(PH₃)₃ReN–AlCl₃ is D _e = 43.7 kcal/mol, which is nearly as high as the bond energy of Me₃N–AlCl₃. The donor-acceptor bonds of the other Cl₂(PH₃)₃ReN–AY₃ complexes are also very strong. Even stronger N–X bonds are predicted for most of the nitrido-chalcogen complexes, which exhibit the trend X = O ≫ S > Se > Te. Analysis of the electronic structure shows that the parent compound Cl₂(PH₃)₃ReN has a Re–N triple bond. The Re–N σ bond is clearly polarized towards nitrogen, while the two π bonds are nearly nonpolar. The Re–N σ and π bonds become more polarized toward nitrogen when a Lewis acid or a chalcogen atom is attached. Bonding in AY₃ complexes should be described as Cl₂(PH₃)₃ReE≡N→AY₃, while the chalcogen complexes should be written with double bonds Cl₂(PH₃)₃Re=N=X. The charge-decomposition analysis indicates that the nitrogen-chalcogen bonds of the heavier chalcogen complexes with X = S, Se, Te can also be interpreted as donor-acceptor bonds between the nitrido complex acting as a Lewis base and the chalcogen atom with an empty p(σ) orbital acting as a Lewis acid. The nitrido oxo complex Cl₂(PH₃)₃ Re=N=O has a covalent N–O double bond. Received: 27 July 1998 / Accepted: 26 October 1998 / Published online: 16 March 1999     

Halogen bonding in mono- and dihydrated halobenzene

Density functional theory calculations were performed on halogen-bonded and hydrogen-bonded systems consisting of a halobenzene (XPh; X = F, Cl, Br, I, and At) and one or two water molecules, using the M06-2X density functional with the 6-31+G(d) (for C, H, F, Cl, and Br) and aug-cc-pVDZ-PP (for I, At) basis sets. The counterpoise procedure was performed to counteract the effect of basis set superposition error. The results show halogen bonds form in the XPh-H₂O system when X > Cl. There is a trend toward stronger halogen bonding as the halogen group is descended, as assessed by interaction energy and X•••O_w internuclear separation (where O_w is the water oxygen). For all XPh-H₂O systems hydrogen-bonded systems exist, containing a combination of CH•••O_w and O_wH_w•••X hydrogen bonds. For all systems except X = At the X•••H_w hydrogen-bonding interaction is stronger than the X•••O_w halogen bond. In the XPh-(H₂O)₂ system halogen bonds form only for X > Br. The two water molecules prefer to form a water dimer, either located around the C H bond (for X = Br, At, and I) or located above the benzene ring (for all halogens). Thus, even in the absence of competing strong interactions, halogen bonds may not form for the lighter halogens due to (1) competition from cooperative weak interactions such as C H•••O and OH•••X hydrogen bonds, or (2) if the formation of the halogen bond would preclude the formation of a water dimer. © 2018 Wiley Periodicals, Inc.     

Weak interactions between hypohalous acids and dimethylchalcogens

The complexes formed between dimethylchalcogens X(CH(3))(2) (X = S, Se, and Te) and hypohalous acids YOH (Y = F, Cl, Br, and I) have been studied at the MP2/aug'-cc-pVTZ computational level, five minima structures being located. Two of them correspond to hydrogen bonds (HB), another two to halogen bonds (XB) with the chalcogen acting as an electron donor, the last one showing a C-H···O contact. The most stable complexes of IOH and BrOH acids present halogen···chalcogen interactions with interaction energies, E(i), up to -49 kJ mol(-1). In the case of the ClOH and FOH molecules, the hydrogen bonded complexes are more stable with interaction energies between -27 and -34 kJ mol(-1). Linear correlations between the molecular electrostatic potential (MEP) stationary points at the van der Waals surface and the interaction energy have been found. The contribution of the different energy terms to the total interaction energy was analyzed by means of the DFT-SAPT theory finding that the electrostatic attractive term is dominant in the complexes with HB and XB, excepting a few cases in which the dispersion and induction terms become more important than the electrostatic one.     

Multinuclear NMR studies on CF3 substituted sulfur,selenium and tellurium compounds

A systematic investigation of thirty-four CF₃Se(II, IV) and eight CF₃Te(II, IV) compounds by ¹³C, ¹⁹F, ⁷⁷Se and ¹²⁵Te NMR spectroscopy resulted in some general features for chemical shifts and coupling constants which agree with the trends of reported ¹⁹F and new ¹³C NMR data of CF₃S(II, IV) compounds. Moreover, the NMR spectra of molecules of the type E=CXY (E = chalcogen, X, Y = halogen) and substances containing a CSe double bond have been studied. From the comparison of these NMR data with those of CF₃ substituted chalcogen compounds, a partial double bond character of the carbon-fluorine and carbon-chalcogen bond in CF₃ substituted chalcogen compounds can be derived:

     

σ-Hole bond tunability in YO<Subscript>2</Subscript>X<Subscript>2</Subscript>:NH<Subscript>3</Subscript> and YO<Subscript>2</Subscript>X<Subscript>2</Subscript>:H<Subscript>2</Subscript>O complexes (X = F,Cl, Br; Y = S,Se): trends and theoretical aspects

A σ-hole is defined as an electron-deficient region on the extension of a covalently bonded group IV–VII atoms. If the electronic density in the σ-hole is sufficiently low, then this region will have a positive electrostatic potential, which allows attractive noncovalent interactions with negative sites. SO₂X₂ and SeO₂X₂ (X = F, Cl and Br) have three Lewis acid sites of σ-hole located in the outermost of chalcogen atom and X end, participating in the chalcogen and halogen bonds with NH₃ and H₂O, respectively. MP2/aug-cc-pVTZ and M06-2X/aug-cc-pVTZ calculations reveal that for a given halogen atom, SeO₂X₂ forms stronger chalcogen bond interactions than SO₂X₂ counterpart. Almost a perfect linear relationship is evident between the interaction energies and the magnitudes of the product of most positive and negative electrostatic potentials. The interaction energies calculated by M06-2X and MP2 methods are almost consistent with each other.     

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.     

Intermolecular Weak Interactions in HTeXH Dimers (X=O,S, Se,Te): Hydrogen Bonds,Chalcogen–Chalcogen Contacts and Chiral Discrimination

A theoretical study of the HTeXH (X=O, S, Se and Te) monomers and homodimers was carried out by means of second‐order Møller‐Plesset perturbation theory (MP2) computational methods. In the case of monomers, the isomerization energy from HTeXH to H₂Te=X and H₂X=Te (X=O, S, Se, and Te) and the rotational transition‐state barriers were obtained. Due to the chiral nature of these compounds, homo and heterochiral dimers were found. The electron density of the complexes was characterized with the atoms‐in‐molecules (AIM) methodology, finding a large variety of interactions. The charge transfer within the dimers was analyzed by means of natural bond orbitals (NBO). The density functional theory‐symmetry adapted perturbation theory (DFT‐SAPT) method was used to compute the components of the interaction energies. Hydrogen bonds and chalcogen–chalcogen interactions were characterized and their influence analyzed concerning the stability and chiral discrimination of the dimers.     

An ab initio investigation of chalcogen–hydride interactions involving HXeH as a chalcogen bond acceptor

This work presents an ab initio study on chalcogen–hydride interactions in several binary complexes of chalcogen-containing molecules with HXeH. The geometries, H–Xe stretching frequencies and interaction energies of XCY···HXeH binary complexes are investigated at MP2/aug-cc-pVTZ and CCSD(T)/aug-cc-pVTZ levels of theory, where X = O, S, Se, Te and Y = S, Se, Te. For each XCY···HXeH complex, a chalcogen–hydride bond is formed between the negatively charged hydrogen atom of the HXeH molecule and the most positive electrostatic potential region (σ-hole) on the surface of the interacting atom Y. Upon complex formation, a notable blue shift is found for the H–Xe stretch vibration. This result reveals that there is a stronger H^?(XeH)⁺ ion-pair character in XCY···HXeH complexes than in free HXeH molecule. In order to shed light on the origin of the chalcogen–hydride interactions, molecular electrostatic potential, quantum theory of atoms in molecules and interaction energy decomposition analyses are performed. Cooperative effects between a conventional chalcogen bond and the chalcogen–hydride interaction in OCY···OCY···HXeH complexes are also investigated.