Interpreting the variations in the kinetic and potential energies in the formation of a covalent bond.

We address the long-standing controversy as to the physical origin of covalent bonding, whether it involves a lowering of the potential energy or a lowering of the kinetic energy. We conclude that both of these do occur and contribute to the formation of the bond. The analysis is in terms of the virial theorem and the variations in the potential energy and the kinetic energy as the atoms approach each other. At large separations, the change in kinetic energy relative to the separated atoms is negative and stabilizing, while the corresponding potential energy change is positive and destabilizing. However, as the atoms approach their equilibrium separation, these rapidly reverse; the kinetic energy increases and the potential energy decreases, so that at equilibrium the net kinetic energy is positive and the net potential energy negative. At equilibrium, the bonding is due solely to the potential energy and is electrostatic.

Can Counter-Intuitive Halogen Bonding Be Coulombic?

We use the term "counter-intuitive" to describe an intermolecular interaction in which the electrostatic potentials of the interacting regions of the ground-state molecules have the same sign, both positive or both negative. In the present work, we consider counter-intuitive halogen bonding with nitrogen bases, in which both the halogen σ-hole and the nitrogen lone pair have negative potentials on their molecular surfaces. We show that these interactions can be treated as Coulombic despite the apparent repulsion between the ground-state molecules, provided that both electrostatics and polarization are explicitly taken into account. We demonstrate first that the energies of 20 counter-intuitive interactions with four nitrogen bases can be expressed very well in terms of just two molecular properties: the electrostatic potential of the halogen σ-hole and the average polarizability of the nitrogen base. Then we show that the same two properties can also represent the energies of an expanded data base that includes the 20 counter-intuitive plus an additional 20 weak and moderately-strong intuitive halogen bonding interactions (in which the σ-hole potentials are now positive).

The π-hole revisited.

It follows from the Schrödinger equation that the forces operating within molecules and molecular complexes are Coulombic, which necessarily entails both electrostatics and polarization. A common and important class of molecular complexes is due to π-holes. These are molecular regions of low electronic density that are perpendicular to planar portions of the molecular frameworks. π-Holes often have positive electrostatic potentials associated with them, which result in mutually polarizing attractive forces with negative sites such as lone pairs, π electrons or anions. In many molecules, π-holes correspond to a flattening of the electronic density surface but in benzene derivatives and in polyazines the π-holes are craters above and below the rings. The interaction energies of π-hole complexes can be expressed quite well in terms of regression relationships that account for both the electrostatics and the polarization. There is a marked gradation in the interaction energies, from quite weak (about -2 kcal mol-1) to relatively strong (about -40 kcal mol-1). Gradations are also evident in the ratios of the intermolecular separations to the sums of the respective van der Waals radii and in the gradual transition of the π-hole atoms from trigonal to quasi-tetrahedral configurations. These trends are consistent with the concept that chemical interactions form a continuum, from very weak to very strong.

The Neglected Nuclei.

Since the nuclei in a molecule are treated as stationary, it is perhaps natural that interpretations of molecular properties and reactivity have focused primarily upon the electronic density distribution. The role of the nuclei has generally received little explicit consideration. Our objective has been to at least partially redress this imbalance in emphasis. We discuss a number of examples in which the nuclei play the determining role with respect to molecular properties and reactive behavior. It follows that conventional interpretations based solely upon electronic densities and donating or withdrawing tendencies should be made with caution.

Electrostatics and Polarization in σ- and π-Hole Noncovalent Interactions: An Overview.

The energetics of σ- and π-hole interactions can be described very well in terms of electrostatics and polarization, consistent with their Coulombic natures. When both of these components are taken into account, very good correlations with quantum-chemically computed interaction energies are obtained. If polarization is only minor, as when the interactions are quite weak, then electrostatics can suffice, as represented by the most positive electrostatic potential associated with the σ- or π-hole. For stronger interactions, the combination of electrostatics plus polarization is very effective even for interaction energies considerably greater in magnitude than what is normally considered noncovalent bonding. Several procedures for treating polarization are summarized, including the use of point charges and the direct inclusion of electric fields.

