Density-functional theory-symmetry-adapted intermolecular perturbation theory with density fitting: A new efficient method to study intermolecular interaction energies

Physisorption of Water on Graphene: Subchemical Accuracy from Many-Body Electronic Structure Methods

Wet carbon interfaces are ubiquitous in the natural world and exhibit anomalous properties, which could be exploited by emerging technologies. However, progress is limited by lack of understanding at the molecular level. Remarkably, even for the most fundamental system (a single water molecule interacting with graphene), there is no consensus on the nature of the interaction. We tackle this by performing an extensive set of complementary state-of-the-art computer simulations on some of the world's largest supercomputers. From this effort a consensus on the water-graphene interaction strength has been obtained. Our results have significant impact for the physical understanding, as they indicate that the interaction is weaker than predicted previously. They also pave the way for more accurate and reliable studies of liquid water at carbon interfaces.

Halogen bonds and σ-holes

The models behind simple bonding theory and the origins of some components often proposed to be involved in weak intermolecular bonds are described with special reference to σ-hole bonding, of which halogen bonds are a subset. A protocol for the analysis of weak intermolecular interactions is proposed on the basis of sound physical principles. This protocol uses three different levels of interaction; "permanent" Coulomb interactions between unperturbed monomers, relaxed Coulomb interactions and dispersion. Of the three, only dispersion is not a real, measurable quantity. It is, however, included in order to describe interactions that cannot be treated entirely by the first two levels.

Adhesion, forces and the stability of interfaces

Weak molecular interactions (WMI) are responsible for processes such as physisorption; they are essential for the structure and stability of interfaces, and for bulk properties of liquids and molecular crystals. The dispersion interaction is one of the four basic interactions types – electrostatics, induction, dispersion and exchange repulsion – of which all WMIs are composed. The fact that each class of basic interactions covers a wide range explains the large variety of WMIs. To some of them, special names are assigned, such as hydrogen bonding or hydrophobic interactions. In chemistry, these WMIs are frequently used as if they were basic interaction types. For a long time, dispersion was largely ignored in chemistry, attractive intermolecular interactions were nearly exclusively attributed to electrostatic interactions. We discuss the importance of dispersion interactions for the stabilization in systems that are traditionally explained in terms of the “special interactions” mentioned above. System stabilization can be explained by using interaction energies, or by attractive forces between the interacting subsystems; in the case of stabilizing WMIs, one frequently speaks of adhesion energies and adhesive forces. We show that the description of system stability using maximum adhesive forces and the description using adhesion energies are not equivalent. The systems discussed are polyaromatic molecules adsorbed to graphene and carbon nanotubes; dimers of alcohols and amines; cellulose crystals; and alcohols adsorbed onto cellulose surfaces.

Understanding non-covalent interactions in larger molecular complexes from first principles

Non-covalent interactions pervade all matter and play a fundamental role in layered materials, biological systems, and large molecular complexes. Despite this, our accumulated understanding of non-covalent interactions to date has been mainly developed in the tens-of-atoms molecular regime. This falls considerably short of the scales at which we would like to understand energy trends, structural properties, and temperature dependencies in materials where non-covalent interactions have an appreciable role. However, as more reference information is obtained beyond moderately sized molecular systems, our understanding is improving and we stand to gain pertinent insights by tackling more complex systems, such as supramolecular complexes, molecular crystals, and other soft materials. In addition, accurate reference information is needed to provide the drive for extending the predictive power of more efficient workhorse methods, such as density functional approximations that also approximate van der Waals dispersion interactions. In this perspective, we discuss the first-principles approaches that have been used to obtain reference interaction energies for beyond modestly sized molecular complexes. The methods include quantum Monte Carlo, symmetry-adapted perturbation theory, non-canonical coupled cluster theory, and approaches based on the random-phase approximation. By considering the approximations that underpin each method, the most accurate theoretical references for supramolecular complexes and molecular crystals to date are ascertained. With these, we also assess a handful of widely used exchange-correlation functionals in density functional theory. The discussion culminates in a framework for putting into perspective the accuracy of high-level wavefunction-based methods and identifying future challenges.

Spectroscopy of a rotating hydrogen molecule in carbon nanotubes

Computing the energy levels of molecular hydrogen rotating in carbon nanotubes of increasing size.

Molecular dimers of methane clathrates: ab initio potential energy surfaces and variational vibrational states

