Journal of Chemistry
Volume 2017 (2017), Article ID 9031494, 8 pages
https://doi.org/10.1155/2017/9031494

Research Article

Description of Weak Halogen Bonding Using Various Levels of Symmetry-Adapted Perturbation Theory Combined with Effective Core Potentials

Piotr Matczak

Department of Theoretical and Structural Chemistry, Faculty of Chemistry, University of Łódź, Pomorska 163/165, 90-236 Lodz, Poland

Correspondence should be addressed to Piotr Matczak; p.a.matczak@gmail.com

Received 5 September 2016; Revised 8 December 2016; Accepted 10 January 2017; Published 9 February 2017

Academic Editor: Franck Rabilloud

Copyright © 2017 Piotr Matczak. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The present work starts with providing a description of the halogen bonding (XB) interaction between the halogen atom of MH₃X (where M = C–Pb and X = I, At) and the N atom of HCN. This interaction leads to the formation of stable yet very weakly bound MH₃X⋯NCH complexes for which the interaction energy () between MH₃X and HCN is calculated using various symmetry-adapted perturbation theory (SAPT) methods combined with the def2-QZVPP basis set and midbond functions. This basis set assigns effective core potentials (ECPs) not only to the I or At atom directly participating in the XB interaction with HCN but also to the M atom when substituted with Sn or Pb. Twelve SAPT methods (or levels) are taken into consideration. According to the SAPT analysis of , the XB interaction in the complexes shows mixed electrostatic-dispersion nature. Next, the accuracy of SAPT is evaluated by comparing with CCSD(T) reference data. This comparison reveals that high-order SAPT2+ method and the much less computationally demanding SAPT(DFT) method perform very well in describing of the complexes. However, the accuracy of these methods decreases dramatically if they are combined with the so-called Hartree-Fock correction.

1. Introduction

Weak intermolecular interactions play an important role in establishing the structure and stability of a broad range of chemical systems, from simple molecular complexes to macromolecular assemblies [1–3]. Computational methods based on electronic structure calculations are useful in determining the magnitude of weak intermolecular interactions as well as in understanding their nature [4–6]. However, the ability to describe weak intermolecular interactions accurately and efficiently differs significantly between various families of computational electronic structure methods. Among these families, the symmetry-adapted perturbation theory (SAPT) [7–11] is capable of providing weak interaction energies accurate even up to 2–4% for benchmark systems of rare-gas dimers [9] but the computational cost of the corresponding calculations in principle scales like with the system size. This makes such calculations very time-consuming, especially if they are carried out for large chemical systems. The high computational cost of SAPT calculations is associated, to a great extent, with the inclusion of intramolecular electron correlation effects. It is particularly true for the traditional formulation of SAPT which uses many-body electron correlation wave functions for molecules forming a complex [7]. An alternative approach to handling the intramolecular electron correlation within SAPT is to utilize the density functional theory (DFT) for the molecular wave functions [12–17]. The resulting SAPT(DFT) formulation scales better than the traditional formulation, and, therefore, the former can be applied to systems with hundreds of atoms [18]. In order to adjust the computational demands of SAPT calculations to the size of a given system, various methods (or levels) have been introduced for the traditional SAPT formulation [10, 19]. In general, these levels differ in their scaling behavior through the omission of some energy correction terms from the SAPT expansion of intermolecular interaction energy.

Weak intermolecular interactions sometimes occur between molecular sites, of which one contains an atom of element with high nuclear charge. Halogen bonding (XB) [20] involving iodine constitutes a good example of such interactions (historically speaking, the XB formation in H₃N⋯I₂ was probably the very first observation of XB interaction [21]). From the computational point of view, describing weak intermolecular interactions involving heavier elements is a real challenge for electronic structure methods for two main reasons. Firstly, the large number of electrons in atoms of heavier elements makes electronic structure calculations very costly. Secondly, the incorporation of relativistic effects in the calculations is essential to obtain accurate energies of weak interactions involving heavier elements. These difficulties can be largely overcome by the use of effective core potentials (ECPs) for atoms of heavier elements (usually with the nuclear charge ) [22]. Within the ECP approximation the explicit treatment of core electrons of heavier atoms is replaced by introducing potential operators fitted to the results of relativistic all-electron calculations for such atoms. Therefore, the ECP approximation is able to implicitly account for leading relativistic effects at low computational cost. The use of ECPs has also been implemented for SAPT calculations of interaction energies [23].

In this work a systematic examination of the accuracy of a combination of SAPT and ECP for the description of weak XB is provided relative to the level of SAPT. The accuracy of a variety of SAPT methods is established across a series of 10 heterodimers containing MH₃X (where M = C–Pb and X = I, At) complexed with HCN. These complexes exemplify the X⋯N type of XB in which the halogen atom possesses a very high nuclear charge (53 and 85 for iodine and astatine, resp.). Thus, the ECP treatment of I- and At-atom core electrons is employed here for the SAPT description of the XB interaction in the MH₃X⋯NCH complexes.

2. Methods

