Next Article in Journal
A Strategy for Selecting “Q-Markers” of Chinese Medical Preparation via Components Transfer Process Analysis with Application to the Quality Control of Shengmai Injection
Previous Article in Journal
Neutrophil Immunomodulatory Activity of Natural Organosulfur Compounds
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Intermolecular Interaction in Methylene Halide (CH2F2, CH2Cl2, CH2Br2 and CH2I2) Dimers

1
State Key Laboratory of Environment-friendly Energy Materials, Southwest University of Science and Technology, Mianyang 621010, China
2
Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, Konkoly Thege út 29-33, 1121 Budapest, Hungary
3
Molecular and Biomolecular Physics Department, National Institute for Research and Development of Isotopic and Molecular Technologies, Donat Street, No. 67-103, Ro-400293 Cluj-Napoca, Romania
*
Authors to whom correspondence should be addressed.
Molecules 2019, 24(9), 1810; https://doi.org/10.3390/molecules24091810
Submission received: 27 April 2019 / Revised: 27 April 2019 / Accepted: 1 May 2019 / Published: 10 May 2019
(This article belongs to the Section Molecular Structure)

Abstract

:
The intermolecular interaction in difluoromethane, dichloromethane, dibromomethane, and diiodomethane dimers has been investigated using high level quantum chemical methods. The potential energy curve of intermolecular interaction along the C⋯C bond distance obtained using the coupled-cluster theory with singles, doubles, and perturbative triples excitations CCSD(T) were compared with values given by the same method, but applying the local (LCCSD(T)) and the explicitly correlated (CCSD(T)-F12) approximations. The accuracy of other theoretical methods—Hartree–Fock (HF), second order Møller–Plesset perturbation (MP2), and dispersion corrected DFT theory—were also presented. In the case of MP2 level, the canonical and the local-correlation cases combined with the density-fitting technique (DF-LMP2)theories were considered, while for the dispersion-corrected DFT, the empirically-corrected BLYP-D and the M06-2Xexchange-correlation functionals were applied. In all cases, the aug-cc-pVTZ basis set was used, and the results were corrected for the basis set superposition error (BSSE) using the counterpoise method. For each molecular system, several dimer geometries were found, and their mutual orientations were compared with the nearest neighbor orientations obtained in recent neutron scattering studies. The nature of the intermolecular interaction energy was discussed.

1. Introduction

Weak intermolecular interactions play important roles in a wide range of chemical and biological processes at the supramolecular level. These supramolecular systems are generally governed by different types of intermolecular interactions, like hydrogen bonds (H-bonds) [1], weak van der Waals (vdW) forces [2], or charge-transfer complexes [3]. Attractive interactions between aromatic π systems are one of the most studied noncovalent vdW forces responsible for many supramolecular organization and recognition processes. They are one of the most important interactions in the vertical base stacking of DNA [4,5] and also could influence the tertiary structure of proteins [6]. However, less is said about other weak vdW forces that have almost the same nature of attractive interaction, but they are formed between saturated hydrocarbons [7,8] or molecular systems that contain halogen atoms [9,10,11,12]. Studies of these interactions have become increasingly important, and many authors point out their special role in mediation of protein–protein, protein –nucleic acid, and receptor–ligand recognition and binding of molecular systems containing halogen atoms (for examples, see [13,14] and the references therein).
Interactions between molecules containing halogen atoms are generally called halogen bonds, and they were the subject of many theoretical [15,16,17,18,19] and experimental [20,21,22,23,24] studies and are discussed in many review articles [25,26,27,28]. Initially, the halogen bonding [29,30] was considered as a charge transfer effect from an electron-donating negatively-charged atom to the neighboring molecule to the σ * orbital of the covalent bond (X–Y) defined by a halogen atom (X) and the “holder” atom (Y, e.g., a carbon atom) [31,32]. However, an increasing number of theoretical investigations have highlighted the electrostatic character of these halogen bonds [18,33,34], particularly the role of the atomic quadrupolar effect [35] in the intermolecular electrostatic interaction. On the other hand, the behavior of the various halogen atoms (F, Cl, Br, or I) in different electron donor-acceptor or charge-transfer supramolecular complexes is still under debate [36,37], as well as the physical origin of these molecular interactions. For example, Osuna et al. [38] found that the relative orientation of the C–F⋯F bond affects the bond strength, while Price et al. [39] proposed an anisotropic model built by a combination of the repulsion, dispersion, and electrostatic forces in the case of –Cl⋯Cl– interaction. A comparative study on the nature and strength of weak hydrogen bonding between the C(sp 3 )–H, C(sp 2 )–H, and C(sp)–H donor bonds and F–C(sp 3 ) acceptors was presented by Grimme’s group [40]. One of their relevant conclusions was that double-zeta quality was not appropriate for the investigation of these weakly-bonded systems, but well-balanced basis sets of at least TZVPP quality are needed. Another important observation by them was that in most of the studied molecular structures, the dispersion interaction term dominated the entire attraction. Joining together all these findings about halogen–halogen and halogen–hydrogen intermolecular interactions, one can conclude that the nature of these interaction shows a complicated picture where many different physical effects, like electrostatic, exchange repulsion, permanent dipole-dipole, or induced dipole-dipole contribute to the final magnitude of the intermolecular interaction energy (IIE) [41,42,43,44].
With two hydrogen and two halogen atoms attached to the same carbon atom, methylene halides ( C H 2 X 2 , X = F , C l , B r , or I) can act either as proton donors or proton acceptors, and thus, they can easily form hydrogen-bonded dimers. In the classical H-bond system (X–H⋯Y), the X–H (proton donor) bond length increases, and the ν X H stretching band undergoes a red shift upon formation of a hydrogen bond. However, studies of Hobza et al. revealed an interesting new class of H-bonds when a halogen atom is present in the molecule. This is called an “improper H-bond”, and one of the main characteristic of them is the blue shift of the ν X H stretching frequencies [45,46,47,48]. The nature of the C–H⋯F–C bond in different fluoromethanes was studied by Kryachko et al. [49] and Ebrahimi et al. [50]; these studies demonstrated the improper H-bond character of the C–H proton donor induced by the presence of the fluorine atom. A systematic analysis of these “improper H-bonds”, for the molecular class of C H 2 X 2 , which would include also chlorine, bromine, and iodine atoms, could not be found in the literature. However, such analysis was performed by Zierkiewicz et al. for molecular families of C X 3 H and X H , (X=F, Cl, Br, or I) in interaction with water [51]. They found that the characteristics of bonding in the hydrogen halide complexes correspond to the standard H-bonding (an elongation of the X H bond and red shift of its stretch frequency), whereas those in the C X 3 H O H 2 complexes (X = F, Cl) are typical of improper blue-shifting H-bonding (a contraction of the CH bond and blue shift of the respective stretching frequency). This finding suggests to us that the nature of the halogen bond depends basically on the particularity of the donor and acceptor molecules.
In this paper, we study the series of methylene halide dimers with various theoretical methods and also relate their interactions to the closest neighbor structures in the liquid phase, as obtained from previous diffraction experiments. After a general description of the used theoretical methods, we present the possible dimer configurations for each methylene halide. This is followed by the comparison of the potential energy curves considering different high-level electron correlation methods for the strongest dimer configuration of each methylene halide with the C s common symmetry group. We discuss the nature of the intermolecular interaction in methylene halides taking into account the H⋯X, H⋯H, and X⋯X pair interactions. In the final section, we compare the computed dimer structures with the most frequent closest neighbor configurations observed in the liquid phase.

2. Computational Methods

For the computation of intermolecular interactions, local (L) electron correlation methods [52,53,54] at the second order perturbation theory level have been proven to give values that are very close to the standard Møller–Plesset perturbation theory (MP2) results, and by construction, they are virtually free of the basis set superposition error (BSSE) [53,54]. Linear scaling of the computational cost as a function of the system size [55] makes it possible to treat larger systems or to use larger basis sets. Using the density fitting (DF) approximation of the electron repulsion integrals [56,57,58], one can reduce again the computation time by about one order of magnitude, applying it both in the Hartree–Fock (HF) and LMP2 cases (DF-HF and DF-LMP2). In particular, the efficiency of DF-LMP2 method was clearly demonstrated by Hill et al. [59] in describing the π –stacked intermolecular interaction dominated by the dispersion forces in the case of the stacking benzene dimer structure. In this way, the computational cost was reduced to O ( N ) O ( N 2 ) without losing much in accuracy compared with the case of the conventional second-order Møller–Plesset perturbation theory (MP2), which scales formally with the order of O ( N 5 ) . In spite of the fact that the local correlation treatment combined with the density-fitting technique can, in general, provide calculations with much lower computational costs, the deficiency of the standard MP2 theory (for example, in π –stacked systems, the LMP2 method overestimates the dispersion forces) remains also an attendant of the LMP2 method. Hill et al. [59] performed a detailed theoretical investigation for different dimer configurations of benzene, and they found that the LMP2 results are quite far from the counterpoise corrected CCSD(T)values. At the same time, when the so-called spin-component scaled (SCS) MP2 theory [60] was applied (both for canonical and localized orbitals), the discrepancy was small compared with the CCSD(T) results. In a later work [61], they have shown that in case of stacked nucleic acid systems, the SCS-LMP2 method fails to describe the IIE correctly, where its mean deviation from the best estimated values of the S22set [62] is −1.62 kcal/mol (usually, it overestimates the interaction energy). Accordingly, instead of the default scaling factor of 6/5 for antiparallel spins and 1/3 for parallel spins, they completely neglected the contribution from antiparallel-spin electron pairs to the MP2 energy and scaled the parallel contribution by 1.76. In this way, they obtained a mean deviation of −0.04 kcal/mol. The method is called spin-component scaled LMP2 for nucleobases (SCSN-LMP2). Choosing the spin-component scaling factors empirically, Hill et al. [59] emphasized that this methodology can no longer be considered as a truly ab initio method, since it provides only a substantial correction to the MP2 overestimation of the dispersion energy at no extra cost. To exclude the BSSE effects in the case of the conventional MP2 method, the Boys–Bernardi counterpoise correction scheme was used [63].
In the past few years, important efforts have been made to develop new, efficient approximation techniques in order to reduce the computational cost of the high electron correlation methods. Good results were obtained by introducing the explicitly-correlated R 12 theories in the electron correlation [64,65,66]. These theories bypass the slow convergence of conventional methods, by augmenting the traditional orbital expansions with a small number of terms that depend explicitly on the interelectronic distance r 12 . Modern R 12 (or F 12 ) methods [67] can deliver MP2 energies that are converged to chemical accuracy (1 kcal/mol) in triple- or even double-zeta basis sets and can also provide high accuracy much faster than the conventional methods. We used the F12bcomputing scheme [67] for the CCSD(T)-F12 method, considering the fixed amplitude ansatz ( T i j i j , 1 = 1/2 and T i j i j , 1 = 1/4) and Slater-type frozen geminals [68]. The aug-cc-pVTZ was considered as the auxiliary basis set for the density-fitting and RI expansion.
The supramolecular ab initio calculations have become very popular tools for the investigation of intermolecular interactions, especially when using the above-presented new approximation techniques. However, these methods do not provide information on the character of these interactions. To get more insight into the physical nature of the intermolecular interaction, we have used the interaction energy decomposition scheme as the symmetry-adapted perturbation-theory (SAPT) method [69] with energy components up to the second order in V. A detailed presentation of the SAPT theory is described in [69], and a brief description of the SAPT theory was also given in our previous work [70].
Using the DF-LMP2 method implemented in the Molpro program package suite [71], we have performed geometry optimization for different dimer conformations of methylene halide ( C H 2 F 2 , C H 2 C l 2 , C H 2 B r 2 , and C H 2 I 2 ) dimers considering the aug-cc-pVTZ basis set [72,73,74,75,76]. Taking the program parameter descriptions as presented in [59], we used the following input settings: (i) we considered the Pipek–Mezey (PM) localization procedure [77]; (ii) in order to solve an occurrent poor orbital localization in the PM technique, when the larger diffuse basis set were used, we eliminated the contribution of the diffuse basis functions to the localization criteria by setting the corresponding rows and columns of the overlap matrix used in the PM localization to zero. Considering the local character of the occupied and virtual orbitals in the local correlation treatment, one can easily obtain also the dispersion part (an intermolecular effect) of the correlation contribution [78].
The intermolecular potential energy curves for methylene halide dimers for several C⋯C bond distances were computed using different theoretical models. Accordingly, the Truhlar’s M06-2x and the dispersion-corrected BLYP-D results were compared with the density-fitting HF, with the density-fitting local correlation (DF-LMP2 and DF-LCCSD(T)), as well as with the explicitly-correlation (CCSD(T)-F12) methods, where the M06, M06-2X [79], and DFT-D(considering the BLYP-D XCfunctional with empirical dispersion corrections [80,81]) energy values were obtained using the NWChem program package [82]. It was shown by Bauzá et al. [83] that most of the hybrid and pure DFT functionals largely overestimate the interaction energies in halogen bonding complexes, and only the M06-2X DFT functional gives reasonably good results compared with the CCSD(T)/aug-cc-pVTZ reference method. The potential energy curves are compared also with two reference methods considering the canonical MP2 and the coupled-cluster with singles, doubles, and perturbative triples excitation CCSD(T) theories, as well as with SCSN-MP2 method, all of them implemented in the same Molpro program package suite [71]. The intermolecular energy decomposition was performed using the SAPT theory implemented in the Molpro package [84]. Molecular structures were visualized and analyzed using the open source Gabedit molecular graphics program [85].