Motivated by the energetic and environmental relevance of methane clathrates, highly accurate ab initio potential energy surfaces (PESs) have been developed for the three possible dimers of the methane and water molecules: (H₂O)₂, CH₄·H₂O, and (CH₄)₂.

Noble gas dimers confined inside C70

The potential energy surfaces for the interior rotation of a series of pairs of noble gas atoms encapsulated in the C₇₀ cavity have been explored.

Quantum-Mechanical Relation between Atomic Dipole Polarizability and the van der Waals Radius

The atomic dipole polarizability α and the van der Waals (vdW) radius R_{vdW} are two key quantities to describe vdW interactions between atoms in molecules and materials. Until now, they have been determined independently and separately from each other. Here, we derive the quantum-mechanical relation R_{vdW}=const×α^{1/7}, which is markedly different from the common assumption R_{vdW}∝α^{1/3} based on a classical picture of hard-sphere atoms. As shown for 72 chemical elements between hydrogen and uranium, the obtained formula can be used as a unified definition of the vdW radius solely in terms of the atomic polarizability. For vdW-bonded heteronuclear dimers consisting of atoms A and B, the combination rule α=(α_{A}+α_{B})/2 provides a remarkably accurate way to calculate their equilibrium interatomic distance. The revealed scaling law allows us to reduce the empiricism and improve the accuracy of interatomic vdW potentials, at the same time suggesting the existence of a nontrivial relation between length and volume in quantum systems.

Do Semilocal Density-Functional Approximations Recover Dispersion Energies at Small Intermonomer Separations?

The methods that add dispersion energies to interaction energies computed using density-functional theory (DFT), known as DFT+D methods, taper off the dispersion energies at distances near van der Waals minima and smaller based on an assumption that DFT starts to reproduce the dispersion energies there. We show that this assumption is not correct as the alleged contribution behaves unphysically and originates to a large extent from nonexchange-correlation terms. Thus, dispersion functions correct DFT in this region for deficiencies unrelated to dispersion interactions.

On the Nature of σ-σ, σ-π, and π-π Stacking in Extended Systems

Stacking interactions have been evaluated, employing computational methods, in dimers formed by analogous aliphatic and aromatic species of increasing size. Changes in stability as the systems become larger are mostly controlled by the balance of increasing repulsion and dispersion contributions, while electrostatics plays a secondary but relevant role. The interaction energy increases as the size of the system grows, but it does much faster in π-π dimers than in σ-π complexes and more remarkably than in σ-σ dimers. The main factor behind the larger stability of aromatic dimers compared to complexes containing aliphatic molecules is related to changes in the properties of the aromatic systems due to electron delocalization leading to larger dispersion contributions. Besides, an extra stabilization in π-π complexes is due to the softening of the repulsive wall in aromatic species that allows the molecules to come closer.

Dispersion interactions between neighboring Bi atoms in (BiH3 )2 and Te(BiR2 )2

Triggered by the observation of a short Bi⋯Bi distance and a BiTeBi bond angle of only 86.6° in the crystal structure of bis(diethylbismuthanyl)tellurane quantum chemical computations on interactions between neighboring Bi atoms in Te(BiR₂ )₂ molecules (R = H, Me, Et) and in (BiH₃ )₂ were undertaken. Bi⋯Bi distances atoms were found to significantly shorten upon inclusion of the d shells of the heavy metal atoms into the electron correlation treatment, and it was confirmed that interaction energies from spin component-scaled second-order Møller-Plesset theory (SCS-MP2) agree well with coupled-cluster singles and doubles theory including perturbative triples (CCSD(T)). Density functional theory-based symmetry-adapted perturbation theory (DFT-SAPT) was used to study the anisotropy of the interplay of dispersion attraction and steric repulsion between the Bi atoms. Finally, geometries and relative stabilities of syn-syn and syn-anti conformers of Te(BiR₂ )₂ (R = H, Me, Et) and interconversion barriers between them were computed. © 2018 Wiley Periodicals, Inc.

Ab initio design of light absorption through silver atomic cluster decoration of TiO2

