London Dispersion Enables the Shortest Intermolecular Hydrocarbon H···H Contact

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.

Long but Strong C-C Single Bonds: Challenges for Theory

Theoretical challenges in describing molecules with anomalously long single C-C bonds are analyzed in terms of the relative contributions of stabilizing and destabilizing intramolecular interactions. Diamondoid dimers that are stable despite the presence of C-C bonds up to 1.7 Å long, as well as other bulky molecules stabilized due to intramolecular noncovalent interactions (London dispersions) are discussed. The unexpected stability of highly crowded molecules, such as diamondoid dimers and tert-butyl-substituted hexaphenylethanes, calls for reconsideration of the "steric effect" traditionally thought to destabilize the molecule. Alternatively, "steric attraction" helps to understand bonding in sterically overloaded molecules, whose structural and energetic analysis requires a proper theoretical description of noncovalent interactions.

Molecular Wind-Up Meter for the Quantification of London Dispersion Interactions

ConspectusThe experimental quantification of interactions on the molecular level provides the necessary basis for the design of functional materials and chemical processes. The interplay of multiple parameters and the small quantity of individual interactions pose a special challenge for such endeavors. The common method is the use of molecular balances, which can exist in two different states. Thereby, a stabilizing interaction can occur in one of the states, favoring its formation and thus affecting the thermodynamic equilibrium of the system. One challenge is determining the change in this equilibrium since various analytical methods could not be applied to fast-changing equilibria. A new and promising method for quantifying molecular interactions is the use of Molecular Wind-up Meters (MWM) in which the change in kinetics, rather than the effect on thermodynamics, is investigated. An MWM is transformed with an energy input (e.g. irradiation) into a metastable state. Then, the rate of thermal transformation back to the ground state is measured. The strength of interactions present in the metastable state controls the kinetics of the back reactions, allowing direct correlation. The advantage of this approach lies in the high sensitivity (energy differences can be larger by 1 order of magnitude) and, in general, allows the use of a broader range of solvents and analytical methods. An Azobenzene-based MWM has been established as a powerful tool to quantify London dispersion interactions. London dispersion (LD) represents the attractive part of the van der Waals potential. Although neglected in the past due to its weak character, it has been shown that the influence of LD on the structure, stability, and reactivity of matter can be decisive. Especially in larger molecules, its energy contribution increases overproportionately with the number of atoms, which has sparked increasing interest in the use of so-called dispersion energy donors (DED) as a new structural element. Application of the azobenzene-based MWM not only allowed the differentiation of bulkiness, but also systematically addressed the influence of the length of n-alkyl chains. Additionally, the solvent influence on LD was studied. Based on the azobenzene MWM, an increment system has been proposed, allowing a rough estimate of the effect of a specific DED.

Influence of Steric and Dispersion Interactions on the Thermochemistry of Crowded (Fluoro)alkyl Compounds

ConspectusAlkanes play a pivotal role in industrial, environmental, and biological processes. They are characterized by their carbon-carbon single-bond structure, remarkable stability, and conformational diversity. Fluorination of such compounds imparts unique physicochemical properties that often enhance pharmacokinetic profiles, metabolic stability, and receptor interactions while keeping beneficial properties. However, such per- and polyfluoroalkyl substances (PFAS) show a persistent presence in the environment and potential adverse health effects, which propelled them to the forefront of global environmental and health discussions. Alkyl compounds are also prototypical for stereoelectronic (SE) effects that are widely applied in chemistry. Substituents are typically described as electron-density-donating/withdrawing and/or responsible for sterically interacting with reagents or strategic groups in the molecule. That alkane branching can result in higher stability compared to less-branched isomers has been investigated in detail also by testing quantum chemical methods, in particular density functional theory (DFT). Alkane branching results in spatially compact structures with close intramolecular contacts so that at a specific size the detailed balance of attractive London dispersion and covalent versus repulsive Pauli exchange interactions shifts to new, chemically unfragile situations. This may lead to dissociation at room temperature and opens the central question: what is the smallest crowed alkane that cannot be made synthetically? In this Account, we try to shed light on the interplay among the various (free) energy components for crowded (fluoro)alkane dissociation. In this context, homolytic cleavage of the central C-C bond in a series of model alkanes of increasing size with tert-butyl (tBu), adamantyl (Ad), and [1.1.1]propellanyl (Prop) substituents is investigated. Reference energies are calculated at the PNO-LCCSD(T)-F12b level and used to benchmark the performance of contemporary DFT functionals. In line with previous conclusions, the application of dispersion corrections to density functionals is mandatory. For crowed structures, the accurate description of the midrange correlation effects, specifically repulsive van der Waals interactions, is crucial, and we observed that the density-dependent VV10 correction is superior to D4 in this context, although the asymptotic region is better described by the latter. The best available dispersion-inclusive functionals show systematic and reasonably small residual errors and can be safely applied to large systems (>100 atoms), for which coupled cluster methods with large basis sets are not computationally feasible anymore. For qualitatively correct predictions of synthetic accessibility under equilibrium conditions (free energy), the inclusion of thermostatistical (entropy) contributions is also essential. According to our results, tetra-tert-butylmethane (C17tBu) is the largest and most crowded system with a positive dissociation free energy and should be synthesizable. The difference between hydrogenated and perfluorinated systems originates from the increase in the steric repulsion of spatially close substituents, which is not compensated to the same extent by attractive orbital and dispersion interactions. A sometimes-assumed similar steric demand for fluorine and hydrogen atoms is not corroborated by our investigations on crowded systems. Perfluorination is found to substantially decrease thermal stability, rendering perfluorinated hexamethylethane (C8tBuF) the last potentially stable representative.

