Application of polarizable ellipsoidal force field model to pnicogen bonds

Noncovalent interactions, such as hydrogen bonds and halogen bonds, are frequently used in drug designing and crystal engineering. Recently, a novel noncovalent pnicogen bonds have been identified as an important driving force in crystal structures with similar bonding mechanisms as hydrogen bond and halogen bond. Although the pnicogen bond is highly anisotropic, the pnicogen bond angles range from 160° to 180° due to the complicated substituent effects. To understand the anisotropic characters of pnicogen bond, a modification of the polarizable ellipsoidal force field (PEff) model previously used to define halogen bonds was proposed in this work. The potential energy surfaces (PESs) of mono‐ and polysubstituted PH₃–NH₃ complexes were calculated at CCSD(T), MP2, and density functional theory levels and were used to examine the modified PEff model. The results indicate that the modified PEff model can precisely characterize pnicogen bond. The root mean squared error of PES obtained with PEff model is less than 0.5 kcal/mol, compared with MP2 results. In addition, the modified PEff model may be applied to other noncovalent bond interactions, which is important to understand the role of intermolecular interactions in the self‐assembly structures. © 2015 Wiley Periodicals, Inc.     

Ab initio study of the complexes of halogen-containing molecules RX (X=Cl, Br, and I) and NH3: towards understanding the nature of halogen bonding and the electron-accepting propensities of covalently bonded halogen atoms

Ab initio calculations have been performed on a series of complexes formed between halogen-containing molecules and ammonia to gain a deeper insight into the nature of halogen bonding. It appears that the dihalogen molecules form the strongest halogen-bonded complexes with ammonia, followed by HOX; the charge-transfer-type contribution has been demonstrated to dominate the halogen bonding in these complexes. For the complexes involving carbon-bound halogen molecules, our calculations clearly indicate that electrostatic interactions are mainly responsible for their binding energies. Whereas the halogen-bond strength is significantly enhanced by progressive fluorine substitution, the substitution of a hydrogen atom by a methyl group in the CH(3)X...NH(3) complex weakened the halogen bonding. Moreover, remote substituent effects have also been noted in the complexes of halobenzenes with different para substituents. The influence of the hybridization state of the carbon atom bonded to the halogen atom has also been examined and the results reveal that halogen-bond strengths decrease in the order HC triple bond CX > H(2)C=CHX approximately O=CHX approximately C(6)H(5)X > CH(3)X. In addition, several excellent linear correlations have been established between the interaction energies and both the amount of charge transfer and the electrostatic potentials corresponding to an electron density of 0.002 au along the R-X axis; these correlations provide good models with which to evaluate the electron-accepting abilities of the covalently bonded halogen atoms. Finally, some positively charged halogen-bonded systems have been investigated and the effect of the charge has been discussed.     

Complexes of HArF and AuX (X = F,Cl, Br,I). Comparison of H-bonds,halogen bonds,F-shared bonds and covalent bonds

The various sorts of complexes in which HArF and AuX (X = F, Cl, Br, I) can engage are probed by MP2/aug-cc-pVTZ calculations. The most weakly bound are those containing a halogen bond (XB) of the AuX⋯FArH sort, with binding energies less than 8 kcal/mol. H-bonded dimers FArH⋯XAu are a little stronger, held together by some 12 kcal/mol. Being the most strongly bound places the F atom of HArF roughly midway between Ar and Au in an F-shaped structure, bound by some 43–54 kcal/mol. The last sort of product involves atomic rearrangements wherein the H atom migrates from Ar to Au, followed by formation of a covalent Ar–Au bond. The resulting molecular unit is stabilized by 30–40 kcal/mol relative to the original HArF and AuX reactants. The H-bonded dimers are held together by an unusually large polarization component, surpassing electrostatic attraction, while dispersion predominates for the halogen bonds. Perturbations of the geometries and stretching frequencies offer a ready means of distinguishing the different types of complexes by spectroscopic techniques.     

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.     

Theoretical investigations into the blue-shifting hydrogen bond in benzene complexes.

The benzene...X complexes (X=benzene, antracene, ovalene) were optimised at the MP2/6-31G** level with the C2v symmetry of the complex and planarity of the proton acceptor being preserved. The resulting stabilisation energies amount to 1.2, 2.3 and 2.9 kcal mol(-1), and the C-H bond of the proton donor is contracted by 0.0035, 0.0052 and 0.0055 A, respectively. The contraction is connected with a blue-shift of the C-H stretch vibration frequency. A two-dimensional anharmonic vibration treatment based on a MP2/6-31G** potential energy surface yields the following blue shifts for the complexes studied: 28, 42 and 43 cm(-1). The dominant attraction in the complexes is London dispersion, while the attractive contribution from electrostatic quadrupole-quadrupole interactions is considerably smaller.     

