The origins of the directionality of noncovalent intermolecular interactions#

Osme Bond: Geometric and Energetic Features in the Adducts between OsO4 and Lewis Bases

Adducts between OsO4 and Lewis bases exert a role in important oxidation processes such as epoxidation and dihydroxylation. It has been shown that the attractive interaction driving the formation of these adducts is a σ-hole bond involving the metal as the electrophilic species; the term Osme Bond (OmB) was proposed for designating it. Here some new adducts between OsO4 and various bases have been characterized through single crystal x-ray diffraction (XRD) and computational studies (density functional theory, DFT), confirming the existence of a robust correlation between σ-hole interaction energy and deformation of the tetrahedral geometry of OsO4. Also, some adducts formed by RuO4 with nucleophiles were investigated computationally.

Uncovering the mechanisms of cyclic peptide self-assembly in membranes with the chirality-aware MA(R/S)TINI forcefield

Cyclic peptides (CPs) formed by alternation of D- and L-amino acids (D,L-CPs) can self-assemble into nanotubes (SCPNs) by parallel or/and antiparallel stacking. Different applications have been attributed to these nanotubes, including the disruption of lipid bilayers of specific compositions and the selective transport of ions throughout membranes. Molecular dynamics (MD) simulations have significantly contributed to understand the interaction between CPs, including the structural, dynamic and transport properties of their supramolecular aggregates. The high computational cost of atomic resolution forcefields makes them impractical for simulating the self-assembly of macromolecules, so coarse-grained (CG) models might represent a more feasible solution for this purpose. However, general CG models used for the simulation of biomolecules such as the MARTINI forcefield do not explicitly consider the non-covalent interactions leading to the formation of secondary structure patterns in proteins. This becomes particularly important in the case of CPs due to the D- and L-chirality alternation in their sequence, leading to opposite orientations of the backbone polar groups on both sides of the cyclic ring plane. In order to overcome this limitation, we have extended the MARTINI forcefield to introduce chirality in each residue of the CPs. The new parametrization, which we have called MA(R/S)TINI, reproduces the expected self-assembly patterns for several CP sequences in the presence of different membrane models, explicitly considering the chirality of the CPs and with no significant extra computational cost. Our simulations provide new mechanistic information of how these systems self-assemble in presence of different lipid scenarios, showing that the CP-CP and CP-membrane interactions are sensitive to the peptide sequence chirality. This opens the door to design new bioactive CPs based on CG-MD simulations. A web-based tool for the automatic parameterization of new CP sequences using MA(R/S)TINI, among other functionalities, is under construction (see http://cyclopep.com).

Intermolecular Covalent Interactions: Nature and Directionality

Quantum chemical methods were employed to analyze the nature and the origin of the directionality of pnictogen (PnB), chalcogen (ChB), and halogen bonds (XB) in archetypal F_m Z⋅⋅⋅F^- complexes (Z=Pn, Ch, X), using relativistic density functional theory (DFT) at ZORA-M06/QZ4P. Quantitative Kohn-Sham MO and energy decomposition analyses (EDA) show that all these intermolecular interactions have in common that covalence, that is, HOMO-LUMO interactions, provide a crucial contribution to the bond energy, besides electrostatic attraction. Strikingly, all these bonds are directional (i.e., F-Z⋅⋅⋅F^- is approximately linear) despite, and not because of, the electrostatic interactions which, in fact, favor bending. This constitutes a breakdown of the σ-hole model. It was shown how the σ-hole model fails by neglecting both, the essential physics behind the electrostatic interaction and that behind the directionality of electron-rich intermolecular interactions. Our findings are general and extend to the neutral, weaker ClI⋅⋅⋅NH₃ , HClTe⋅⋅⋅NH₃ , and H₂ ClSb⋅⋅⋅NH₃ complexes.

Hydrogen Bonding versus Halogen Bonding: Spectroscopic Investigation of Gas-Phase Complexes Involving Bromide and Chloromethanes

Hydrogen bonding and halogen bonding are important non-covalent interactions that are known to occur in large molecular systems, such as in proteins and crystal structures. Although these interactions are important on a large scale, studying hydrogen and halogen bonding in small, gas-phase chemical species allows for the binding strengths to be determined and compared at a fundamental level. In this study, anion photoelectron spectra are presented for the gas-phase complexes involving bromide and the four chloromethanes, CH₃ Cl, CH₂ Cl₂ , CHCl₃ , and CCl₄ . The stabilisation energy and electron binding energy associated with each complex are determined experimentally, and the spectra are rationalised by high-level CCSD(T) calculations to determine the non-covalent interactions binding the complexes. These calculations involve nucleophilic bromide and electrophilic bromine interactions with chloromethanes, where the binding motifs, dissociation energies and vertical detachment energies are compared in terms of hydrogen bonding and halogen bonding.

The roles of charge transfer and polarization in non-covalent interactions: a perspective from ab initio valence bond methods

Noncovalent interactions are ubiquitous and have been well recognized in chemistry, biology and material science. Yet, there are still recurring controversies over their natures, due to the wide range of noncovalent interaction terms. In this Essay, we employed the Valence Bond (VB) methods to address two types of interactions which recently have drawn intensive attention, i.e., the halogen bonding and the CH‧‧‧HC dihydrogen bonding. The VB methods have the advantage of interpreting molecular structures and properties in the term of electron-localized Lewis (resonance) states (structures), which thereby shed specific light on the alteration of the bonding patterns. Due to the electron localization nature of Lewis states, it is possible to define individually and measure both polarization and charge transfer effects which have different physical origins. We demonstrated that both the ab initio VB method and the block-localized wavefunction (BLW) method can provide consistent pictures for halogen bonding systems, where strong Lewis bases NH₃, H₂O and NMe₃ partake as the halogen bond acceptors, and the halogen bond donors include dihalogen molecules and XNO₂ (X = Cl, Br, I). Based on the structural, spectral, and energetic changes, we confirm the remarkable roles of charge transfer in these halogen bonding complexes. Although the weak C-H∙∙∙H-C interactions in alkane dimers and graphene sheets are thought to involve dispersion only, we show that this term embeds delicate yet important charge transfer, bond reorganization and polarization interactions.

Nature of hydrogen-bond-enhanced halogen bonding viewed through electron density changes

Elucidating the mechanism of how we can achieve fine tuning of intermolecular interaction strength will be helpful for designing functionally important molecules. In the present study, a theoretical analysis is conducted, by examining the electron density changes, for two halogen-bonding iodinated systems whose halogen-bond strengths have been considered to be enhanced by the presence of a hydrogen-bond donating group (termed hydrogen-bond-enhanced halogen bonding). It is shown that, contrary to the expectation obtained from the enhancement of electrostatic potential along the line extended from the C-I bond, the anisotropy of electron distribution on the iodine atom remains nearly the same. This means that the hydrogen bond and halogen bond contribute almost independently and additively to the enhancement of electrostatic potential, indicating the nature of this enhancement and, in a more general sense, the relationship between the strength and the extent of directionality of halogen bonding.

The use of electrostatic potential at nuclei in the analysis of halogen bonding

Molecular electrostatic potential data at interacting nuclei provide strong evidence of bond formation in many intermolecular halogen bonded complexes.

Are Heavy Pnictogen-π Interactions Really "π Interactions"?

The noncovalent interactions of heavy pnictogens with π-arenes play a fundamental role in fields like crystal engineering or catalysis. The strength of such bonds is based on an interplay between dispersion and donor/acceptor interactions, and is generally attributed to the presence of π-arenes. Computational studies of the interaction between the heavy pnictogens As, Sb and Bi and cyclohexane, in comparison with previous studies on the interaction between heavy pnictogens and benzene, show that this concept probably has to be revised. A thorough analysis of all the different energetic components that play a role in these systems, carried out with state-of-the-art computational methods, sheds light on how they influence one another and the effect that their interplay has on the overall system. Furthermore, the analysis of such interactions leads us to the unexpected finding that the presence of the pnictogen compounds strongly affects the conformational equilibrium of cyclohexane, reversing the relative stability of the chair and boat-twist conformers, and thus suggesting a possible application of tuneable dispersion energy donors to stabilise the desired conformation.

Hydrogen and Halogen Bonding in Homogeneous External Electric Fields: Modulating the Bond Strengths

Recent years have witnessed various fascinating phenomena arising from the interactions of noncovalent bonds with homogeneous external electric fields (EEFs). Here we performed a computational study to interpret the sensitivity of intrinsic bond strengths to EEFs in terms of steric effect and orbital interactions. The block-localized wavefunction (BLW) method, which combines the advantages of both ab initio valence bond (VB) theory and molecular orbital (MO) theory, and the subsequent energy decomposition (BLW-ED) approach were adopted. The sensitivity was monitored and analyzed using the induced energy term, which is the variation in each energy component along the EEF strength. Systems with single or multiple hydrogen (H) or halogen (X) bond(s) were also examined. It was found that the X-bond strength change to EEFs mainly stems from the covalency change, while generally the steric effect rules the response of H-bonds to EEFs. Furthermore, X-bonds are more sensitive to EEFs, with the key difference between H- and X-bonds lying in the charge transfer interaction. Since phenylboronic acid has been experimentally used as a smart linker in EEFs, switchable sensitivity was scrutinized with the example of the phenylboronic acid dimer, which exhibits two conformations with either antiparallel or parallel H-bonds, thereby, opposite or consistent responses to EEFs. Among the studied systems, the quadruple X-bonds in molecular capsules exhibit remarkable sensitivity, with its interaction energy increased by -95.2 kJ mol^-1 at the EEF strength 0.005 a.u.

Yet another perspective on hole interactions

Hole interactions are known by different names depending on the key atom of the bond (halogen bond, chalcogen bond, hydrogen bond, etc.), and the geometry of the interaction (σ if in line, π if perpendicular to the Lewis acid plane). However, its origin starts with the creation of a Lewis acid by an underlying covalent bond, which forms an electrostatic depletion and a virtual antibonding orbital, which can create non-covalent interactions with Lewis bases. In this (maybe subjective) perspective, we will claim that hole interactions must be defined via the molecular orbital origin of the molecule. Under this premise we can better explore the richness of such bonding patterns. For that, we will study old, recent and new systems, trying to pinpoint some misinterpretations that are often associated with them. We will use as exemplars the triel bonds, a couple of metal complexes, a discussion on convergent σ-holes, and many cases of anti-electrostatic hole interactions.

Nonspecific Binding-Fundamental Concepts and Consequences for Biosensing Applications

Nature achieves differentiation of specific and nonspecific binding in molecular interactions through precise control of biomolecules in space and time. Artificial systems such as biosensors that rely on distinguishing specific molecular binding events in a sea of nonspecific interactions have struggled to overcome this issue. Despite the numerous technological advancements in biosensor technologies, nonspecific binding has remained a critical bottleneck due to the lack of a fundamental understanding of the phenomenon. To date, the identity, cause, and influence of nonspecific binding remain topics of debate within the scientific community. In this review, we discuss the evolution of the concept of nonspecific binding over the past five decades based upon the thermodynamic, intermolecular, and structural perspectives to provide classification frameworks for biomolecular interactions. Further, we introduce various theoretical models that predict the expected behavior of biosensors in physiologically relevant environments to calculate the theoretical detection limit and to optimize sensor performance. We conclude by discussing existing practical approaches to tackle the nonspecific binding challenge in vitro for biosensing platforms and how we can both address and harness nonspecific interactions for in vivo systems.

Role of Charge Transfer in Halogen Bonding

Halogen bonding has received intensive attention recently for its applications in the construction of supramolecular assemblies and crystal engineering and its implications and potentials in chemical and biological processes and rational drug design. Peculiarly, in intermolecular interactions, halogen atoms are known as electron-donating groups carrying partial negative charges in molecules due to its high electronegativity, but they can counterintuitively act as Lewis acids and bind with Lewis bases in the form of a halogen bond. The unsettling issue regarding the nature of the halogen bonding is whether the electrostatics or charge transfer interaction dominates. The recently proposed σ-hole concept nicely reinforces the role of electrostatic attraction. Also, good correlations between the halogen bonding strength and the interaction energy from the simple point-charge model have been found. This leads to the claim that there is no need to invoke the charge transfer concept in the halogen bond. But there is alternative evidence supporting the importance of charge transfer interaction. Here, we visited a series of prominent halogen bonded complexes of the types Y₃C-X···Z (X = Br, I; Y = F, Cl, Br; Z = F^-, Cl^-, Br^-, I^-, NMe₃) with the block-localized wave function (BLW) method at the M06-2X-D3/6-311+G(d,p) (def2-SVP for iodine) level of theory. As the simplest variant of ab initio valence bond (VB) theory, the BLW method is unique in the strict localization of electrons within interacting moieties, allowing for quantitative evaluation of the charge transfer effect on geometries, spectral properties, and energetics in halogen bonding complexes. By comparing the halogen bonding complexes with and without the charge transfer interaction, we proved that the charge transfer interaction significantly shortens the X···Z bonding distance and stretches the C-X bonds. But the shortening of the halogen bonding results in the less favorable steric effect, which is composed of Pauli repulsion, electrostatics, and electron correlation. There are approximate linear correlations between the charge transfer effect and binding energy and between bonding distance and binding energy. These correlations may lead to the illusion that the charge transfer interaction is unimportant or irrelevant, but further analyses showed that the inclusion of charge transfer is critical for the proper description of the halogen bonding, as considering only electrostatics and polarization leads to only about 45-60% of the binding strengths and much elongated bonding distances.

Not Only Hydrogen Bonds: Other Noncovalent Interactions

In this review, we provide a consistent description of noncovalent interactions, covering most groups of the Periodic Table. Different types of bonds are discussed using their trivial names. Moreover, the new name “Spodium bonds” is proposed for group 12 since noncovalent interactions involving this group of elements as electron acceptors have not yet been named. Excluding hydrogen bonds, the following noncovalent interactions will be discussed: alkali, alkaline earth, regium, spodium, triel, tetrel, pnictogen, chalcogen, halogen, and aerogen, which almost covers the Periodic Table entirely. Other interactions, such as orthogonal interactions and π-π stacking, will also be considered. Research and applications of σ-hole and π-hole interactions involving the p-block element is growing exponentially. The important applications include supramolecular chemistry, crystal engineering, catalysis, enzymatic chemistry molecular machines, membrane ion transport, etc. Despite the fact that this review is not intended to be comprehensive, a number of representative works for each type of interaction is provided. The possibility of modeling the dissociation energies of the complexes using different models (HSAB, ECW, Alkorta-Legon) was analyzed. Finally, the extension of Cahn-Ingold-Prelog priority rules to noncovalent is proposed.

On the Interplay between Charge-Shift Bonding and Halogen Bonding

The nature of halogen-bond interactions has been analysed from the perspective of the astatine element, which is potentially the strongest halogen-bond donor. Relativistic quantum calculations on complexes formed between halide anions and a series of Y₃ C-X (Y=F to X, X=I, At) halogen-bond donors disclosed unexpected trends, e. g., At₃ C-At revealing a weaker donating ability than I₃ C-I despite a stronger polarizability. All the observed peculiarities have their origin in a specific component of C-Y bonds: the charge-shift bonding. Descriptors of the Quantum Chemical Topology show that the halogen-bond strength can be quantitatively anticipated from the magnitude of charge-shift bonding operating in Y₃ C-X. The charge-shift mechanism weakens the ability of the halogen atom X to engage in halogen bonds. This outcome provides rationales for outlier halogen-bond complexes, which are at variance with the consensus that the halogen-bond strength scales with the polarizability of the halogen atom.

The Nature of Hydrogen Bonds: A Delineation of the Role of Different Energy Components on Hydrogen Bond Strengths and Lengths

Hydrogen bonds are a complex interplay between different energy components, and their nature is still subject of an ongoing debate. In this minireview, we therefore provide an overview of the different perspectives on hydrogen bonding. This will be done by discussing the following individual energy components: 1) electrostatic interactions, 2) charge-transfer interactions, 3) π-resonance assistance, 4) steric repulsion, 5) cooperative effects, 6) dispersion interactions and 7) secondary electrostatic interactions. We demonstrate how these energetic factors are essential in a correct description of the hydrogen bond, and discuss several examples of systems whose energetic and geometrical features are not captured by easy-to-use predictive models.