The Decay of Dispersion Interaction and Its Remarkable Effects on the Kinetics of Activation Reactions Involving Alkyl Chains

The importance of dispersion interactions in many chemical processes is well recognized. It is known that the dispersion strength would decay with the increasing separation between the interacting groups; however, its effects on chemical reactivity have not been well understood. Here we reveal the decay law of dispersion interactions along the n-alkyl chain, its effective interaction ranges for common functional groups, and their remarkable effects on the kinetics of activation reactions involving alkyl chains. This is achieved by DLPNO-CCSD(T) calculations and the local energy decomposition analysis and is supported by experimental findings. In particular, our calculations indicate that the lifetime of alkyl-substituted cis-azobenzenes increases with the alkyl chain length but reaches a steady value when alkyl chains are longer than butyl groups, which is in satisfactory accordance with experimental measurements. We also propose a concise expression to describe the dispersion decay, which shows excellent agreement with our computed results.

An NMR Spectroscopy View on London Dispersion in Catalysis: Detection, Quantification, and Application in Ion Pair and Transition Metal Catalysis

ConspectusThe energetic contribution of London dispersion (LD) can cover a broad range from very few to hundreds of kJ mol-1 for extended interaction interfaces due to its pairwise additivity. However, for a designed and successful application of LD in chemical catalysis, there are still many obstacles and questions that remain. In principle, LD can be regarded as the attractive part of the van der Waals potential. Thus, considering the whole van der Waals potential, including the repulsive part (steric repulsion), the ideal solution to the problem in catalysis would be to design compatible interaction interfaces at exactly the correct distance. In the case of a self-assembled, flexible structure arrangement, entropic contributions and solvent interactions might be detrimental. In the case of a rigid catalyst pocket, steric hindrance might not allow for large substituents that are usually applied as dispersion energy donors (DEDs). For a working catalytic system, the following question arises: how is it possible to dissect the complex interaction interfaces in terms of energetic contributions? Usually, the energetic contribution of LD to catalysis is addressed by using calculations. However, adequately computing the correct energetic contributions can be extremely challenging for a vast conformational space with all kinds of intermolecular interactions. Thus, experimental data are essential for comparison or benchmarking.Therefore, in this Account, we describe our quest for detailed experimental data obtained via NMR spectroscopy to experimentally dissect and quantify LD in catalytic systems. In addition, we address the question of whether bulky substituents used as DEDs can be used in confined catalytic pockets. With the example of Pd phosphoramidite complexes, we show how it is possible to experimentally dissect and quantify the contribution of individual interaction areas in complicated transition metal complexes. Furthermore, a correlation between conformational rigidity and heterodimer preference clearly reveals that LD can only unfold its full potential in cases where entropic contributions are minimized. This finding can also explain the small contribution of LD in flexible and solvent-exposed molecular balances. In the field of Brønsted acid catalysis, we demonstrated that LD has a strong influence on the structures, stability, and populations of confined catalytic intermediates. LD is key for populating higher aggregates such as dimers. In addition, offsets between the experimental and computational results were observed and attributed to solvent-solute dispersion interactions. We studied the delicate interplay of attractive and repulsive interactions by adding bulky DED substituents onto a substrate, which can function as a molecular balance system. Intriguingly, the effect of LD on the free substrate was straightforwardly transferred onto the highly confined intermediates. Furthermore, this effect could even be read out in the enantioselectivities of the underlying reaction. This conceptualized a general approach regarding how LD can be used beneficially in catalysis to convert from moderate/good to excellent stereoselectivities. It showcased that bulky groups such as tert-butyl must not only be regarded as occupied volumes.

