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.     

Contra-thermodynamic behavior in intermolecular hydrogen transfer of alkylperoxy radicals.

Quantum chemical investigation of bimolecular hydrogen transfer involving alkylperoxy radicals, a key reaction family in the free-radical oxidation of hydrocarbons, was performed to establish structure-reactivity relationships. Eight different reactions were investigated featuring four different alkane substrates (methane, ethane, propane and isobutane) and two different alkylperoxy radicals (methylperoxy and iso-propylperoxy). Including forward and reverse pairs, sixteen different activation energies and enthalpies of reaction were used to formulate structure-reactivity relationships to describe this chemistry. We observed that the enthalpy of formation of loosely bound intermediate states has a strong inverse correlation with the overall heat of reaction and that this results in unique contra-thermodynamic behavior such that more exothermic reactions have higher activation barriers. A new structure-reactivity relationship was proposed that fits the calculated data extremely well: E(A)=E(o)+alphaDeltaH(rxn) where alpha=-0.10 for DeltaH(rxn)<0, and alpha=1.10 for DeltaH(rxn)>0 and E(o)=3.05 kcal mol(-1).     

Properties of halogen atoms for representing intermolecular electrostatic interactions related to halogen bonding and their substituent effects

The electrostatic properties of halogen atoms are studied theoretically in relation to their ability of halogen bonding, which is an attractive intermolecular interaction of a covalently bonded halogen atom with a negatively charged atom of a neighboring molecule. The electric quadrupole (of electronic origin) with a positive zz component Θ_zz of a covalently bonded halogen atom, where the z axis is taken along the covalent bond involving the halogen atom, is mainly responsible for the attractive electrostatic interaction with a negatively charged atom. This positive Θ_zz is an intrinsic property of halogen atoms with the p_x²p_y²p_z configuration of the valence electronic shell, as shown by ab initio molecular orbital calculations for isolated halogen atoms with this electronic configuration, and increases in the order of F < Cl < Br < I, in parallel with the known general sequence of the strength of halogen bonding. For halogen‐containing aromatic compounds, the substituent effects on the electrostatic properties are also studied. It is shown that the magnitude of Θ_zz and the electric field originating from it are rather insensitive to the substituent effect, whereas the electric field originating from atomic partial charges has a large substituent effect. The latter electric field tends to partially cancel the former. The extent of this partial cancellation is reduced in the order of Cl < Br < I and is also reducible by proper substitution on or within the six‐membered ring of halobenzene. Perspectives on the development of potential function parameters applicable to halogen‐bonding systems are also briefly discussed. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

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.     

Theoretical study of the intermolecular hydrogen bond interaction for furan-HCl and furan-CHCl_3 complexes

The importance of intermolecular interactions in biology and material science has prompted chemists to explore the nature of the variety of such interactions. The strongest of these interac-tions are the hydrogen bonds, which play an important role in determining the molecular confor-mation, crystal packing, and the structure of biological systems such as nucleic acids. Extensive experimental and theoretical efforts[1—5] have been devoted to the studies of this type of interac-tions, such as …     

Substituent effects on noncovalent halogen/π interactions: Theoretical study

Noncovalent halogen/π interactions of FCl with substituted benzenes have been investigated using ab initio calculations. It was shown that the predicted maximum interaction energy gap between the substituted and unsubstituted systems amounts to 1.14 kcal/mol, and therefore substituents on benzene have a pronounced effect on the strength of halogen/π interactions. While the presence of electron‐donating groups (NH₂, CH₃, and OH) on benzene enhances the interaction energy appreciably, an opposite effect is observed for electron‐accepting groups (NO₂, CN, Br, Cl, and F). The large gain of the attraction by electron correlation illustrates that the stabilities of the systems considered arise primarily from the dispersion interaction. Beside the dispersion interaction, the charge‐transfer interaction also plays an important role in halogen/π interactions, as a charge density analysis suggested. To provide more insight into the nature of halogen/π interactions, topological analysis of the electron density distribution and properties of bond critical points were determined in terms of the atoms in molecules (AIM) theory. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Single‐electron halogen bond: Ab initio study