Oriented External Electric Fields: Tweezers and Catalysts for Reactivity in Halogen-Bond Complexes

This theoretical study establishes ways of controlling and enabling an uncommon chemical reaction, the displacement reaction, B:---(X-Y) → (B-X)⁺ + :Y^-, which is nascent from a B:---(X-Y) halogen bond (XB) by nucleophilic attack of the base, B:, on the halogen, X. In most of the 14 cases examined, these reactions possess high barriers either in the gas phase (where the X-Y bond dissociates to radicals) or in solvents such as CH₂Cl₂ and CH₃CN (which lead to endothermic processes). Thus, generally, the XB species are trapped in deep minima, and their reactions are not allowed without catalysis. However, when an oriented-external electric field (OEEF) is directed along the B---X---Y reaction axis, the field acts as electric tweezers that orient the XB along the field's axis, and intensely catalyze the process, by tens of kcal/mol, thus rendering the reaction allowed. Flipping the OEEF along the reaction axis inhibits the reaction and weakens the interaction of the XB. Furthermore, at a critical OEEF, each XB undergoes spontaneous and barrier-free reaction. As such, OEEF achieves quite tight control of the structure and reactivity of XB species. Valence bond modeling is used to elucidate the means whereby OEEFs exert their control.

A Continuum from Halogen Bonds to Covalent Bonds: Where Do λ3 Iodanes Fit?