3. Intermolecular Interactions and Dimer Structures

Throughout the following section, we present a detailed theoretical analysis for methylene halide ( C H 2 F 2 , C H 2 C l 2 , C H 2 B r 2 , and C H 2 I 2 ) dimers describing their intermolecular interactions and dimer configurations. A suite of a selected common configuration for methylene halide dimers considering their atomic volumes is presented in Figure 1.
The geometry optimizations were done using the DF-LMP2 level of theory considering the aug-cc-pVTZ basis set. We have found several geometry configurations for each methylene halide dimer: three bounded dimers for C H 2 F 2 and five bounded structures for C H 2 C l 2 , C H 2 B r 2 , and C H 2 I 2 . Since these dimer structures are very similar for different methylene halides, we denote them with A, B, C, D, and E. Configuration A is totally asymmetric ( C 1 ); B belongs to the C s point symmetry group with a single symmetry plane; C and D belong to the C 2 h with a C 2 rotation axis and a symmetry plane with inversion center; while E belongs to C 2 v symmetry group having a simple C 2 rotation axis with two symmetry planes. The graphical view of the five conformations is shown in Figure 2. Since, for the C H 2 F 2 case, the A and C dimer configurations are missing (they converged to B), we computed and compared the potential energy curves for the B-type dimer configuration for all methylene halide systems.

3.1. Methylene Chloride ( C H 2 C l 2 )

In the case of methylene chloride (or dichloromethane), all five dimer configurations were found. The strongest bounded dimer was A, while the weakest one was the E configuration. For the five methylene chloride dimers, the interaction energy diminishes in alphabetic order. The interaction energies at the HF, LMP2, and LCCSD(T) levels of theory, as well as the C⋯C bond distances and dipole moments were presented in our previous work [86]. Here, we continue to characterize the intermolecular interaction of the B dimer considering further theoretical methods such as the coupled-cluster theory with singles, doubles, and perturbative triples excitation CCSD(T) and its explicitly-correlated approximations (CCSD(T)-F12), the canonical MP2 with BSSE correction, the spin-component scaled MP2 for nucleobases, and the dispersion corrected DFT (empirically-corrected BLYP-D and M06-2X exchange-correlation functionals) methods and compare the results obtained by them. The CCSD(T) theory was taken as a reference for comparison with the other theoretical methods.
The potential energy curves obtained with these methods are shown in Figure 3. Generally, one can conclude that all applied methods gave bounded structures for the B dimer geometry of C H 2 C l 2 . In the first step, the coupled-cluster methods were compared, where, in addition, two different approximations of explicitly- (F12) and local-correlation (L) were also considered. Since the canonical CCSD(T) calculation together with the aug-cc-pVTZ basis set requires large computational effort, the IIEs were calculated only for five representative C⋯C intermolecular bond distances (3.7, 3.8, 3.9, 4.0, and 4.1 Å) around the minimum. In the right graph of Figure 3, one can observe that the LCCSD(T) energies showed an excellent agreement with the reference curve of the CCSD(T) near the energy minimum, with an average deviation of 0.057 kcal/mol. For the CCSD(T)-F12 method, the discrepancies in the energy values were somewhat larger, and we found on average a deviation of 0.258 kcal/mol; the CCSD(T)-F12 method overestimated a bit the conventional CCSD(T) results. This difference between the CCSD(T) and the CCSD(T)-F12 energies may be due to the basis set incompleteness error, which is still present in the energy values obtained with the well-known CCSD(T)/aug-cc-pVTZ level of theory, but counted by the CCSD(T)-F12 through the F12 technique. On the other hand, one should not forget that CCSD(T)-F12 uses the density-fitting and resolution of identity techniques for speeding up the electron integral calculations, the errors of which were also present in the energy difference between the two methods. Accordingly, the real effect of the basis set incompleteness error can be estimated considering the difference between the LCCSD(T) and CCSD(T)-F12 energies (0.2 kcal/mol). If one compares the conventional MP2 results with our reference method (BSSE corrected in both cases), the well-known deficiency of binding energy overestimation for the MP2 method is confirmed again. The mean deviation for MP2 was 0.415 kcal/mol, being the largest among the methods applied for the B dimer structure of C H 2 C l 2 .
As the next step, we considered theoretical methods at the second order perturbation level including different approximation techniques. The best agreement with the reference method was obtained for SCSN-MP2 with a mean deviation of 0.123 kcal/mol, followed by the LMP2 method with a value of 0.175 kcal/mol. Regarding the DFT methods with different XC functionals, they did not match so well as was found for the perturbation and coupled-cluster methods presented before. Although the discrepancy in the binding energy was not so high, the C⋯C intermolecular bond distance was much shorter than for the earlier investigated cases. At the same time, the long range behavior of the potential curves drawn by the DFT methods was quite different from those obtained with perturbational or coupled-cluster theories.
Summarizing the results presented up to now, one can conclude that the LCCSD(T), SCSN-MP2, and LMP2 results were in a very good agreement with the CCSD(T) reference method. Furthermore, the CCSD(T)-F12 method can also be considered as a solution to reproduce the high level electron correlation effect with an acceptable accuracy. At the same time, the conventional MP2 and the two DFT results were a bit far from the correct description of the full amount of the electron correlation.
As already mentioned, in our previous work [86], we compared the IIEs for the five dimer geometries obtained at different levels of theory. Here, we complete this investigation by introducing the intermolecular energy decomposition scheme defined by the SAPT, as well as by the local electron correlation methods. The complete set of energy values for intermolecular interactions, their dispersion components, and the conformational energy difference of different dimers related to the total energy of structure A are presented in Table 1. Schütz et al. [78] pointed out that the dispersion energy derived from the local electron correlation theory was not directly comparable with the SAPT counterpart E d i s p 20 (the second-order dispersion energy between two monomers expanded in terms of HF molecular orbitals). They found that the E d i s p 20 values were substantially more negative than the dispersion energy defined in the framework of local electron correlation theory. Indeed, comparing the two different dispersion energy results, we can observe on average a shift of 0.6 kcal/mol in favor of the SAPT-type. Therefore, we will not compare these values between them, only show them in the same table. In contrast, the total amounts of the IIEs obtained with the SAPT and local perturbation theories could be compared. Therefore, if one considers the SAPT IIE as the reference energy, the mean absolute deviation (MAD) of the LCCSD(T) energy values was 0.14 kcal/mol, while the similar value for the LMP2 results was 0.16 kcal/mol. Based on this comparison, we can conclude that the intermolecular energy results obtained with different theories were in a good agreement. We performed another comparison for the total dimer energies of the five dimer configurations. In the last row of Table 1, the configuration energy differences are presented. The energy value of the most strongly-bounded dimer structure (Conformation A) was chosen as the reference energy. Analyzing the energy differences, one can see that the A and B structures were very close to each other, with only a 0.08-kcal/mol energy difference. The third structure of C was also a strongly-bound dimer; the difference was 0.37 kcal/mol, which is well below the thermal energy (TE) at room temperature (298 K, E k B T 0.6 kcal/mol). The last two structures of D and E were energetically far from the first two configurations, showing more than a 1-kcal/mol weaker intermolecular bonding.

3.2. Methylene Bromide ( C H 2 B r 2 )

