Differences in the Sublimation Energy of Benzene and Hexahalogenbenzenes Are Caused by Dispersion Energy

Quantifying the Impact of Halogenation on Intermolecular Interactions and Binding Modes of Aromatic Molecules

Halogenation of aromatic molecules is frequently used to modulate intermolecular interactions with ramifications for optoelectronic and mechanical properties. In this work, we accurately quantify and understand the nature of intermolecular interactions in perhalogenated benzene (PHB) clusters. Using benchmark binding energies from the fixed-node diffusion Monte Carlo (FN-DMC) method, we show that generalized Kohn-Sham semicanonical projected random phase approximation (GKS-spRPA) plus approximate exchange kernel (AKX) provides reliable interaction energies with mean absolute error (MAE) of 0.23 kcal/mol. Using the GKS-spRPA+AXK method, we quantify the interaction energies of several binding modes of PHB clusters ((C₆X₆)_n; X = F, Cl, Br, I; n = 2, 3). For a given binding mode, the interaction energies increase 3-4 times from X = F to X = I; the X-X binding modes have energies in the range of 2-4 kcal/mol, while the π-π binding mode has interaction energies in the range of 4-12 kcal/mol. SAPT-DFT-based energy decomposition analysis is then used to show that the equilibrium geometries are dictated primarily by the dispersion and exchange interactions. Finally, we test the accuracy of several dispersion-corrected density functional approximations and show that only the r2SCAN-D4 method has a low MAE and correct long-range behavior, which makes it suitable for large-scale simulations and for developing structure-function relationships of halogenated aromatic systems.

Three types of noncovalent interactions studied between pyrazine and XF

Three types noncovalent interactions (type I, II and III) between pyrazine (C₄H₄N₂) and XF (X = F, Cl, Br, and I) have been discovered at the MP2/aug-cc-pVTZ level. TypeI is σ-hole interaction between the positive site on the halogen X of XF and the negative site on one of the pyrazine nitrogens. Type II is counterintuitive σ-hole interaction driven by polarization between the positive site on the halogen X of XF and a portion of the pyrazine ring. Type III is an interaction between the lateral regions of the halogen X of XF and the position of the pyrazine ring. Through comparing the calculated interaction energy, we can know that the type II and type III interactions are weaker than the corresponding type I interactions, and type III interactions are weaker than the corresponding type II interactions in C₄H₄N₂-XF complexes. SAPT analysis shows that the electrostatic energy are the major source of the attraction for the type I (σ-hole) interactions while the type III interactions are mainly dispersion energy. For the type II (counterintuitive σ-hole) interactions in C₄H₄N₂-XF (X = F and Cl) complexes, electrostatic energy are the major source of the attraction, while in C₄H₄N₂-XF (X = Br and I) complexes, the electrostatic term, induction and dispersion play equally important role in the total attractive interaction. NBO analysis, AIM theory, and conceptual DFT are also being utilized.

Separation of halogenated benzenes enabled by investigation of halogen-π interactions with carbon materials

We report the existence of bimodal interactions, the π–π and halogen–π interactions, between the halogenated benzenes and aromatic materials.

The halogen–π (X–π) interaction is an intermolecular interaction between the electron-poor region of bonded halogen atoms and aromatic rings. We report an experimental evaluation of the halogen–π (X–π) interaction using liquid chromatography with carbon-material coated columns providing strong π interactions in the normal phase mode. A C₇₀-fullerene (C70)-coated column showed higher retentions for halogenated benzenes as the number of halogen substitutions increased as a result of X–π interactions. In addition, the strength of the X–π interaction increased in the order of F < Cl < Br < I. Changes to the UV absorption of C70 and the brominated benzenes suggested that the intermolecular interaction changed from the π–π interaction to X–π interaction as the number of bromo substitutions increased. Computer simulations also showed that the difference in dipole moments among structural isomers affected the strength of the π–π interaction. Furthermore, we concluded from small peak shifts in ¹H NMR and from computer simulations that the orbital interaction contributes to the X–π interactions. Finally, we succeeded in the one-pot separation of all isomers of brominated benzenes using the C70-coated column by optimizing the mobile phase conditions.

Comparison of the DFT-SAPT and Canonical EDA Schemes for the Energy Decomposition of Various Types of Noncovalent Interactions

Interaction energies computed with density functional theory can be divided into physically meaningful components by symmetry-adapted perturbation theory (DFT-SAPT) or the canonical energy decomposition analysis (EDA). In this work, the decomposition results obtained by these schemes were compared for more than 200 hydrogen-, halogen-, and pnicogen-bonded, dispersion-bound, and mixed complexes to investigate their similarity in the evaluation of the nature of noncovalent interactions. BLYP functional with D3(BJ) correction was used for the EDA scheme, whereas asymptotically corrected PBE0 functional for DFT-SAPT provided some of the best combinations for description of noncovalent interactions. Both schemes provide similar results concerning total interaction energies and insight into the individual energy components. For most complexes, the dominant energetic term was identified equally by both decomposition schemes. Because the canonical EDA is computationally less demanding than the DFT-SAPT, the former can be especially used in cases where the systems investigated are very large.