The intrinsic bonding nature of λ 3 -iodanes was investigated to determine where its hypervalent bonds fit along the spectrum between halogen bonding and covalent bonding. Density functional theory with an augmented Dunning valence triple zeta basis set ( ω B97X-D/aug-cc-pVTZ) coupled with vibrational spectroscopy was utilized to study a diverse set of 34 hypervalent iodine compounds. This level of theory was rationalized by comparing computational and experimental data for a small set of closely-related and well-studied iodine molecules and by a comparison with CCSD(T)/aug-cc-pVTZ results for a subset of the investigated iodine compounds. Axial bonds in λ 3 -iodanes fit between the three-center four-electron bond, as observed for the trihalide species IF 2 − and the covalent FI molecule. The equatorial bonds in λ 3 -iodanes are of a covalent nature. We explored how the equatorial ligand and axial substituents affect the chemical properties of λ 3 -iodanes by analyzing natural bond orbital charges, local vibrational modes, the covalent/electrostatic character, and the three-center four-electron bonding character. In summary, our results show for the first time that there is a smooth transition between halogen bonding → 3c–4e bonding in trihalides → 3c–4e bonding in hypervalent iodine compounds → covalent bonding, opening a manifold of new avenues for the design of hypervalent iodine compounds with specific properties.

Tetrel bonding interaction: an analysis with the block-localized wavefunction (BLW) approach