Explicit Inclusion of Polarizing Electric Fields in σ- and π-Hole Interactions.

The interactions between a wide variety of molecules having σ-holes or π-holes and several nitrogen bases have been analyzed computationally. The σ- and π-hole atoms span groups III-VII of the periodic table. The interaction energies range from quite weak, typical of non-covalent bonding, to unusually strong: from -4.6 to -22.0 kcal/mol for σ-hole bonding and from -4.0 to -42.4 kcal/mol for π-hole bonding. The markedly greater strengths of some bonds does not imply that any new factors or types of bonding are involved; they simply reflect higher degrees of the polarization that is part of any Coulombic interaction. To explain the stronger bonding, this polarization must be explicitly taken into account. We show that the interaction energies can be related quite well to (a) the maximum positive electrostatic potentials associated with the σ- or π-holes on their molecular surfaces, (b) the polarizabilities of the nitrogen bases, and especially (c) the polarizing electric fields of the σ- or π-hole molecules at the positions of the nitrogens.

σ-holes and π-holes: Similarities and differences.

σ-Holes and π-holes are regions of molecules with electronic densities lower than their surroundings. There are often positive electrostatic potentials associated with them. Through these potentials, the molecule can interact attractively with negative sites, such as lone pairs, π electrons, and anions. Such noncovalent interactions, "σ-hole bonding" and "π-hole bonding," are increasingly recognized as being important in a number of different areas. In this article, we discuss and compare the natures and characteristics of σ-holes and π-holes, and factors that influence the strengths and locations of the resulting electrostatic potentials. © 2017 Wiley Periodicals, Inc.

The σ-Hole Coulombic Interpretation of Trihalide Anion Formation.

It is shown that the interactions of dihalogen molecules XY with halide anions Z- to form trihalide anions (XYZ)- can be satisfactorily described as Coulombic, involving the σ-holes on the atoms Y, but only if polarization is taken into account. We have approximated the polarizing effect of the halide anion Z- by means of a unit negative point charge. The CCSD/aug-cc-pVTZ computed interaction energies ΔE correlate well with the most positive electrostatic potentials associated with the induced σ-holes over a ΔE range of -12 to -63 kcal mol-1 . The (XYZ)- anions are more stable when the central atom is the largest, as has been observed, because the central atom is then the most polarizable, making the electrostatic potential associated with its σ-hole more positive.

A perspective on quantum mechanics and chemical concepts in describing noncovalent interactions.

Since quantum mechanical calculations do not typically lend themselves to chemical interpretation, analyses of bonding interactions depend largely upon models (the octet rule, resonance theory, charge transfer, etc.). This sometimes leads to a blurring of the distinction between mathematical modelling and physical reality. The issue of polarization vs. charge transfer is an example; energy decomposition analysis is another. The Hellmann-Feynman theorem at least partially bridges the gap between quantum mechanics and conceptual chemistry. It proceeds rigorously from the Schrödinger equation to demonstrating that the forces exerted upon the nuclei in molecules, complexes, etc., are entirely classically coulombic attractions with the electrons and repulsions with the other nuclei. In this paper, we discuss these issues in the context of noncovalent interactions. These can be fully explained in coulombic terms, electrostatics and polarization (which include electronic correlation and dispersion).

Close contacts and noncovalent interactions in crystals.

Close contacts, defined as interatomic separations less than the sum of the respective van der Waals radii, are commonly invoked to identify attractive nonbonded interactions in crystal lattices. While this is often effective, it can also be misleading because (a) there are significant uncertainties associated with van der Waals radii, and (b) it may not be valid to attribute the interactions solely to specific pairs of atoms. The interactions within crystal lattices are Coulombic, and the strongest positive and/or negative regions do not always correspond to the positions of atoms; they are sometimes located between atoms. Examples of both types are given and discussed, focusing in particular upon σ-hole interactions.

The σ-hole revisited.

