London Dispersion Interactions Research Articles

Antiferromagnetic ordering based on intermolecular London dispersion interactions in amphiphilic TEMPO ammonium salts.

Antiferromagnetic coupling in TEMPO-based radicals can be enhanced via self-assembly through London dispersion interactions in amphiphilic solids. The synthesis, magnetic characterization, and three crystal structures of the solid radical ion salts (R-DMAT-n)X with various counterions X and alkyl chain lengths n are reported. Magnetic susceptibility and absolute EPR signal intensity measurements show singlet-triplet transitions in a number of cases, which is discussed in relation to the crystal structures. Antiferromagnetic ordering effects are sensitive to both the length of the alkyl chain and the counter anion.

Molecular interactions in nanocellulose assembly.

The contribution of hydrogen bonds and the London dispersion force in the cohesion of cellulose is discussed in the light of the structure, spectroscopic data, empirical molecular-modelling parameters and thermodynamics data of analogue molecules. The hydrogen bond of cellulose is mainly electrostatic, and the stabilization energy in cellulose for each hydrogen bond is estimated to be between 17 and 30 kJ mol-1 On average, hydroxyl groups of cellulose form hydrogen bonds comparable to those of other simple alcohols. The London dispersion interaction may be estimated from empirical attraction terms in molecular modelling by simple integration over all components. Although this interaction extends to relatively large distances in colloidal systems, the short-range interaction is dominant for the cohesion of cellulose and is equivalent to a compression of 3 GPa. Trends of heat of vaporization of alkyl alcohols and alkanes suggests a stabilization by such hydroxyl group hydrogen bonding to be of the order of 24 kJ mol-1, whereas the London dispersion force contributes about 0.41 kJ mol-1 Da-1 The simple arithmetic sum of the energy is consistent with the experimental enthalpy of sublimation of small sugars, where the main part of the cohesive energy comes from hydrogen bonds. For cellulose, because of the reduced number of hydroxyl groups, the London dispersion force provides the main contribution to intermolecular cohesion.This article is part of a discussion meeting issue 'New horizons for cellulose nanotechnology'.

Can Coordination-Driven Supramolecular Self-Assembly Reactions Be Conducted from Fully Aliphatic Linkers?

The reaction between a preassembled CuI bimetallic molecular clip with a short intermetallic distance and a series of fully aliphatic cyano-capped ditopic linkers with increasing lengths was investigated. It is shown that, depending on the length of the ditopic linkers, the rational design of unprecedented supramolecular compact metallacycles containing fully aliphatic walls is possible. The specific preorganized molecular arrangement of the molecular clip used favors stabilizing interlinker London dispersion interactions, which allow, as the length of the linkers increases, the selective formation of discrete compact metallacycles at the expense of 1D coordination polymers. The generalizability of this approach was demonstrated by the reaction of fully aliphatic cyano-capped linkers with two other types of preassembled CuI bimetallic molecular clips that also had short intermetallic distances.

Evaluating Force-Field London Dispersion Coefficients Using the Exchange-Hole Dipole Moment Model

London dispersion interactions play an integral role in materials science and biophysics. Force fields for atomistic molecular simulations typically represent dispersion interactions by the 12-6 Lennard-Jones potential using empirically determined parameters. These parameters are generally underdetermined, and there is no straightforward way to test if they are physically realistic. Alternatively, the exchange-hole dipole moment (XDM) model from density-functional theory predicts atomic and molecular London dispersion coefficients from first principles, providing an innovative strategy to validate the dispersion terms of molecular-mechanical force fields. In this work, the XDM model was used to obtain the London dispersion coefficients of 88 organic molecules relevant to biochemistry and pharmaceutical chemistry and the values compared with those derived from the Lennard-Jones parameters of the CGenFF, GAFF, OPLS, and Drude polarizable force fields. The molecular dispersion coefficients for the CGenFF, GAFF, and OPLS models are systematically higher than the XDM-calculated values by a factor of roughly 1.5, likely due to neglect of higher order dispersion terms and premature truncation of the dispersion-energy summation. The XDM dispersion coefficients span a large range for some molecular-mechanical atom types, suggesting an unrecognized source of error in force-field models, which assume that atoms of the same type have the same dispersion interactions. Agreement with the XDM dispersion coefficients is even poorer for the Drude polarizable force field. Popular water models were also examined, and TIP3P was found to have dispersion coefficients similar to the experimental and XDM references, although other models employ anomalously high values. Finally, XDM-derived dispersion coefficients were used to parametrize molecular-mechanical force fields for five liquids-benzene, toluene, cyclohexane, n-pentane, and n-hexane-which resulted in improved accuracy in the computed enthalpies of vaporization despite only having to evaluate a much smaller section of the parameter space.

