Advances and Prospects in Understanding London Dispersion Interactions in Molecular Chemistry

London dispersion (LD) interactions are the main contribution of the attractive part of the van der Waals potential. Even though LD effects are the driving force for molecular aggregation and recognition, the role of these omnipresent interactions in structure and reactivity had been largely underappreciated over decades. However, in the recent years considerable efforts have been made to thoroughly study LD interactions and their potential as a chemical design element for structures and catalysis. This was made possible through a fruitful interplay of theory and experiment. This review highlights recent results and advances in utilizing LD interactions as a structural motif to understand and utilize intra- and intermolecularly LD-stabilized systems. Additionally, we focus on the quantification of LD interactions and their fundamental role in chemical reactions.

Conformational energies of biomolecules in solution: Extending the MPCONF196 benchmark with explicit water molecules

A prerequisite for the computational prediction of molecular properties like conformational energies of biomolecules is a reliable, robust, and computationally affordable method usually selected according to its performance for relevant benchmark sets. However, most of these sets comprise molecules in the gas phase and do not cover interactions with a solvent, even though biomolecules typically occur in aqueous solution. To address this issue, we introduce a with explicit water molecules solvated version of a gas-phase benchmark set containing 196 conformers of 13 peptides and other relevant macrocycles, namely MPCONF196 [J. Řezáč et al., JCTC 2018, 14, 1254-1266], and provide very accurate PNO-LCCSD(T)-F12b/AVQZ' reference values. The novel solvMPCONF196 benchmark set features two additional challenges beyond the description of conformers in the gas phase: conformer-water and water-water interactions. The overall best performing method for this set is the double hybrid revDSDPBEP86-D4/def2-QZVPP yielding conformational energies of almost coupled cluster quality. Furthermore, some (meta-)GGAs and hybrid functionals like B97M-V and ω B97M-D with a large basis set reproduce the coupled cluster reference with an MAD below 1 kcal mol- 1  . If more efficient methods are required, the composite DFT-method r2 SCAN-3c (MAD of 1.2 kcal mol- 1 ) is a good alternative, and when conformational energies of polypeptides or macrocycles with more than 500-1000 atoms are in the focus, the semi-empirical GFN2-xTB or the MMFF94 force field (for very large systems) are recommended.

Affinity of Telluronium Chalcogen Bond Donors for Lewis Bases in Solution: A Critical Experimental-Theoretical Joint Study

Telluronium salts [Ar2 MeTe]X were synthesized, and their Lewis acidic properties towards a number of Lewis bases were addressed in solution by physical and theoretical means. Structural X-ray diffraction analysis of 21 different salts revealed the electrophilicity of the Te centers in their interactions with anions. Telluroniums' propensity to form Lewis pairs was investigated with OPPh3 . Diffusion-ordered NMR spectroscopy suggested that telluroniums can bind up to three OPPh3 molecules. Isotherm titration calorimetry showed that the related heats of association in 1,2-dichloroethane depend on the electronic properties of the substituents of the aryl moiety and on the nature of the counterion. The enthalpies of first association of OPPh3 span -0.5 to -5 kcal mol-1 . Study of the affinity of telluroniums for OPPh3 by state-of-the-art DFT and ab-initio methods revealed the dominant Coulombic and dispersion interactions as well as an entropic effect favoring association in solution. Intermolecular orbital interactions between [Ar2 MeTe]+ cations and OPPh3 are deemed insufficient on their own to ensure the cohesion of [Ar2 MeTe ⋅ Bn ]+ complexes in solution (B=Lewis base). Comparison of Grimme's and Tkatchenko's DFT-D4/MBD-vdW thermodynamics of formation of higher [Ar2 MeTe ⋅ Bn ]+ complexes revealed significant molecular size-dependent divergence of the two methodologies, with MBD yielding better agreement with experiment.

Boosting Energy-Transfer Processes via Dispersion Interactions