A first-principles study of the stability and optical response of subnanometer silver clusters Agn (n ≤ 5) on a TiO2(110) surface is presented. First, the adequacy of the vdW-corrected DFT-D3 approach is assessed using the domain-based pair natural orbital correlation DLPNO-CCSD(T) calculations along with the Symmetry-Adapted Perturbation Theory [SAPT(DFT)] applied to a cluster model. Next, using the DFT-D3 treatment with a periodic slab model, we analyze the interaction energies of the atomic silver clusters with the TiO2(110) surface. Finally, the hybrid HSE06 functional and a reduced density matrix treatment are applied to obtain the projected electronic density of states and photo-absorption spectra of the TiO2(110) surface, with and without adsorbed silver clusters. Our results show the stability of the supported clusters, the enhanced light absorbance intensity of the material upon their deposition, and the appearance of intense secondary broad peaks in the near-infrared and the visible regions of the spectrum, with positions depending on the size and shape of the supported clusters. The secondary peaks arise from the photo-induced transfer of electrons from intra-band valence 5s orbitals of the noble-metal cluster to 3d Ti band states of the supporting material.

Comparison of the DFT-SAPT and Canonical EDA Schemes for the Energy Decomposition of Various Types of Noncovalent Interactions

Interaction energies computed with density functional theory can be divided into physically meaningful components by symmetry-adapted perturbation theory (DFT-SAPT) or the canonical energy decomposition analysis (EDA). In this work, the decomposition results obtained by these schemes were compared for more than 200 hydrogen-, halogen-, and pnicogen-bonded, dispersion-bound, and mixed complexes to investigate their similarity in the evaluation of the nature of noncovalent interactions. BLYP functional with D3(BJ) correction was used for the EDA scheme, whereas asymptotically corrected PBE0 functional for DFT-SAPT provided some of the best combinations for description of noncovalent interactions. Both schemes provide similar results concerning total interaction energies and insight into the individual energy components. For most complexes, the dominant energetic term was identified equally by both decomposition schemes. Because the canonical EDA is computationally less demanding than the DFT-SAPT, the former can be especially used in cases where the systems investigated are very large.

The Nature of Chemical Bonding in Lewis Adducts as Reflected by 27Al NMR Quadrupolar Coupling Constant: Combined Solid-State NMR and Quantum Chemical Approach

Lewis acids and Lewis adducts are widely used in the chemical industry because of their high catalytic activity. Their precise geometrical description and understanding of their electronic structure are a crucial step for targeted synthesis and specific use. Herein, we present an experimental/computational strategy based on a solid-state NMR crystallographic approach allowing for detailed structural characterization of a wide range of organoaluminum compounds considerably differing in their chemical constitution. In particular, we focus on the precise measurement and subsequent quantum-chemical analysis of many different ²⁷Al NMR resonances in the extremely broad range of quadrupolar coupling constants from 1 to 50 MHz. In this regard, we have optimized an experimental strategy combining a range of static as well as magic angle spinning experiments allowing reliable detection of the entire set of aluminum sites present in trimesitylaluminum (AlMes₃) reaction products. In this way, we have spectroscopically resolved six different products in the resulting polycrystalline mixture. All ²⁷Al NMR resonances are precisely recorded and comprehensively analyzed by a quantum-chemical approach. Interestingly, in some cases the recorded ²⁷Al solid-state NMR spectra show unexpected quadrupolar coupling constant values reaching up to ca. 30 MHz, which are attributed to tetra-coordinated aluminum species (Lewis adducts with trigonal pyramidal geometry). The cause of this unusual behavior is explored by analyzing the natural bond orbitals and complexation energies. The linear correlation between the quadrupolar coupling constant value and the nature of bonds in the Lewis adducts is revealed. Moreover, the ²⁷Al NMR data are shown to be sensitive to the geometry of the tetra-coordinated organoaluminum species. Our findings thus provide a viable approach for the direct identification of Lewis acids and Lewis adducts, not only in the investigated multicomponent organoaluminum compounds but also in inorganic zeolites featuring catalytically active trigonal (Al^III) and strongly perturbed Al^IV sites.

First-order symmetry-adapted perturbation theory for multiplet splittings

We present a symmetry-adapted perturbation theory (SAPT) for the interaction of two high-spin open-shell molecules (described by their restricted open-shell Hartree-Fock determinants) resulting in low-spin states of the complex. The previously available SAPT formalisms, except for some system-specific studies for few-electron complexes, were restricted to the high-spin state of the interacting system. Thus, the new approach provides, for the first time, a SAPT-based estimate of the splittings between different spin states of the complex. We have derived and implemented the lowest-order SAPT term responsible for these splittings, that is, the first-order exchange energy. We show that within the so-called S² approximation commonly used in SAPT (neglecting effects that vanish as fourth or higher powers of intermolecular overlap integrals), the first-order exchange energies for all multiplets are linear combinations of two matrix elements: a diagonal exchange term that determines the spin-averaged effect and a spin-flip term responsible for the splittings between the states. The numerical factors in this linear combination are determined solely by the Clebsch-Gordan coefficients: accordingly, the S² approximation implies a Heisenberg Hamiltonian picture with a single coupling strength parameter determining all the splittings. The new approach is cast into both molecular-orbital and atomic-orbital expressions: the latter enable an efficient density-fitted implementation. We test the newly developed formalism on several open-shell complexes ranging from diatomic systems (Li⋯H, Mn⋯Mn, …) to the phenalenyl dimer.

