Motion of Br2 molecules in clathrate cages. A computational study of the dynamic effects on its spectroscopic behavior

The anomalous halogen bonding interactions between chlorine and bromine with water in clathrate hydrates

Clathrate hydrate phases of Cl₂ and Br₂ guest molecules have been known for about 200 years. The crystal structure of these phases was recently re-determined with high accuracy by single crystal X-ray diffraction. In these structures, the water oxygen–halogen atom distances are determined to be shorter than the sum of the van der Waals radii, which indicates the action of some type of non-covalent interaction between the dihalogens and water molecules. Given that in the hydrate phases both lone pairs of each water oxygen atom are engaged in hydrogen bonding with other water molecules of the lattice, the nature of the oxygen–halogen interactions may not be the standard halogen bonds characterized recently in the solid state materials and enzyme–substrate compounds. The nature of the halogen–water interactions for the Cl₂ and Br₂ molecules in two isolated clathrate hydrate cages has recently been studied with ab initio calculations and Natural Bond Order analysis (Ochoa-Resendiz et al. J. Chem. Phys. 2016, 145, 161104). Here we present the results of ab initio calculations and natural localized molecular orbital analysis for Cl₂ and Br₂ guests in all cage types observed in the cubic structure I and tetragonal structure I clathrate hydrates to characterize the orbital interactions between the dihalogen guests and water. Calculations with isolated cages and cages with one shell of coordinating molecules are considered. The computational analysis is used to understand the nature of the halogen bonding in these materials and to interpret the guest positions in the hydrate cages obtained from the X-ray crystal structures.

Molecular simulations and density functional theory calculations of bromine in clathrate hydrate phases

Bromine forms a tetragonal clathrate hydrate structure (TS-I) very rarely observed in clathrate hydrates of other guest substances. The detailed structure, energetics, and dynamics of Br2 and Cl2 in TS-I and cubic structure I (CS-I) clathrate hydrates are studied in this work using molecular dynamics and quantum chemical calculations. X-ray diffraction studies show that the halogen-water-oxygen distances in the cages of these structures are shorter than the sum of the van der Waals radii of halogen and oxygen atoms. This suggests that the stabilizing effects of halogen bonding or other non-covalent interactions (NCIs) may contribute to the formation of the unique tetragonal bromine hydrate structure. We performed molecular dynamics simulations of Br2 and Cl2 clathrate hydrates using our previously developed five-site charge models for the dihalogen molecules [Dureckova et al. Can. J. Chem. 93, 864 (2015)] which reproduce the computed electrostatic potentials of the dihalogens and account for the electropositive σ-hole of the halogen bond donor (the dihalogen). Analysis of the radial distribution functions, enthalpies of encapsulation, velocity and orientation autocorrelation functions, and polar angle distributions are carried out for Br2 and Cl2 guests in various cages to contrast the properties of these guests in the TS-I and CS-I phases. Quantum chemical partial geometry optimizations of Br2 and Cl2 guests in the hydrate cages using the M06-2X functional give short halogen-water distances compatible with values observed in X-ray diffraction experiments. NCI plots of guest-cage structures are generated to qualitatively show the relative strength of the non-bonding interactions between dihalogens and water molecules. The differences between behaviors of Br2 and Cl2 guests in the hydrate cages may explain why bromine forms the unique TS-I phase.

Performance of local correlation methods for halogen bonding: The case of Br2-(H2O)n,n = 4,5 clusters and Br2@5(12)6(2) clathrate cage

The performance of local correlation methods is examined for the interactions present in clusters of bromine with water where the combined effect of hydrogen bonding (HB), halogen bonding (XB), and hydrogen-halogen (HX) interactions lead to many interesting properties. Local methods reproduce all the subtleties involved such as many-body effects and dispersion contributions provided that specific methodological steps are followed. Additionally, they predict optimized geometries that are nearly free of basis set superposition error that lead to improved estimates of spectroscopic properties. Taking advantage of the local correlation energy partitioning scheme, we compare the different interaction environments present in small clusters and those inside the 5(12)6(2) clathrate cage. This analysis allows a clear identification of the reasons supporting the use of local methods for large systems where non-covalent interactions play a key role.

Molecular dynamics simulation of halogen bonding in Cl2, BrCl, and mixed Cl2/Br2 clathrate hydrates

Clathrate hydrate phases of dihalogen molecules have properties that differ from those of other guest molecules of similar size. The water oxygen–chlorine distances in the structure I (sI) Cl2 hydrate are smaller than the sum of the van der Waals radii of oxygen and chlorine. Bromine hydrate forms a unique clathrate hydrate structure that is not seen in other guest substances. In mixed Cl2/Br2 structure I hydrate, the water oxygen–bromine distances are also smaller than the sum of the oxygen and bromine van der Waals radii. We previously studied the structure of three dihalogen clathrate hydrates using single crystal X-ray diffraction and described these structural features in terms of halogen bonding between the dihalogen and water molecules. In this work, we perform molecular dynamics simulations of cubic sI Cl2, mixed Cl2/Br2, and BrCl clathrate hydrate phases. We perform quantum chemical computations on the dihalogen molecules to determine the nature of σ-hole near the halogen atoms. We fit the electrostatic potential of the molecules to point charge models including dummy atoms that represent σ-holes adjacent to the halogen molecules. Molecular dynamics simulations are used to determine the lattice constants, radial distribution functions, and guest dynamics in these phases. We determine the effect of guest size and difference in halogen bonding on the properties of the clathrate hydrate phase. Simulations for the Cl2, BrCl, and mixed Cl2/Br2 hydrates are performed with small cages of the sI clathrate hydrate phases completely full or filled with experimental occupancies with Cl2 guests.