The generation of open-shell intermediates under mild conditions has opened broad synthetic opportunities during this century. However, these reactive species often require a case specific and tailored tuning of experimental parameters in order to efficiently convert substrates into products. We report a general approach that can overcome these ubiquitous limitations for several visible-light promoted energy-transfer processes. The use of either naphthalene (5-20 equiv.) or simple binaphthyl derivatives (10-30 mol %) greatly increases their efficiency, giving rise to a new strategy for catalysis. The trend is consistent among different media, photocatalysts, light sources and substrates, allowing one to improve existing methods, to more easily optimize conditions for new ones, and, moreover, to disclose otherwise inaccessible reaction pathways.

Contributions of London Dispersion Forces to Solution-Phase Association Processes

ConspectusDespite their ubiquity and early discovery, London dispersion forces are often overlooked. This is due, in part, to the difficulty in assessing their contributions to molecular and polymeric structure, stability, properties, and reactivities. However, recent advances in modeling have revealed that dispersion interactions play an important role in many important chemical and biological processes. Experimental confirmation of their impact in solution has been challenging, leading to controversies about their relative importance.In the course of studying noncovalent interactions using molecular devices, our understanding and appreciation for the importance of dispersion interactions have evolved. This Account follows this intellectual journey by using examples from the literature. The goals are twofold: to describe recent advances in understanding the interaction and to provide guidance to researchers studying weak noncovalent interactions. However, first, the experimental methods for measuring the effects of dispersion interactions and the strategies for isolating their influence are described. These include the design of molecular devices to measure these weak noncovalent interactions and the strategies to disentangle the solvation, solvophobic, and dispersion components of the resulting equilibria.The literature examples are organized around five fundamental questions. (1) Do dispersion interactions have a measurable effect on solution equilibria? (2) To what extent do solvents attenuate or compensate for dispersion interactions? (3) To what extent do the solvation and solvophobic terms influence the dispersion equilibria? (4) Can we predict whether a system will form attractive dispersion or repulsive steric interactions? (5) Can the dispersion term be isolated and interrogated? We were often surprised by the answers to these questions. In each case, we describe how the systems were designed to address these questions and discuss possible interpretations of the results.While dispersion interactions in solution were weak (usually <1 kcal/mol), their influence on complexation and conformational equilibria can be observed and measured. This underscores the significance of these interactions in molecular recognition, coordination chemistry, reaction design, and catalysis. The solvent components of the dispersion equilibria can also be significant. Therefore, the isolation of the dispersion contributions from the solvation and solvophobic effects represents an ongoing challenge. The experimental studies also provide important benchmarks and offer valuable insights to help refine the next generation of computational solvent models.

Dispersion Interactions in Condensed Phases and inside Molecular Containers

ConspectusThe past decade has seen significant progress in the understanding and appreciation of the importance of London dispersion interactions (LDIs) in supramolecular systems and solutions. The Slater-Kirkwood formula relates LDIs to the molecular polarizabilities of the two interacting molecular species (α) and their interaction distance (a dependence of R-6). When advancing arguments related to intermolecular interactions, it is frequently assumed that molecules with larger molecular polarizabilities are more amenable to larger LDIs. However, arguments related to molecular polarizabilities are not always transferable to the condensed phase. In fact, the underlying bulk and molecular polarizabilities of common solvents show opposing trends. The intuitive concept that aromatic molecules are more polarizable than saturated hydrocarbons and that perfluorinated molecules are less polarizable than saturated hydrocarbons applies to the condensed phase only. When treating association phenomena in solution, where LDIs are generally very attenuated, the use of bulk polarizabilities is recommended, which are experimentally accessible through either refractive index measurements or suitable solvatochromic probes. Such probes can also be used to assess polarizabilities inside molecular container compounds, such as cucurbit[n]urils (CBn), cyclodextrins, calixarenes, and hemicarcerands. These macrocyclic cavities can have extreme microenvironments. For example, the inner concave phase of CB7 has been shown to be weakly polarizable, falling in between the gas phase and perfluorohexane; those of β-cyclodextrin and p-sulfonatocalix[4]arene have been found to be similarly polarizable as water and alkanes, respectively, and the inside of hemicarcerands displays a very large bulk polarizability, exceeding that of diiodomethane. CBn compounds are privileged molecular container compounds, which we exemplify in this Account through case studies. (1) CBn macrocycles are prime water-soluble receptors for hydrocarbons, allowing for the reduction of the binding free energies to two components: the hydrophobic effect and dispersion interactions. To understand hydrocarbon binding, we initiated the HYDROPHOBE challenge, which revealed the shortcomings of both quantum-chemical and molecular dynamics approaches. (2) The smallest CBn receptor, CB5, is uniquely suited to bind the entire noble gas series, where hydrophobic effects and dispersion interactions operate in opposite directions. CB5 was revaled to be a unique synthetic receptor for noble gases, with the dominant driving force being the recovery of the cavitation energies for the hydration of noble gases in aqueous solution. Computational methods that encounter challenges in predicting hydrocarbon affinities and trends for CB6 and CB7 perform well for noble gases binding to CB5. (3) The larger homologue, CB8, allows one to set up intermolecular interaction chambers by the encapsulation of a (first) aromatic guest, thereby tuning LDIs inside the receptor cavity. In this manner, CB8 can be modulated to preferentially bind unsaturated and aromatic rather than saturated hydrocarbons, while the unmodified cavities of the smaller macrocycles CB6 and CB7 show selective binding of saturated hydrocarbons. (4) The (charged) host-guest complexes of CBn hosts are sufficiently stable in the gas phase, allowing for the study of the influence of LDIs on inner-phase chemical reactions. These studies are particularly interesting for the theoretical analysis of isolated host-guest LDIs, as experimental and computational data are directly comparable in the gas phase due to the absence of the solvation effect.