Ab initio calculations have been performed on single‐electron halogen bonds between methyl radical and bromine‐containing molecules to gain a deeper insight into the nature of such noncovalent interactions. Bader's atoms in molecules (AIM) theory have also been applied to the analysis of the linking of the single‐electron halogen bond. Various characteristics of the R? Br…CH₃ interaction, i.e., binding energies, geometrical parameters and topological properties of the electron density have been determined. The presence of the bond critical points (BCPs) between the bromine atom and methyl radical and the values of electron density and Laplacian of electron density at these BCPs indicate the closed‐shell interactions in the complexes. The single‐electron halogen bonds, which are significantly weaker than the normal halogen bonds, exhibit equally bond strength as compared to the single‐electron hydrogen bond. It has been also found that plotting of the binding energies versus topological properties of the electron density at the BCPs gives two straight lines. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

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.     

四氯化苯醌-二苯撑体系的光诱导电子转移研究

用从头算方法对四氯化苯醌-二苯撑体系分子间相互作用进行了理论计算研究.用MP2/6-31G^**方法,分别优化电子给体二苯撑,受体四氯化苯醌的稳定构型,用同样的方法优化配合物的层间距得到其最稳定构型,并计算了BSSE校正后的电子给受体配合物的稳定化能.用CIS/6-31++G^**方法,计算了给体、受体及配合物的电子激发态.理论计算验证了给体和受体间能形成稳定的电子给受体配合物,该配合物受光激发能直接产生电荷转移态.在球孔穴近似和点偶极近似下,对电荷转移吸收的理论计算结果进行了非平衡溶剂化能校正.经非平衡溶剂化能校正的电荷转移跃迁能与实验值符合较好.     

Ab initio calculations on halogen-bonded complexes and comparison with density functional methods

A systematic theoretical investigation on a series of dimeric complexes formed between some halocarbon molecules and electron donors has been carried out by employing both ab initio and density functional methods. Full geometry optimizations are performed at the Moller-Plesset second-order perturbation (MP2) level of theory with the Dunning's correlation-consistent basis set, aug-cc-pVDZ. Binding energies are extrapolated to the complete basis set (CBS) limit by means of two most commonly used extrapolation methods and the aug-cc-pVXZ (X = D, T, Q) basis sets series. The coupled cluster with single, double, and noniterative triple excitations [CCSD(T)] correction term, determined as a difference between CCSD(T) and MP2 binding energies, is estimated with the aug-cc-pVDZ basis set. In general, the inclusion of higher-order electron correlation effects leads to a repulsive correction with respect to those predicted at the MP2 level. The calculations described herein have shown that the CCSD(T) CBS limits yield binding energies with a range of -0.89 to -4.38 kcal/mol for the halogen-bonded complexes under study. The performance of several density functional theory (DFT) methods has been evaluated comparing the results with those obtained from MP2 and CCSD(T). It is shown that PBEKCIS, B97-1, and MPWLYP functionals provide accuracies close to the computationally very expensive ab initio methods.     

Further insights in the ability of classical nonadditive potentials to model actinide ion–water interactions

Pursuing our efforts on the development of accurate classical models to simulate radionuclides in complex environments (Réal et al., J. Phys. Chem. A 2010, 114, 15913; Trumm et al. J. Chem. Phys. 2012, 136, 044509), this article places a large emphasis on the discussion of the influence of models/parameters uncertainties on the computed structural, dynamical, and temporal properties. Two actinide test cases, trivalent curium and tetravalent thorium, have been studied with three different potential energy functions, which allow us to account for the polarization and charge‐transfer effects occurring in hydrated actinide ion systems. The first type of models considers only an additive energy term for modeling ion/water charge‐transfer effects, whereas the other two treat cooperative charge‐transfer interactions with two different analytical expressions. Model parameters are assigned to reproduce high‐level ab initio data concerning only hydrated ion species in gas phase. For the two types of cooperative charge‐transfer models, we define two sets of parameters allowing or not to cancel out possible errors inherent to the force field used to model water/water interactions at the ion vicinity. We define thus five different models to characterize the solvation of each ion. For both ions, our cooperative charge‐transfer models lead to close results in terms of structure in solution: the coordination number is included within 8 and 9, and the mean ion/water oxygen distances are 2.45 and 2.49 Å, respectively, for Th(IV) and Cm(III). © 2012 Wiley Periodicals, Inc.     

On the origin of red and blue shifts of X-H and C-H stretching vibrations in formic acid (formate ion) and proton donor complexes.