In this study, fifty-one iconic tetrel bonding complexes were studied using the block localized wave function (BLW) method which can derive the self-consistent wavefunction for an electron-localized (diabatic) state where charge transfer is strictly deactivated.

High-Level Ab Initio Calculations of Intermolecular Interactions: Heavy Main-Group Element π-Interactions

This work reports high-level ab initio calculations and a detailed analysis on the nature of intermolecular interactions of heavy main-group element compounds and π systems. For this purpose we have chosen a set of benchmark molecules of the form MR₃ , in which M=As, Sb, or Bi, and R=CH₃ , OCH₃ , or Cl. Several methods for the description of weak intermolecular interactions are benchmarked including DFT-D, DFT-SAPT, MP2, and high-level coupled cluster methods in the DLPNO-CCSD(T) approximation. Using local energy decomposition (LED) and an analysis of the electron density, details of the nature of this interaction are unraveled. The results yield insight into the nature of dispersion and donor-acceptor interactions in this type of system, including systematic trends in the periodic table, and also provide a benchmark for dispersion interactions in heavy main-group element compounds.

On the Computational Characterization of Charge-Transfer Effects in Noncovalently Bound Molecular Complexes

Charge-transfer (CT) is an important binding force in the formation of intermolecular complexes, and there have been a variety of theoretical models proposed to quantify this effect. These approaches, which typically rely on a definition of a "CT-free" state based on a partition of the system, sometimes yield significantly different results for a given intermolecular complex. Two widely used definitions of the "CT-free" state, the absolutely localized molecular orbitals (ALMO) method (where only on-fragment orbital mixings are permitted) and the constrained density functional theory (CDFT) approach (where fragment electron populations are fixed), are carefully examined in this work. Natural bond orbital (NBO) and the regularized symmetry-adapted perturbation theory (SAPT) are also briefly considered. Results for the ALMO and CDFT definitions of CT are compared on a broad range of model systems, including hydrogen-bonding systems, borane complexes, metal-carbonyl complexes, and complexes formed by water and metal cations. For most of these systems, CDFT yields a much smaller equilibrium CT energy compared to that given by the ALMO-based definition. This is mainly because the CDFT population constraint does not fully inhibit CT, which means that the CDFT "CT-free" state is in fact CT-contaminated. Examples of this contamination include (i) matching forward and backward donation (e.g., formic acid dimer) and (ii) unidirectional CT without changing fragment populations. The magnitude of the latter effect is quantified in systems such as the water dimer by employing a 3-space density constraint in addition to the orbital constraint. Furthermore, by means of the adiabatic EDA, it is shown that several observable effects of CT, such as the "pyramidalization" of the planar BH₃ molecule upon the complexation with Lewis bases, already appear on the "CT-free" CDFT surface. These results reveal the essential distinctions between the ALMO and CDFT definitions of CT and suggest that the former is more consistent with accepted understanding of the role of CT in intermolecular binding.