Dispersion Energy-Stabilized Boron and Phosphorus Lewis Pairs

An isostructural series of boron/phosphorus Lewis pairs was systematically investigated. The association constants of the Lewis pairs were determined at variable temperatures, enabling the extraction of thermodynamic parameters. The stabilization of the Lewis adduct increased with increasing size of the dispersion energy donor groups, although the donor and acceptor properties of the Lewis pairs remained largely unchanged. This data was utilized to challenge state-of-the-art quantum chemical methods, which finally led to an enhanced workflow for the determination of thermochemical properties of weakly bound Lewis pairs within an accuracy of 0.6 to 1.0 kcal mol-1 for computed association free energies.

Silyl Groups Are Strong Dispersion Energy Donors

We present an experimental and computational study to investigate noncovalent interactions between silyl groups that are often employed as "innocent" protecting groups. We chose an extended cyclooctatetraene (COT)-based molecular balance comprising unfolded (1,4-disubstituted) and folded (1,6-disubstituted) valance bond isomers that typically display remote and close silyl group contacts, respectively. The thermodynamic equilibria were determined using nuclear magnetic resonance measurements. Additionally, we utilized Boltzmann weighted symmetry-adapted perturbation theory (SAPT) at the sSAPT0/aug-cc-pVDZ level of theory to dissect and quantify noncovalent interactions. Apart from the extremely bulky tris(trimethylsilyl)silyl "supersilyl" group, there is a preference for the folded 1,6-COT valence isomer, with London dispersion interactions being the main stabilizing factor. This makes silyl groups excellent dispersion energy donors, a finding that needs to be taken into account in synthesis planning.

Model Compounds for Intermediates and Transition States in Sonogashira and Negishi Coupling: d8-d10 Bonds in Large Heterobimetallic Complexes Are Weaker than Computational Chemistry Predicts

We report an experimental study, with accompanying DFT calculations, on a series of heterobimetallic complexes with Pd(II) and Cu(I), Ag(I), or Au(I). The isolable pincer complexes are models for the intermediates and transition state for the transmetalation step in Sonogashira and Negishi coupling reactions, among which, according to the DFT calculations, only the transition state has the two metal centers within bonding distance. Furthermore, we report a substituted version of an analogous heterobimetallic complex in which a competing dissociation sets an upper-bound on the strength of the d⁸-d¹⁰ metal-metal bond. Analysis of the structures in the solid state and in solution, and the competitive dissociation experiment in the gas phase, indicate that the dispersion-corrected DFT method used in the study appears to overestimate the strength of the metal-metal interaction, thus distorting the shape of the computed potential energy surface systematically for transmetalation.

Increasing Complexity in a Conformer Space Step-by-Step: Weighing London Dispersion against Cation-π Interactions