Direction-dependent secondary bonds and their stepwise melting in a uracil-based molecular crystal studied by infrared spectroscopy and theoretical modeling

Three types of supramolecular interactions are identified in the three crystallographic directions in crystals of 1,4-bis[(1-hexylurac-6-yl) ethynyl]benzene, a uracil-based molecule with a linear backbone. These three interactions, characterized by their strongest component, are: intermolecular double H-bonds along the molecular axis, London dispersion interaction of hexyl chains connecting these linear assemblies, and π–π stacking of the aromatic rings perpendicular to the molecular planes. On heating, two transitions happen, disordering of hexyl chains at 473 K, followed by H-bond melting at 534 K. The nature of the bonds and transitions was established by matrix-isolation and temperature-dependent infrared spectroscopy and supported by theoretical computations.

London Dispersion Directs On-Surface Self-Assembly of [121]Tetramantane Molecules.

London dispersion (LD) acts between all atoms and molecules in nature, but the role of LD interactions in the self-assembly of molecular layers is still poorly understood. In this study, direct visualization of single molecules using atomic force microscopy with CO-functionalized tips revealed the exact adsorption structures of bulky and highly polarizable [121]tetramantane molecules on Au(111) and Cu(111) surfaces. We determined the absolute molecular orientations of the completely sp3-hybridized tetramantanes on metal surfaces. Moreover, we demonstrate how LD drives this on-surface self-assembly of [121]tetramantane hydrocarbons, resulting in the formation of a highly ordered 2D lattice. Our experimental findings were underpinned by a systematic computational study, which allowed us to quantify the energies associated with LD interactions and to analyze intermolecular close contacts and attractions in detail.

Physico‐chemical and Spectroscopic Approach to Analyse the Behaviour of Surface‐Active Ionic Liquid and Conventional Surfactant in Aqueous Glycine

AbstractThe present work reports a comparison of direct or indirect (water‐mediated) interactions of surface active ionic liquid SAIL, tetrabutylammonium dodecylsulphate (TBADS) in aqueous solutions of glycine with those of sodium dodecylsulphate (SDS) in same solvent system, by employing techniques like conductivity, fluorescence, UV–Vis probe studies and clouding phenomenon. High hydrophobicity and penetrating ability of TBA+ counterions and most significant salting‐out effect due to the preferential hydration of glycine molecules accounts for the enhanced micellization and better surface activity of TBADS as compared to its conventional analogue SDS. Further, driving force for the micellization in both the cases is the hydrophobic interactions, whereas in the case of SDS, London dispersion interactions as alternative forces also contribute considerably at higher temperatures. In addition, the decrease in the cloud point () of aqueous TBADS solutions in the presence of glycine is attributed to the dehydration of hydrophilic head groups of micelles, resulting in their collapse and hence phase separation.

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

Neutron diffraction of tri(3,5-tert-butylphenyl)methane at 20 K reveals an intermolecular C-H···H-C distance of only 1.566(5) Å, which is the shortest reported to date. The compound crystallizes as a C3-symmetric dimer in an unusual head-to-head fashion. Quantum chemical computations of the solid state at the HSE-3c level of theory reproduce the structure and the close contact well (1.555 Å at 0 K) and emphasize the significance of packing effects; the gas-phase dimer structure at the same level shows a 1.634 Å C-H···H-C distance. Intermolecular London dispersion interactions between contacting tert-butyl substituents surrounding the central contact deliver the decisive energetic contributions to enable this remarkable bonding situation.