The borderless world of chemical bonding across the van der Waals crust and the valence region

The definition of the van der Waals crust as the spherical section between the atomic radius and the van der Waals radius of an element is discussed and a survey of the application of the penetration index between two interacting atoms in a wide variety of covalent, polar, coordinative or noncovalent bonding situations is presented. It is shown that this newly defined parameter permits the comparison of bonding between pairs of atoms in structural and computational studies independently of the atom sizes.

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.

London Dispersion Effects on the Stability of Heavy Tetrel Molecules

London dispersion (LD) interactions, which stem from long-range electron correlations arising from instantaneously induced dipoles can occur between neighboring atoms or molecules, for example, between H atoms within ligand C-H groups. These interactions are currently of interest as a new method of stabilizing long bonds and species with unusual oxidation states. They can also limit reactivity by installing LD enhanced groups into organic frameworks or ligand substituents. Here, we address the most recent advances in the design of LD enhanced ligands, the sterically counterintuitive structures that can be generated and the consequences that these interactions can have on the structures and reactivity of sterically crowded heavy group 14 species.

London Dispersion Interactions Imitate Pressure for Molecular Crystals

The packing of molecular crystals, in which the constituent molecular units have no directional forces, is primarily controlled by weak London dispersion (LD) forces. These forces assist in stabilizing the system by bringing the molecular units into the proximity of each other. In this paper, the same effect is shown to be externally induced by pressure. The minimal pressure required to correctly describe the crystal structure without LD interactions (P_LD) provides a quantifiable measure for the weak intermolecular interactions. LD forces are shown to be essential for an accurate description of the pressure-induced phase transitions across examples of linear, trigonal-planar, square-planar, tetrahedral, trigonal bipyramidal, and octahedral molecules.

Exploring the Limits of Intramolecular London Dispersion Stabilization with Bulky Dispersion Energy Donors in Alkane Solution

We present an experimental study of a cyclooctatetraene-based molecular balance disubstituted with increasingly bulky tert-butyl (tBu), adamantyl (Ad), and diamantyl (Dia) substituents in the 1,4-/1,6-positions for which we determined the valence-bond shift equilibrium in n-hexane (hex), n-octane (oct), and n-dodecane (dod). Computations including implicit and explicit solvation support our temperature-dependent NMR equilibrium measurements indicating that the more sterically crowded 1,6-isomer is always favored, irrespective of solvent, and that the free energy is quite insensitive to substituent size.

Probing the Size Limit of Dispersion Energy Donors with a Bifluorenylidene Balance: Magic Cyclohexyl

We report the synthesis of 14 2,2'-disubstituted 9,9'-bifluorenylidenes as molecular balances for the quantification of London dispersion interactions between various dispersion energy donors. For all balances, we measured ΔG_Z/E at 333 K using ¹H NMR in seven organic solvents. For various alkyl and aryl substituents, we generally observe a preference for the "folded" Z-isomer due to attractive London dispersion interactions. The cyclohexyl-substituted system shows the largest Z-preference in this study with ΔG_Z/E = -0.6 ± 0.05 kcal mol^-1 in all solvents, owing to the rotational freedom of cyclohexyl groups paired with their large polarizability that maximizes London dispersion interactions. On the other hand, rigid and sterically more demanding substituents like tert-butyl unexpectedly favor the unfolded E-isomer. This is a result of the close relative position in which the functional groups are positioned in this molecular balance. This close proximity is the reason for the increase of Pauli repulsion in the Z-isomers with large rigid substituents (tert-butyl, adamantyl, and diamantyl) which leads to an equilibrium shift toward the unfolded E-form. While we were able to reproduce most of our experimental trends qualitatively using contemporary computational chemistry methods, quantitative accuracy of the employed methods still needs further improvement.

Synthesis, Structural Characterization, Hirshfeld Surface Analysis, and Antibacterial Study of Pd(II) and Ni(II) Schiff Base Complexes Derived from Aliphatic Diamine