Structure-Based Analysis of Boronic Acids as Inhibitors of Acinetobacter-Derived Cephalosporinase-7, a Unique Class C β-Lactamase

Acinetobacter baumannii is a multidrug resistant pathogen that infects more than 12 000 patients each year in the US. Much of the resistance to β-lactam antibiotics in Acinetobacter spp. is mediated by class C β-lactamases known as Acinetobacter-derived cephalosporinases (ADCs). ADCs are unaffected by clinically used β-lactam-based β-lactamase inhibitors. In this study, five boronic acid transition state analog inhibitors (BATSIs) were evaluated for inhibition of the class C cephalosporinase ADC-7. Our goal was to explore the properties of BATSIs designed to probe the R1 binding site. K_i values ranged from low micromolar to subnanomolar, and circular dichroism (CD) demonstrated that each inhibitor stabilizes the β-lactamase-inhibitor complexes. Additionally, X-ray crystal structures of ADC-7 in complex with five inhibitors were determined (resolutions from 1.80 to 2.09 Å). In the ADC-7/CR192 complex, the BATSI with the lowest K_i (0.45 nM) and greatest Δ T_m (+9 °C), a trifluoromethyl substituent, interacts with Arg340. Arg340 is unique to ADCs and may play an important role in the inhibition of ADC-7. The ADC-7/BATSI complexes determined in this study shed light into the unique recognition sites in ADC enzymes and also offer insight into further structure-based optimization of these inhibitors.