Manganese(I)-Catalyzed Dispersion-Enabled C-H/C-C Activation.

C-H/C-C Functionalizations were achieved with the aid of a versatile manganese(I) catalyst. Thus, an organometallic manganese-catalyzed C-H activation set the stage for silver-free C-H/C-C transformations with ample substrate scope and excellent levels of chemo-, site-, and diastereo-selectivities. The robust nature of the manganese(I) catalysis regime was reflected by the first C-H/C-C functionalization on amino acids under racemization-free reaction conditions. Detailed experimental and computational mechanistic studies provided strong evidence for a facile C-H activation and a rate-determining C-C cleavage, with considerable contribution from London dispersion interactions.

The role of London dispersion interactions in strong and moderate intermolecular hydrogen bonds in the crystal and in the gas phase

Two variants of density functional theory computations have been applied to characterization of hydrogen bonds of the 1-(2-hydroxylethyl)-3-methylimidazolium acetate ([C2OHmim][OAc]), i.e. with and without inclusion of dispersion interactions. A comparison of the results demonstrates that London dispersion interactions have very little impact on the energetical, geometrical, infrared spectroscopic and electron density parameters of charge-assisted intermolecular hydrogen bonds functioning both in the crystal of the [C2OHmim][OAc] and in the isolated [C2OHmim]+ [OAc]− ion pairs.

Computational insights into the oxidation of mono- and 1,4 disubstituted arenes by the Toluene Dioxygenase enzymatic complex

Toluene Dioxygenase (TDO) enzymatic complex has been widely used for preparation of enantiopure cis-cyclohexadienediols for synthetic applications. Along the last 30 years, a variety of mono- and di-substituted arenes have been studied as substrates for this enzyme, and some interesting observations have been reported regarding the yield and selectivity of the biotransformation. Nevertheless, none of them has been explained considering the active site’ structural and electronic features. In this work we present the first computational model of TDO’s active site, with a description of its architecture and interactions with the substrate to understand and predict substrate orientation. Our findings indicate that in the O2-free TDO, the iron(II) is stabilized by the coordination of an oxygen atom from the neighboring Gln215 residue. Besides, the active site is comprised by four pockets with different relative affinities for the substrates’ substituents. Monosubstituted arenes adopt a pose in which the alkyl chains maximize the London dispersion interactions with minimal steric clashes, giving an explanation for the observed trend in the benzylic hydroxylation yields. Finally, the computational results allowed us to rationalize the enantiomeric excess of 1,4-disubstituted arenes and the regioselectivity of the dihydroxylation of fluorobenzene. These data were used to develop and validate an affinity index (AI), as a quantitative indicator of the orientation preference for these substrates. This practical and easy-to-use tool can be applied to successfully predict the orientation of the para-disubstituted benzenes into the TDO active site.

Benchmark tests of a strongly constrained semilocal functional with a long-range dispersion correction

The strongly constrained and appropriately normed (SCAN) semilocal density functional [J. Sun, A. Ruzsinszky, J. P. Perdew \textit{Phys. Rev. Lett.} {\bf 115}, 036402 (2015)] obeys all 17 known exact constraints for meta-generalized-gradient approximations (meta-GGA) and includes some medium range correlation effects. Long-range London dispersion interactions are still missing, but can be accounted for via an appropriate correction scheme. In this study, we combine SCAN with an efficient London dispersion correction and show that lattice energies of simple organic crystals can be improved with the applied correction by 50\%. The London-dispersion corrected SCAN meta-GGA outperforms all other tested London-dispersion corrected meta-GGAs for molecular geometries. Our new method delivers mean absolute deviations (MADs) for main group bond lengths that are consistently below 1\,pm, rotational constants with MADs of 0.2\%, and noncovalent distances with MADs below 1\%. For a large database of general main group thermochemistry and kinetics, it also delivers a weighted mean absolute deviation below 4 kcal/mol, one of the lowest for long-range corrected meta-GGA functionals. We also discuss some consequences of numerical sensitivity encountered for meta-GGAs.