A Schiff base bearing two methyl substituents, namely, 6,6'-((1E,1'E)-((2,2-dimethylpropane-1,3-diyl) bis(azanylylidene)) bis(methanylylidene)) bis(2-methylphenol) [H₂AD1Me] was synthesized and characterized through physicochemical and spectroscopic analyses. Then, the Schiff base was complexed with Pd(II) and Ni(II) to form [Pd(AD1Me)] and [Ni(AD1Me)], respectively. Both metal complexes were successfully obtained and characterized through several analyses, viz., melting point, elemental analysis, molar conductivity, magnetic susceptibility, FTIR, ¹H NMR, UV-vis, and single crystal X-ray diffraction. A quantitative analysis of the intermolecular interactions in the crystal structures has been performed using Hirshfeld surface analysis. Both metal complexes were crystallized in a monoclinic crystal system with the space group of P2₁/c. Additionally, the deprotonated phenolic oxygen atom (O1/O2) and azomethine nitrogen atom (N1/N2) of the ligand chelate the Pd(II) and Ni(II) ions, forming a slightly distorted square-planar complex containing three six-membered rings encircling the metal core with dsp² hybridization. The shift of ν(C=N) to a higher frequency in FTIR by 26-28 cm^-1 indicated that the complexation to Pd(II) and Ni(II) through the azomethine N was established. It was further supported through the shifting of the azomethine proton signal to higher or lower chemical shifts with Δδ = 0.43-1.15 ppm in ¹H NMR. In addition, the shifting of the n-π*(C=N) band in UV-vis spectra with Δλ = 24-40 nm indicated the involvement of azomethine nitrogen in the complexation. All the compounds showed no significant antibacterial activity against three bacterial strains, namely, Staphylococcus aureus subsp. aureus Rosenbach (ATCC 6538), Streptococcus mutans Clarke (ATCC 700,610), and Proteus vulgaris (ATCC 6380), as the percent growth inhibition calculated was less than 90%.

Fabrication of Functional bioMOF-100 Prototype as Drug Delivery System for Breast Cancer Therapy

Breast cancer is the most frequent cause of cancer death in women, representing the fifth leading cause of cancer death overall. Therefore, the growing search for the development of new treatments for breast cancer has been developed lately as well as drug delivery systems such as biocompatible metal-organic Frameworks (bio-MOFs). These may be promising and attractive for drug incorporation and release. The present study aims to develop a drug carrier system RCA (bioMOF-100 submitted to the activation process) containing incorporated curcumin (CCM), whose material surface is coated with folic acid molecules (FA) to promote the targeting of drug carrier systems to the tumor region. They were synthesized and characterized using several characterization techniques. The materials were submitted to drug encapsulation tests, whose encapsulation efficiency was 32.80% for CCM@RCA-1D. Using the ¹H nuclear magnetic resonance (NMR) spectroscopy technique, it was possible to verify the appearance of signals referring to folic acid, suggesting success in the functionalization of these matrices. In vitro tests such as cell viability and type of cell death were evaluated in both series of compounds (CCM@RCA-1D, CCM@RCA-1D/FA) in breast tumor lines. The results revealed low toxicity of the materials and cell death by late apoptosis. Thus, these results indicate that the matrices studied can be promising carriers in the treatment of breast cancer.

Tilting the Balance: London Dispersion Systematically Enhances Enantioselectivities in Brønsted Acid Catalyzed Transfer Hydrogenation of Imines

London dispersion (LD) is attracting more and more attention in catalysis since LD is ubiquitously present and cumulative. Since dispersion is hard to grasp, recent research has concentrated mainly on the effect of LD in individual catalytic complexes or on the impact of dispersion energy donors (DEDs) on balance systems. The systematic transfer of LD effects onto confined and more complex systems in catalysis is still in its infancy, and no general approach for using DED residues in catalysis has emerged so far. Thus, on the example of asymmetric Brønsted acid catalyzed transfer hydrogenation of imines, we translated the findings of previously isolated balance systems onto confined catalytic intermediates, resulting in a systematic enhancement of stereoselectivity when employing DED-substituted substrates. As the imine substrate is present as Z- and E-isomers, which can, respectively, be converted to R- and S-product enantiomers, implementing tert-butyl groups as DED residues led to an additional stabilization of the Z-imine by up to 4.5 kJ/mol. NMR studies revealed that this effect is transferred onto catalyst/imine and catalyst/imine/nucleophile intermediates and that the underlying reaction mechanism is not affected. A clear correlation between ee and LD stabilization was demonstrated for 3 substrates and 10 catalysts, allowing to convert moderate-good to good-excellent enantioselectivities. Our findings conceptualize a general approach on how to beneficially employ DED residues in catalysis: they clearly showcase that bulky alkyl residues such as tert-butyl groups must be considered regarding not only their repulsive steric bulk but also their attractive properties even in catalytic complexes.

Best-Practice DFT Protocols for Basic Molecular Computational Chemistry