Complexes between formic acid or formate anion and various proton donors (HF, H(2)O, NH(3), and CH(4)) are studied by the MP2 and B3LYP methods with the 6-311++G(3df,3pd) basis set. Formation of a complex is characterized by electron-density transfer from electron donor to ligands. This transfer is much larger with the formate anion, for which it exceeds 0.1 e. Electron-density transfer from electron lone pairs of the electron donor is directed into sigma* antibonding orbitals of X--H bonds of the electron acceptor and leads to elongation of the bond and a red shift of the X--H stretching frequency (standard H-bonding). However, pronounced electron-density transfer from electron lone pairs of the electron donor also leads to reorganization of the electron density in the electron donor, which results in changes in geometry and vibrational frequency. These changes are largest for the C--H bonds of formic acid and formate anion, which do not participate in H-bonding. The resulting blue shift of this stretching frequency is substantial and amounts to almost 35 and 170 cm(-1), respectively.     

Disentanglement of orthogonal hydrogen and halogen bonds via natural orbital for chemical valence: A charge displacement analysis

As known, the electron density of covalently bound halogen atoms is anisotropically distributed, making them potentially able to establish many weak interactions, acting at the same time as halogen bond donors and hydrogen bond acceptors. Indeed, there are many examples in which the halogen and hydrogen bond coexist in the same structure and, if a correct bond analysis is required, their separation is mandatory. Here, the advantages and limitations of coupling the charge displacement analysis with natural orbital for chemical valence method (NOCV-CD) to separately analyze orthogonal weak interactions are shown, for both symmetric and asymmetric adducts. The methodology gives optimal results with intermolecular adducts but, in the presence of an organometallic complex, also intramolecular interactions can be correctly analyzed. Beyond the methodological aspects, it is shown that correctly separate and quantify the interactions can give interesting chemical insights about the systems.     

Competition between hydrogen and halogen bonding in halogenated 1‐methyluracil: Water systems

The competition between hydrogen‐ and halogen‐bonding interactions in complexes of 5‐halogenated 1‐methyluracil (XmU; X = F, Cl, Br, I, or At) with one or two water molecules in the binding region between C5‐X and C4?O4 is investigated with M06‐2X/6‐31+G(d). In the singly‐hydrated systems, the water molecule forms a hydrogen bond with C4?O4 for all halogens, whereas structures with a halogen bond between the water oxygen and C5‐X exist only for X = Br, I, and At. Structures with two waters forming a bridge between C4?O and C5‐X (through hydrogen‐ and halogen‐bonding interactions) exist for all halogens except F. The absence of a halogen‐bonded structure in singly‐hydrated ClmU is therefore attributed to the competing hydrogen‐bonding interaction with C4?O4. The halogen‐bond angle in the doubly‐hydrated structures (150–160°) is far from the expected linearity of halogen bonds, indicating that significantly non‐linear halogen bonds may exist in complex environments with competing interactions. © 2016 Wiley Periodicals, Inc.     

Do intramolecular halogen bonds exist? Ab initio calculations and crystal structures’ evidences

Cl···O intramolecular contacts are analyzed here. Hence, the Cambridge Structural Database search was performed to find five-membered Cl–C=C–C=O pseudo-rings closed through Cl···O intramolecular contacts. Such rings were analyzed and it was found that some of Cl···O distances are smaller than the corresponding sum of van der Waals radii. Additionally ab initio MP2/6-311++G(d,p) calculations and AIM analyses were carried out here on 3-chloropropenal and its fluoro derivatives to deepen the nature of intramolecular Cl···O interactions, some of findings may indicate that they are attractive and stabilizing ones. This article is dedicated to Prof. T. M. Krygowski on the occasion of his 70th birthday.     

Do intramolecular halogen bonds exist? Ab initio calculations and crystal structures’ evidences

Cl···O intramolecular contacts are analysed here. Hence the Cambridge Structural Database (CSD) search was performed to find five-membered Cl–C=C–C=O pseudo-rings closed through Cl···O intramolecular contacts. Such rings were analysed and it was found that some of Cl···O distances are smaller than the corresponding sum of van der Waals radii. Additionally ab initio MP2/6-311++G(d,p) calculations and AIM analyses were carried out here on 3-chloropropenal and its fluoro derivatives to deepen the nature of intramolecular Cl···O interactions, some of findings may indicate that they are attractive and stabilizing ones. This paper is dedicated to Prof. T. M. Krygowski on the occasion of his 70th birthday