Energy decomposition analysis in an adiabatic picture

Energy decomposition analysis (EDA) of electronic structure calculations has facilitated quantitative understanding of diverse intermolecular interactions. Nevertheless, such analyses are usually performed at a single geometry and thus decompose a "single-point" interaction energy. As a result, the influence of the physically meaningful EDA components on the molecular structure and other properties are not directly obtained. To address this gap, the absolutely localized molecular orbital (ALMO)-EDA is reformulated in an adiabatic picture, where the frozen, polarization, and charge transfer energy contributions are defined as energy differences between the stationary points on different potential energy surfaces (PESs), which are accessed by geometry optimizations at the frozen, polarized and fully relaxed levels of density functional theory (DFT). Other molecular properties such as vibrational frequencies can thus be obtained at the stationary points on each PES. We apply the adiabatic ALMO-EDA to different configurations of the water dimer, the water-Cl^- and water-Mg²⁺/Ca²⁺ complexes, metallocenes (Fe²⁺, Ni²⁺, Cu²⁺, Zn²⁺), and the ammonia-borane complex. This method appears to be very useful for unraveling how physical effects such as polarization and charge transfer modulate changes in molecular properties induced by intermolecular interactions. As an example of the insight obtained, we find that a linear hydrogen bond geometry for the water dimer is preferred even without the presence of polarization and charge transfer, while the red shift in the OH stretch frequency is primarily a charge transfer effect; by contrast, a near-linear geometry for the water-chloride hydrogen bond is achieved only when charge transfer is allowed.

The intrinsic strength of the halogen bond: electrostatic and covalent contributions described by coupled cluster theory

36 halogen-bonded complexes YXAR_m (X: F, Cl, Br; Y: donor group; AR_m acceptor group) have been investigated at the CCSD(T)/aug-cc-pVTZ level of theory. Binding energies, geometries, NBO charges, charge transfer, dipole moments, electrostatic potential, electron and energy density distributions, difference density distributions, vibrational frequencies, local stretching and bending force constants, and relative bond strength orders n have been calculated and used to order the halogen bonds according to their intrinsic strength. Halogen bonding is found to arise from electrostatic and strong covalent contributions. It can be strengthened by H-bonding or lone pair delocalization. The covalent character of a halogen bond increases in the way 3c-4e (three-center-four-electron) bonding becomes possible. One can characterize halogen bonds by their percentage of 3c-4e bonding. FCl-phosphine complexes can form relatively strong halogen bonds provided electronegative substituents increase the covalent contributions in form of 3c-4e halogen bonding. Binding energies between 1 and 45 kcal mol^-1 are calculated, which reflects the large variety in halogen bonding.

Site specificity of halogen bonding involving aromatic acceptors

Based on Cambridge structural database survey and quantum chemical studies, aromatic halogen bond (XB) acceptors are found to have unique pattern of XB binding sites and rim specificity.

Systematic Coupled Cluster Study of Noncovalent Interactions Involving Halogens, Chalcogens, and Pnicogens