The next investigated molecular system was the dibromomethane dimer. Here, we followed the same geometry optimization procedure as in case of the dichloromethane dimer. The equilibrium dimer geometries showed the same spatial symmetries and rank of energies. The IIEs and their dispersion components for the five dimer geometries obtained at different levels of theory, as well as conformational energy differences are presented in Table 2. The HF contribution was either close to zero or repulsive. The energy stability of the dimers was due to the electron correlations. Similar findings were obtained in the case of alkane chains [87], where the HF energy part showed a repulsive behavior, while the electron correlation gave a large attractive contribution. The strongest bounded molecular complexes were the A and B dimers. Depending on the method used for obtaining the intermolecular energy, the A dimer showed a slightly stronger interaction ( Δ E A L M P 2 = −4.15 kcal/mol and Δ E B L M P 2 = −3.98 kcal/mol, as well as Δ E A L C C S D ( T ) = −4.00 kcal/mol and Δ E B L C C S D ( T ) = −3.86 kcal/mol) in the case of LMP2 or LCCSD(T), while for the SAPT results, the B and A forms presented almost the same strength ( Δ E B S A P T = −3.85 kcal/mol and Δ E A S A P T = −4.33 kcal/mol). The nature of interaction is however different: the A dimer showed larger repulsion and stronger dispersion attraction than the B form. The other three complexes (C, D, and E) showed weaker intermolecular forces, but all their conformational energy differences were close to TE. The d(C⋯C) intermolecular bond distance varied between 3.7 and 4.4 Å; the shortest value was taken by the C dimer and the largest one by the D complex.
In order to compare the performance of the applied theoretical methods, we have computed several potential energy curves along the C⋯C line. In Figure 4, we present the potential energy curves obtained at the HF, CCSD(T), LCCSD(T), CCSD(T)-F12, MP2, LMP2, SCSN-MP2, M06-2X, and BLYP-D levels of theory using the aug-cc-pVTZ basis set. In the right side of Figure 4, the curves near the energy minimum are shown in more detail. The CCSD(T) energy was taken again as the reference potential energy curve. First, we compared methods based on the coupled-cluster theory where different theoretical approximations were included. In this case, one can see that energies obtained using localized orbitals (using the LCCSD(T) method) were in a very good agreement with the reference energies. At the energy minimum points, the difference was only 0.07 kcal/mol. The other coupled-cluster method with explicitly-correlated (-F12) approximation overestimated the energy minimum with about 0.27 kcal/mol. Again, this overestimation mostly came from the basis set incompleteness error effects, which was quite similar with the C H 2 C l 2 case (0.2 kcal/mol by taking the energy difference between the CCSD(T)-Fi2 and the LCCSD(T) energies). Almost the same amount of energy deviation was obtained in the previous case of the B configuration of the C H 2 C l 2 dimer. The second group of theoretical methods was defined in the framework of MP2 theory. For this case, the canonical MP2 and its local and spin-component scaling approximations were considered. Comparing them with the reference curve, we can see that the canonical MP2 strongly underestimated the reference energy curve, giving energy differences of more than 0.6 kcal/mol. The LMP2 and SCSN-MP2 methods were closer to the reference curve, with a difference of less than 0.2 kcal/mol. Energy curves drawn using the meta-hybrid GGAand the dispersion corrected DFT XC functionals were also in a relatively good agreement with the reference curve. This is especially true for the M06-2X XC functional. However, unlike the previous cases of coupled-cluster and MP2 theories, the DFT minima were shifted by more than 0.1 Å, showing shorter equilibrium bond distances.

3.3. Methylene Iodide ( C H 2 I 2 )