Synthesis, characterization and crystal structure of 2-chloroethyl(methylsulfonyl)methanesulfonate

This is a study of structure – reactivity relationship of clomesone.

Theoretical study of intermolecular interactions in crystalline arene–perhaloarene adducts in terms of the electron density

The intermolecular interactions in C₆X₆–arene crystals (X = F, Cl) and the halogen substitution effect can be quantified by the electron density.

Noncovalent π⋅⋅⋅π interaction between graphene and aromatic molecule: structure, energy, and nature

Noncovalent π⋅⋅⋅π interactions between graphene and aromatic molecules have been studied by using density functional theory with empirical dispersion correction (ωB97X-D) combined with zeroth-order symmetry adapted perturbation theory (SAPT0). Excellent agreement of the interaction energies computed by means of ωB97X-D and spin component scaled (SCS) SAPT0 methods, respectively, shows great promise for the two methods in the study of the adsorption of aromatic molecules on graphene. The other important finding in this study is that, according to SCS-SAPT0 analyses, π⋅⋅⋅π interactions between graphene and aromatic molecules are largely dependent on both dispersion and electrostatic type interactions. It is also noticed that π⋅⋅⋅π interactions become stronger and more dispersive (less electrostatic) upon substitution of the very electronegative fluorine atoms onto the aromatic molecules.

Tuning band gaps of BN nanosheets and nanoribbons via interfacial dihalogen bonding and external electric field

Density functional theory computations with dispersion corrections (DFT-D) were performed to investigate the dihalogen interactions and their effect on the electronic band structures of halogenated (fluorinated and chlorinated) BN bilayers and aligned halogen-passivated zigzag BN nanoribbons (BNNRs). Our results reveal the presence of considerable homo-halogen (FF and ClCl) interactions in bilayer fluoro (chloro)-BN sheets and the aligned F (Cl)-ZBNNRs, as well as substantial hetero-halogen (FCl) interactions in hybrid fluoro-BN/chloro-BN bilayer and F-Cl-ZBNNRs. The existence of interfacial dihalogen interactions leads to significant band-gap modifications for the studied BN nanosystems. Compared with the individual fluoro (chloro)-BN monolayers or pristine BNNRs, the gap reduction in bilayer fluoro-BN (B-FF-N array), hybrid fluoro-BN/chloro-BN bilayer (N-FCl-N array), aligned Cl-ZBNNRs (B-ClCl-N alignment), and hybrid F-Cl-ZBNNRs (B-FCl-N alignment) is mainly due to interfacial polarizations, while the gap narrowing in bilayer chloro-BN (N-ClCl-N array) is ascribed to the interfacial nearly-free-electron states. Moreover, the binding strengths and electronic properties of the interactive BN nanosheets and nanoribbons can be controlled by applying an external electric field, and extensive modulation from large-gap to medium-gap semiconductors, or even metals can be realized by adjusting the direction and strength of the applied electric field. This interesting strategy for band gap control based on weak interactions offers unique opportunities for developing BN nanoscale electronic devices.

Selective induced polarization through electron transfer in acetone and pyrazole ester derivatives via C–H⋯OC interaction

A set of organic compounds (pyrazole ester derivatives,viz.5-[3-(substituted)-propoxy]-3-methyl-1-phenyl-1H-pyrazole-4-carboxylic acid methyl ester and 5-[2-(substituted)-ethoxy]-3-methyl-1-phenyl-1H-pyrazole-4-carboxylic acid methyl ester) was synthesized and their affinity and stability towards the acetone molecule were tested by NMR.

Cl···Cl interactions in molecular crystals: insights from the theoretical charge density analysis

The structure, IR harmonic frequencies and intensities of normal vibrations of 20 molecular crystals with the X-Cl···Cl-X contacts of different types, where X = C, Cl, and F and the Cl···Cl distance varying from ~3.0 to ~4.0 Å, are computed using the solid-state DFT method. The obtained crystalline wave functions have been further used to define and describe quantitatively the Cl···Cl interactions via the electron-density features at the Cl···Cl bond critical points. We found that the electron-density at the bond critical point is almost independent of the particular type of the contact or hybridization of the ipso carbon atom. The energy of Cl···Cl interactions, E(int), is evaluated from the linking E(int) and local electronic kinetic energy density at the Cl···Cl bond critical points. E(int) varies from 2 to 12 kJ/mol. The applicability of the geometrical criterion for the detection of the Cl···Cl interactions in crystals with two or more intermolecular Cl···Cl contacts for the unique chlorine atom is not straightforward. The detection of these interactions in such crystals may be done by the quantum-topological analysis of the periodic electron density.