Nowadays, many chemical investigations are supported by routine calculations of molecular structures, reaction energies, barrier heights, and spectroscopic properties. The lion's share of these quantum-chemical calculations applies density functional theory (DFT) evaluated in atomic-orbital basis sets. This work provides best-practice guidance on the numerous methodological and technical aspects of DFT calculations in three parts: Firstly, we set the stage and introduce a step-by-step decision tree to choose a computational protocol that models the experiment as closely as possible. Secondly, we present a recommendation matrix to guide the choice of functional and basis set depending on the task at hand. A particular focus is on achieving an optimal balance between accuracy, robustness, and efficiency through multi-level approaches. Finally, we discuss selected representative examples to illustrate the recommended protocols and the effect of methodological choices.

Towards Hexagonal Planar Nickel: A Dispersion-Stabilised Tri-Lithium Nickelate

Advancing the understanding of lithum nickelate complexes, here we report a family of homoleptic organonickelate complexes obtained by reacting Ni(COD)₂ and lithium aryl-acetylides in the presence of the bidentate donor TMEDA. These compounds represent rare examples of low-valent transition-metals supported solely by organolithium ligands. Whilst the solid-state structures indicate a hexagonal planar geometry around Ni⁰ with Ni-Li bonds, bonding analysis via QTAIM, NCI, NBO and ELI methods reveals that the Ni-Li interactions are repulsive in nature, characterising these complexes as tri-coordinated. London dispersion forces between TMEDA and the organic substituents on nickel are found to play a crucial role in the stabilisation and thus isolation of these complexes. Preliminary reactivity studies demonstrate that the homoleptic lithium nickelates undergo stoichiometric cross-coupling with PhI to give dinickel clusters containing both anionic acetylide and neutral alkyne ligands.

Inhibition of Alkali Metal Reduction of 1-Adamantanol by London Dispersion Effects

A series of alkali metal 1-adamantoxide (OAd¹ ) complexes of formula [M(OAd¹ )(HOAd¹ )₂ ], where M=Li, Na or K, were synthesised by reduction of 1-adamantanol with excess of the alkali metal. The syntheses indicated that only one out of every three HOAd¹ molecules was reduced. An X-ray diffraction study of the sodium derivative shows that the complex features two unreduced HOAd¹ donors as well as the reduced alkoxide (OAd¹ ), with the Ad¹ fragments clustered together on the same side of the NaO₃ plane, contrary to steric considerations. This is the first example of an alkali metal reduction of an alcohol that is inhibited from completion due to the formation of the [M(OAd¹ )(HOAd¹ )₂ ] complexes, stabilized by London dispersion effects. NMR spectroscopic studies revealed similar structures for the lithium and potassium derivatives. Computational analyses indicate that decisive London dispersion effects in the molecular structure are a consequence of the many C-H⋅⋅⋅H-C interactions between the OAd¹ groups.

Synthetic Doping of Diamondoids through Skeletal Editing

We present a strategy for the skeletal editing of diamondoid structures to selectively displace methylene for heteroatom moieties in the carbon framework. This constitutes a synthetic approach to doping diamond-like structures with electron donor dopants (O, N, and S). The key steps involve two subsequent retro-Barbier fragmentations followed by cage reconstruction in the presence of a dopant. Remarkably, the incorporation of n-dopants reduces the strain of the diamondoid cage as shown through homodesmotic equations.

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.

Open-Shell Variant of the London Dispersion-Corrected Hartree-Fock Method (HFLD) for the Quantification and Analysis of Noncovalent Interaction Energies

The London dispersion (LD)-corrected Hartree–Fock (HF) method (HFLD) is an ab initio approach for the quantification and analysis of noncovalent interactions (NCIs) in large systems that is based on the domain-based local pair natural orbital coupled-cluster (DLPNO-CC) theory. In the original HFLD paper, we discussed the implementation, accuracy, and efficiency of its closed-shell variant. Herein, an extension of this method to open-shell molecular systems is presented. Its accuracy is tested on challenging benchmark sets for NCIs, using CCSD(T) energies at the estimated complete basis set limit as reference. The HFLD scheme was found to be as accurate as the best-performing dispersion-corrected exchange-correlation functionals, while being nonempirical and equally efficient. In addition, it can be combined with the well-established local energy decomposition (LED) for the analysis of NCIs, thus yielding additional physical insights.

London Dispersion Helps Refine Steric A-Values: Dispersion Energy Donor Scales