For the C H 2 I 2 dimer, besides the basis set for the iodine atom, the effective core potential in the case of core electrons was also used [76]. The same geometry optimization procedure was applied (see Section 3.1 and Section 3.2). The optimized dimer geometries showed the same spatial symmetries and rank of energies as for the previous two binary systems. The IIEs, at HF and correlation levels, as well as the C⋯C intermolecular distance and the dispersion part of the IIE for the five dimer geometries of dibromomethane obtained at the DF-LMP2 level of theory are presented in Table 3. The C⋯C distances followed a similar trend as the one found for the C H 2 C l 2 and C H 2 B r 2 cases (Table 3). The energy stability of different C H 2 I 2 dimer structures was ensured exclusively by the attractive forces at the correlation because HF components for the all five geometries were repulsive. The highest value was obtained for the C configuration ( Δ E C H F = 2.86 kcal/mol), where also the C⋯C distance was the shortest one, while the E dimer geometry showed the “weakest repulsion” ( Δ E C H F = 0.81 kcal/mol) and the longest C⋯C distance. This repulsion was outweighed by the dispersion forces, which showed up at the electron correlation level. The most strongly-bound structure was the A geometry ( Δ E A L M P 2 = −4.82 kcal/mol), followed by the B ( Δ E B L M P 2 = −4.54 kcal/mol), C ( Δ E C L M P 2 = −3.90 kcal/mol), D ( Δ E D L M P 2 = −3.85 kcal/mol), and E ( Δ E E L M P 2 = −3.31 kcal/mol) geometries. The corresponding values obtained with the SAPT theory were: A Δ E A S A P T = −4.83 kcal/mol, B Δ E B S A P T = −4.72 kcal/mol, C Δ E C S A P T = −3.46 kcal/mol, D Δ E D S A P T = −3.89 kcal/mol, and E Δ E E S A P T = −3.61 kcal/mol. It can be observed that in the case of A and D geometries, we had an excellent agreement between the LMP2 and SAPT values, but also for B and E structures, this difference remained below 0.2 kcal/mol. Only in the C case, the difference was larger, 0.44 kcal/mol. The A and B structures were relatively close to each other ( E B A C o n f = 0.25 kcal/mol, while the C, D, and E dimers showed more than 0.9-kcal/mol higher total energy values.
Also in this case, we computed potential energy curves along the C⋯C line considering different theoretical models similarly to the previous two cases (see Figure 5). Unfortunately, the BLYP-D calculation could not be performed because of the lack of atomic parameters for iodine atoms in the corresponding quantum chemistry computer code. As seen in Figure 5, left, the HF curve lied very far from the curves obtained with methods where the electron correlation was taken into account, showing that the electron correlation plays an important role in the intermolecular interaction in the C H 2 I 2 dimer structure. Considering methods only at the CCSD(T) level of theory, but with different approximation techniques, we can see that the CC method that used localized orbitals (LCCSD(T)) was in a very good agreement with the canonical CCSD(T) results. The difference between the energy minimums was only 0.15 kcal/mol. For the CCSD(T)-F12 method, the deviation was larger, around 0.55 kcal/mol. Accordingly, the basis set incompleteness error can be estimated as the energy difference between the CCSD(T)-F12 and LCCSD(T) energies (∼0.4 kcal/mol), which shows that in the case of the C H 2 I 2 , the aug-cc-pVTZ basis set was not large enough to cover most of the finite basis set effects. In the case of MP2 quality methods, the best performance was obtained for the LMP2: it gave 0.32 kcal/mol lower energy compared to the reference curve. This was followed by the SCSN-MP2 result, with 0.46 kcal/mol stronger bonding, while the worst approximation was obtained by the canonical MP2 method. It underestimated the reference energy by more than 1.0 kcal/mol. The M06-2X DFT functional gave accurate binding energy, but the equilibrium C⋯C intermolecular distance was shorter by 0.1 Å , and its long-range behavior was different from the others.

3.4. Methylene Fluoride ( C H 2 F 2 )

The intermolecular interaction in the methylene fluoride showed a different picture compared with the other three methylene halides. Only three different dimer structures were obtained, as during the optimization, both the A and C starting geometries turned over to B. The shortest C⋯C intermolecular distance was obtained for the B geometry, while the largest one for the D structure. The nature of the intermolecular interaction was also different from that in the previous three methylene halide dimers, the E structure being energetically more favorable than the D one. In this case, we already had a significant binding at the HF level, which was further increased by the electron correlation effects. The HF energies, presented in Table 4, were: Δ E B H F = −1.75 kcal/mol, Δ E D H F = −1.02 kcal/mol, and Δ E E H F = −1.70 kcal/mol, while the total IIE at the correlation level for the LMP2 case were: Δ E B L M P 2 = −2.62 kcal/mol, Δ E D L M P 2 = −1.67 kcal/mol, and Δ E E L M P 2 = −2.08 kcal/mol, and for SAPT, they were: Δ E B S A P T = −3.47 kcal/mol, Δ E D S A P T = −1.99 kcal/mol, and Δ E E S A P T = −2.99 kcal/mol. This binding energy increase was mainly due to the dispersion-type electron correlation effects, both in the LMP2 and SAPT cases. Comparing the LMP2 and SAPT energies, one can observe that the discrepancy between these methods was much larger than for the previous three methylene halides. In order to give an explanation for this discrepancy, one should compare the results with the CCSD(T) potential energy curve around the equilibrium C⋯C intermolecular position (see Figure 6b). The energy minimum obtained with the CCSD(T) method near the equilibrium geometry was 3.0 kcal/mol, which was lower than the LMP2, but higher than the SAPT results (see Figure 6a). Comparison of the potential energy curves with the reference CCSD(T) showed that the MP2, LMP2, SCSN-MP2, and LCCSD(T) methods underestimated, while the CCSD(T)-F12 method overestimated the reference values. It seems that for the case of C H 2 F 2 , both the basis set incompleteness errors and those induced by the density-fitting and RI approximations were important, but mainly they were canceling each other. Among these methods, the best approximation of the reference energy curves was obtained by the SCSN-MP2 and CCSD(T)-F12 methods, with energy differences near the equilibrium geometry (d(C⋯C) = 3.5Å) of +0.17 kcal/mol and −0.20 kcal/mol, respectively. For the DFT methods, the M06-2X overestimated by ∼3.0 kcal/mol, while the BLYP XC functional with empirical dispersion corrections underestimated by ∼0.24 kcal/mol the minimum energy. In both cases, the d(C⋯C) equilibrium distance was shifted by ±0.1Å.

3.5. Intermolecular Interaction in Methylene Halide Dimers

The range of the intermolecular d(C⋯C) distances is given by the relative orientation of the halogen or the hydrogen atoms to each other and thus by the balance of different repulsive and attractive type forces between the halogen and the hydrogen atoms. It was shown by Riley et al. [88] that the electrostatic interaction between the halogen σ -hole and the electronegative halogen bond donor is responsible for the high degree of directionality exhibited by halogen bonds. In our case, we do not have a real halogen σ -hole; there are electrostatic interactions between the electronegative halogen atoms and the electropositive protons. The final electrostatic contribution is the sum of the attractive H⋯X and the repulsive H⋯H and X⋯X pair interactions. On the other hand, these noncovalent interactions have also a strong dispersion component, and therefore, it is important to choose the computational method to treat halogen bonding systems correctly.
Each dimer configuration can be coded through the relative positions of the hydrogen and halogen atoms, considering the different numbers of H⋯X pair interactions. In this way, the A configuration is given by 2H⋯3X, the B by 3H⋯3X, the C by 2H⋯4X, the D by 2H⋯2X, and the E by 4H⋯2X. Of course, as mentioned before, not only the number of H⋯X pair interactions is important, but also the number of X⋯X and H⋯H pair interactions. Comparing the d(C⋯C) pair distances in the five dimer structures of dichloro-, dibromo-, and diiodomethane systems, we see that the most compact dimer packing was obtained for the C dimer, followed by the B, A, D, and E configurations. It is interesting that the conformational and intermolecular energy ranking order did not follow the order of the d(CC) distances. In the case of difluoromethane, the configurations A and C were missing; therefore, the most compact packing was obtained for the B form, followed by the E and D configurations. Analyzing the different energy contributions to the total intermolecular interaction energy, like electrostatic-, exchange-repulsion-, induction- and dispersion-energy terms, we observed that the total IIEs were obtained as the energy balance of the attractive electrostatic- and dispersion-type and the repulsive exchange-repulsion-type interactions (see Table 5 and Figure 7). The largest contributions were obtained for the exchange-repulsion and dispersion parts, having different signs, followed by the electrostatic contribution and the relatively small induction part. The energy difference between the exchange-repulsion and dispersion parts was in general between −0.5 and +0.6 kcal/mol, mostly canceling each other. There were only two cases—the conformation E of dichloromethane and conformation C of diiodomethane—where these values were somewhat larger (0.69 kcal/mol and −0.94 kcal/mol, respectively). These results tell us that the strength of the total IIE is determined mainly by the electrostatic contribution, namely the balance of the electrostatic attraction between hetero-atoms (H and halogens) and the electrostatic repulsion between the same type of atoms [89,90]. The relatively small induction contribution added to this electrostatic part can slightly influence the final total interaction energy values. Comparing the total IIEs for all four methylene halides obtained at the local CCSD(T) level of theory, we can conclude that in the gas phase, the most likely conformations for dichloro-, dibromo-, and diiodo-methane dimers are the A and B forms and for the difluoromethane the B form.

4. Orientational Correlations in the Liquid Phase

The energy minimum configurations of molecular clusters, obtained by quantum chemical methods, are often useful approximations for the arrangement of the neighboring molecules in the liquid phase. In cases were there are strong specific interactions, such as hydrogen bonding, hydrogen bonded dimers or larger associates are abundant in the liquid phase, as evidenced by diffraction experiments and molecular dynamics simulations. Our molecules lacked such strong characteristic interactions, and the pairwise interaction energies were relatively weak; therefore, the main contribution can be expected to come from vdW forces, and the mutual arrangement and orientation of the close neighbors is likely to be determined by the steric effects. It is interesting to see to what extent the pair interactions determine the mutual arrangement of the closest neighbors in the liquid and if it is possible to learn something about, or predict the mutual orientation, using only ab initio methods for pairs of molecules. Here, we make an attempt to compare the ab initio geometries of the dimers, with the most detailed information available to date from diffraction experiments on the corresponding molecular liquids.
Traditionally, the orientational ordering in liquid phase is described by taking into account the angle between two characteristic vectors of the neighboring molecules, and mapping this angle as a function of distance. Recent studies by Pothoczki et al. [91,92] on liquid methylene chloride, bromide, and iodide revealed more detailed features of the mutual orientation of these molecules in the liquid state, not known before. This advance has been achieved by applying a novel geometrical method of analysis of the atomic configurations obtained in force-field computer simulations or neutron and X-ray diffraction experiments combined with reverse Monte Carlo (RMC) simulations.
The most common description of the orientational order for asymmetric molecules, having a unique axis (such as dipole moment), is the variation of the angle between the dipole moments of two molecules with the distance between them. As shown first by Rey [93], substantially more detailed information can be extracted from a three-dimensional atomic model of a liquid, introducing some classification of the mutual orientations of the neighboring molecules. Depending on the shape of the molecule, certain groups can be defined, and the number of molecule pairs belonging to the given group can be calculated. This recipe is not universal, but specific for the geometry of the molecules. For the mutual orientation of tetrahedral molecules, Rey defined six groups, distinguished by the orientation of the closest corner atoms of the pair of molecules, and gave a recipe for calculation using the coordinates of the corner atoms of the tetrahedron [93]. In brief, given a pair of tetrahedral molecules, two parallel planes containing the molecule centers are constructed, and the molecular pairs are classified by the number of the corners belonging to the two molecules, laying between the planes, i.e., looking towards each other. This scheme results in six groups (1:1, corner-to-corner, 1:2, corner-to-edge, and so on.) For molecules with distorted tetrahedral shape, such as methylene halides, or haloforms, the original classification of Rey had been expanded by Pothoczki et al., by introducing several subgroups in which not only the number of atoms between the planes was counted, but also their type was taken into account. For methylene halides, this classification gave 28 subgroups within the six Rey groups, according to the closeness of the given type of atom on one molecule to a distinct atom, side, or edge of the other one [92].
Calculation of the number of pairs of neighbors belonging to a certain group in the simulation boxes can show the most probable mutual orientations in the simulated liquid model.

5. Comparison of Nearest Neighbor Orientations in the Liquid and in the Gas Phase

The structure of the difluoromethane dimer in the gas phase has been probed by molecular beam microwave adsorption spectroscopy [20,94]. In both studies, the dimer rotational spectra were found to correspond to the strongest, B (face-to-face) dimer configuration, having triple hydrogen bonding structure, with short, improper, or anti-hydrogen bonding. Studies of the heavier methylene halides have not been reported to date. In a recent work, difluoromethane-difluoroethane heterodimers have been studied by pulsed-jet Fourier transform rotational spectroscopy [95]. Interestingly, three dimer configurations could be identified in the spectra, with relative populations following the order of interaction energies calculated at the MP2/6-311++G(3df,3pd) level of theory. The most abundant dimer is a linear one, the C H F 2 group of difluoroethane oriented to difluoromethane in the same way as the B configuration of the difluoromethane dimer. Larger clusters of three and four C H 2 F 2 molecules have also been measured [96,97], and for each case, one computed configuration could be identified in the corresponding rotational spectra [98]. In the trimer, two molecules are oriented to each other like a distorted D configuration of the dimers.
As already mentioned, the computed dimers of the three methylene halides look rather similar, and they can be grouped according to the classification of Pothoczki et al. [92], into the following subgroups:
A: Y,H - Y,Y,H (edge-to-face)
B: H,H,Y - Y,Y,H (face-to-face)
C: Y,Y,H - Y,Y,H face-to-face
D: Y,Y - H,H edge-to-edge
E: H,H,Y - H,H,Y face-to-face
At first glance, it is seen that there are no corner-to-edge or corner-to-face configurations, and that three out of five are face-to-face configurations. It is interesting to look for whether such orientations can be found in the liquid state as well. According to the diffraction data, the first neighbor mutual orientations in the three liquids are dominantly edge-to-face and edge-to-edge, and they account for about 80% of the contact pairs [92]. The lack of corner orientations amongst the dimers is due to the highly unfavorable energy of such configurations, as they would counteract the vdW attraction, as well as to the optimal dipole-dipole orientation. Furthermore, no corner-to-corner neighbors are present in the liquid. The edge-to-edge orientation can be regarded as the intermediate level of “compactness”, situated between the closest face-to-edge/face and the loosest corner-to-corner/edge orientations. Such an arrangement of the molecules in the liquid, in which most of the contacts are of the intermediate compactness, is therefore highly plausible, since the closest orientations for the majority of molecules cannot be simultaneously achieved.
From the analysis of X-ray and neutron scattering measurements of methylene halide molecular liquids [91], Pothoczki et al. summarized the results in a few statements [92]. We cite here (with minor text editing) those that concern the closest neighbor molecular pairs and relate them to the ab initio results for methylene halide dimers.
  • “Dipolar effects, although they are visible, do not have a decisive role in forming pairwise molecular arrangements in C H 2 X 2 liquids. Steric effects, on the other hand, are more important.” Comparing the five calculated configurations, we see that fully-antiparallel dipole orientations appear in Structures C and E, while a parallel dipole orientation is only seen for Structure D. In the remaining two structures with strongest binding, the dipole vectors are nearly perpendicular. Thus, it can be seen that also for the dimers in a vacuum, the dipolar forces are weak, and the attraction of the molecules is rather due to vdW interactions (which are often referred to as steric forces when discussing RMC computer modeling). The weak contribution of the electrostatic forces in molecular liquids C B r C l 3 and C H C l 3 was also confirmed in molecular dynamics simulations with Coulomb interactions turned off [99,100].
  • “The most frequent orientation of molecules is of the 2:2 (edge-to-edge) type over the entire distance range in each liquid. Within the 2:2 original group, the H,X-H,X subgroup is the most prominent, apart from the short range orientations in C H 2 C l 2 where the H,H-H,Cl arrangement is the most frequent.” As mentioned above, we can attribute the dominance of the edge-to-edge orientations to the optimal intermediately-compact arrangement of the molecules. The only 2:2-type mutual orientation occurs in the weakly-bound D-type dimer, in which the molecules are aligned with parallel dipole moments, the corresponding subgroup being H,H-X,X. This orientation however does not appear in the liquid, indicating that the dipole interaction is too weak to orient the neighboring molecules.
  • “The structure of liquid methylene chloride appears to be different from the structure of the other two materials. The origin of structural differences is the significant size difference between C H 2 C l 2 and C H 2 B r 2 / C H 2 I 2 molecules”.
From the comparisons described above, one can conclude that the mutual orientations of the closest molecules in the liquid rarely adopt one of the five calculated energy minimum configurations. For the studied methylene halide liquids, this is the consequence of the roughly spherical shape of the molecules, for which the interaction energy in the preferred orientations is not much lower than the energies of other mutual orientation. The orientation analysis of the diffraction data also revealed that only a small difference is seen between the actual liquids, as compared to a reference model liquid having tetrahedral molecules with steric repulsion, but no atomic charges.
The situation can by quite different for other molecular liquids, in which the dipole interaction combined with the highly asymmetrical shape of the molecules appreciably determines the mutual orientation of the closest neighbors, such as nitromethane [101] or 1,3,5-trifluorobenzene [102]; in these studies, the diffraction data revealed the existence of the preferred arrangement of the molecular pairs similar to those predicted by ab initio calculations. In hydrogen bonded liquids also, such as water and aqueous solutions, the structure seen by diffraction and simulations usually corresponds rather well to the force-field model predictions.
At low temperature and/or high pressures, methylene halides form crystalline phases [103,104]. Analyzing the crystal packing and cell configurations only, the D dimer configuration (head-to-tail) occurs in solid methylene bromide and iodide, while for C H 2 C l 2 , none of the A-E dimer configurations can be recognized in the crystal.

6. Conclusions

In this paper, the series of methylene halide ( C H 2 X 2 , X = F, Cl, Br, or I) dimers was studied using various theoretical methods. The conformational energy analysis has shown that in the case of the dichloro-, dibromo-, and diiodo-methane systems, there are two, energetically almost identical conformations (A and B) with the lowest total energy value. Several dimer conformations, e.g., Conformations D and E of dichloro- and dibromo-methane or C and D of the diiodomethane system have quite similar total energy, close to the thermal energy ( k B T at room temperature ≈ 0.59 kcal/mol). The largest conformational energy difference of 1.41 kcal/mol indicates that, in the liquid phase at room temperature, all dimer configurations A, B, C, D, and E can exist with lower or higher probabilities, and none of them can be excluded based on energy considerations. The intermolecular interaction energies of the most stable dimer conformations of different methylene halide dimers are between −3.6 and −4.6 kcal/mol (obtained with the LCCSD(T) method) and are strongly based on electron correlation effects. On the other hand, these electron correlation effects can be surprisingly well covered by the pair-correlation contributions, because the higher order electron correlations give a relatively small contribution to the intermolecular interaction energy. Similar to the saturated hydrocarbons, this intermolecular interaction energy is dominated mainly by the competition of the repulsive exchange and the attractive dispersion energy components, with the difference that in the methylene halide dimers, the electrostatic component has a larger contribution than that found for the saturated hydrocarbons.
The dimer orientations were compared to the closest neighbor orientations observed in the liquid phase in diffraction experiments. Both in liquid, and in vacuo, the dipole-dipole forces were not dominating the mutual interactions. The most frequent orientation observed in the liquid was not very probable in the gas phase, and the strongest dimer was relatively rare in the liquid phase, indicating that for these systems, the molecule pair interactions cannot be simply used to predict the structures expected in the liquid phase. These examples show that for molecular liquids consisting of quasi spherical molecules, more effective theoretical methods such ab initio molecular dynamics simulations or improved ab initio force-fields [105,106] are necessary.

Author Contributions

Conceptualization, L.A. and A.B.; methodology, A.B.; investigation, A.B.; writing, A.B. and L.A.; funding acquisition, A.B.

Funding

This research was funded by the Romanian National Authority for Scientific Research and Innovation (ANCSI) through the Core Program IZO-MOL-EA, Grant Number PN-19-35.

Acknowledgments

A.B. thanks the support of the Data Center of NIRDIMTCluj-Napoca for providing the computational infrastructure and technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rajput, L.; Biradha, K. Robust hydrogen bonding synthon in one-dimensional and two-dimensional coordination polymers of pyridine-appended reverse amides and amides. CrystEngComm 2009, 11, 1220–1222. [Google Scholar] [CrossRef]
  2. Verdan, S.; Melich, X.; Bernardinelli, G.; Williams, A.F. Molecular bricklaying II. CrystEngComm 2009, 11, 1416–1426. [Google Scholar] [CrossRef]
  3. Jérome, D.; Shulz, H.J. Organic conductors and superconductors. Adv. Phys. 1982, 31, 299–490. [Google Scholar] [CrossRef]
  4. Hunter, C.A.; Singh, J.; Thornton, J.M. ππ interactions: The geometry and energetics of phenylalanine-phenylalanine interactions in proteins. J. Mol. Biol. 1991, 218, 837–846. [Google Scholar] [CrossRef]
  5. Saenger, W. Principles of Nucleic Acid Structure; Springer: New York, NY, USA, 1984. [Google Scholar]
  6. Hobza, P. Theoretical studies of hydrogen bonding. Annu. Rep. Prog. Chem. Sect. C Phys. Chem. 2004, 100, 3–27. [Google Scholar] [CrossRef]
  7. Grimme, S. Do special noncovalent ππ stacking interactions really exist? Angew. Chem. Int. Ed. 2008, 47, 3430–3434. [Google Scholar] [CrossRef]
  8. Tsuzuki, S.; Honda, K.; Uchimaru, T.; Mikami, M. Estimated MP2 and CCSD(T) interaction energies of n-alkane dimers at the basis set limit: Comparison of the methods of Helgaker et al. and Feller. J. Chem. Phys. 2006, 124, 114304. [Google Scholar] [CrossRef]
  9. Li, Q.; Xu, X.; Liu, T.; Jing, B.; Li, W.; Cheng, J.; Gong, B.; Sun, J. Competition between hydrogen bond and halogen bond in complexes of formaldehyde with hypohalous acids. Phys. Chem. Chem. Phys. 2010, 12, 6837–6843. [Google Scholar] [CrossRef] [PubMed]
  10. Zhang, X.; Zeng, Y.; Li, X.; Meng, L.; Zheng, S. Comparison in the complexes of oxygen-containing σ-electron donor with hydrogen halide and dihalogen molecules. J. Mol. Struct. Theochem. 2010, 950, 27–35. [Google Scholar] [CrossRef]
  11. Cukiernik, F.D.; Zelcer, A.; Garland, M.T.; Baggio, R. Halogen bonding in 1,2-dibromo-4,5-dimethoxybenzene and 1,2-diiodo-4,5-dimethoxybenzene. Acta Cryst. 2008, C64, o604–o608. [Google Scholar] [CrossRef]
  12. Cukiernik, F.D.; Cecchi, F.; Baggio, R. Comparison of halogen bonding and van der Waals and ππ interactions in 4,5-dibromo-2-hexylphenol. Acta Cryst. 2009, C65, o233–o236. [Google Scholar]
  13. Zhou, P.; Lv, J.; Zou, J.; Tian, F.; Shang, Z. Halogen–water–hydrogen bridges in biomolecules. J. Struct Biol. 2010, 169, 172–182. [Google Scholar] [CrossRef] [PubMed]
  14. Lu, Y.; Wang, Y.; Zhu, W. Nonbonding interactions of organic halogens in biological systems: Implications for drug discovery and biomolecular design. Phys. Chem. Chem. Phys. 2010, 12, 4543–4551. [Google Scholar] [CrossRef] [PubMed]
  15. Nguyen, H.M.T.; Peeters, J.; Zeegers-Huyskens, T. Theoretical study of the blue-shifting hydrogen bonds between CH2X2 and CHX3 (X = F, Cl, Br) and hydrogen peroxide. J. Mol. Struct. 2006, 792–793, 16–22. [Google Scholar] [CrossRef]
  16. Clark, T.; Hennemann, M.; Murray, J.S.; Politzer, P. Halogen bonding: The σ-hole. J. Mol. Model. 2007, 13, 291–296. [Google Scholar] [CrossRef] [PubMed]
  17. Politzer, P.; Lane, P.; Concha, M.C.; Ma, Y.; Murray, J.S. An overview of halogen bonding. J. Mol. Model. 2007, 13, 305–311. [Google Scholar] [CrossRef] [PubMed]
  18. Torii, H.; Yoshida, M. Properties of halogen atoms for representing intermolecular electrostatic interactions related to halogen bonding and their substituent effects. J. Comput. Chem. 2010, 31, 107–116. [Google Scholar] [CrossRef]
  19. Politzer, P.; Murray, J.S. Halogen Bonding: An Interim Discussion. ChemPhysChem 2013, 14, 278–294. [Google Scholar] [CrossRef] [PubMed]
  20. Blanco, S.; López, J.C.; Lesarri, A.; Alonso, J.L. A molecular-beam Fourier transform microwave study of difluoromethane dimer. J. Mol. Struct. 2002, 612, 255–260. [Google Scholar] [CrossRef]
  21. Lago, A.F.; Kercher, J.P.; Böldi, A.; Sztáray, B.; Miller, B.; Wurzelmann, D.; Baer, T. Dissociative Photoionization and Thermochemistry of Dihalomethane Compounds Studied by Threshold Photoelectron Photoion Coincidence Spectroscopy. J. Phys. Chem. A 2005, 109, 1802–1809. [Google Scholar] [CrossRef] [PubMed]
  22. Makogon, O.; Flyunt, R.; Tobien, T.; Naumov, S.; Bonifačić, M. Dimethylselenide as a Probe for Reactions of Halogenated Alkoxyl Radicals in Aqueous Solution. Degradation of Dichloro- and Dibromomethane. J. Phys. Chem. A 2008, 112, 5908–5916. [Google Scholar] [CrossRef]
  23. Hauchecorne, F.; Herrebout, W.A. Experimental Characterization of C–X⋯Y–C (X = Br, I; Y = F, Cl) Halogen–Halogen Bonds. J. Phys. Chem. A 2013, 117, 11548–11557. [Google Scholar] [CrossRef] [PubMed]
  24. Shimizu, K.; da Silva, J.F. Halogen and Hydrogen Bonding Interplay in the Crystal Packing of Halometallocenes. Molecules 2018, 23, 2959. [Google Scholar] [CrossRef] [PubMed]
  25. Gilday, L.C.; Robinson, S.W.; Barendt, T.A.; Langton, M.J.; Mullaney, B.R.; Beer, P.D. Halogen Bonding in Supramolecular Chemistry. Chem. Rev. 2015, 115, 7118–7195. [Google Scholar] [CrossRef] [PubMed]
  26. Pothoczki, S.; Temleitner, L.; Pusztai, L. Structure of Neat Liquids Consisting of (Perfect and Nearly) Tetrahedral Molecules. Chem. Rev. 2015, 115, 13308–13361. [Google Scholar] [CrossRef]
  27. Cavallo, G.; Metrangolo, P.; Milani, R.; Pilati, T.; Priimagi, A.; Resnati, G.; Terraneo, G. The Halogen Bond. Chem. Rev. 2016, 116, 2478–2601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Wang, H.; Wang, W.; Jin, W.J. σ-Hole Bond vs π-Hole Bond: A Comparison Based on Halogen Bond. Chem. Rev. 2016, 116, 5072–5104. [Google Scholar] [CrossRef]
  29. Desiraju, G.R.; Ho, P.S.; Kloo, L.; Legon, A.C.; Marquardt, R.; Metrangolo, P.; Politzer, P.; Resnati, G.; Rissanen, K. Definition of the halogen bond (IUPAC Recommendations 2013). Pure Appl. Chem. 2013, 85, 1711–1713. [Google Scholar] [CrossRef] [Green Version]
  30. Legon, A.C. Prereactive Complexes of Dihalogens XY with Lewis Bases B in the Gas Phase: A Systematic Case for the Halogen Analogue B⋯XY of the Hydrogen Bond B⋯HX. Angew. Chem. Int. Ed. 1999, 38, 2686–2714. [Google Scholar] [CrossRef]
  31. Karpfen, A. The intermolecular interaction in the charge-transfer complexes between amines and halogens. Theor. Chem. Acc. 2003, 110, 1–10. [Google Scholar] [CrossRef]
  32. Rosokha, S.V.; Neretin, I.S.; Rosokha, T.Y.; Hecht, J.; Kochi, J.K. Charge-Transfer Character of Halogen Bonding: Molecular Structures and Electronic Spectroscopy of Carbon Tetrabromide and Bromoform Complexes with Organic σ- and π-Donors. Heteroatom Chem. 2006, 17, 449–459. [Google Scholar] [CrossRef]
  33. Awwadi, F.F.; Willett, R.D.; Peterson, K.A.; Twamley, B. The Nature of Halogen⋯Halide Synthons: Theoretical and Crystallographic Studies. J. Phys. Chem. A 2007, 111, 2319–2328. [Google Scholar] [CrossRef] [PubMed]
  34. Wang, W.; Hobza, P. Origin of the X-Hal (Hal = Cl, Br) Bond-Length Change in the Halogen-Bonded Complexes. J. Phys. Chem. A 2008, 112, 4114–4119. [Google Scholar] [CrossRef] [PubMed]
  35. Torii, H. Atomic quadrupolar effect in intermolecular electrostatic interactions of chloroalkanes: the cases of chloroform and dichloromethane. J. Mol. Liq. 2005, 119, 31–39. [Google Scholar] [CrossRef]
  36. Metrangolo, P.; Resnati, G. Type II halogen⋯halogen contacts are halogen bonds. IUCrJ 2014, 1, 5–7. [Google Scholar] [CrossRef] [PubMed]
  37. Nguyen, H.L.; Horton, P.N.; Hursthouse, M.B.; Legon, A.C.; Bruce, D.W. Halogen bonding: A new interaction for liquid crystal formation. J. Am. Chem. Soc. 2003, 126, 16–17. [Google Scholar] [CrossRef] [PubMed]
  38. Osuna, R.M.; Hernández, V.; Navarrete, J.T.L.; D’Oria, E.; Novoa, J.J. Theoretical evaluation of the nature and strength of the F⋯F intermolecular interactions present in fluorinated hydrocarbons. Theor. Chem. Acc. 2011, 128, 541–553. [Google Scholar] [CrossRef]
  39. Price, S.L.; Stone, A.J.; Lucas, J.; Rowland, R.S.; Thornley, A.E. The Nature of -Cl⋯Cl- Intermolecular Interactions. J. Am. Chem. Soc. 1994, 116, 4910–4918. [Google Scholar] [CrossRef]
  40. Hyla-Kryspin, I.; Haufe, G.; Grimme, S. Weak Hydrogen Bridges: A Systematic Theoretical Study on the Nature and Strength of C–H⋯F–C Interactions. Chem. Eur. J. 2004, 10, 3411–3422. [Google Scholar] [CrossRef] [PubMed]
  41. Rednic, M.I.; Varga, R.A.; Bende, A.; Grosu, I.G.; Miclăuş, M.; Hădade, N.D.; Terec, A.; Bogdan, B.; Grosu, I. Supramolecular anion recognition by β-HCH. Chem. Commun. 2016, 52, 12322–12325. [Google Scholar] [CrossRef]
  42. Grosu, I.G.; Rednic, M.I.; Miclăuş, M.; Grosu, I.; Bende, A. The nature of intermolecular interactions in pyridinium–anion–β-hexachlorocyclohexane molecular crystals. Phys. Chem. Chem. Phys. 2017, 19, 20691–20698. [Google Scholar] [CrossRef]
  43. Novák, M.; Marek, R.; Foroutan-Nejad, C. Anti-Electrostatic CH–Ion Bonding in Decorated Graphanes. Chem. Eur. J. 2017, 23, 14931–14936. [Google Scholar]
  44. Varadwaj, A.; Marques, H.M.; Varadwaj, P.R. Is the Fluorine in Molecules Dispersive? Is Molecular Electrostatic Potential a Valid Property to Explore Fluorine-Centered Non-Covalent Interactions? Molecules 2019, 24, 379. [Google Scholar] [CrossRef] [PubMed]
  45. Hobza, P.; Špirko, V.; Selzle, H.L.; Schlag, E.W. Anti-Hydrogen Bond in the Benzene Dimer and Other Carbon Proton Donor Complexes. J. Phys. Chem. A 1998, 102, 2501–2504. [Google Scholar] [CrossRef]
  46. Cubero, E.; Orozco, M.; Hobza, P.; Luque, F.J. Hydrogen Bond versus Anti-Hydrogen Bond: A Comparative Analysis. J. Phys. Chem. A 1999, 103, 6394–6401. [Google Scholar] [CrossRef]
  47. Van der Veken, B.; Herrebout, W.A.; Szostak, R.; Shchepkin, D.N.; Havlas, Z.; Hobza, P. The Nature of Improper, Blue-Shifting Hydrogen Bonding Verified Experimentally. J. Am. Chem. Soc. 2001, 123, 12290–12293. [Google Scholar] [CrossRef]
  48. Hobza, P.; Havlas, Z. Improper, blue-shifting hydrogen bond. Theor. Chem. Acc. 2002, 108, 325–334. [Google Scholar] [CrossRef]
  49. Kryachko, E.; Scheiner, S. CH⋯F Hydrogen Bonds. Dimers of Fluoromethanes. J. Phys. Chem. A 2004, 108, 2527–2535. [Google Scholar] [CrossRef]
  50. Ebrahimi, A.; Roohi, H.; Habibi, S.M. The characterization of stationary points in the potential energy surface of difluoromethane dimer. J. Mol. Struct. (THEOCHEM) 2004, 684, 87–93. [Google Scholar] [CrossRef]
  51. Zierkiewicz, W.; Michalska, D.; Havlas, Z.; Hobza, P. Study of the Nature of Improper Blue-Shifting Hydrogen Bonding and Standard Hydrogen Bonding in the X3CH⋯OH2 and XH⋯OH2 Complexes (X = F, Cl, Br, I): A Correlated Ab Initio Study. ChemPhysChem 2002, 3, 511–518. [Google Scholar] [CrossRef]
  52. Pulay, P. Localizability of dynamic electron correlation. Chem. Phys. Lett. 1983, 100, 151–154. [Google Scholar] [CrossRef]
  53. Saebø, S.; Pulay, P. Local treatment of electron correlation. Annu. Rev. Phys. Chem. 1993, 44, 213–236. [Google Scholar] [CrossRef]
  54. Hampel, C.; Werner, H.-J. Local treatment of electron correlation in coupled cluster theory. J. Chem. Phys. 1996, 104, 6286–6297. [Google Scholar] [CrossRef]
  55. Hetzer, G.; Schütz, M.; Stoll, H.; Werner, H.-J. Low-order scaling local correlation methods II: Splitting the Coulomb operator in linear scaling local second-order Møller–Plesset perturbation theory. J. Chem. Phys. 2000, 113, 9443–9455. [Google Scholar] [CrossRef]
  56. Schütz, M. Low-order scaling local electron correlation methods. III. Linear scaling local perturbative triples correction (T). J. Chem. Phys. 2000, 113, 9986–10001. [Google Scholar] [CrossRef]
  57. Schütz, M.; Werner, H.-J. Low-order scaling local electron correlation methods. IV. Linear scaling local coupled-cluster (LCCSD). J. Chem. Phys. 2001, 114, 661–681. [Google Scholar] [CrossRef]
  58. Vahtras, O.; Almlöf, J.; Feyereisen, M.W. Integral approximations for LCAO-SCF calculations. Chem. Phys. Lett. 1993, 213, 514–518. [Google Scholar] [CrossRef]
  59. Hill, J.G.; Platts, J.A.; Werner, H.-J. Calculation of intermolecular interactions in the benzene dimer using coupled-cluster and local electron correlation methods. Phys. Chem. Chem. Phys. 2006, 8, 4072–4078. [Google Scholar] [CrossRef]
  60. Grimme, S. Improved second-order Møller–Plesset perturbation theory by separate scaling of parallel- and antiparallel-spin pair correlation energies. J. Chem. Phys. 2003, 118, 9095–9102. [Google Scholar] [CrossRef]
  61. Hill, J.G.; Platts, J.A. Spin-Component Scaling Methods for Weak and Stacking Interactions. J. Chem. Theory Comput. 2007, 3, 80–85. [Google Scholar] [CrossRef]
  62. Jurečka, P.; Šponer, J.; Cerný, J.; Hobza, P. Benchmark database of accurate (MP2 and CCSD(T) complete basis set limit) interaction energies of small model complexes, DNA base pairs, and amino acid pairs. Phys. Chem. Chem. Phys. 2006, 8, 1985–1993. [Google Scholar] [CrossRef]
  63. Boys, S.F.; Bernardi, F. The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors. Mol. Phys. 1970, 19, 553–566. [Google Scholar] [CrossRef]
  64. Kutzelnigg, W. r12-Dependent terms in the wave function as closed sums of partial wave amplitudes for large l. Theor. Chim. Acta 1985, 68, 445–469. [Google Scholar] [CrossRef]
  65. Klopper, W.; Manby, F.R.; Ten-No, S.; Valeev, E.F. R12 methods in explicitly correlated molecular electronic structure theory. Int. Rev. Phys. Chem. 2006, 25, 427–468. [Google Scholar] [CrossRef]
  66. Noga, J.; Kedžuch, S.; Šimunek, J.; Ten-No, S. Explicitly correlated coupled cluster F12 theory with single and double excitations. J. Chem. Phys. 2008, 128, 174103. [Google Scholar] [CrossRef]
  67. Adler, T.B.; Knizia, G.; Werner, H.-J. A simple and efficient CCSD(T)-F12 approximation. J. Chem. Phys. 2007, 127, 221106. [Google Scholar] [CrossRef]
  68. May, A.J.; Manby, F.R. An explicitly correlated second order Møller–Plesset theory using a frozen Gaussian geminal. J. Chem. Phys. 2004, 121, 4479–4485. [Google Scholar] [CrossRef]
  69. Jeziorski, B.; Moszyński, R.; Szalewicz, K. Perturbation Theory Approach to Intermolecular Potential Energy Surfaces of van der Waals Complexes. Chem. Rev. 1994, 94, 1887–1930. [Google Scholar] [CrossRef]
  70. Bende, A.; Almásy, L. Weak intermolecular bonding in N,N’-dimethylethyleneurea dimers and N,N’-dimethylethyleneurea–water systems: The role of the dispersion effects in intermolecular interaction. Chem. Phys. 2008, 354, 202–210. [Google Scholar] [CrossRef]
  71. Werner, H.-J.; Knowles, P.J.; Lindh, R.; Manby, F.R.; Schütz, M.; Celani, P.; Korona, T.; Mitrushenkov, A.; Rauhut, G.; Adler, T.B.; et al. MOLPRO, Version 2010.1, a Package of ab Initio Programs. Available online: http://www.molpro.net (accessed on 9 May 2019).
  72. Dunning, T.H., Jr. Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. J. Chem. Phys. 1989, 90, 1007–1023. [Google Scholar] [CrossRef]
  73. Kendall, R.A.; Dunning, T.H., Jr.; Harrison, R.J. Electron affinities of the first-row atoms revisited. Systematic basis sets and wave functions. J. Chem. Phys. 1992, 96, 6796–6806. [Google Scholar] [CrossRef]
  74. Woon, D.E.; Dunning, T.H., Jr. Gaussian basis sets for use in correlated molecular calculations. III. The atoms aluminum through argon. J. Chem. Phys. 1993, 98, 1358–1371. [Google Scholar] [CrossRef]
  75. Wilson, A.K.; Woon, D.E.; Peterson, K.A.; Dunning, T.H., Jr. Gaussian basis sets for use in correlated molecular calculations. IX. The atoms gallium through krypton. J. Chem. Phys. 1999, 110, 7667–7676. [Google Scholar] [CrossRef]
  76. Peterson, K.A.; Shepler, B.C.; Figgen, D.; Stoll, H. On the Spectroscopic and Thermochemical Properties of ClO, BrO, IO, and Their Anions. J. Phys. Chem. A 2006, 110, 13877–13883. [Google Scholar] [CrossRef] [PubMed]
  77. Pipek, J.; Mezey, P.G. A fast intrinsic localization procedure applicable for ab initio and semiempirical linear combination of atomic orbital wave functions. J. Chem. Phys. 1989, 90, 4916–4926. [Google Scholar] [CrossRef]
  78. Schütz, M.; Rauhut, G.; Werner, H.-J. Local Treatment of Electron Correlation in Molecular Clusters: Structures and Stabilities of (H2O)n, n = 2–4. J. Phys. Chem. A 1998, 102, 5997–6003. [Google Scholar] [CrossRef]
  79. Zhao, Y.; Truhlar, D.G. The M06 Suite of Density Functionals for Main Group Thermochemistry, Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06 Functionals and Twelve Other Functionals. Theo. Chem. Acc. 2008, 120, 215–241. [Google Scholar] [CrossRef]
  80. Grimme, S. Accurate description of van der Waals complexes by density functional theory including empirical corrections. J. Comput. Chem. 2004, 25, 1463–1473. [Google Scholar] [CrossRef]
  81. Grimme, S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 2006, 27, 1787–1799. [Google Scholar] [CrossRef] [PubMed]
  82. Bylaska, E.J.; de Jong, W.A.; Govind, N.; Kowalski, K.; Straatsma, T.P.; Valiev, M.; Wang, D.; Apra, E.; Windus, T.L.; Hammond, J.; et al. NWChem, A Computational Chemistry Package for Parallel Computers, Version 6.1; Pacific Northwest National Laboratory: Richland, WA, USA, 2009.
  83. Bauzá, A.; Alkorta, I.; Frontera, A.; Elguero, J. On the Reliability of Pure and Hybrid DFT Methods for the Evaluation of Halogen, Chalcogen, and Pnicogen Bonds Involving Anionic and Neutral Electron Donors. J. Chem. Theory Comput. 2013, 9, 5201–5210. [Google Scholar] [CrossRef] [Green Version]
  84. Heßelmann, A.; Jansen, G. The helium dimer potential from a combined density functional theory and symmetry-adapted perturbation theory approach using an exact exchange–correlation potential. Phys. Chem. Chem. Phys. 2003, 5, 5010–5014. [Google Scholar] [CrossRef]
  85. Allouche, A.-R. Gabedit—A graphical user interface for computational chemistry softwares. J. Comput. Chem. 2011, 32, 174–182. [Google Scholar] [CrossRef]
  86. Almásy, L.; Bende, A. Ab initio structures of interacting methylene chloride molecules with comparison to the liquid phase. J. Mol. Liq. 2011, 158, 205–207. [Google Scholar] [CrossRef] [Green Version]
  87. Bende, A.; Grosu, I.; Turcu, I. Molecular Modeling of Phenothiazine Derivatives: Self-Assembling Properties. J. Phys. Chem. A 2010, 114, 12479–12489. [Google Scholar] [CrossRef]
  88. Riley, K.E.; Hobza, P. The relative roles of electrostatics and dispersion in the stabilization of halogen bonds. Phys. Chem. Chem. Phys. 2013, 15, 17742–17751. [Google Scholar] [CrossRef]
  89. Stone, A.J. Are Halogen Bonded Structures Electrostatically Driven? J. Am. Chem. Soc. 2013, 135, 7005–7009. [Google Scholar] [CrossRef]
  90. Hill, J.G.; Legon, A.C. On the directionality and non-linearity of halogen and hydrogen bonds. Phys. Chem. Chem. Phys. 2015, 17, 858–867. [Google Scholar] [CrossRef] [Green Version]
  91. Pothoczki, S.; Kohara, S.; Pusztai, L. Partial radial distribution functions of methylene halide molecular liquids. J. Mol. Liq. 2010, 153, 112–116. [Google Scholar] [CrossRef]
  92. Pothoczki, S.; Temleitner, L.; Pusztai, L. Extended orientational correlation study for molecular liquids containing distorted tetrahedral molecules: Application to methylene halides. J. Chem. Phys. 2010, 132, 164511. [Google Scholar] [CrossRef]
  93. Rey, R. Quantitative characterization of orientational order in liquid carbon tetrachloride. J. Chem. Phys. 2007, 126, 164506. [Google Scholar] [CrossRef]
  94. Caminati, W.; Melandri, S.; Moreschini, P.; Favero, P.G. The C-F⋯H-C “anti-hydrogen bond” in the gas phase: Microwave structure of the difluoromethane dimer. Angew. Chem. Int. Ed. 1999, 38, 2924–2925. [Google Scholar] [CrossRef]
  95. Lu, T.; Chen, J.; Zhang, J.; Gou, Q.; Xia, Z.; Feng, G. Conformational landscape of the weakly bound difluoromethane-1,1-difluoroethane dimer explored by rotational spectroscopy and quantum chemical calculations. J. Mol. Spectr. 2019, 357, 32–37. [Google Scholar] [CrossRef]
  96. Blanco, S.; Melandri, S.; Ottaviani, P.; Caminati, W. Shapes and Noncovalent Interactions of Oligomers: The Rotational Spectrum of the Difluoromethane Trimer. J. Am. Chem. Soc. 2007, 129, 2700–2703. [Google Scholar] [CrossRef] [PubMed]
  97. Feng, G.; Evangelisti, L.; Cacelli, I.; Carbonaro, L.; Prampolini, G.; Caminati, W. Oligomers based on weak hydrogen bond networks: A rotational study of the tetramer of difluoromethane. ChemComm 2014, 50, 171–173. [Google Scholar] [CrossRef]
  98. Prampolini, G.; Carbonaro, L.; Feng, G.; Evangelisti, L.; Caminati, W.; Cacelli, I. Computational Screening of Weak Hydrogen Bond Networks: Predicting Stable Structures for Difluoromethane Oligomers. J. Chem. Theor. Comp. 2014, 10, 2204–2211. [Google Scholar] [CrossRef] [PubMed]
  99. Pothoczki, S.; Ottochian, A.; Rovira-Esteva, M.; Pardo, L.C.; Tamarit, J.L.; Cuello, G.J. Role of steric and electrostatic effects in the short-range order of quasitetrahedral molecular liquids. Phys. Rev. B 2012, 85, 014202. [Google Scholar] [CrossRef]
  100. Karnes, J.J.; Benjamin, I. On the local intermolecular ordering and dynamics of liquid chloroform. J. Mol. Liq. 2017, 248, 121–126. [Google Scholar] [CrossRef]
  101. Megyes, T.; Bálint, S.; Grósz, T.; Radnai, T.; Bakó, I.; Almásy, L. Structure of liquid nitromethane: Comparison of simulation and diffraction studies. J. Chem. Phys. 2007, 126, 164507. [Google Scholar] [CrossRef] [Green Version]
  102. Cabaco, M.I.; Danten, Y.; Besnard, M.; Guissani, Y.; Guillot, B. Evidence of dimer formation in neat liquid 1,3,5-trifluorobenzene. Chem. Phys. Lett. 1996, 262, 120–124. [Google Scholar] [CrossRef]
  103. Podsiadło, M.; Dziubek, K.; Katrusiak, A. In situ high-pressure crystallization and compression of halogen contacts in dichloromethane. Acta Cryst. 2005, B61, 595–600. [Google Scholar] [CrossRef]
  104. Podsiadło, M.; Dziubek, K.; Szafranski, M.; Katrusiak, A. Molecular interactions in crystalline dibromomethane and diiodomethane, and the stabilities of their high-pressure and low-temperature phases. Acta Cryst. 2006, B62, 1090–1098. [Google Scholar] [CrossRef]
  105. Li, A.H.-T.; Huang, S.-C.; Chao, S.D. Molecular dynamics simulation of liquid carbon tetrachloride using ab initio force field. J. Chem. Phys. 2010, 132, 024506. [Google Scholar] [CrossRef] [PubMed]
  106. Yin, C.-C.; Li, A.H.-T.; Chao, S.D. Liquid chloroform structure from computer simulation with a full ab initio intermolecular interaction potential. J. Chem. Phys. 2013, 139, 194501. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Selected configurations for methylene halide dimers considering their atomic volumes.
Figure 1. Selected configurations for methylene halide dimers considering their atomic volumes.
Molecules 24 01810 g001
Figure 2. Relative orientation of the methylene halide monomers in the five stable dimer configurations (A—asymmetric; B—plane-symmetric; C and D—with C 2 h point group symmetry and E—with C 2 v point group symmetry).
Figure 2. Relative orientation of the methylene halide monomers in the five stable dimer configurations (A—asymmetric; B—plane-symmetric; C and D—with C 2 h point group symmetry and E—with C 2 v point group symmetry).
Molecules 24 01810 g002
Figure 3. Potential energy curves for the B form of the dichloromethane dimer obtained at the HF, CCSD(T), LCCSD(T), CCSD(T)-F12,Møller–Plesset perturbation (MP2), LMP2, SCSN-MP2, M06-2X, and BLYP-D levels of theory using the aug-cc-pVTZ basis set; (a) for larger intermolcular C⋯C bond separation (3–6 Å) and (b) around the equilibrium C⋯C distance (3.5–4.5 Å).
Figure 3. Potential energy curves for the B form of the dichloromethane dimer obtained at the HF, CCSD(T), LCCSD(T), CCSD(T)-F12,Møller–Plesset perturbation (MP2), LMP2, SCSN-MP2, M06-2X, and BLYP-D levels of theory using the aug-cc-pVTZ basis set; (a) for larger intermolcular C⋯C bond separation (3–6 Å) and (b) around the equilibrium C⋯C distance (3.5–4.5 Å).
Molecules 24 01810 g003
Figure 4. Potential energy curves for the B form of the dibromomethane dimer obtained at the HF, CCSD(T), LCCSD(T), CCSD(T)-F12, MP2, LMP2, SCSN-MP2, M06-2X, and BLYP-D levels of theory using the aug-cc-pVTZ basis set; (a) for larger intermolcular C⋯C bond separation (2.75–7.25 Å) and (b) around the equilibrium C⋯C distance (3.6–4.5 Å).
Figure 4. Potential energy curves for the B form of the dibromomethane dimer obtained at the HF, CCSD(T), LCCSD(T), CCSD(T)-F12, MP2, LMP2, SCSN-MP2, M06-2X, and BLYP-D levels of theory using the aug-cc-pVTZ basis set; (a) for larger intermolcular C⋯C bond separation (2.75–7.25 Å) and (b) around the equilibrium C⋯C distance (3.6–4.5 Å).
Molecules 24 01810 g004
Figure 5. Potential energy curves for the B form of the diiodomethane dimer obtained at the HF, CCSD(T), LCCSD(T), CCSD(T)-F12, MP2, LMP2, SCSN-MP2, M06-2X, and BLYP-D levels of theory using the aug-cc-pVTZ basis set. (a) for larger intermolcular C···C bond separation (3.25–7 Å) and (b) around the equilibrium C···C distance (3.6–5 Å).
Figure 5. Potential energy curves for the B form of the diiodomethane dimer obtained at the HF, CCSD(T), LCCSD(T), CCSD(T)-F12, MP2, LMP2, SCSN-MP2, M06-2X, and BLYP-D levels of theory using the aug-cc-pVTZ basis set. (a) for larger intermolcular C···C bond separation (3.25–7 Å) and (b) around the equilibrium C···C distance (3.6–5 Å).
Molecules 24 01810 g005
Figure 6. Potential energy curves for the B form of the difluoromethane dimer obtained at the HF, CCSD(T), LCCSD(T), CCSD(T)-F12, MP2, LMP2, SCSN-MP2, M06-2X, and BLYP-D levels of theory using the aug-cc-pVTZ basis set; (a) for larger intermolcular C⋯C bond separation (3–6 Å) and (b) around the equilibrium C⋯C distance (3–4.25 Å).
Figure 6. Potential energy curves for the B form of the difluoromethane dimer obtained at the HF, CCSD(T), LCCSD(T), CCSD(T)-F12, MP2, LMP2, SCSN-MP2, M06-2X, and BLYP-D levels of theory using the aug-cc-pVTZ basis set; (a) for larger intermolcular C⋯C bond separation (3–6 Å) and (b) around the equilibrium C⋯C distance (3–4.25 Å).
Molecules 24 01810 g006
Figure 7. The relative contribution of different energy components (electrostatic, exchange-repulsion, induction, and dispersion) for the four methylene halide dimers (for Conf. A-E see Figure 2).
Figure 7. The relative contribution of different energy components (electrostatic, exchange-repulsion, induction, and dispersion) for the four methylene halide dimers (for Conf. A-E see Figure 2).
Molecules 24 01810 g007
Table 1. The intermolecular interaction energies and their dispersion components for the five dimer geometries of dichloromethane obtained at the HF, LMP2, LCCSD(T), and SAPT levels of theory, as well as the conformational energy differences calculated by the LMP2 method (i = B, C, D, E).
Table 1. The intermolecular interaction energies and their dispersion components for the five dimer geometries of dichloromethane obtained at the HF, LMP2, LCCSD(T), and SAPT levels of theory, as well as the conformational energy differences calculated by the LMP2 method (i = B, C, D, E).
DimerUnitABCDE
d(C⋯C)Å3.913.853.654.094.28
Δ E H F kcal/mol0.57−0.041.27−0.70−1.02
Δ E L M P 2 kcal/mol−3.72−3.64−3.34−2.81−2.55
E L M P 2 D i s p kcal/mol−3.86−3.35−4.00−2.53−2.56
Δ E L C C S D ( T ) kcal/mol−3.58−3.51−3.18−2.72−2.47
Δ E S A P T kcal/mol−3.89−3.98−3.23−3.16−2.73
E S A P T D i s p kcal/mol−4.56−4.05−4.60−3.13−2.96
E i A C o n f kcal/mol-0.080.371.151.39
Table 2. The intermolecular interaction energies and their dispersion components for the five dimer geometries of dibromomethane obtained at the HF, LMP2, LCCSD(T), and SAPT levels of theory, as well as the conformational energy differences calculated by the LMP2 method (i = B, C, D, E).
Table 2. The intermolecular interaction energies and their dispersion components for the five dimer geometries of dibromomethane obtained at the HF, LMP2, LCCSD(T), and SAPT levels of theory, as well as the conformational energy differences calculated by the LMP2 method (i = B, C, D, E).
DimerUnitABCDE
d(C⋯C)Å4.043.963.714.394.25
Δ E H F kcal/mol1.170.401.960.620.17
Δ E L M P 2 kcal/mol−4.15−3.98−3.57−3.06−3.01
E L M P 2 D i s p kcal/mol−4.60−3.90−4.63−3.27−2.92
Δ E L C C S D ( T ) kcal/mol−4.00−3.86−3.40−2.97−2.93
Δ E S A P T kcal/mol−4.33−3.85−3.69−3.20−3.17
E S A P T D i s p kcal/mol−6.26−4.22−5.98−4.36−4.12
E i A C o n f kcal/mol-0.140.591.011.07
Table 3. The intermolecular interaction energies and their dispersion components for the five dimer geometries of diiodomethane obtained at the HF, LMP2, LCCSD(T), and SAPT levels of theory, as well as the conformational energy differences calculated by the LMP2 method (i = B, C, D, E).
Table 3. The intermolecular interaction energies and their dispersion components for the five dimer geometries of diiodomethane obtained at the HF, LMP2, LCCSD(T), and SAPT levels of theory, as well as the conformational energy differences calculated by the LMP2 method (i = B, C, D, E).
DimerUnitABCDE
d(C⋯C)Å4.254.163.804.534.56
Δ E H F kcal/mol2.141.362.861.640.81
Δ E L M P 2 kcal/mol−4.82−4.54−3.90−3.85−3.31
E L M P 2 D i s p kcal/mol−5.82−4.97−5.54−4.59−3.59
Δ E L C C S D ( T ) kcal/mol−4.61−4.34−3.68−3.67−3.17
Δ E S A P T kcal/mol−4.83−4.72−3.46−3.89−3.61
E S A P T D i s p kcal/mol−7.01−6.12−6.37−5.56−4.46
E i A C o n f kcal/mol-0.250.920.901.41
Table 4. The intermolecular interaction energies and their dispersion components for the five dimer geometries of difluoromethane obtained at the HF, LMP2, LCCSD(T), and SAPT levels of theory, as well as the conformational energy differences using the aug-cc-pVTZ basis set (i = B, C, D, E).
Table 4. The intermolecular interaction energies and their dispersion components for the five dimer geometries of difluoromethane obtained at the HF, LMP2, LCCSD(T), and SAPT levels of theory, as well as the conformational energy differences using the aug-cc-pVTZ basis set (i = B, C, D, E).
DimerUnitBDE
d(C⋯C)Å−3.54−3.98−3.61
Δ E H F kcal/mol−1.75−1.02−1.50
Δ E L M P 2 kcal/mol−2.62−1.67−2.08
E L M P 2 D i s p kcal/mol−1.58−1.37−0.90
Δ E L C C S D ( T ) kcal/mol−2.53−1.63−1.99
Δ E S A P T kcal/mol−3.47−1.99−2.99
E S A P T D i s p kcal/mol−2.48−1.35−2.21
E i A C o n f kcal/mol-0.550.87
Table 5. SAPT/aug-cc-pVTZ electrostatic (E E S ), exchange-repulsion (E E R ), induction (E I ), and dispersion (E D i s p ) energy terms for methylene halide ( C H 2 F 2 , C H 2 C l 2 , C H 2 B r 2 , and C H 2 I 2 ) dimers. Δ E T o t is the total SAPT intermolecular interaction energy and Δ E H F the total intermolecular interaction energy at the Hartree–Fock level. All energy values are in kcal/mol.
Table 5. SAPT/aug-cc-pVTZ electrostatic (E E S ), exchange-repulsion (E E R ), induction (E I ), and dispersion (E D i s p ) energy terms for methylene halide ( C H 2 F 2 , C H 2 C l 2 , C H 2 B r 2 , and C H 2 I 2 ) dimers. Δ E T o t is the total SAPT intermolecular interaction energy and Δ E H F the total intermolecular interaction energy at the Hartree–Fock level. All energy values are in kcal/mol.
ConfigurationComplexE ES E ER E I E Disp Δ E HF Δ E Tot
ACCl 2 H 2 −3.125.13−0.90−4.560.57−3.89
CBr 2 H 2 −4.746.340.11−6.261.17−4.33
CI 2 H 2 −4.387.14−1.75−7.012.79−4.83
BCF 2 H 2 −2.992.60−0.44−2.48−1.75−3.47
CCl 2 H 2 −3.274.69−0.84−4.05−0.04−3.98
CBr 2 H 2 −4.515.58−0.08−5.430.40−4.22
CI 2 H 2 −4.316.55−1.23−6.122.27−4.72
CCCl 2 H 2 −2.224.63−0.81−4.601.27−3.23
CBr 2 H 2 −3.445.460.18−5.981.96−3.69
CI 2 H 2 −2.135.43−1.62−6.373.32−3.46
DCF 2 H 2 −1.491.12−0.19−1.35−1.02−1.99
CCl 2 H 2 −2.083.06−0.46−3.130.23−3.16
CBr 2 H 2 −2.744.20−0.44−4.360.62−3.20
CI 2 H 2 −3.645.82−0.96−5.562.22−3.89
ECF 2 H 2 −2.862.50−0.42−2.21−1.50−2.99
CCl 2 H 2 −2.673.65−0.58−2.96−0.07−2.73
CBr 2 H 2 −2.744.27−0.79−4.120.17−3.17
CI 2 H 2 −3.514.88−0.90−4.461.39−3.61

Share and Cite

MDPI and ACS Style

Almásy, L.; Bende, A. Intermolecular Interaction in Methylene Halide (CH2F2, CH2Cl2, CH2Br2 and CH2I2) Dimers. Molecules 2019, 24, 1810. https://doi.org/10.3390/molecules24091810

AMA Style

Almásy L, Bende A. Intermolecular Interaction in Methylene Halide (CH2F2, CH2Cl2, CH2Br2 and CH2I2) Dimers. Molecules. 2019; 24(9):1810. https://doi.org/10.3390/molecules24091810

Chicago/Turabian Style

Almásy, László, and Attila Bende. 2019. "Intermolecular Interaction in Methylene Halide (CH2F2, CH2Cl2, CH2Br2 and CH2I2) Dimers" Molecules 24, no. 9: 1810. https://doi.org/10.3390/molecules24091810

APA Style

Almásy, L., & Bende, A. (2019). Intermolecular Interaction in Methylene Halide (CH2F2, CH2Cl2, CH2Br2 and CH2I2) Dimers. Molecules, 24(9), 1810. https://doi.org/10.3390/molecules24091810

Article Metrics

Back to TopTop