Perspective: Found in translation: Quantum chemical tools for grasping non-covalent interactions

Today's quantum chemistry methods are extremely powerful but rely upon complex quantities such as the massively multidimensional wavefunction or even the simpler electron density. Consequently, chemical insight and a chemist's intuition are often lost in this complexity leaving the results obtained difficult to rationalize. To handle this overabundance of information, computational chemists have developed tools and methodologies that assist in composing a more intuitive picture that permits better understanding of the intricacies of chemical behavior. In particular, the fundamental comprehension of phenomena governed by non-covalent interactions is not easily achieved in terms of either the total wavefunction or the total electron density, but can be accomplished using more informative quantities. This perspective provides an overview of these tools and methods that have been specifically developed or used to analyze, identify, quantify, and visualize non-covalent interactions. These include the quantitative energy decomposition analysis schemes and the more qualitative class of approaches such as the Non-covalent Interaction index, the Density Overlap Region Indicator, or quantum theory of atoms in molecules. Aside from the enhanced knowledge gained from these schemes, their strengths, limitations, as well as a roadmap for expanding their capabilities are emphasized.

Correlation effects and many-body interactions in water clusters

Background: The quantum-chemical description of the interactions in water clusters is an essential basis for deriving accurate and physically sound models of the interaction potential for water to be used in molecular simulations. In particular, the role of many-body interactions beyond the two-body interactions, which are often not explicitly taken into account by empirical force fields, can be accurately described by quantum chemistry methods on an adequate level, e.g., random-phase approximation electron correlation methods. The relative magnitudes of the different interaction energy contributions obtained by accurate ab initio calculations can therefore provide useful insights that can be exploited to develop enhanced force field methods. Results: In line with earlier theoretical studies of the interactions in water clusters, it has been found that the main contribution to the many-body interactions in clusters with a size of up to N = 13 molecules are higher-order polarisation interaction terms. Compared to this, many-body dispersion interactions are practically negligible for all studied sytems. The two-body dispersion interaction, however, plays a significant role in the formation of the structures of the water clusters and their stability, since it leads to a distinct compression of the cluster sizes compared to the structures optimized on an uncorrelated level. Overall, the many-body interactions amount to about 13% of the total interaction energy, irrespective of the cluster size. The electron correlation contribution to these, however, amounts to only about 30% to the total many-body interactions for the largest clusters studied and is repulsive for all structures considered in this work. Conclusion: While this shows that three- and higher-body interactions can not be neglected in the description of water complexes, the electron correlation contributions to these are much smaller in comparison to the two-body electron correlation effects. Efficient quantum chemistry approaches for describing intermolecular interactions between water molecules may therefore describe higher-body interactions on an uncorrelated Hartree-Fock level without a serious loss in accuracy.

Ab Initio Force Fields for Organic Anions: Properties of [BMIM][TFSI], [BMIM][FSI], and [BMIM][OTf] Ionic Liquids

Room-temperature ionic liquids (ILs) composed of organic anions bis(trifluoromethanesulfonyl)imide (TFSI), bis(fluorosulfonyl)imide (FSI), and trifluoromethanesulfonate (OTf) exhibit interesting physical properties and are important for many electrochemical applications. TFSI and FSI form "hydrophobic" ILs, immiscible with water but miscible with many organic solvents and polymers; for computer simulation studies, it is thus essential to develop force fields for these anions that are transferable among this wide variety of chemical environments. In this work, we develop entirely ab initio force fields for the TFSI, FSI, and OTf anions and predict the properties of corresponding 1-butyl-3-methylimidazolium ILs. We discuss important subtleties in the force field development related to accurately modeling conformational flexibility, that is, relaxed torsional profiles and intramolecular electrostatic interactions. The TFSI anions have notable conformational flexibility in the IL, and we predict approximately 70% cisoid and 20% transoid conformations, which is largely driven by cation/anion ion-pair interactions and is opposite to the trend expected from the anion ab initio potential energy surface. The favorable interactions between the cation and cisoid TFSI conformations result in a shoulder in the cation/anion radial distribution function at short distances, whereas interconversion between cisoid and transoid conformations occurs on a commensurate time scale as ion diffusion processes. In addition to this physical insight on anion effects, we expect that these force fields will have important applications for studying a variety of complex electrolyte systems.