We report an evaluation of the importance of London dispersion in moderately large (up to 36 heavy atoms) organic molecules by means of a molecular torsion balance whose conformations "weigh" one interaction against another in the absence of solvents. The experimental study, with gas-phase cryogenic ion vibrational predissociation (CIVP) spectroscopy, solid-state Fourier transfer infrared (FT-IR), and single-crystal X-ray crystallography, is accompanied by density functional theory calculations, including an extensive search and analysis of accessible conformations. We begin with the unsubstituted molecular torsion balance, and then step up the complexity systematically by adding alkyl groups incrementally as dispersion energy donors (DEDs) to achieve a degree of chemical complexity comparable to what is typically found in transition states for many regio- and stereoselective reactions in organic and organometallic chemistry. We find clear evidence for the small attractive contribution by DEDs, as had been reported in other studies, but we also find that small individual contributions by London dispersion, when they operate in opposition to other weak noncovalent interactions, produce composite effects on the structure that are difficult to predict intuitively, or by modern quantum chemical calculations. The experimentally observed structures, together with a reasonable value for a reference cation-π interaction, indicate that the pairwise interaction between two tert-butyl groups, in the best case, is modest. Moreover, the visualization of the conformational space, and comparison to spectroscopic indicators of the structure, as one steps up the complexity of the manifold of noncovalent interactions, makes clear that in silico predictive ability for the structure of moderately large, flexible, organic molecules falters sooner than one might have expected.

Beyond Steric Crowding: Dispersion Energy Donor Effects in Large Hydrocarbon Ligands