We suggest a scale of dispersion energy donors (DEDs) that allows for direct comparisons with steric effects. This scale is based on the classic A-values and allows groups to reorient to minimize strain, thereby providing an advantage over raw group polarizabilities. The A-value can no longer be considered purely a steric factor. Even for groups that do not participate in charge transfer or electrostatic interactions, the A-value includes Pauli repulsion (steric hindrance) and attractive London dispersion (LD) interactions. Although the common assumption is that, at the distances found in monosubstituted cyclohexanes, steric demands are the key factors influencing conformer preferences, we show in this computational study that there is a non-negligible LD part. We use this system to build a DED scale and a complementary steric scale. These scales are quantitatively comparable, as they are based on the same system, and allow for comparison of the two competing interactions in experimentally relevant settings. In addition, we show that LD interactions can be used to explain puzzling data regarding relative group sizes.

Could London Dispersion Force Control Regioselective (2 + 2) Cyclodimerizations of Benzynes? YES: Application to the Synthesis of Helical Biphenylenes

In recent years, London dispersion interactions, which are the attractive component of the van der Waals potential, have been found to play an important role in controlling the regio- and/or stereoselectivity of various reactions. Particularly, the dispersion interactions between substrates and catalysts (or ligands) are dominant in various selective catalyzes. In contrast, repulsive steric interactions, rather than the attractive dispersion interactions, between bulky substituents are predominant in most of the noncatalytic reactions. Herein, we demonstrate the first example of London dispersion-controlled noncatalytic (2 + 2) cyclodimerization of substituted benzynes to selectively afford proximal biphenylenes in high yields and regioselectivities, depending on the extent of dispersion interactions in the substituents. This method can be applied for the synthesis of novel helical biphenylenes, which would be fascinating for chemists as these compounds are potential skeletons for ligands, catalysts, and medicines.

Understanding the Nature and Properties of Hydrogen-Hydrogen Bonds: The Stability of a Bulky Phosphatetrahedrane as a Case Study

Recently, the first mixed C/P phosphatetrahedranes (^tBuC)₃P and (^tBuCP)₂ were reported. Unlike (^tBuCP)₂, (^tBuC)₃P exhibits remarkable thermal stability, which can be partially attributed to a network of nine hydrogen-hydrogen bonds (HHBs) localized between the tert-butyl substituents. The stabilizing contribution arising from this network of HHBs was obtained from local energy decomposition (LED) analysis calculated at the domain-based local pair natural orbital CCSD(T) (DLPNO-CCSD(T)) level of theory. These calculations suggest that each HHB contributes approximately -0.7 kcal/mol of stabilization; however, the net stabilization energy likely lies between -0.25 and -0.5 kcal/mol because of steric repulsion. Spatial analysis of the London dispersion energy via a dispersion interaction density (DID) plot reveals that the DID surface is localized at key C-H groups involved in HHBs, consistent with London dispersion interactions predominantly arising from HHBs. In addition, we present a computed mechanism that supports a phosphinidenoid species as a key reaction intermediate in the synthesis of (^tBuC)₃P.

Ultrashort Hδ+Hδ- intermolecular distance in a supramolecular system in the solid state

Herein we report experimental evidence for the shortest intermolecular distance reported for two electronically-different hydrogen atoms in the solid state. The Hδ+Hδ- non-covalent interaction was studied using theoretical calculations indicating that electrostatic and dispersion forces are of paramount importance.

Dispersion-Bound Isolated Dimers in the Gas Phase: Observation of the Shortest Intermolecular CH⋅⋅⋅H-C Distance via Stimulated Raman Spectroscopy

The triphenylmethane and all‐meta tert‐butyl triphenylmethane dimers, (TPM)₂ and (T ^tBuPM)₂, respectively, were studied with ionization loss stimulated Raman spectroscopy in molecular beam experiments to resolve structure sensitive vibrations. This answers the question whether the recently reported linear head‐to‐head arrangement in (T ^tBuPM)₂ results from crystal packing or prevails also in the gas phase, and therefore must result from extraordinarily strong London dispersion (LD) interactions. Our study clearly demonstrates that the head‐to‐head arrangement is maintained even under isolated molecular beam conditions in the absence of crystal packing effects. The central Raman‐active aliphatic C−D vibration of appropriately deuterated (T ^tBuPM)₂ associated with an unusually short C−D⋅⋅⋅D−C distance exhibits a strong blue‐shift compared to the undisturbed case. As the LD stabilizing tert‐butyl groups are absent in (TPM)₂, it displays an approximately S₆‐symmetric tail‐to‐tail arrangement.

CHAL336 Benchmark Set: How Well Do Quantum-Chemical Methods Describe Chalcogen-Bonding Interactions?