Physisorption of H2 on Fullerenes and the Solvation of C60 by Hydrogen Clusters at Finite Temperature: A Theoretical Assessment

The interaction between hydrogen and carbonaceous nanostructures is of fundamental interest in various areas of physical chemistry. In this contribution we have revisited the physisorption of hydrogen molecules and H₂ clusters on fullerenes, following a first-principles approach in which the interaction is quantitatively evaluated for the C₂₀ system using high-level electronic structure methods. Relative to coupled cluster data at the level of single, double, and perturbative triple excitations taken as a benchmark, the results for rotationally averaged physisorbed H₂ show a good performance of MP2 variants and symmetry-adapted perturbation theory, but significant deviations and basis set convergence issues are found for dispersion-corrected density functional theory. These electronic structure data are fitted to produce effective coarse-grained potentials for use in larger systems such as C₆₀-H₂. Using path-integral molecular dynamics, the potentials are also applied to parahydrogen clusters solvated around fullerenes, across the regime where the first solvation shell becomes complete and as a function of increasing temperature. For C₆₀ our findings indicate a sensible dependence of the critical solvation size on the underlying potential. As the temperature is increased, a competition is found between the surface and radial expansions of the solvation shell, with one molecule popping away at intermediate temperatures but getting reinserted at even higher temperatures.

Extension of the interacting quantum atoms (IQA) approach to B3LYP level density functional theory (DFT)

An interaction between two atoms, bonded or non-bonded, consists of interatomic contributions: electrostatic energy, exchange energy and electronic correlation energy. Together with the intra-atomic energy of an atom, these contributions are the basic components of the Interacting Quantum Atom (IQA) energy decomposition scheme. Here, we investigate IQA's proper use in conjunction with an explicit implementation of the B3LYP functional. The recovery of the total molecular energy from the IQA components is emphasised, for the first time. A systematic study of three model systems of biological relevance, N-methylacetamide (NMA), the doubly capped tripeptide GlyGlyGly and an alloxan dimer, shows the stabilization effect of B3LYP on most of the interatomic exchange energies (V) compared to their Hartree-Fock values. Diagrams of exchange energies versus interatomic distance show the clustering of interactions, one cluster for each 1,n (n = 1 to 6 where the atoms are separated by n - 1 bonds). The positioning of some V values outside their expected cluster marks interesting interactions.

Stability of melamine-exfoliated graphene in aqueous media: quantum-mechanical insights at the nanoscale

In recent experiments, melamine (1,3,5-triazine-2,4,6-triamine) has been proposed as an effective exfoliating agent to obtain high quality graphene from graphite. After washing out the melamine in excess, small amounts (ppm) are still needed to stabilize the dispersion of graphene flakes in aqueous media. To understand the origin of this behaviour, we investigated the melamine-graphene-water system and the fundamental interactions that determine its structure and energetics. To disentangle the subtle interplay of hydrogen-bonding and dispersive forces we used state-of-the-art ab initio calculations based on density functional theory. First, we focused on the case of water molecules interacting with melamine-graphene assemblies at different melamine coverages. We found that water-melamine interactions provide the driving force for washing off the melamine from graphene. Then, we addressed the interaction of single and double layers of water molecules with the graphene surface in the presence of an adsorbed melamine molecule. We found that this melamine acts as a non-covalent anchor for keeping a number of water molecules conveniently close to the graphene surface, thus helping its stabilization in aqueous media. Our analysis helps understanding how competing weak forces can lead to a stable graphene water suspension thanks to small amounts of adsorbed melamine. From our results, we derive simple indications on how the water-graphene interfacial properties can be tuned via non-covalent adsorption of small functional molecules with H-bond donor/acceptor groups. These new hints can be helpful to prepare stable graphene dispersions in water and so to unlock graphene potential in aqueous environments.

Dissecting the concave-convex π-π interaction in corannulene and sumanene dimers: SAPT(DFT) analysis and performance of DFT dispersion-corrected methods