The noncovalent interactions of 32 complexes involving pnicogens, chalcogens, and halogens atoms were investigated at the CCSD(T)/aug-cc-pVTZ level of theory. Two different types of complexes could be distinguished on the basis of geometric parameters, electron difference densities, and the charge transfer mechanisms associated with each type. In the type I conformation, the monomers adopt a skewed orientation allowing charge to be transfer between both monomers, whereas in the type II conformation the monomers adopt a linear arrangement, maximizing charge transfer in only one direction. Type I complexes involving the interaction between pnicogens and chalcogens cannot be unambiguously defined as chalcogen or pnicogen bonds, they are an admixture of both. The charge transfer dependence on the conformation adopted by the complexes described in this work can serve as a novel conformationally driven design concept for materials.

C-H···O Hydrogen Bonding. The Prototypical Methane-Formaldehyde System: A Critical Assessment

Distinguishing the functionality of C-H···O hydrogen bonds (HBs) remains challenging, because their properties are difficult to quantify reliably. Herein, we present a study of the model methane-formaldehyde complex (MFC). Six stationary points on the MFC potential energy surface (PES) were obtained at the CCSD(T)/ANO2 level. The CCSDT(Q)/CBS interaction energies of the conformers range from only -1.12 kcal mol^-1 to -0.33 kcal mol^-1, denoting a very flat PES. Notably, only the lowest energy stationary point (MFC1) corresponds to a genuine minimum, whereas all other stationary points-including the previously studied ideal case of a_e(C-H···O) = 180°-exhibit some degree of freedom that leads to MFC1. Despite the flat PES, we clearly see that the HB properties of MFC1 align with those of the prototypical water dimer O-H···O HB. Each HB property generally becomes less prominent in the higher-energy conformers. Only the MFC1 conformer prominently exhibits (1) elongated C-H donor bonds, (2) attractive C-H···O═C interactions, (3) n(O) → σ*(C-H) hyperconjugation, (4) critical points in the electron density from Bader's method and from the noncovalent interactions method, (5) positively charged donor hydrogen, and (6) downfield NMR chemical shifts and nonzero ²J(C_M-H_M···O_F) coupling constants. Based on this research, some issues merit further study. The flat PES hinders reliable determinations of the HB-induced shifts of the C-H stretches; a similarly difficult challenge is observed for the experiment. The role of charge transfer in HBs remains an intriguing open question, although our BLW and NBO computations suggest that it is relevant to the C-H···O HB geometries. These issues notwithstanding, the prominence of the HB properties in MFC1 serves as clear evidence that the MFC is predominantly bound by a C-H···O HB.

The Many Facets of Chalcogen Bonding: Described by Vibrational Spectroscopy

A diverse set of 100 chalcogen-bonded complexes comprising neutral, cationic, anionic, divalent, and double bonded chalcogens has been investigated using ωB97X-D/aug-cc-pVTZ to determine geometries, binding energies, electron and energy density distributions, difference density distributions, vibrational frequencies, local stretching force constants, and associated bond strength orders. The accuracy of ωB97X-D was accessed by CCSD(T)/aug-cc-pVTZ calculations of a subset of 12 complexes and by the CCSD(T)/aug-cc-pVTZ //ωB97X-D binding energies of 95 complexes. Most of the weak chalcogen bonds can be rationalized on the basis of electrostatic contributions, but as the bond becomes stronger, covalent contributions can assume a primary role in the strength and geometry of the complexes. Covalency in chalcogen bonds involves the charge transfer from a lone pair orbital of a Lewis base into the σ* orbital of a divalent chalcogen or a π* orbital of a double bonded chalcogen. We describe for the first time a symmetric chalcogen-bonded homodimer stabilized by a charge transfer from a lone pair orbital into a π* orbital. New polymeric materials based on chalcogen bonds should take advantage of the extra stabilization granted by multiple chalcogen bonds, as is shown for 1,2,5-telluradiazole dimers.

Effect of Fluorination on the Competition of Halogen Bonding and Hydrogen Bonding: Complexes of Fluoroiodomethane with Dimethyl Ether and Trimethylamine

To further rationalize the competition between halogen and hydrogen bonding, a combined experimental and theoretical study on the weakly bound molecular complexes formed between the combined halogen bond/hydrogen bond donor fluoroiodomethane and the Lewis bases dimethyl ether and trimethylamine (in standard and fully deuterated form) is presented. The experimental data are obtained by recording infrared and Raman spectra of mixtures of the compounds in liquid krypton, at temperatures between 120 and 156 K. The experiments are supported by ab initio calculations at the MP2/aug-cc-pVDZ-PP level, statistical thermodynamics and Monte Carlo free energy perturbation calculations. For the mixtures containing fluoroiodomethane and dimethyl ether a hydrogen-bonded complex with an experimental complexation enthalpy of -7.0(2) kJ mol^-1 is identified. Only a single weak spectral feature is observed which can be tentatively assigned to the halogen-bonded complex. For the mixtures involving trimethylamine, both halogen- and hydrogen-bonded complexes are observed, the experimental complexation enthalpies being -12.5(1) and -9.6(2) kJ mol^-1 respectively. To evaluate the influence of fluorination on the competition between halogen and hydrogen bonding, the results obtained for fluoroiodomethane are compared with those of a previous study involving difluoroiodomethane.