We present the CHAL336 benchmark set-the most comprehensive database for the assessment of chalcogen-bonding (CB) interactions. After careful selection of suitable systems and identification of three high-level reference methods, the set comprises 336 dimers each consisting of up to 49 atoms and covers both σ- and π-hole interactions across four categories: chalcogen-chalcogen, chalcogen-π, chalcogen-halogen, and chalcogen-nitrogen interactions. In a subsequent study of DFT methods, we re-emphasize the need for using proper London dispersion corrections when treating noncovalent interactions. We also point out that the deterioration of results and systematic overestimation of interaction energies for some dispersion-corrected DFT methods does not hint at problems with the chosen dispersion correction but is a consequence of large density-driven errors. We conclude this work by performing the most detailed DFT benchmark study for CB interactions to date. We assess 109 variations of dispersion-corrected and dispersion-uncorrected DFT methods and carry out a detailed analysis of 80 of them. Double-hybrid functionals are the most reliable approaches for CB interactions, and they should be used whenever computationally feasible. The best three double hybrids are SOS0-PBE0-2-D3(BJ), revDSD-PBEP86-D3(BJ), and B2NCPLYP-D3(BJ). The best hybrids in this study are ωB97M-V, PW6B95-D3(0), and PW6B95-D3(BJ). We do not recommend using the popular B3LYP functional nor the MP2 approach, which have both been frequently used to describe CB interactions in the past. We hope to inspire a change in computational protocols surrounding CB interactions that leads away from the commonly used, popular methods to the more robust and accurate ones recommended herein. We would also like to encourage method developers to use our set for the investigation and reduction of density-driven errors in new density functional approximations.

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.

Assessing the Applicability of the Geometric Counterpoise Correction in B2PLYP/Double-ζ Calculations for Thermochemistry, Kinetics, and Noncovalent Interactions

We present a proof-of-concept study of the suitability of Kruse and Grimme’s geometric counterpoise correction (gCP) for basis set superposition errors (BSSEs) in double-hybrid density functional calculations with a double-ζ basis set. The gCP approach only requires geometrical information as an input and no orbital/density information is needed. Therefore, this correction is practically free of any additional cost. gCP is trained against the Boys and Bernardi counterpoise correction across a set of 528 noncovalently bound dimers. We investigate the suitability of the approach for the B2PLYP/def2-SVP level of theory, and reveal error compensation effects—missing London dispersion and the BSSE—associated with B2PLYP/def2-SVP calculations, and present B2PLYP-gCP-D3(BJ)/def2-SVP with the reparametrised DFT-D3(BJ) and gCP corrections as a more balanced alternative. Benchmarking results on the S66x8 benchmark set for noncovalent interactions and the GMTKN55 database for main-group thermochemistry, kinetics, and noncovalent interactions show a statistical improvement of the B2PLYP-gCP-D3(BJ) scheme over plain B2PLYP and B2PLYP-D3(BJ). B2PLYP-D3(BJ) shows significant overestimation of interaction energies, barrier heights with larger deviations from the reference values, and wrong relative stabilities in conformers, all of which can be associated with BSSE. We find that the gCP-corrected method represents a significant improvement over B2PLYP-D3(BJ), particularly for intramolecular noncovalent interactions. These findings encourage future developments of efficient double-hybrid DFT strategies that can be applied when double-hybrid calculations with large basis sets are not feasible due to system size.

Cost-effective composite methods for large-scale solid-state calculations

Following the development in recent years of progressively more accurate approximations to the exchange-correlation functional, the use of density functional theory (DFT) methods to examine increasingly large and complex systems has grown, in particular for solids and other condensed matter systems. However the cost of these calculations is high, often requiring the use of specialist HPC facilities. As such, for the purpose of large-scale high-throughput screening of material properties, a hierarchy of simplified DFT methods has been proposed that allows rapid electronic structure calculation of large systems, and we have recently extended this scheme to the solid state (sol-3c). Here, we analyze the applicability and scaling of the new sol-3c DFT methods to molecules and crystals composed of light-elements, such as small proteins and model DNA-helices. Furthermore, the calculation of the electronic structure of large to very large porous systems, such as metal-organic frameworks and inorganic nanoparticles, is discussed. The new composite methods have been implemented in the CRYSTAL17 code, which efficiently implements hybrid functionals and enables routine application of the new methods to large-scale calculations of such materials with excellent performance, even with small-scale computing resources.

Molecular designs for expanding the limits of ultralong C-C bonds and ultrashort HH non-bonded contacts