The characteristics of the concave-convex π-π interactions are evaluated in 32 buckybowl dimers formed by corannulene, sumanene, and two substituted sumanenes (with S and CO groups), using symmetry-adapted perturbation theory [SAPT(DFT)] and density functional theory (DFT). According to our results, the main stabilizing contribution is dispersion, followed by electrostatics. Regarding the ability of DFT methods to reproduce the results obtained with the most expensive and rigorous methods, TPSS-D seems to be the best option overall, although its results slightly tend to underestimate the interaction energies and to overestimate the equilibrium distances. The other two tested DFT-D methods, B97-D2 and B3LYP-D, supply rather reasonable results as well. M06-2X, although it is a good option from a geometrical point of view, leads to too weak interactions, with differences with respect to the reference values amounting to about 4 kcal/mol (25% of the total interaction energy). © 2017 Wiley Periodicals, Inc.

Noncovalent interactions between cisplatin and graphene prototypes

Cisplatin (CP) has been widely used as an anticancer drug for more than 30 years despite severe side effects due to its low bioavailability and poor specificity. For this reason, it is paramount to study and design novel nanomaterials to be used as vectors capable to effectively deliver the drug to the biological target. The CP square-planar geometry, together with its low water solubility, suggests that it could be possibly easily adsorbed on 2D graphene nanostructures through the interaction with the related highly conjugated π-electron system. In this work, pyrene has been first selected as the minimum approximation to the graphene plane, which allows to properly study the noncovalent interactions determining the CP adsorption. In particular, electronic structure calculations at the MP2C and DFT-SAPT levels of theory have allowed to obtain benchmark interaction energies for some limiting configurations of the CP-pyrene complex, as well as to assess the role of the different contributions to the total interaction: it has been found that the parallel configurations of the aggregate are mainly stabilized around the minimum region by dispersion, in a similar way as for complexes bonded through π-π interactions. Then, the benchmark interaction energies have been used to test corresponding estimations obtained within the less expensive DFT to validate an optimal exchange-correlation functional which includes corrections to take properly into account for the dispersion contribution. Reliable DFT interaction energies have been therefore obtained for CP adsorbed on graphene prototypes of increasing size, ranging from coronene, ovalene, and up to C₁₅₀ H₃₀ . Finally, DFT geometry optimizations and frequency calculations have also allowed a reliable estimation of the adsorption enthalpy of CP on graphene, which is found particularly favorable (about -20 kcal/mol at 298 K and 1 bar) being twice that estimated for the corresponding benzene adsorption. © 2017 Wiley Periodicals, Inc.

The coulombic σ-hole model describes bonding in CX3I⋯Y− complexes completely

Contrary to recent reports, the σ-hole interaction energies of complexes between the carbon tetrahalides CX₃I (X = F, Cl, Br, I) and halide anions Y⁻ (Y = F, Cl, Br, I) are described very well by the simple Coulombic σ-hole concept if it is applied properly.

The role of the σ-holes in stability of non-bonded chalcogenide⋯benzene interactions: the ground and excited states

Different characters of molecular electrostatic potential (MEP) in the ground and excited states of chalcogenides are responsible for changes in conformer stability of T-shape and stacked non-bonded chalcogenide–benzene complexes upon electronic excitation.

Tetracyclo(9-methyl-2,7-carbazole) as a promising nanohoop for gas trapping: a multiscale study

H₂S, CS₂, NO₂, Br₂, HF, and C₂H₆ are the ideal adsorbates within the TCC host from their respective congeners.

Symmetry-adapted perturbation theory potential for the adenine dimer

A new ab initio intermolecular interaction potential for the adenine dimer has been developed.

Intramolecular interactions in sterically crowded hydrocarbon molecules

A molecular fragmentation method has been used to analyze the intramolecular interactions in the three molecules coupled diamantane, hexaphenylethane, and all-meta-tert-butyl substituted hexaphenylethane. The significance of these systems lies in the fact, that steric crowding effects enable a stabilization of the central carbon bond that possesses an extended length (1.6 to 1.7 Å) beyond conventional carbon-carbon bonds due to the steric repulsion of the attached hydrocarbon groups. The total stability of these molecules therefore depends on a delicate balance between attractive interaction forces on the one hand and on repulsive forces on the other hand. We have quantified the different interaction energy contributions using symmetry-adapted perturbation theory based on a density functional theory description of the monomers. It has been found that the attractive dispersion interactions increase more strongly with the level of crowding in the systems than the counteracting exchange interactions. This shows that steric crowding effects can have a significant impact on the structure and stability of large and branched molecules. © 2017 Wiley Periodicals, Inc.