Interactions between sterically crowded hydrocarbon-substituted ligands are widely considered to be repulsive because of the intrusion of the electron clouds of the ligand atoms into each other's space, which results in Pauli repulsion. Nonetheless, there is another interaction between the ligands which is less widely publicized but is always present. This is the London dispersion (LD) interaction which can occur between atoms or molecules in which dipoles can be induced instantaneously, for example, between the H atoms from the ligand C-H groups.These LD interactions are always attractive, but their effects are not as widely recognized as those of the Pauli repulsion despite their central role in the formation of condensed matter. Their relatively poor recognition is probably due to the relative weakness (ca. 1 kcal mol^-1) of individual H···H interactions owing to their especially strong distance dependence. In contrast, where there are numerous H···H interactions, a collective LD energy equaling several tens of kcal mol^-1 may ensue. As a result, in some molecules the latent importance of the LD attraction energies emerges and assumes a prominence that can overshadow the Pauli effects (e.g., in the stabilization of high-oxidation-state transition-metal alkyls, inducing disproportionation reactions, or in the stabilization of otherwise unstable bonds).Despite being known for over a century, the accurate quantification of individual H···H LD effects in molecular species is a relatively recent phenomenon and at present is based mainly on modified DFT calculations. A few leading reviews summarized these earlier studies of the C-H···H-C LD interactions in organic molecules, and their effects on the structures and stabilities were described. LD effects in sterically crowded inorganic and organometallic molecules have been recognized.The author's interest in these LD effects arose fortuitously over a decade ago during research on sterically crowded heavier main-group element carbene analogues and two-coordinate, open-shell (d¹-d⁹) transition-metal complexes where counterintuitive steric effects were observed. More detailed explanations of these effects were provided by dispersion-corrected DFT calculations in collaboration with the groups of Tuononen and Nagase (see below).This Account describes our development of these initial results for other inorganic molecular classes. More recently, the work has led us to move to the planned inclusion of dispersion effects in ligands to stabilize new molecular types with theoretical input from the groups of Vasko and Grimme (see below). Our approach sought to use what Grimme has described as dispersion effect donor (DED) groups (i.e., spatially close-lying, densely packed substituents either as ligands (e.g., -C₆H₂-2,4,6-Cy₃, Cy = cyclohexyl) or as parts of ligands (e.g., a Cy substituent) that produce relatively large dispersion energies to stabilize these new compounds.We predict that the future design of sterically crowding hydrocarbon ligands will include the consideration and incorporation of LD effects as a standard methodology for directed use in the attainment of new synthetic targets.

Assessing the persistence of chalcogen bonds in solution with neural network potentials

Non-covalent bonding patterns are commonly harvested as a design principle in the field of catalysis, supramolecular chemistry and functional materials to name a few. Yet, their computational description generally neglects finite temperature and environment effects, which promote competing interactions and alter their static gas-phase properties. Recently, neural network potentials (NNPs) trained on Density Functional Theory (DFT) data have become increasingly popular to simulate molecular phenomena in condensed phase with an accuracy comparable to ab initio methods. To date, most applications have centered on solid-state materials or fairly simple molecules made of a limited number of elements. Herein, we focus on the persistence and strength of chalcogen bonds involving a benzotelluradiazole in condensed phase. While the tellurium-containing heteroaromatic molecules are known to exhibit pronounced interactions with anions and lone pairs of different atoms, the relevance of competing intermolecular interactions, notably with the solvent, is complicated to monitor experimentally but also challenging to model at an accurate electronic structure level. Here, we train direct and baselined NNPs to reproduce hybrid DFT energies and forces in order to identify what are the most prevalent non-covalent interactions occurring in a solute-Cl$^-$-THF mixture. The simulations in explicit solvent highlight competition with chalcogen bonds formed with the solvent and the short-range directionality of the interaction with direct consequences for the molecular properties in the solution. The comparison with other potentials (e.g., AMOEBA, direct NNP and continuum solvent model) also demonstrates that baselined NNPs offer a reliable picture of the non-covalent interaction interplay occurring in solution.

An Incremental System To Predict the Effect of Different London Dispersion Donors in All‐ meta ‐Substituted Azobenzenes

Predictive models based on incremental systems exist for many chemical phenomena, thus allowing easy estimates. Despite their low magnitude in isolated systems London dispersion interactions are ubiquitous in manifold situations ranging from solvation to catalysis or in biological systems. Based on our azobenzene system, we systematically determined the London dispersion donor strength of the alkyl substituents Me, Et, iPr up to tBu. Based on this data, we were able to implement an incremental system for London dispersion for the azobenzene scheme. We propose an equation that allows the prediction of the effect of change of substituents on London dispersion interactions in azobenzenes, which has to be validated in similar molecular arrangements in the future.

Actual Symmetry of Symmetric Molecular Adducts in the Gas Phase, Solution and in the Solid State

This review discusses molecular adducts, whose composition allows a symmetric structure. Such adducts are popular model systems, as they are useful for analyzing the effect of structure on the property selected for study since they allow one to reduce the number of parameters. The main objectives of this discussion are to evaluate the influence of the surroundings on the symmetry of these adducts, steric hindrances within the adducts, competition between different noncovalent interactions responsible for stabilizing the adducts, and experimental methods that can be used to study the symmetry at different time scales. This review considers the following central binding units: hydrogen (proton), halogen (anion), metal (cation), water (hydrogen peroxide).

Intramolecular London Dispersion Interactions Do Not Cancel in Solution

We present a comprehensive experimental study of a di-t-butyl-substituted cyclooctatetraene-based molecular balance to measure the effect of 16 different solvents on the equilibrium of folded versus unfolded isomers. In the folded 1,6-isomer, the two t-butyl groups are in close proximity (H···H distance ≈ 2.5 Å), but they are far apart in the unfolded 1,4-isomer (H···H distance ≈ 7 Å). We determined the relative strengths of these noncovalent intramolecular σ-σ interactions via temperature-dependent nuclear magnetic resonance measurements. The origins of the interactions were elucidated with energy decomposition analysis at the density functional and ab initio levels of theory, pinpointing the predominance of London dispersion interactions enthalpically favoring the folded state in any solvent measured.

London Dispersion in Alkane Solvents

The importance of London dispersion interactions in solution is an ongoing debate. Although the significance of dispersion for structure and stability is widely accepted, the degree of its attenuation in solution is still not properly understood. Quantitative evaluations are derived mostly from computations. Experimental data provide guidelines to include London dispersion in solution phase design. Herein, dispersive interactions were examined with an azobenzene probe. Alkyl substituents in meta positions of the azobenzene core were systematically varied and the effect on the half-lives for the thermally induced Z to E isomerization in several alkane solvents was determined. The results show that intramolecular dispersion is only marginally influenced. In solvents with low surface tension, reduced destabilizing solvent-solvent interactions increase the half-life up to 20 %. Specific individual interactions between alkyl chains on the azobenzene and those of the solvent lead to additional fluctuations of the half-lives. These presumably result from structural changes of the conformer ensemble.

Quantifying Interactions and Solvent Effects Using Molecular Balances and Model Complexes

Where the basic units of molecular chemistry are the bonds within molecules, supramolecular chemistry is based on the interactions that occur between molecules. Understanding the "how" and "why" of the processes that govern molecular self-assembly remains an open challenge to the supramolecular community. While many interactions are readily examined in silico through electronic structure calculations, such insights may not be directly applicable to experimentalists. The practical limitations of computationally accounting for solvation is perhaps the largest bottleneck in this regard, with implicit solvation models failing to comprehensively account for the specific nature of solvent effects and explicit models incurring a prohibitively high computational cost. Since molecular recognition processes usually occur in solution, insight into the nature and effect of solvation is imperative not only for understanding these phenomena but also for the rational design of systems that exploit them.Molecular balances and supramolecular complexes have emerged as useful tools for the experimental dissection of the physicochemical basis of various noncovalent interactions, but they have historically been underexploited as a platform for the evaluation of solvent effects. Contrasting with large biological complexes, smaller synthetic model systems enable combined experimental and computational analyses, often facilitating theoretical analyses that can work in concert with experiment.Our research has focused on the development of supramolecular systems to evaluate the role of solvents in molecular recognition, and further characterize the underlying mechanisms by which molecules associate. In particular, the use of molecular balances has provided a framework to measure the magnitude of solvent effects and to examine the accuracy of solvent models. Such approaches have revealed how solvation can modulate the electronic landscape of a molecule and how competitive solvation and solvent cohesion can provide thermodynamic driving forces for association. Moreover, the use of simple model systems facilitates the interrogation and further dissection of the physicochemical origins of molecular recognition. This tandem experimental/computational approach has married less common computational techniques, like symmetry adapted perturbation theory (SAPT) and natural bonding orbital (NBO) analysis, with experimental observations to elucidate the influence of effects that are difficult to resolve experimentally (e.g., London dispersion and electron delocalization).

Adduct under Field-A Qualitative Approach to Account for Solvent Effect on Hydrogen Bonding

The location of a mobile proton in acid-base complexes in aprotic solvents can be predicted using a simplified Adduct under Field (AuF) approach, where solute–solvent effects on the geometry of hydrogen bond are simulated using a fictitious external electric field. The parameters of the field have been estimated using experimental data on acid-base complexes in CDF₃/CDClF₂. With some limitations, they can be applied to the chemically similar CHCl₃ and CH₂Cl₂. The obtained data indicate that the solute–solvent effects are critically important regardless of the type of complexes. The temperature dependences of the strength and fluctuation rate of the field explain the behavior of experimentally measured parameters.

Hydrogen peroxide adducts of triarylphosphine oxides

Five new hydrogen peroxide adducts of phosphine oxides (p-Tol₃PO·H₂O₂)₂ (1), (o-Tol₃PO·H₂O₂)₂ (2), (o-Tol₂PhPO·H₂O₂)₂ (3), (p-Tol₃PO)₂·H₂O₂ (4), and (o-TolPh₂PO)₂·H₂O₂ (5), and the water adduct (o-Tol₂PhPO·H₂O)₂ (6) have been synthesized and fully characterized. Their single crystal X-ray structures have been determined and analyzed. The IR and ³¹P NMR data are in accordance with strong hydrogen bonding of the hydrogen peroxide. The mono- versus dimeric nature of the adduct assemblies has been investigated by DOSY NMR experiments. Raman spectroscopy of the symmetric adducts and the ν(O-O) stretching bands confirm the presence of hydrogen-bonded hydrogen peroxide in the solid materials. The solubilities in organic solvents have been quantified. Due to the high solubilities of 1-6 in organic solvents their ¹⁷O NMR spectra could be recorded in natural abundance, providing well-resolved signals for the P[double bond, length as m-dash]O and O-O groups. The adducts 1-5 have been probed regarding their stability in solution at 105 °C. The decomposition of the adduct 1 takes place by loss of the active oxygen atoms in two steps.