Computational study of the process of hydrogen bond breaking: the case of the formamide-formic acid complex

MP2/6-311++G(d,p) and B3LYP/6-311++G(d,p) quantum calculations are used to study the formamide-formic acid complex (FFAC), a system bound by two hydrogen bonds, N--H...O and O--H...O, that forms a bond ring at equilibrium. When the intermolecular separation between monomers R increases, this ring opens at a distance for which the weaker N--H...O bond breaks remaining the stronger O--H...O bond. The computational study characterizes that process addressing changes of interaction energy DeltaE, structure and properties of the electron density rho(r) as well as spatial distributions of rho(r), the electrostatic potential U(r), and the electron localization function eta(r). It is shown that the spatial derivatives of DeltaE, the topology of rho(r), and qualitative changes noticed in U(r) = 0 isocontours allow to identify a precise distance R for which one can say the N--H...O hydrogen bond has broken. Both levels of theory predict essentially the same changes of structure and electron properties associated to the process of breaking and virtually identical distances at which it takes place.     

Comparison between Hydrogen and Halogen Bonds in Complexes of 6-OX-Fulvene with Pnicogen and Chalcogen Electron Donors

Quantum chemical calculations are applied to complexes of 6-OX-fulvene (X=H, Cl, Br, I) with ZH₃/H₂Y (Z=N, P, As, Sb; Y=O, S, Se, Te) to study the competition between the hydrogen bond and the halogen bond. The H-bond weakens as the base atom grows in size and the associated negative electrostatic potential on the Lewis base atom diminishes. The pattern for the halogen bonds is more complicated. In most cases, the halogen bond is stronger for the heavier halogen atom, and pnicogen electron donors are more strongly bound than chalcogen. Halogen bonds to chalcogen atoms strengthen in the order O<S<Se<Te, whereas the pattern is murkier for the pnicogen donors. In terms of competition, most halogen bonds to pnicogen donors are stronger than their H-bond analogues, but there is no clear pattern with respect to chalcogen donors. O prefers a H-bond, while halogen bonds are favored by Te. For S and Se, I-bonds are strongest, followed Br, H, and Cl-bonds in that order.     

The role of the halogen bond in iodothyronine deiodinase: Dependence on chalcogen substitution in naphthyl-based mimetics

The effects on the activity of thyroxine (T4) due to the chalcogen replacement in a series of peri-substituted naphthalenes mimicking the catalytic function of deiodinase enzymes are computationally examined using density functional theory. In particular, T4 inner-ring deiodination pathways assisted by naphthyl-based models bearing two tellurols and a tellurol-thiol pair in peri-position are explored and compared with the analogous energy profiles for the naphthalene mimic having two selenols. The presence of a halogen bond (XB) in the intermediate formed in the first step and involved in the rate-determining step of the reaction is assumed to facilitate the process increasing the rate of the reaction. The rate-determining step calculated energy barrier heights allow rationalizing the experimentally observed superior catalytic activity of tellurium containing mimics. Charge displacement analysis is used to ascertain the presence and the role of the electron density charge transfer occurring in the rate-determining step of the reaction, suggesting the incipient formation or presence of a XB interaction. © 2019 Wiley Periodicals, Inc.     

Ab initio studies of electron acceptor-donor interactions with blue- and red-shifted hydrogen bonds.

The features of blue- and red-shifted electron acceptor-donor (ACH/B) hydrogen bonds have been compared by using quantum chemical calculations. The geometry, the interaction energy and the vibrational frequencies of both blue- (ACH=F3CH, Cl3CH with B=FCD3) and red-shifted (ACH=F3CH, Cl3CH with B=NH3 and ACH=CH3CCH with B=FCD3, NH3) complexes were obtained by using ab initio MP2(Full)/6-31+G(d,p) calculations with the a priori basis-set superposition error (BSSE) correction method. One-dimensional potential energy and dipole moment functions of the dimensionless normal coordinate Q1, corresponding to the CH stretching mode of ACH, have been compared for both types of complexes. Contributions of separate components of the interaction energy to the frequency shift and the effect of electron charge transfer were examined for a set of intermolecular distances by using the symmetry-adapted perturbation theory (SAPT) approach and natural bond orbitals (NBO) population analysis.     

Accurate binding energies of hydrogen,halogen, and dihydrogen bonded complexes and cation enhanced binding strengths