A covalently-bonded atom typically has a region of lower electronic density, a "σ-hole," on the side of the atom opposite to the bond, along its extension. There is frequently a positive electrostatic potential associated with this region, through which the atom can interact attractively but noncovalently with negative sites. This positive potential reflects not only the lower electronic density of the σ-hole but also contributions from other portions of the molecule. These can significantly influence both the value and also the angular position of the positive potential, causing it to deviate from the extension of the covalent bond. We have surveyed these effects, and their consequences for the directionalities of subsequent noncovalent intermolecular interactions, for atoms of Groups IV-VII. The overall trends are that larger deviations of the positive potential result in less linear intermolecular interactions, while smaller deviations lead to more linear interactions. We find that the deviations of the positive potentials and the nonlinearities of the noncovalent interactions tend to be greatest for atoms of Groups V and VI. We also present arguments supporting the use of the 0.001 a.u. contour of the electronic density as the molecular surface on which to compute the electrostatic potential.

Hydrogen Bonding between Metal-Ion Complexes and Noncoordinated Water: Electrostatic Potentials and Interaction Energies.

The hydrogen bonding of noncoordinated water molecules to each other and to water molecules that are coordinated to metal-ion complexes has been investigated by means of a search of the Cambridge Structural Database (CSD) and through quantum chemical calculations. Tetrahedral and octahedral complexes that were both charged and neutral were studied. A general conclusion is that hydrogen bonds between noncoordinated water and coordinated water are much stronger than those between noncoordinated waters, whereas hydrogen bonds of water molecule in tetrahedral complexes are stronger than in octahedral complexes. We examined the possibility of correlating the computed interaction energies with the most positive electrostatic potentials on the interacting hydrogen atoms prior to interaction and obtained very good correlation. This study illustrates the fact that electrostatic potentials computed for ground-state molecules, prior to interaction, can provide considerable insight into the interactions.

σ-Hole bonding: a physical interpretation.

The anisotropic electronic densities of covalently-bonded Group IV-VII atoms frequently give rise to regions of positive electrostatic potential on the extensions of covalent bonds to these atoms. Through such positive "σ-holes," the atoms can interact attractively and highly directionally with negative sites such as the lone pairs of Lewis bases, anions, π electrons, etc. In the case of Group VII this is called "halogen bonding." Hydrogen bonding can be viewed as a less directional subset of σ-hole interactions. Since positive σ-holes often exist in conjunction with regions of negative potential, the atoms can also interact favorably with positive sites. In accordance with the Hellmann-Feynman theorem, all of these interactions are purely Coulombic in nature (which encompasses polarization and dispersion). The strength of σ-hole bonding increases with the magnitudes of the potentials of the positive σ-hole and the negative site; their polarizabilities must sometimes also be taken explicitly into account.

Complementarity of reaction force and electron localization function analyses of asynchronicity in bond formation in Diels-Alder reactions.

We have computationally compared three Diels-Alder cycloadditions involving cyclopentadiene and substituted ethylenes; one of the reactions is synchronous, while the others are slightly or highly asynchronous. Synchronicity and weak asynchronicity are characterized by the reaction force constant κ(ξ) having just a single minimum in the transition region along the intrinsic reaction coordinate ξ, while for high asynchronicity κ(ξ) has a negative maximum with minima on both sides. The electron localization function (ELF) shows that the features of κ(ξ) can be directly related to the formation of the new C-C bonds between the diene and the dienophile. There is thus a striking complementarity between κ(ξ) and ELF; κ(ξ) identifies the key points along ξ and ELF describes what is happening at those points.

Halogen bonding: an interim discussion.

Halogen bonding is a noncovalent interaction that is receiving rapidly increasing attention because of its significance in biological systems and its importance in the design of new materials in a variety of areas, for example, electronics, nonlinear optical activity, and pharmaceuticals. The interactions can be understood in terms of electrostatics/polarization and dispersion; they involve a region of positive electrostatic potential on a covalently bonded halogen and a negative site, such as the lone pair of a Lewis base. The positive potential, labeled a σ hole, is on the extension of the covalent bond to the halogen, which accounts for the characteristic near-linearity of halogen bonding. In many instances, the lateral sides of the halogen have negative electrostatic potentials, allowing it to also interact favorably with positive sites. In this discussion, after looking at some of the experimental observations of halogen bonding, we address the origins of σ holes, the factors that govern the magnitudes of their electrostatic potentials, and the properties of the resulting complexes with negative sites. The relationship of halogen and hydrogen bonding is examined. We also point out that σ-hole interactions are not limited to halogens, but can also involve covalently bonded atoms of Groups IV-VI. Examples of applications in biological/medicinal chemistry and in crystal engineering are mentioned, taking note that halogen bonding can be "tuned" to fit various requirements, that is, strength of interaction, steric factors, and so forth.