The role of dispersion type metal⋯π interaction in the enantiotropic phase transition of two polymorphs of tris-(thienyl)bismuthine

Bi(2-C₄H₃S)₃ shows an enantiotropic phase transition that is dominated by London dispersion forces. DFT calculations on model compounds were carried out in order to investigate the competition between Bi⋯S and Bi⋯π heteroarene interaction.

Quantitative Assessment of Halogen Bonding Utilizing Vibrational Spectroscopy

A total of 202 halogen-bonded complexes have been studied using a dual-level approach: ωB97XD/aug-cc-pVTZ was used to determine geometries, natural bond order charges, charge transfer, dipole moments, electron and energy density distributions, vibrational frequencies, local stretching force constants, and relative bond strength orders n. The accuracy of these calculations was checked for a subset of complexes at the CCSD(T)/aug-cc-pVTZ level of theory. Apart from this, all binding energies were verified at the CCSD(T) level. A total of 10 different electronic effects have been identified that contribute to halogen bonding and explain the variation in its intrinsic strength. Strong halogen bonds are found for systems with three-center-four-electron (3c-4e) bonding such as chlorine donors in interaction with substituted phosphines. If halogen bonding is supported by hydrogen bonding, genuine 3c-4e bonding can be realized. Perfluorinated diiodobenzenes form relatively strong halogen bonds with alkylamines as they gain stability due to increased electrostatic interactions.

How to Twist, Split and Warp a σ-Hole with Hypervalent Halogens

Halogen bonds (XB) are no longer newcomers in the chemistry family. However, XB in hypervalent halogens has not been thoroughly studied. We provide a molecular orbital explanation of the shape and strength of XBs in hypervalent halogens and other species, focusing on the charge transfer and electrostatic aspects of these bonds. Our results show that σ-holes (and subsequently the XBs associated with them) can be easily divided and bent by the influence of equatorial substituents. The inductive effect of both the equatorial and axial groups can affect these distortions, but also the angle between the equatorial ligands has a large influence on the shape of the σ-holes and the molecular orbitals acting as electron acceptor. Although the observation of these warped XB can be hindered by other noncovalent interactions, they may be ubiquitous in crystal structures of hypervalent species, where multiple XB can appear as secondary interactions on each halogen. We propose what can be considered the archetypal hypervalent halogen donor (a pincer type iodosodilactone) and a Lewis dot structure that includes the σ-holes.

Expanding the applicability of electrostatic potentials to the realm of transition states

Central to any reaction mechanism study, and sometimes a challenging job, is tracing a transition state in a reaction path. For the first time, electrostatic potentials (ESP) of the reactants were used as guiding tactics to predict whether there is a possibility of any transition state in a reaction surface. The main motive behind this strategy is to see whether the directionality nature of the transition state has something to do with the anisotropic natures of the ESP with their embedded directionalities. Strategically, some atmospherically important, but simple, reactions have been chosen for this study, which heretofore were believed to be barrierless. By carefully analysing the ESP maps of the reactants, regions of possible interactions were located. Using the bilinear interpolation of the 2D grids of the ESP surfaces, search co-ordinates were fine-tuned for a local gradient based approach for the search of a transition state. Out of the three reactions studied in this work, we were able to successfully locate transition states, for the first time, in two cases and the third one still proved to be barrierless. This gives a clear indication that though ESP maps can qualitatively predict the possibility of a transition state; it is not always true that there should definitely be a transition state, as some of the reaction surfaces may genuinely be barrierless. But, nevertheless this strategy definitely has credential to be tested for many more reactions, either new or already established, and may be applied to create the initial search co-ordinates for any well-established transition state search method. Moreover, we have observed that the analysis of the ESP maps of the reactants were very much useful in explaining the nature of interactions existing in those observed transition states and we hope the same can also be extended to any transition state in an electrostatically driven reaction potential energy surface.

Selective detection of non-aromatic pesticides via cyclodextrin-promoted fluorescence modulation

The introduction of pesticides to a cyclodextrin-fluorophore solution leads to highly specific fluorescence changes that are used for pesticide detection.