Organic crystal polymorphism: a benchmark for dispersion corrected mean field electronic structure methods

We analyze the energy landscape of the sixth crystal structure prediction blind test targets with various first principles and semi-empirical quantum chemical methodologies. A new benchmark set of 59 crystal structures (termed POLY59) for testing quantum chemical methods based on the blind test target crystals is presented. We focus on different means to include London dispersion interactions within the density functional theory (DFT) framework. We show the impact of pairwise dispersion corrections like the semi-empirical D2 scheme, the Tkatchenko–Scheffler (TS) method, and the density-dependent dispersion correction dDsC. Recent methodological progress includes higher-order contributions in both the many-body and multipole expansions. We use the D3 correction with Axilrod–Teller–Muto type three-body contribution, the TS based many-body dispersion (MBD), and the nonlocal van der Waals density functional (vdW-DF2). The density functionals with D3 and MBD correction provide an energy ranking of the blind test polymorphs in excellent agreement with the experimentally found structures. As a computationally less demanding method, we test our recently presented minimal basis Hartree–Fock method (HF-3c) and a density functional tight-binding Hamiltonian (DFTB). Considering the speed-up of three to four orders of magnitudes, the energy ranking provided by the low-cost methods is very reasonable. We compare the computed geometries with the corresponding X-ray data where TPSS-D3 performs best. The importance of zero-point vibrational energy and thermal effects on crystal densities is highlighted.

Organic crystal polymorphism: a benchmark for dispersion-corrected mean-field electronic structure methods.

We analyze the energy landscape of the sixth crystal structure prediction blind test targets with various first principles and semi-empirical quantum chemical methodologies. A new benchmark set of 59 crystal structures (termed POLY59) for testing quantum chemical methods based on the blind test target crystals is presented. We focus on different means to include London dispersion interactions within the density functional theory (DFT) framework. We show the impact of pairwise dispersion corrections like the semi-empirical D2 scheme, the Tkatchenko-Scheffler (TS) method, and the density-dependent dispersion correction dDsC. Recent methodological progress includes higher-order contributions in both the many-body and multipole expansions. We use the D3 correction with Axilrod-Teller-Muto type three-body contribution, the TS based many-body dispersion (MBD), and the nonlocal van der Waals density functional (vdW-DF2). The density functionals with D3 and MBD correction provide an energy ranking of the blind test polymorphs in excellent agreement with the experimentally found structures. As a computationally less demanding method, we test our recently presented minimal basis Hartree-Fock method (HF-3c) and a density functional tight-binding Hamiltonian (DFTB). Considering the speed-up of three to four orders of magnitudes, the energy ranking provided by the low-cost methods is very reasonable. We compare the computed geometries with the corresponding X-ray data where TPSS-D3 performs best. The importance of zero-point vibrational energy and thermal effects on crystal densities is highlighted.

ChemInform Abstract: Small Atomic Orbital Basis Set First‐Principles Quantum Chemical Methods for Large Molecular and Periodic Systems: A Critical Analysis of Error Sources

AbstractReview: 123 refs.

Dispersion-Corrected Mean-Field Electronic Structure Methods.

Mean-field electronic structure methods like Hartree-Fock, semilocal density functional approximations, or semiempirical molecular orbital (MO) theories do not account for long-range electron correlation (London dispersion interaction). Inclusion of these effects is mandatory for realistic calculations on large or condensed chemical systems and for various intramolecular phenomena (thermochemistry). This Review describes the recent developments (including some historical aspects) of dispersion corrections with an emphasis on methods that can be employed routinely with reasonable accuracy in large-scale applications. The most prominent correction schemes are classified into three groups: (i) nonlocal, density-based functionals, (ii) semiclassical C6-based, and (iii) one-electron effective potentials. The properties as well as pros and cons of these methods are critically discussed, and typical examples and benchmarks on molecular complexes and crystals are provided. Although there are some areas for further improvement (robustness, many-body and short-range effects), the situation regarding the overall accuracy is clear. Various approaches yield long-range dispersion energies with a typical relative error of 5%. For many chemical problems, this accuracy is higher compared to that of the underlying mean-field method (i.e., a typical semilocal (hybrid) functional like B3LYP).