Halogen bonding and other σ-hole interactions: a perspective.

A σ-hole bond is a noncovalent interaction between a covalently-bonded atom of Groups IV-VII and a negative site, e.g. a lone pair of a Lewis base or an anion. It involves a region of positive electrostatic potential, labeled a σ-hole, on the extension of one of the covalent bonds to the atom. The σ-hole is due to the anisotropy of the atom's charge distribution. Halogen bonding is a subset of σ-hole interactions. Their features and properties can be fully explained in terms of electrostatics and polarization plus dispersion. The strengths of the interactions generally correlate well with the magnitudes of the positive and negative electrostatic potentials of the σ-hole and the negative site. In certain instances, however, polarizabilities must be taken into account explicitly, as the polarization of the negative site reaches a level that can be viewed as a degree of dative sharing (coordinate covalence). In the gas phase, σ-hole interactions with neutral bases are often thermodynamically unfavorable due to the relatively large entropy loss upon complex formation.

The reaction force constant as an indicator of synchronicity/nonsynchronicity in [4+2] cycloaddition processes.

A variety of experimental and computational analyses support the concept that a chemical reaction has a transition region, in which the system changes from distorted states of the reactants to distorted states of the products. The boundaries of this region along the intrinsic reaction coordinate ξ, which includes the traditional transition state, are defined unambiguously by the minimum and maximum of the reaction force F(ξ), which is the negative gradient of the potential energy V(ξ). The transition region is characterized by the reaction force constant κ(ξ), the second derivative of V(ξ), being negative throughout. It has recently been demonstrated that the profile of κ(ξ) in the transition region is a sensitive indicator of the degree of synchronicity of a concerted reaction: a single κ(ξ) minimum is associated with full or nearly full synchronicity, while a κ(ξ) maximum (negative) between two minima is a sign of considerable nonsynchronicity, i.e. a two-stage concerted process. We have now applied reaction force analysis to the Diels-Alder cycloadditions of the various cyanoethylenes to cyclopentadiene. We examine the relative energy requirements of the structurally- and electronically-intensive phases of the activation processes. We demonstrate that the variation of κ(ξ) in the transition region is again indicative of the level of synchronicity. The fully synchronous cycloadditions are those in which the cyanoethylenes are symmetrically substituted. Unsymmetrical substitution leads to minor nonsynchronicity for monocyanoethylene but much more - i.e. two stages - for 1,1-dicyano- and 1,1,2-tricyanoethylene. We also show that the κ(ξ) tend to become less negative as the activation energies decrease.

The reaction force constant: an indicator of the synchronicity in double proton transfer reactions.

Earlier work, both experimental and computational, has drawn attention to the transition region in a chemical reaction, which includes the traditional transition state but extends along the intrinsic reaction coordinate ξ from perturbed forms of the reactants to perturbed forms of the products. The boundaries of this region are defined by the reaction force F(ξ), which is the negative gradient of the potential energy V(ξ) of the system along ξ. The reaction force constant κ(ξ), the second derivative of V(ξ), is negative throughout the transition region. We have now demonstrated, for a series of twelve double proton transfer processes, that the profile of κ(ξ) in the transition region is an indicator of the synchronicity of the two proton migrations in each case. When they are fully or nearly fully synchronous, κ(ξ) has a single minimum in the transition region. When the migrations are considerably nonsynchronous, κ(ξ) has two minima separated by a local maximum. Such an assessment of the degree of synchronicity cannot readily be made from an examination of the transition state alone, nor it is easily detected in the profiles of V(ξ) and F(ξ).