skip to main content

Title: Characterizing the interplay of Pauli repulsion, electrostatics, dispersion and charge transfer in halogen bonding with energy decomposition analysis
The halogen bond is a class of non-covalent interaction that has attracted considerable attention recently. A widespread theory for describing them is the σ-hole concept, which predicts that the strength of the interaction is proportional to the size of the σ-hole, a region of positive electrostatic potential opposite a σ bond. Previous work shows that in the case of CX 3 I, with X equal to F, Cl, Br, and I, the σ-hole trend is exactly opposite to the trend in binding energy with common electron pair donors. Using energy decomposition analysis (EDA) applied to a potential energy scan as well as the recent adiabatic EDA technique, we show that the observed trend is a result of charge transfer. Therefore a picture of the halogen bond that excludes charge transfer cannot be complete, and permanent and induced electrostatics do not always provide the dominant stabilizing contributions to halogen bonds. Overall, three universally attractive factors, polarization, dispersion and charge transfer, together with permanent electrostatics, which is usually attractive, drive halogen bonding, against Pauli repulsion.
Authors:
; ; ;
Award ID(s):
1665315 1363342
Publication Date:
NSF-PAR ID:
10066954
Journal Name:
Physical Chemistry Chemical Physics
Volume:
20
Issue:
2
Page Range or eLocation-ID:
905 to 915
ISSN:
1463-9076
Sponsoring Org:
National Science Foundation
More Like this
  1. Quantum chemistry in the form of density functional theory (DFT) calculations is a powerful numerical experiment for predicting intermolecular interaction energies. However, no chemical insight is gained in this way beyond predictions of observables. Energy decomposition analysis (EDA) can quantitatively bridge this gap by providing values for the chemical drivers of the interactions, such as permanent electrostatics, Pauli repulsion, dispersion, and charge transfer. These energetic contributions are identified by performing DFT calculations with constraints that disable components of the interaction. This review describes the second-generation version of the absolutely localized molecular orbital EDA (ALMO-EDA-II). The effects of different physical contributions on changes in observables such as structure and vibrational frequencies upon complex formation are characterized via the adiabatic EDA. Example applications include red- versus blue-shifting hydrogen bonds; the bonding and frequency shifts of CO, N 2 , and BF bound to a [Ru(II)(NH 3 ) 5 ] 2 + moiety; and the nature of the strongly bound complexes between pyridine and the benzene and naphthalene radical cations. Additionally, the use of ALMO-EDA-II to benchmark and guide the development of advanced force fields for molecular simulation is illustrated with the recent, very promising, MB-UCB potential.
  2. In order to explore how σ-hole potentials, as evaluated by molecular electrostatic potential (MEP) calculations, affect the ability of halogen atoms to engage in structure-directing intermolecular interactions, we synthesized four series of ethynyl halogen-substituted amide containing pyridines (activated targets); ( N -(pyridin-2-yl)benzamides (Bz-act-X), N -(pyridin-2-yl)picolinamides (2act-X), N -(pyridin-2-yl)nicotinamides (3act-X) and N -(pyridin-2-yl) isonicotinamides (4act-X), where X = Cl/Br/I. The molecules are deliberately equipped with three distinctly different halogen-bond acceptor sites, π, N(pyr), and OC, to determine binding site preferences of different halogen-bond donors. Crystallographic data for ten (out of a possible twelve) new compounds were thus analyzed and compared with data for the corresponding unactivated species. The calculated MEPs of all the halogen atoms were higher in the activated targets in comparison to the unactivated targets and were in the order of iodine ≈ chloroethynyl < bromoethynyl < iodoethynyl. This increased positive σ-hole potential led to a subsequent increase in propensity for halogen-bond formation. Two of the four chloroethynyl structures showed halogen bonding, and all three of the structurally characterized bromoethynyl species engaged in halogen bonding. The analogous unactived species showed no halogen bonds. Each chloroethynyl donor selected a π-cloud as acceptor and the bromoethynyl halogen-bond donors opted for eithermore »π or N(pyr) sites, whereas all halogen bonds involving an iodoethynyl halogen-bond donor (including both polymorphs of Bz-act-I ) engaged exclusively with a N(pyr) acceptor site.« less
  3. Molecules of the type XYT = Ch (T = C, Si, Ge; Ch = S, Se; X,Y = H, CH3, Cl, Br, I) contain a σ-hole along the T = Ch bond extension. This hole can engage with the N lone pair of NCH and NCCH3 so as to form a chalcogen bond. In the case of T = C, these bonds are rather weak, less than 3 kcal/mol, and are slightly weakened in acetone or water. They owe their stability to attractive electrostatic energy, supplemented by dispersion, and a much smaller polarization term. Immersion in solvent reverses the electrostatic interaction to repulsive, while amplifying the polarization energy. The σ-holes are smaller for T = Si and Ge, even negative in many cases. These Lewis acids can nonetheless engage in a weak chalcogen bond. This bond owes its stability to dispersion in the gas phase, but it is polarization that dominates in solution.
  4. To facilitate computational investigation of intermolecular interactions in the solution phase, we report the development of ALMO-EDA(solv), a scheme that allows the application of continuum solvent models within the framework of energy decomposition analysis (EDA) based on absolutely localized molecular orbitals (ALMOs). In this scheme, all the quantum mechanical states involved in the variational EDA procedure are computed with the presence of solvent environment so that solvation effects are incorporated in the evaluation of all its energy components. After validation on several model complexes, we employ ALMO-EDA(solv) to investigate substituent effects on two classes of complexes that are related to molecular CO 2 reduction catalysis. For [FeTPP(CO 2 -κC)] 2− (TPP = tetraphenylporphyrin), we reveal that two ortho substituents which yield most favorable CO 2 binding, –N(CH 3 ) 3 + (TMA) and –OH, stabilize the complex via through-structure and through-space mechanisms, respectively. The coulombic interaction between the positively charged TMA group and activated CO 2 is found to be largely attenuated by the polar solvent. Furthermore, we also provide computational support for the design strategy of utilizing bulky, flexible ligands to stabilize activated CO 2 via long-range Coulomb interactions, which creates biomimetic solvent-inaccessible “pockets” in that electrostatics is unscreened.more »For the reactant and product complexes associated with the electron transfer from the p -terphenyl radical anion to CO 2 , we demonstrate that the double terminal substitution of p -terphenyl by electron-withdrawing groups considerably strengthens the binding in the product state while moderately weakens that in the reactant state, which are both dominated by the substituent tuning of the electrostatics component. These applications illustrate that this new extension of ALMO-EDA provides a valuable means to unravel the nature of intermolecular interactions and quantify their impacts on chemical reactivity in solution.« less
  5. As known, small HCl–water nanoclusters display a particular dissociation behaviour, whereby at least four water molecules are required for the ionic dissociation of HCl. In this work, we examine how intermolecular interactions promote the ionic dissociation of such nanoclusters. To this end, a set of 45 HCl–water nanoclusters with up to four water molecules is introduced. Energy decomposition analysis based on absolutely localized molecular orbitals (ALMO-EDA) is employed in order to study the importance of frozen interaction, dispersion, polarization, and charge-transfer for the dissociation. The vertical ALMO-EDA scheme is applied to HCl–water clusters along a proton-transfer coordinate varying the amount of spectator water molecules. The corresponding ALMO-EDA results show a clear preference for the dissociated cluster only in the presence of four water molecules. Our analysis of adiabatic ALMO-EDA results reveals a push–pull mechanism for the destabilization of the HCl bond based on the synergy between forward and backward charge-transfer.