First-principles Study of Electronic and Magnetic Properties of Manganese Decorated Graphene

<p>The functionalization of graphene by the addition of Manganese (Mn) atom to its surface has been investigated computationally by using density-functional theory (DFT) based first-principles method within DFT-D2 level of approximations. The calculations have been computed employing the Quantum ESPRESSO codes. The stability, geometrical structures, electronic and magnetic properties of pure and Mn adatom graphene systems have been studied. From the information of adsorption energies of Mn atom on the different sites of graphene, the top site is found to be the most favorable one for its adsorption. Present study finds that the London dispersion interaction plays a major role in the weak binding of Mn on graphene. The study of electronic and magnetic properties of Mn decorated graphene shows that the conduction and valence band are overlapped with finite density of states (DOS) at Fermi level. The dissimilar DOS for up and down spin calculations quantify magnetic moment as 5.48 μ<sub>B</sub> which is consistent with the previous study</p><p>Journal of Nepal Physical Society Vol.3(1) 2015: 24-34</p>

Small Atomic Orbital Basis Set First-Principles Quantum Chemical Methods for Large Molecular and Periodic Systems: A Critical Analysis of Error Sources.

In quantum chemical computations the combination of Hartree–Fock or a density functional theory (DFT) approximation with relatively small atomic orbital basis sets of double‐zeta quality is still widely used, for example, in the popular B3LYP/6‐31G* approach. In this Review, we critically analyze the two main sources of error in such computations, that is, the basis set superposition error on the one hand and the missing London dispersion interactions on the other. We review various strategies to correct those errors and present exemplary calculations on mainly noncovalently bound systems of widely varying size. Energies and geometries of small dimers, large supramolecular complexes, and molecular crystals are covered. We conclude that it is not justified to rely on fortunate error compensation, as the main inconsistencies can be cured by modern correction schemes which clearly outperform the plain mean‐field methods.

Treating London-Dispersion Effects with the Latest Minnesota Density Functionals: Problems and Possible Solutions.

It is shown that the latest Minnesota density functionals (SOGGA11, M11-L, N12, MN12-L, SOGGA11-X, M11, N12-SX, and MN12-SX) do not properly describe London-dispersion interactions. Grimme's DFT-D3 correction can solve this problem partially; however, double-counting of medium-range electron correlation can occur. For the related M06-L functional, the alternative VV10 van der Waals kernel is tested, but it experiences similar double-counting. Most functionals give unphysical dissociation curves for the argon dimer, an indication for method-inherent problems, and further investigation is recommended. These results are further evidence that the London-dispersion problem in density functional theory approximations is unlikely to be solved by mere empirical optimization of functional parameters, unless the functionals contain components that ensure the correct asymptotic long-range behavior. London dispersion is ubiquitous, which is why the reported findings are not only important for theoreticians but also a reminder to the general chemist to carefully consider their choice of method before undertaking computational studies.

Graphene oxide as an optimal candidate material for methane storage.

Methane, the primary constituent of natural gas, binds too weakly to nanostructured carbons to meet the targets set for on-board vehicular storage to be viable. We show, using density functional theory calculations, that replacing graphene by graphene oxide increases the adsorption energy of methane by 50%. This enhancement is sufficient to achieve the optimal binding strength. In order to gain insight into the sources of this increased binding, that could also be used to formulate design principles for novel storage materials, we consider a sequence of model systems that progressively take us from graphene to graphene oxide. A careful analysis of the various contributions to the weak binding between the methane molecule and the graphene oxide shows that the enhancement has important contributions from London dispersion interactions as well as electrostatic interactions such as Debye interactions, aided by geometric curvature induced primarily by the presence of epoxy groups.