Recent experiments have reported the formation of very long C-C bonds (d_C-C > 1.80 Å) and very short HH non-bonded contacts (d_HH < 1.5 Å) in several sets of molecules. Both these rare phenomena arise due to specific donor-acceptor interactions and London dispersion interactions respectively. Favorable negative hyperconjugation, namely H₂N(lone-pair) →σ*(C-C), creates an ultralong C-C bond in diamino-o-carborane with d_C-C > 1.829 Å and a planar amine reminiscent of a transition-state like structure for ammonia inversion. The small and narrow barrier favours rapid inversion through quantum mechanical tunnelling (QMT) and produces a translationally averaged planar amine as observed in the experiments. On the other hand, designing specific confined molecular cavities or chambers like in,in-bis(hydrosilane) or its germanane analogs furnishes an ultrashort HH distance = 1.47 Å and 1.38 Å respectively. The predisposition of such closely placed HH contacts arises from the rather effective attractive dispersion interactions between them. Controlling the strength of the dispersion interactions provides a rich landscape for realizing such close HH distances. Molecular design ably assisted by computational modeling to further tune these interactions provides new avenues to break the glass-ceilings of ultralong C-C bonds or ultrashort HH contacts. Dispersion-corrected DFT calculations and ab initio molecular dynamics simulations generate a large library of such unique features in a diverse class of molecules. This feature article highlights the design principles to realize hitherto longest C-C bonds/shortest HH contacts.

Liquid water contains the building blocks of diverse ice phases

Water molecules can arrange into a liquid with complex hydrogen-bond networks and at least 17 experimentally confirmed ice phases with enormous structural diversity. It remains a puzzle how or whether this multitude of arrangements in different phases of water are related. Here we investigate the structural similarities between liquid water and a comprehensive set of 54 ice phases in simulations, by directly comparing their local environments using general atomic descriptors, and also by demonstrating that a machine-learning potential trained on liquid water alone can predict the densities, lattice energies, and vibrational properties of the ices. The finding that the local environments characterising the different ice phases are found in water sheds light on the phase behavior of water, and rationalizes the transferability of water models between different phases.

Molecular understanding of water is challenging due to the structural complexity of liquid water and the large number of ice phases. Here the authors use a machine-learning potential trained on liquid water to demonstrate the structural similarity of liquid water and that of 54 real and hypothetical ice phases.

Selective Activation of Aromatic Aldehydes Promoted by Dispersion Interactions: Steric and Electronic Factors of a π-Pocket within Cage-Shaped Borates for Molecular Recognition

Selective bond formations are one of the most important reactions in organic synthesis. In the Lewis acid mediated electrophile reactions of carbonyls, the selective formation of a carbonyl-acid complex plays a critical role in determining selectivity, which is based on the difference in the coordinative interaction between the carbonyl and Lewis acid center. Although this strategy has attained progress in selective bond formations, the discrimination between similarly sized aromatic and aliphatic carbonyls that have no functional anchors to strongly interact with the metal center still remains a challenging issue. Herein, this work focuses on molecular recognition driven by dispersion interactions within some aromatic moieties. A Lewis acid catalyst with a π-space cavity, which is referred to as a π-pocket, as the recognition site for aromatic carbonyls is designed. Cage-shaped borates 1B with various π-pockets demonstrated significant chemoselectivity for aromatic aldehydes 3 b-f over that of aliphatic 3 a in competitive hetero-Diels-Alder reactions. The effectiveness of our catalysts was also evidenced by intramolecular recognition of the aromatic carbonyl within a dicarbonyl substrate. Mechanistic and theoretical studies demonstrated that the selective activation of aromatic substrates was driven by the preorganization step with a larger dispersion interaction, rather than the rate-determining step of the C-C bond formation, and this was likely to contribute to the preferred activation of aromatic substrates over that of aliphatic ones.

A Comprehensive Study on Tetraaryltetrabenzoporphyrins

Tetraaryltetrabenzoporphyrins (TATBPs) show, due to their optoelectronic properties, rising potential as dyes in various fields of physical and biomedical sciences. However, unlike in the case of porphyrins, the potential structural diversity of TATBPs has been explored only to little extent, owed mainly to synthetic hurdles. Herein, we prepared a comprehensive library of 30 TATBPs and investigated their fundamental properties. We elucidated structural properties by X-ray crystallography and found explanations for physical properties such as solubility. Fundamental electronic aspects were studied by optical spectroscopy as well as by electrochemistry and brought in context to the stability of the molecules. Finally, we were able to develop a universal synthetic protocol, utilizing a readily established isoindole synthon, which gives TATBPs in high yields, regardless of the nature of the used arylaldehyde and without meticulous chromatographic purifications steps. This work serves as point of orientation for scientists, that aim to utilize these molecules in materials, nanotechnological, and biomedical applications.