Interaction energy (E_int) values of a variety of hydrogen, halogen, and dihydrogen bonded complexes in the weak, medium, and strong regimes have been computed using W1BD, MP2, M06L density functional theory, and hybrid methods MP4//MP2, MP4//M06L, and CCSD(T)//MP2. W1BD level E_int and CCSD(T) results reported in the literature show very good agreement (mean absolute deviation = 0.19 kcal/mol). MP2 underestimates E_int while M06L shows accurate behavior for all except halogen and charge‐assisted hydrogen bonds. MP4//MP2, MP4//M06L, and CCSD(T)//MP2 yield E_int very close to those obtained from W1BD. The high accuracy energy data at MP4/MP2 is used to study the effect of a cation (Li⁺, NH₄⁺) on the E_int. The cation enhances electron donation from the donor to noncovalent bonding region leading to substantial enhancement in E_int (～141–566% for Li⁺ and ～105–539% for NH₄⁺) and promotes a noncovalent bond in the weak regime to medium regime and that in the medium regime to strong regime. © 2014 Wiley Periodicals, Inc.     

On differences between hydrogen bonding and improper blue-shifting hydrogen bonding.

Twenty two hydrogen-bonded and improper blue-shifting hydrogen-bonded complexes were studied by means of the HF, MP2 and B3LYP methods using the 6-31G(d,p) and 6--311 ++G(d,p) basis sets. In contrast to the standard H bonding, the origin of the improper blue-shifting H bonding is still not fully understood. Contrary to a frequently presented idea, the electric field of the proton acceptor cannot solely explain the different behavior of the H-bonded and improper blue-shifting H-bonded complexes. Compression of the hydrogen bond due to different attractive forces-dispersion or electrostatics--makes an important contribution as well. The symmetry-adapted perturbation theory (SAPT) has been utilized to decompose the total interaction energy into physically meaningful contributions. In the red-shifting complexes, the induction energy is mostly larger than the dispersion energy while, in the case of blue-shifting complexes, the situation is opposite. Dispersion as an attractive force increases the blue shift in the blue-shifting complexes as it compresses the H bond and, therefore, it increases the Pauli repulsion. On the other hand, dispersion in the red-shifting complexes increases their red shift.     

Building a Better Halide Receptor: Optimum Choice of Spacer,Binding Unit,and Halosubstitution

Quantum calculations are used to measure the binding of halides to a number of bipodal dicationic receptors, constructed as a pair of binding units separated by a spacer group. A number of variations are studied. A H atom on each binding unit (imidazolium or triazolium) is replaced by Br or I. Benzene, thiophene, carbazole, and dimethylnaphthalene are considered as spacer groups. Each receptor is paired with halides F^?, Cl^?, Br^?, and I^?. Substitution with I on the binding unit yields a large enhancement of binding, as much as 13 orders of magnitude; a much smaller increase occurs for substitution with Br. Imidazolium is a more effective binding agent than is triazolium. Benzene and dimethylnaphthalene represent the best spacers, followed by thiophene and carbazole. F^? binds much more strongly than do the other halides, which obey the order Cl^?>Br^?>I^?.     

A new method for quick predicting the strength of intermolecular hydrogen bonds

A new method is proposed to quick predict the strength of intermolecular hydrogen bonds. The method is employed to produce the hydrogen-bonding potential energy curves of twenty-nine hydrogen-bonded dimers. The calculation results show that the hydrogen-bonding potential energy curves obtained from this method are in good agreement with those obtained from MP2/6-31+G** calculations by including the BSSE correction, which demonstrate that the method proposed in this work can be used to calculate the hydrogen-bonding interactions in peptides. Supported by the National Natural Science Foundation of China (Grants Nos. 20573049 and 20633050) and the research fund of the Educational Department of Liaoning Province (2004C019, 20060469)     

Tuning the Competition between Hydrogen and Tetrel Bonds by a Magnesium Bond

A computational study of the complexes formed by TF₃OH (T=C, Si, Ge) with three nitrogen-containing bases NCH, NH₃, and imidazole (IM) is carried out at the MP2/aug-cc-pVTZ level. TF₃OH can participate in two different types of noncovalent interactions: a hydrogen bond (HB) involving the hydroxyl proton and a tetrel bond (TB) with the tetel atom T. The strength of the HB is largely unaffected by the identity of T while the TB is enhanced as T grows larger. The HB is preferred over the TB for most systems, with the exception of GeF₃OH with either NH₃ or IM. MgCl₂ engages in a Mg⋅⋅⋅O Magnesium bond (Mg-bond) with the TF₃OH O atom, which cooperatively enhances both the HB and TB. The HB strengthening is particularly large for the NH₃ or IM bases, and especially for CF₃OH, but is slowly reduced as the T atom grows larger. The TB enhancement, on the other hand, behaves in the opposite fashion, accelerating for the larger T atoms. As a bottom line, the Mg-bond generally reinforces and accentuates the preference for the HB or TB that is already present in the dimer. The Mg-bond is also responsible for a proton transfer in the HB configurations with NH₃ and IM.