The geometries of the MH₃X⋯NCH complexes have been optimized at the MP2/def2-QZVPP level of theory [24–27]. These geometries are characterized by the absence of imaginary vibrational frequencies, and, therefore, they correspond to minima on the potential energy surfaces of the complexes. These geometries will also be kept fixed in all subsequent calculations of the interaction energy () between the MH₃X and HCN parts of the complexes. The MP2 method has been combined with the resolution-of-the-identity approximation [28]. The def2-QZVPP basis set makes use of Stuttgart relativistic ECPs [26] to describe the core electrons of Sn, Pb, I, and At. These calculations have been carried out with the TURBOMOLE 6.6 program [29].

The values of between MH₃X and HCN in the complexes have been calculated using such SAPT methods as SAPT0, SAPT2, SAPT2+, SAPT2+, SAPT2+3, and SAPT(DFT) [10, 19]. The first five methods are based on the traditional (or regular) wave function-based formulation of SAPT. They truncate the SAPT expansion of either at the second or at the third order in the perturbation operator of intermolecular interaction, and they differ in the degree of the completeness of included intramolecular electron correlation effects. Details of the truncation of SAPT expansion for individual SAPT methods are shown in the following definitions:

The explanation of SAPT energy correction terms on the right sides of (1) can be found in a number of reviews covering the fundamentals of SAPT [7–11]. Suffice it to say here that SAPT0, SAPT2, and SAPT2+ treat the intermolecular interaction perturbation up to second order, while SAPT2+ and SAPT2+3 additionally include third-order terms in the perturbation operator of intermolecular interaction. Regarding the intramolecular electron correlation, it is completely neglected by SAPT0 and included up to second order for electrostatic, exchange, and induction contributions to in SAPT2 and additionally for the dispersion contribution by SAPT2+, SAPT2+, and SAPT2+3. The PBE0 density functional [30] with the Fermi-Amaldi asymptotic correction [31] and the Tozer-Handy splicing scheme [32] have been used in SAPT(DFT). To compute the asymptotic correction, either available experimental or calculated ionization potentials have been employed (see Supplementary Section  1 in Supplementary Material available online at https://doi.org/10.1155/2017/9031494). Both the wave function-based SAPT methods and SAPT(DFT) can be supplemented with the so-called Hartree-Fock (HF) correction ( or , depending on the SAPT method to which the correction is added). For SAPT0, SAPT2, SAPT2+, SAPT2+, and SAPT(DFT) such a correction is determined as

whereas a different form is employed for the SAPT2+3 method

The term in (2) denotes the supermolecular Hartree-Fock interaction energy in the complex for which the SAPT interaction energy is calculated. The HF correction is a way to capture some higher-order terms not explicitly evaluated by SAPT. All SAPT results have been produced using the dimer-centered plus basis set (DC⁺BS) scheme, with the def2-QZVPP basis set extended with a set of midbond (mb) functions [33]. The resulting basis set will be denoted by def2-QZVPP+mb further in the text. The set of mb functions is of the 3s3p2d2f form, with the exponents 0.9, 0.3, and 0.1 for and p and 0.6 and 0.2 for and . These functions are centered midway between the X and N atoms of the complexes. All SAPT calculations have been done using the SAPT 2012.2 [34] program interfaced with DALTON 2.0 [35]. The SAPT(DFT) calculations have additionally been carried out with the SAPT(DFT) implementation available in the MOLPRO 2012.1 program [36], but these calculations do not end with successful reproduction of the results obtained with SAPT 2012.2.

Apart from the SAPT description of the XB interaction in MH₃X⋯NCH, additional supermolecular estimations of are provided in this work for comparison purposes. A number of popular wave function theory methods, mainly those belonging to the family of Møller-Plesset methods [24], will be taken into account. MP2 [37], SCS-MP2 [38], SCS(MI)-MP2 [39], MP2C [40], MP2.5 [41], MP3 [37], MP4(SDQ) [37], MP4 [37], and CCSD [42] are selected. The resolution-of-the-identity approximation is used wherever possible. The aforementioned methods are combined with the def2-QZVPP+mb basis set in order to calculate the magnitude of the XB interaction in the MH₃X⋯NCH complexes. The resulting values are corrected for the basis set superposition error by utilizing a counterpoise correction [43]. All the supermolecular calculations of have been performed using TURBOMOLE 6.6.

The accuracy of SAPT methods is evaluated by comparing their values calculated for MH₃X⋯NCH with the corresponding reference data. The reference values of have been obtained from the CCSD(T) method [42] combined with the def2-QZVPP+mb basis set. The CCSD(T)/def2-QZVPP+mb values of are corrected for the basis set superposition error [43]. Additionally, the CCSD(T) values extrapolated to the complete basis set (CBS) limit [44–47] will also serve as reference values. The CCSD(T)/CBS level of theory is commonly perceived as the “gold standard” for predicting reliable energies of weak intermolecular interactions [48]. This level of theory was previously used to provide benchmark values in a set of 40 complexes that were stabilized by a variety of weak intermolecular interactions in which halogens participated [49]. Details of the CCSD(T)/CBS calculations performed for MH₃X⋯NCH are given in Supplementary Section  3.

3. Results and Discussion

3.1. Fundamental Properties of the Complexes

The geometry optimization of the complexes with the assumed X⋯N XB interaction leads to structures in which the HCN molecule is oriented precisely along the M–X bond axis (see Figure 1). The values of the distance between the X and N atoms () for the complexes in their optimized geometries are listed in Table 1. These values are rather large ( Å) but in the case of the complexes with X = I they are still shorter than the sum of the I- and N-atom van der Waals radii. This sum is equal to 3.70 Å if the van der Waals radii proposed by Alvarez [50] are taken. The fulfillment of the geometrical criterion formulated as  Å suggests the existence of XB in the MH₃I⋯NCH complexes. Such a geometrical criterion cannot be applied to the complexes with X = At because the van der Waals radius of astatine is unknown [50].

Table 1: Distance between the X and N atoms (; in Å), maximal value of electrostatic potential in the σ-hole region on the X atom (; in kcal/mol), CCSD(T)/CBS and SAPT(DFT) interaction energies between MH₃X and HCN (; in kcal/mol), and SAPT(DFT) principal interaction energy components (, , , and ; in kcal/mol) for the investigated MH₃X⋯NCH complexes. The percentage of each attractive SAPT(DFT) component relative to the total attraction is given in parentheses.

Complex

CCSD(T)/CBS
SAPT(DFT)

CH₃I⋯NCH
3.261
12.2
−1.46
−1.23
−2.23

−0.52

−2.45

3.97

SiH₃I⋯NCH
3.451
8.7
−1.14
−1.08
−1.41

−0.35

−1.83

2.51

GeH₃I⋯NCH
3.426
7.3
−1.01
−0.94
−1.34

−0.38

−1.95

2.73

SnH₃I⋯NCH
3.470
4.4
−0.79
−0.74
−1.02

−0.36

−1.86

2.50

PbH₃I⋯NCH
3.481
1.6
−0.55
−0.49
−0.73

−0.37

−1.87

2.47

CH₃At⋯NCH
3.171
18.4
−2.36
−2.02
−3.95

−0.84

−3.19

5.96

SiH₃At⋯NCH
3.357
13.6
−1.80
−1.72
−2.58

−0.55

−2.41

3.83

GeH₃At⋯NCH
3.322
12.9
−1.74
−1.59
−2.68

−0.60

−2.59

4.28

SnH₃At⋯NCH
3.353
10.8
−1.52
−1.38
−2.37

−0.57

−2.51

4.08

PbH₃At⋯NCH
3.343
8.1
−1.33
−1.15
−2.25

−0.59

−2.61

4.30

Figure 1: Structure of the investigated MH₃X⋯NCH complexes. The optimized geometry of CH₃I⋯NCH is shown as an example. The XB interaction is marked by a dotted line.

The structure of the MH₃X⋯NCH complexes results from the highly directional ability of the X atom, covalently bonded to the M atom, to attract the N-atom lone electron pair of HCN, acting as a nucleophilic species in this case. A widely accepted mechanism of XB formation postulates that XB originates from the interaction between two centers possessing opposite electrostatic potentials [51]. For the MH₃X⋯NCH complexes the center exhibiting a positive electrostatic potential is located on the X atom in the region on the extension of the M–X bond (i.e., the -hole), whereas the negative electrostatic potential characterizes the region around the lone electron pair of the N atom. The maximal value of electrostatic potential on the molecular surface delineated by the electron density contour of 0.001 au () in the -hole region on the X atom is a convenient tool for quantifying the -hole [52]. The calculated values of are presented in Table 1. An increase in values is observed if the X atom changes from I to At, and thus it becomes heavier and more polarizable [52]. Moreover, values increase with the growing electron-attracting character of the MH₃ group [53]. The electron-attracting character can be represented by the so-called group electronegativity [54] which is calculated for the MH₃ groups (for details, see Supplementary Section 4). The electronegativity of the MH₃ groups increases gradually while going up group 14 from M = Pb to M = C, and, therefore, the most positive values are achieved for CH₃X in either of the sequences of the MH₃X⋯NCH complexes with the X atom kept fixed.

The CCSD(T)/CBS values shown in Table 1 are negative for all MH₃X⋯NCH complexes. This implies their energetic stabilization and this also constitutes an energetic criterion for XB occurrence in the complexes. The values of fall in a narrow range from −0.55 to −2.36 kcal/mol. The XB interaction in MH₃X⋯NCH can be classified as weak in view of the fact that the strongest XB interactions reported in literature exhibit interaction energies exceeding −40 kcal/mol [55]. The CH₃I molecule forms weaker XB with HCN than with NH₃ [56] or CH₂O [57]. The weakness of the XB interaction in the MH₃X⋯NCH complexes, as indicated by their values, is connected with the large distances. If pairs of the complexes with the M atom fixed are considered, MH₃At always forms stronger XB than MH₃I. This is reflected in the geometrical relation for each such pair. The XB interaction becomes weaker and weaker when the M atom traverses down group 14 for two sequences of the complexes with the X atom kept fixed and varying M. The clear trend in is only roughly followed by a trend in for these sequences. These trends do not correlate with each other ideally but the monotonous weakening of the XB interaction is often accompanied by the gradual elongation of . It is noteworthy that correlates well with for all MH₃X⋯NCH complexes: the greater is developed in the -hole region, the stronger the XB occurs.

In order to establish the nature of the XB interaction in MH₃X⋯NCH, the energy correction terms from the SAPT(DFT) expansion of have been grouped into four principal components corresponding to electrostatics, induction, dispersion, and exchange. Details of the grouping scheme can be found in Supplementary Section5. The values of the electrostatic (), induction (), dispersion (), and exchange () components of obtained from SAPT(DFT) are listed in Table 1. The most important aspect of these results is the fact that and play the central role in the attraction between MH₃X and HCN. These two components are responsible for at least 88% of the total attraction occurring in the complexes. Thus, the XB interaction in MH₃X⋯NCH shows mixed electrostatic-dispersion nature. This finding is in agreement with the analogical conclusion drawn in an earlier work taking various types of XB into account [58]. It is interesting to examine trends in the attractive components for the sequences of the complexes with the X atom kept fixed. For such sequences the percentage of electrostatics decreases (and simultaneously the percentage of dispersion increases) while descending group 14 from M = C to M = Pb. This leads to the conclusion that the role of dispersion in the stabilization of the complexes is most significant for the weakest XB interaction. component does not change monotonously if M descends group 14. Moreover, the percentage share of in the total attraction between MH₃X and HCN remains practically constant for all 10 complexes.

3.2. Assessment of Various SAPT Methods

In order to evaluate the performance of various SAPT methods in describing the XB interaction in the MH₃X⋯NCH complexes, values predicted by these methods are compared with CCSD(T) reference results in a statistical manner. This statistical comparison is based on three metrics of errors in the SAPT values. The root mean-squared error (RMSE), the percentage relative root mean-squared error (rRMSE  (%)), and the mean signed error (MSE) are calculated using the following formulas:

where is the number of the complexes studied (M = 10). RMSE is used to characterize the overall performance of each SAPT method in predicting . Because values are small themselves and the same applies to the RMSE values, then rRMSE  (%) should more conveniently and reliably illustrate the accuracy of the SAPT values relative to the corresponding CCSD(T) results. Finally, MSE provides information about systematic errors occurring in values obtained from each SAPT method. If MSE is negative (positive) for a given SAPT method, then this method tends to overestimate (underestimate) the strength of the XB interaction in the MH₃X⋯NCH complexes.

The evaluation of SAPT accuracy has been performed for 12 SAPT methods. Table 2 lists these methods, together with the three mean errors determined for their values. Based on the RMSE values one can think that SAPT generally performs very well because RMSE is less than 1 kcal/mol for all methods. This is not necessarily true if rRMSE  (%) is analyzed as well. rRMSE  (%) is high for several SAPT methods, which indicates that the RMSE values are in fact quite large when compared to both the CCSD(T)/def2-QZVPP+mb and CCSD(T)/CBS values. It is justified to rank SAPT methods’ performances relative to the rRMSE  (%) values. From the magnitude of rRMSE  (%) one can conclude that the SAPT2+ method provides the best approximation of the CCSD(T) reference results. As with SAPT2+, the SAPT(DFT) method reproduces perfectly the CCSD(T)/def2-QZVPP+mb energies. A slightly larger rRMSE  (%) value is observed for SAPT2+3 but this value still does not exceed 10%. The values of MSE for the three methods indicate that SAPT2+ tends to slightly overestimate the strength of the XB interaction in MH₃X⋯NCH, while SAPT(DFT) and SAPT2+3 show the opposite tendency (the mean underestimation of reaches a maximum of 0.14 kcal/mol for SAPT(DFT)). Unsurprisingly, two regular high-order SAPT methods SAPT2+ and SAPT2+3 yield with very good accuracy. It is rather puzzling, however, that SAPT2+3 performs slightly worse than SAPT2+ while the reverse ordering of these two methods could actually be expected. This means that the inclusion of all third-order interaction energy correction terms does not lead to any improvement in for the MH₃X⋯NCH complexes. The SAPT2+ accuracy superior to that of SAPT2+3 is mostly due to the fact that the third-order induction and exchange-induction terms are poorly recovered by SAPT combined with ECPs [23] (see also Supplementary Section 2).

Table 2: Errors (RMSE, rRMSE (%), and MSE; in kcal/mol or %) in obtained from various SAPT methods. The errors without parentheses have been calculated relative to the CCSD(T)/def2-QZVPP+mb results, while those put in parentheses are based on the CCSD(T)/CBS reference values of .

Method
RMSE
rRMSE (%)
MSE

SAPT0
0.28
(0.23)
21
(16)
−0.26
(−0.21)

SAPT0 +
0.67
(0.62)
51
(44)
−0.64
(−0.59)

SAPT2
0.39
(0.34)
36
(29)
−0.38
(−0.33)

SAPT2 +
0.78
(0.73)
66
(58)
−0.76
(−0.71)

SAPT2+
0.30
(0.25)
27
(21)
−0.30
(−0.25)

SAPT2+ +  
0.70
(0.65)
58
(50)
−0.67
(−0.62)

SAPT2+
0.06
(0.04)
7
(3)
−0.05
(−0.01)

SAPT2+ +
0.44
(0.39)
37
(30)
−0.43
(−0.38)

SAPT2+3
0.13
(0.16)
9
(9)
0.03
(0.08)

SAPT2+3 +
0.51
(0.46)
41
(34)
−0.49
(−0.44)

SAPT(DFT)
0.12
(0.16)
7
(10)
0.08
(0.14)

SAPT(DFT) +
0.30
(0.25)
26
(20)
−0.29
(−0.24)

From the results in Table 2 the effect of the HF correction on can be established. This correction collects mainly higher-order induction and exchange-induction terms but it may also contain some spurious, unphysical terms. It is evident that the HF correction adversely affects the accuracies of all SAPT methods in describing the weak XB interaction in MH₃X⋯NCH. Six SAPT methods utilizing this correction (the so-called hybrid approach) show many times larger RMSE values than the same methods without or . As illustrated by the negative values of MSE, the majority of SAPT methods tend to predict an overbinding between MH₃X and HCN. Adding the HF correction results in further overestimation of the XB strength in the complexes. A similar overestimation of I⋯O XB strength was also detected in a previous SAPT study of the CH₃I⋯OCH₂ complex [57].

The deterioration of SAPT accuracy due to the incorporation of or into is likely caused by a spurious, unphysical part of the HF correction. It is known that spurious, unphysical effects contribute considerably to the HF correction if a molecular complex is bound primarily by dispersion [59]. This seems to be the case for the MH₃X⋯NCH complexes for which the dispersion component of their is always significant or even dominant (the percentage of ; see Table 1).

The simplest SAPT method, that is, SAPT0, outperforms both SAPT2 and SAPT2+, although the latter methods incorporate more energy correction terms, and, therefore, they are computationally more demanding. This superior accuracy of SAPT0 seems to be due to a fairly good compensation of energy correction terms appearing in the SAPT0 expansion of . It has been demonstrated in another recent work [19] that SAPT0 yields in many cases not much worse interaction energies compared to more advanced SAPT methods which account for intramolecular electron correlation effects. The combination of SAPT0, SAPT2, or SAPT2+ with fails badly in describing the XB interaction in MH₃X⋯NCH; rRMSE  (%) of the resulting hybrid levels exceeds 50% relative to the CCSD(T)/def2-QZVPP+mb reference results.

It is instructive to compare the performance of the SAPT methods with the performance of a series of popular wave function theory methods that compute within the supermolecular approach. Statistical metrics measuring the accuracy of supermolecular values obtained from 10 wave function theory methods against the CCSD(T) reference data are presented in Table 3. Of the 10 methods, SCS(MI)-MP2, MP2C, and MP2.5 perform best. They afford very small RMSE and rRMSE  (%) and these errors are even slightly smaller than those made by SAPT2+ and SAPT2+3. What is particularly important is that the SCS(MI)-MP2 calculations entail much lower computational cost (MP2 formally scales as with the system size). The good performance of SCS(MI)-MP2 has actually been expected because this method has been specially designed for weak intermolecular interactions. Similarly, the MP2C and MP2.5 methods were also developed for providing accurate interaction energies in various molecular complexes. For the MH₃X⋯NCH complexes, these modifications of the original MP2 method do show a great improvement in the reproduction of over their ancestor. The performances of MP4 and CCSD are rather disappointing and these advanced methods yield less accurate values of than the regular high-order SAPT and SAPT(DFT) methods. According to the HF method, the interaction between MH₃X and HCN destabilizes the MH₃X⋯NCH complexes, which disqualifies this method as a computational tool for describing weak XB interactions reliably.

Table 3: Errors (RMSE, rRMSE  (%), and MSE; in kcal/mol or %) in obtained from various wave function theory methods within the supermolecular approach. The errors without parentheses have been calculated relative to the CCSD(T)/def2-QZVPP+mb results, while those put in parentheses are based on the CCSD(T)/CBS reference values of .

Method
RMSE
rRMSE  (%)
MSE

HF
2.19
(2.24)
200
(193)
2.16
(2.22)

MP2
0.46
(0.40)
41
(34)
−0.45
(−0.39)

SCS-MP2
0.11
(0.16)
10
(14)
0.10
(0.15)

SCS(MI)-MP2
0.08
(0.05)
7
(3)
−0.06
(−0.01)

MP2C
0.06
(0.08)
3
(5)
0.01
(0.06)

MP2.5
0.06
(0.11)
4
(8)
0.04
(0.10)

MP3
0.55
(0.60)
48
(50)
0.54
(0.59)

MP4(SDQ)
0.23
(0.28)
23
(26)
0.22
(0.27)

MP4
0.34
(0.29)
31
(25)
−0.33
(−0.28)

CCSD
0.40
(0.45)
36
(38)
0.39
(0.44)

4. Conclusions

A description of the X⋯N XB interaction in a series of 10 MH₃X⋯NCH complexes (where M = C–Pb and X = I, At) has been obtained from various SAPT methods combined with the def2-QZVPP basis set and midbond functions. This basis set assigns ECPs not only to the I or At atom directly participating in the XB interaction with HCN but also to the M atom when substituted with Sn or Pb. The accuracy of SAPT has been evaluated by comparing with CCSD(T) reference data. The main findings made in this work can be summarized as follows:

(1)The CCSD(T)/CBS calculations of reveal that the complexes are stabilized by a weak X⋯N XB interaction whose values fall in the range between −0.55 and −2.36 kcal/mol. Astatine forms stronger XB interaction with HCN than iodine, by no more than 1 kcal/mol in each pair of the complexes with a given M. A relationship between and the kind of group 14 element substituting the M atom has been identified for two sequences of the complexes with the X atom kept fixed. The strength of the XB interaction decreases gradually when the M atom varies down group 14 from C to Pb. This is rationalized in terms of and MH₃-group electronegativity.(2)Grouping the terms from the SAPT(DFT) expansion of into four principal components indicates the mixed electrostatic-dispersion nature of the XB interaction in MH₃X⋯NCH. The XB weakening caused by the substitution of the M and X atoms clearly affects the percentage shares of and . This weakening is associated with the decrease of percentage share and the simultaneous increase of percentage share.(3)Regular high-order SAPT methods such as SAPT2+ and SAPT2+3 provide XB description very close to that obtained at the CCSD(T)/def2-QZVPP+mb and CCSD(T)/CBS levels of theory. The accuracy of obtained from the computationally favorable SAPT(DFT) method is comparable to that found for SAPT2+ interaction energies. The HF correction present in the hybrid approach leads to a significant deterioration in the accuracy of the SAPT description of the XB interaction in the MH₃X⋯NCH complexes.

Competing Interests

The author declares no competing interests.

Acknowledgments

This work was partially supported by PL-Grid Infrastructure.

References

G. A. Jeffrey and W. Saenger, Hydrogen Bonding in Biological Structures, Springer, Berlin, Germany, 1991.

I. Haiduc and F. T. Edelmann, Supramolecular Organometallic Chemistry, Wiley-VCH, Weinheim, Germany, 1999.

A. Karshikoff, Non-Covalent Interactions in Proteins, Imperial College Press, Singapore, 2006.

S. Tsuzuki, “Interactions with aromatic rings,” Structure and Bonding, vol. 115, pp. 149–193, 2005.

I. G. Kaplan, Intermolecular Interactions: Physical Picture, Computational Methods and Model Potentials, John Wiley & Sons, Chichester, UK, 2006.

P. Hobza and K. Müller-Dethlefs, Non-Covalent Interactions. Theory and Experiment, RSC Publishing, Cambridge, UK, 2010.

B. Jeziorski, R. Moszynski, and K. Szalewicz, “Perturbation theory approach to intermolecular potential energy surfaces of van der Waals complexes,” Chemical Reviews, vol. 94, no. 7, pp. 1887–1930, 1994.

K. Szalewicz, K. Patkowski, and B. Jeziorski, “Intermolecular interactions via perturbation theory: from diatoms to biomolecules,” Structure and Bonding, vol. 116, pp. 43–117, 2005.

K. Szalewicz, “Symmetry-adapted perturbation theory of intermolecular forces,” Wiley Interdisciplinary Reviews: Computational Molecular Science, vol. 2, no. 2, pp. 254–272, 2012.

E. G. Hohenstein and C. D. Sherrill, “Wavefunction methods for noncovalent interactions,” Wiley Interdisciplinary Reviews: Computational Molecular Science, vol. 2, no. 2, pp. 304–326, 2012.

G. Jansen, “Symmetry-adapted perturbation theory based on density functional theory for noncovalent interactions,” Wiley Interdisciplinary Reviews: Computational Molecular Science, vol. 4, no. 2, pp. 127–144, 2014.

H. L. Williams and C. F. Chabalowski, “Using Kohn-Sham orbitals in symmetry-adapted perturbation theory to investigate intermolecular interactions,” The Journal of Physical Chemistry A, vol. 105, no. 3, pp. 646–659, 2001.

A. J. Misquitta and K. Szalewicz, “Intermolecular forces from asymptotically corrected density functional description of monomers,” Chemical Physics Letters, vol. 357, no. 3-4, pp. 301–306, 2002.

A. Heßelmann and G. Jansen, “First-order intermolecular interaction energies from Kohn-Sham orbitals,” Chemical Physics Letters, vol. 357, no. 5-6, pp. 464–470, 2002.

A. Heßelmann and G. Jansen, “Intermolecular induction and exchange-induction energies from coupled-perturbed Kohn-Sham density functional theory,” Chemical Physics Letters, vol. 362, no. 3-4, pp. 319–325, 2002.

A. J. Misquitta, B. Jeziorski, and K. Szalewicz, “Dispersion energy from density-functional theory description of monomers,” Physical Review Letters, vol. 91, no. 3, Article ID 033201, 4 pages, 2003.

A. Heßelmann and G. Jansen, “Intermolecular dispersion energies from time-dependent density functional theory,” Chemical Physics Letters, vol. 367, no. 5-6, pp. 778–784, 2003.

R. Podeszwa, W. Cencek, and K. Szalewicz, “Efficient calculations of dispersion energies for nanoscale systems from coupled density response functions,” Journal of Chemical Theory and Computation, vol. 8, no. 6, pp. 1963–1969, 2012.

T. M. Parker, L. A. Burns, R. M. Parrish, A. G. Ryno, and C. D. Sherrill, “Levels of symmetry adapted perturbation theory (SAPT). I. Efficiency and performance for interaction energies,” The Journal of Chemical Physics, vol. 140, no. 9, Article ID 094106, 16 pages, 2014.

P. Metrangolo and G. Resnati, “Halogen bonding: a paradigm in supramolecular chemistry,” Chemistry—A European Journal, vol. 7, no. 12, pp. 2511–2519, 2001.

F. Guthrie, “On the iodide of iodammonium,” Journal of the Chemical Society, vol. 16, pp. 239–244, 1863.

P. Schwerdtfeger, Relativistic Electronic Structure Theory. Part 1. Fundamentals, Elsevier, Amsterdam, The Netherlands, 2002.

K. Patkowski and K. Szalewicz, “Frozen core and effective core potentials in symmetry-adapted perturbation theory,” The Journal of Chemical Physics, vol. 127, no. 16, Article ID 164103, 17 pages, 2007.

C. Møller and M. S. Plesset, “Note on an approximation treatment for many-electron systems,” Physical Review, vol. 46, no. 7, pp. 618–622, 1934.

F. Weigend and R. Ahlrichs, “Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: design and assessment of accuracy,” Physical Chemistry Chemical Physics, vol. 7, no. 18, pp. 3297–3305, 2005.

B. Metz, H. Stoll, and M. Dolg, “Small-core multiconfiguration-Dirac-Hartree-Fock-adjusted pseudopotentials for post-d main group elements: application to PbH and PbO,” The Journal of Chemical Physics, vol. 113, no. 7, pp. 2563–2569, 2000.

C. Hättig, A. Hellweg, and A. Köhn, “Distributed memory parallel implementation of energies and gradients for second-order Møller-Plesset perturbation theory with the resolution-of-the-identity approximation,” Physical Chemistry Chemical Physics, vol. 8, no. 10, pp. 1159–1169, 2006.

F. Weigend and M. Häser, “RI-MP2: first derivatives and global consistency,” Theoretical Chemistry Accounts, vol. 97, no. 1–4, pp. 331–340, 1997.

R. Ahlrichs, M. K. Armbruster, R. A. Bachorz et al., TURBOMOLE 6.6, a development of University of Karlsruhe and Forschungszentrum Karlsruhe GmbH, 1989–2007, TURBOMOLE GmbH, since 2007, 2014, http://www.turbomole.com.

C. Adamo and V. Barone, “Toward reliable density functional methods without adjustable parameters: the PBE0 model,” The Journal of Chemical Physics, vol. 110, no. 13, pp. 6158–6170, 1999.

E. Fermi and E. Amaldi, “Le orbite s degli elementi,” Memorie dell'Accademia d'Italia, vol. 6, pp. 119–149, 1934.

D. J. Tozer and N. C. Handy, “Improving virtual Kohn-Sham orbitals and eigenvalues: application to excitation energies and static polarizabilities,” The Journal of Chemical Physics, vol. 109, no. 23, pp. 10180–10189, 1998.

H. L. Williams, E. M. Mas, K. Szalewicz, and B. Jeziorski, “On the effectiveness of monomer‐, dimer‐, and bond‐centered basis functions in calculations of intermolecular interaction energies,” The Journal of Chemical Physics, vol. 103, no. 17, pp. 7374–7391, 1995.

R. Bukowski, W. Cencek, P. Jankowski et al., “SAPT2012: an ab initio program for symmetry-adapted perturbation theory calculations of intermolecular interaction energies. Sequential and parallel versions,” University of Delaware, Newark, Del, USA; University of Warsaw, Warsaw, Poland, 2013.

C. Angeli, K. L. Bak, V. Bakken et al., “DALTON 2.0, a molecular electronic structure program,” 2005, http://www.kjemi.uio.no/software/dalton/dalton.html.

H.-J. Werner, P. J. Knowles, G. Knizia et al., “MOLPRO 2012.1, a package of ab initio programs,” University College Cardiff Consultants Limited, Cardiff, UK, 2012, http://www.molpro.net.

D. Cremer, “Møller-Plesset perturbation theory,” in Encyclopedia of Computational Chemistry, P. V. R. Schleyer, N. L. Allinger, T. Clark et al., Eds., pp. 1706–1735, John Wiley & Sons, Chichester, UK, 1998.

S. Grimme, “Improved second-order Møller-Plesset perturbation theory by separate scaling of parallel- and antiparallel-spin pair correlation energies,” The Journal of Chemical Physics, vol. 118, no. 20, pp. 9095–9102, 2003.

R. A. Distasio Jr. and M. Head-Gordon, “Optimized spin-component scaled second-order Møller-Plesset perturbation theory for intermolecular interaction energies,” Molecular Physics, vol. 105, no. 8, pp. 1073–1083, 2007.

A. Heßelmann, “Improved supermolecular second order Møller-Plesset intermolecular interaction energies using time-dependent density functional response theory,” The Journal of Chemical Physics, vol. 128, no. 14, Article ID 144112, 9 pages, 2008.

M. Pitoňák, P. Neogrády, J. Černý, S. Grimme, and P. Hobza, “Scaled MP3 non-covalent interaction energies agree closely with accurate CCSD(T) benchmark data,” ChemPhysChem, vol. 10, no. 1, pp. 282–289, 2009.

J. Gauss, “Coupled-cluster theory,” in Encyclopedia of Computational Chemistry, P. V. R. Schleyer, N. L. Allinger, T. Clark et al., Eds., pp. 615–636, John Wiley & Sons, Chichester, UK, 1998.

S. F. Boys and F. Bernardi, “The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors,” Molecular Physics, vol. 19, no. 4, pp. 553–566, 1970.

A. Karton and J. M. L. Martin, “Comment on: ‘Estimating the Hartree-Fock limit from finite basis set calculations’ [Jensen F (2005) Theoretical Chemistry Accounts 113: 267],” Theoretical Chemistry Accounts, vol. 115, no. 4, pp. 330–333, 2005.

A. Halkier, T. Helgaker, P. Jørgensen et al., “Basis-set convergence in correlated calculations on Ne, N₂, and H₂O,” Chemical Physics Letters, vol. 286, no. 3-4, pp. 243–252, 1998.

D. G. Truhlar, “Basis-set extrapolation,” Chemical Physics Letters, vol. 294, no. 1–3, pp. 45–48, 1998.

F. Neese and E. F. Valeev, “Revisiting the atomic natural orbital approach for basis sets: robust systematic basis sets for explicitly correlated and conventional correlated ab initio methods?” Journal of Chemical Theory and Computation, vol. 7, no. 1, pp. 33–43, 2011.

K. E. Riley, M. Pitoňák, P. Jurečka, and P. Hobza, “Stabilization and structure calculations for noncovalent interactions in extended molecular systems based on wave function and density functional theories,” Chemical Reviews, vol. 110, no. 9, pp. 5023–5063, 2010.

J. Řezáč, K. E. Riley, and P. Hobza, “Benchmark calculations of noncovalent interactions of halogenated molecules,” Journal of Chemical Theory and Computation, vol. 8, no. 11, pp. 4285–4292, 2012.

S. Alvarez, “A cartography of the van der Waals territories,” Dalton Transactions, vol. 42, no. 24, pp. 8617–8636, 2013.

T. Clark, M. Hennemann, J. S. Murray, and P. Politzer, “Halogen bonding: the σ-hole,” Journal of Molecular Modeling, vol. 13, no. 2, pp. 291–296, 2007.

P. Politzer, J. S. Murray, and M. C. Concha, “Halogen bonding and the design of new materials: organic bromides, chlorides and perhaps even fluorides as donors,” Journal of Molecular Modeling, vol. 13, no. 6-7, pp. 643–650, 2007.

A. Bundhun, P. Ramasami, J. S. Murray, and P. Politzer, “Trends in σ-hole strengths and interactions of F₃MX molecules (M = C, Si, Ge and X = F, Cl, Br, I),” Journal of Molecular Modeling, vol. 19, no. 7, pp. 2739–2746, 2013.

F. De Proft, W. Langenaeker, and P. Geerlings, “Ab initio determination of substituent constants in a density functional theory formalism: calculation of intrinsic group electronegativity, hardness, and softness,” The Journal of Physical Chemistry, vol. 97, no. 9, pp. 1826–1831, 1993.

P. Metrangolo, H. Neukirch, T. Pilati, and G. Resnati, “Halogen bonding based recognition processes: a world parallel to hydrogen bonding,” Accounts of Chemical Research, vol. 38, no. 5, pp. 386–395, 2005.

J.-W. Zou, Y.-J. Jiang, M. Guo et al., “Ab initio study of the complexes of halogen-containing molecules RX (X=Cl, Br, and I) and NH₃: towards understanding the nature of halogen bonding and the electron-accepting propensities of covalently bonded halogen atoms,” Chemistry—A European Journal, vol. 11, no. 2, pp. 740–751, 2005.

K. E. Riley and P. Hobza, “Investigations into the nature of halogen bonding including symmetry adapted perturbation theory analyses,” Journal of Chemical Theory and Computation, vol. 4, no. 2, pp. 232–242, 2008.

K. E. Riley and P. Hobza, “The relative roles of electrostatics and dispersion in the stabilization of halogen bonds,” Physical Chemistry Chemical Physics, vol. 15, no. 41, pp. 17742–17751, 2013.

K. Patkowski, K. Szalewicz, and B. Jeziorski, “Third-order interactions in symmetry-adapted perturbation theory,” The Journal of Chemical Physics, vol. 125, no. 15, Article ID 154107, 20 pages, 2006.