σ/σ- and π/π-Interactions Are Equally Important: Multilayered Graphanes

Assessing the domain-based local pair natural orbital (DLPNO) approximation for non-covalent interactions in sizable supramolecular complexes

The titular domain-based local pair natural orbital (DLPNO) approximation is the most widely used method for extending correlated wave function models to large molecular systems, yet its fidelity for intermolecular interaction energies in large supramolecular complexes has not been thoroughly vetted. Non-covalent interactions are sensitive to tails of the electron density and involve nonlocal dispersion that is discarded or approximated if the screening of pair natural orbitals (PNOs) is too aggressive. Meanwhile, the accuracy of the DLPNO approximation is known to deteriorate as molecular size increases. Here, we test the DLPNO approximation at the level of second-order Møller-Plesset perturbation theory (MP2) and coupled-cluster theory with singles, doubles, and perturbative triples [CCSD(T)] for a variety of large supramolecular complexes. DLPNO-MP2 interaction energies are within 3% of canonical values for small dimers with ≲10 heavy atoms, but for larger systems, the DLPNO approximation is often quite poor unless the results are extrapolated to the canonical limit where the threshold for discarding PNOs is taken to zero. Counterpoise correction proves to be essential in reducing errors with respect to canonical results. For a sequence of nanoscale graphene dimers up to (C96H24)2, extrapolated DLPNO-MP2 interaction energies agree with canonical values to within 1%, independent of system size, provided that the basis set does not contain diffuse functions; these cause the DLPNO approximation to behave erratically, such that results cannot be extrapolated in a meaningful way. DLPNO-CCSD(T) calculations are typically performed using looser PNO thresholds as compared to DLPNO-MP2, but this significantly impacts accuracy for large supramolecular complexes. Standard DLPNO-CCSD(T) settings afford errors of 2-6 kcal/mol for dimers involving coronene (C24H12) and circumcoronene (C54H18), even at the DLPNO-CCSD(T1) level.

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.

σ-σ Stacked supramolecular junctions

Intermolecular charge transport plays an essential role in organic electronic materials and biological systems. To date, experimental investigations of intermolecular charge transport in molecular materials and electronic devices have been restricted to conjugated systems in which π-π stacking interactions are involved. Herein we demonstrate that the σ-σ stacking interactions between neighbouring non-conjugated molecules offer an efficient pathway for charge transport through supramolecular junctions. The conductance of σ-σ stacked molecular junctions formed between two non-conjugated cyclohexanethiol or single-anchored adamantane molecules is comparable to that of π-π stacked molecular junctions formed between π-conjugated benzene rings. The current-voltage characteristics and flicker noise analysis demonstrate the existence of stacked molecular junctions formed between the electrode pairs and exhibit the characteristics of through-space charge transport. Density functional theory calculations combined with the non-equilibrium Green's function method reveal that efficient charge transport occurs between two molecules configured with σ-σ stacking interactions.

London Dispersion Favors Sterically Hindered Diarylthiourea Conformers in Solution

We present an experimental and computational study on the conformers of N,N'-diphenylthiourea substituted with different dispersion energy donor (DED) groups. While the unfolded anti-anti conformer is the most relevant for thiourea catalysis, intramolecular noncovalent interactions counterintuitively favor the folded syn-syn conformer, as evident from a combination of low-temperature nuclear magnetic resonance measurements and computations. In order to quantify the noncovalent interactions, we utilized local energy decomposition analysis and symmetry-adapted perturbation theory at the DLPNO-CCSD(T)/def2-TZVPP and sSAPT0/6-311G(d,p) levels of theory. Additionally, we applied a double-mutant cycle to experimentally study the effects of bulky substituents on the equilibria. We determined London dispersion as the key interaction that shifts the equilibria towards the syn-syn conformers. This preference is likely a factor why such thiourea derivatives can be poor catalysts.

A Minimum Quantum Chemistry CCSD(T)/CBS Data Set of Dimeric Interaction Energies for Small Organic Functional Groups: Heterodimers

We extend our previous quantum chemistry calculations of interaction energies for 31 homodimers of small organic functional groups (the SOFG-31 data set) by including 239 heterodimers with monomers selected within the SOFG-31 data set, thus resulting in the SOFG-31+239 data set. The minimum-level theoretical scheme contains (1) the basis set superposition error corrected supermolecule (BSSE-SM) approach for intermolecular interactions; (2) the second-order Møller-Plesset perturbation theory (MP2) with the Dunning's aug-cc-pVXZ (X = D, T, Q) basis sets for the geometry optimization and correlation energy calculations; and (3) the single-point energy calculations with the coupled cluster with single, double, and perturbative triple excitations method at the complete basis set limit [CCSD(T)/CBS] using the well-tested extrapolation methods for the MP2 energy calibrations. In addition, we have performed a parallel series of energy decomposition calculations based on the symmetry adapted perturbation theory (SAPT) in order to gain chemical insights. That the above procedure cannot be further reduced has been proven to be very crucial for constructing reliable data sets of interaction energies. The calculated CCSD(T)/CBS interaction energy data can serve as a benchmark for testing or training less accurate but more efficient calculation methods, such as the electronic density functional theory. As an application, we employ a segmental SAPT model previously developed for the SOFG-31 data set to predict binding energies of large heterodimer complexes. These model energy "quanta" can be used in coarse-grained molecular dynamics simulations by avoiding large-scale calculations.

3× Axial vs 3× Equatorial: The ΔGGA Value Is a Robust Computational Measure of Substituent Steric Effects

We define ΔG_GA as the free energy change for the formal equilibrium: [13]G-H + 1-X-adamantane → [13]G-X + adamantane, where [13]G-H is the C₁₃H₂₂ fragment of all-trans graphane with 3-fold symmetry. This compares with a situation where the group X is equatorial to three cyclohexane rings with one where it is axial to three rings. ΔG_GA values vary from 2.9 (CN) to 145.7 kJ mol^-1 (CCl₃), and this wide range means that ΔG can be calculated with confidence. ΔG_GA values for Me, Et, i-Pr, and t-Bu form a regular series, 34.9, 63.3, 101.6, and 142.0, and clearly reflect the steric size of the groups. We propose a model where the six axial hydrogens surrounding X on [13]G-X provide a nearly circular constriction on the substituent close to its point of attachment but which does not extend far above this. We compare these results with A values and with calculations on 2- and 7-substituted [1(2,3)4]pentamantanes. We show that electronic effects on ΔG_GA values are negligible but that they correlate well with computed cone and solid angles subtended by the substituent.

Mechanistic Insights into Activation of Carbon Monoxide, Carbon Dioxide, and Nitrous Oxide by Acyclic Silylene

Owing to an empty p orbital and a lone pair of electrons on the Si center, silylene exhibits reactivity similar to a transition-metal system capable of activating H₂/C-H bonds and small molecules. In this work, with the aid of density functional theory calculations, we systematically investigated the reactions of an acyclic silylene with CO, CO₂, and N₂O. The detailed mechanisms obtained lead to an in-depth understanding of the silylene single-site ambiphilic reactivity.

Intramolecular London Dispersion Interactions Do Not Cancel in Solution

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

A minimum quantum chemistry CCSD(T)/CBS dataset of dimeric interaction energies for small organic functional groups

We have performed a quantum chemistry study on the bonding patterns and interaction energies for 31 dimers of small organic functional groups (dubbed the SOFG-31 dataset), including the alkane-alkene-alkyne (6 + 4 + 4 = 14, AAA) groups, alcohol-aldehyde-ketone (4 + 4 + 3 = 11, AAK) groups, and carboxylic acid-amide (3 + 3 = 6, CAA) groups. The basis set superposition error corrected super-molecule approach using the second order Møller-Plesset perturbation theory (MP2) with the Dunning's aug-cc-pVXZ (X = D, T, Q) basis sets has been employed in the geometry optimization and energy calculations. To calibrate the MP2 calculated interaction energies for these dimeric complexes, we perform single-point calculations with the coupled cluster with single, double, and perturbative triple excitations method at the complete basis set limit [CCSD(T)/CBS] using the well-tested extrapolation methods. In order to gain more physical insights, we also perform a parallel series of energy decomposition calculations based on the symmetry adapted perturbation theory (SAPT). The collection of these CCSD(T)/CBS interaction energy values can serve as a minimum quantum chemistry dataset for testing or training less accurate but more efficient calculation methods. As an application, we further propose a segmental SAPT model based on chemically recognizable segments in a specific functional group. These model interactions can be used to construct coarse-grained force fields for larger molecular systems.

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.

Carbocatalytic Acetylene Cyclotrimerization: A Key Role of Unpaired Electron Delocalization

Development of sustainable catalysts for synthetic transformations is one of the most challenging and demanding goals. The high prices of precious metals and the unavoidable leaching of toxic metal species leading to environmental contamination make the transition metal-free catalytic systems especially important. Here we demonstrate that carbene active centers localized on carbon atoms at the zigzag edge of graphene represent an alternative platform for efficient catalytic carbon-carbon bond formation in the synthesis of benzene. The studied acetylene trimerization reaction is an efficient atom-economic route to build an aromatic ring-a step ubiquitously important in organic synthesis and industrial applications. Computational modeling of the reaction mechanism reveals a principal role of the reversible spin density oscillations that govern the overall catalytic cycle, facilitate the product formation, and regenerate the catalytically active centers. Dynamic π-electron interactions in 2D carbon systems open new opportunities in the field of carbocatalysis, unachievable by means of transition metal-catalyzed transformations. The theoretical findings are confirmed experimentally by generating key moieties of the carbon catalyst and performing the acetylene conversion to benzene.

Intermolecular London Dispersion Interactions of Azobenzene Switches for Tuning Molecular Solar Thermal Energy Storage Systems

The performance of molecular solar thermal energy storage systems (MOST) depends amongst others on the amount of energy stored. Azobenzenes have been investigated as high-potential materials for MOST applications. In the present study it could be shown that intermolecular attractive London dispersion interactions stabilize the (E)-isomer in bisazobenzene that is linked by different alkyl bridges. Differential scanning calorimetry (DSC) measurements revealed, that this interaction leads to an increased storage energy per azo-unit of more than 3 kcal/mol compared to the parent azobenzene. The origin of this effect has been supported by computation as well as X-ray analysis. In the solid state structure attractive London dispersion interactions between the C-H of the alkyl bridge and the π-system of the azobenzene could be clearly assigned. This concept will be highly useful in designing more effective MOST systems in the future.

Prototypical cyclohexane dimers: spectroscopic evidence for σ stacking at low temperatures

In this work, the first unambiguous spectroscopic evidence for the existence of σ stacking interactions in cyclohexane dimers has been provided using matrix isolation infrared spectroscopy.

Ligation state of nickel during C-O bond activation with monodentate phosphines

The oxidative addition of phenolic electrophiles at Ni(0) in the presence of monodentate phosphine ligands was studied with both dispersion-free and dispersion-containing DFT methods. With the popular bulky ligand PCy3, consideration of dispersion has a striking effect on the predicted ligation state of nickel during oxidative addition of aryl sulfamates. Dispersioncontaining methods such as M06L indicate a clear preference for a bis-phosphine ligated transition state (TS), while dispersion free methods like B3LYP strongly favor a monophosphine ligated TS. This discrepancy in predicted ligation state is also found with small phosphines (PMe3) in combination with some aryl electrophiles (carbamates, acetates, pivalates, chlorides), but a bis-PMe3-ligated TS is predicted regardless of dispersion for other electrophiles (sulfamates, mesylates, tosylates). DFT calculations that include dispersion also offer a possible explanation for the observed poor efficacy of P t Bu3 as a ligand in Ni-catalyzed cross-coupling reactions.

Steric "attraction": not by dispersion alone

Non-covalent interactions between neutral, sterically hindered organic molecules generally involve a strong stabilizing contribution from dispersion forces that in many systems turns the 'steric repulsion' into a 'steric attraction'. In addition to London dispersion, such systems benefit from electrostatic stabilization, which arises from a short-range effect of charge penetration and gets bigger with increasing steric bulk. In the present work, we quantify this contribution for a diverse set of molecular cores, ranging from unsubstituted benzene and cyclohexane to their derivatives carrying tert-butyl, phenyl, cyclohexyl and adamantyl substituents. While the importance of electrostatic interactions in the dimers of sp2-rich (e.g., π-conjugated) cores is well appreciated, less polarizable assemblies of sp3-rich systems with multiple short-range CH···HC contacts between the bulky cyclohexyl and adamantyl moieties are also significantly influenced by electrostatics. Charge penetration is drastically larger in absolute terms for the sp2-rich cores, but still has a non-negligible effect on the sp3-rich dimers, investigated herein, both in terms of their energetics and equilibrium interaction distances. These results emphasize the importance of this electrostatic effect, which has so far been less recognized in aliphatic systems compared to London dispersion, and are therefore likely to have implications for the development of force fields and methods for crystal structure prediction.

Ultra-low friction mechanism of highly sp3-hybridized amorphous carbon controlled by interfacial molecule adsorption

The friction behaviors of highly sp³-hybridized carbon films, including ultra-nanocrystalline diamond and diamond-like carbon materials, strongly depend on atmospheres.

The silane-methane dimer revisited: more than a dispersion-bound system?

We present here a comprehensive computational and theoretical analysis of the silane-methane dimer with the goal of understanding the origin of the interactions that hold it together and the factors that affect its strength. Several interaction topologies have been analysed and the associated interaction energies have been evaluated at the CCSD(T)/aug-cc-pVTZ level of theory. Next, substitution effects have been studied on several silane and methane derivatives. The molecular electrostatic potential (MEP) maps of the molecules involved in the interactions have been built to try to correlate the interaction energies with the maximum/minimum EP values (V_s). Furthermore, we have performed an energy decomposition analysis to gain deeper insight into the physical nature of the interactions and to unravel whether dispersion is the primary component of the attraction. Finally, we complete the theoretical analysis with the study of several experimental crystal structures in which there are silylmethyl short contacts.

Implications of Charge Penetration for Heteroatom-Containing Organic Semiconductors

The noncovalent interactions of neutral π-conjugated cores, pertinent to organic semiconductor materials, are intimately related to their charge transport properties and involve a subtle interplay of dispersion, Pauli repulsion, and electrostatic contributions. Realizing structural arrangements that are both energetically preferred and sufficiently conductive is a challenge. We tackle this problem by means of charge penetration contribution to the interaction energy, boosted in systems containing large heteroatoms (e.g., sulfur, selenium, phosphorus, silicon, and arsenic). We find that in both the model and "realistic" dimers of such heteroatom-containing cores dispersion is balanced out by the exchange and interaction energy is instead governed by substantial charge penetration. These systems also feature stronger electronic couplings compared to the dispersion-driven dimers of oligoacenes and/or the herringbone assemblies. Thus, charge penetration, enhanced in the π-conjugated cores comprising larger heteroatoms, arises as an attractive strategy toward potentially more stable and efficient organic electronic materials.

London dispersion in molecular chemistry--reconsidering steric effects

London dispersion, which constitutes the attractive part of the famous van der Waals potential, has long been underappreciated in molecular chemistry as an important element of structural stability, and thus affects chemical reactivity and catalysis. This negligence is due to the common notion that dispersion is weak, which is only true for one pair of interacting atoms. For increasingly larger structures, the overall dispersion contribution grows rapidly and can amount to tens of kcal mol(-1) . This Review collects and emphasizes the importance of inter- and intramolecular dispersion for molecules consisting mostly of first row atoms. The synergy of experiment and theory has now reached a stage where dispersion effects can be examined in fine detail. This forces us to reconsider our perception of steric hindrance and stereoelectronic effects. The quantitation of dispersion energy donors will improve our ability to design sophisticated molecular structures and much better catalysts.

A DFT-D study of hydrogen adsorption on functionalized graphene

In this paper, we use density functional theory with dispersion correction functional (DFT-D) as implemented in the Vienna ab initio simulation package in order to investigate hydrogen adsorption on graphane (GH) and fluorographene (GF).

Graphane versus graphene: a computational investigation of the interaction of nucleobases, aminoacids, heterocycles, small molecules (CO2, H2O, NH3, CH4, H2), metal ions and onium ions

We compared the binding affinity of graphane and graphene with various molecules and ions.

Saturated vs. unsaturated hydrocarbon interactions with carbon nanostructures

The interactions of various acyclic and cyclic hydrocarbons in both saturated and unsaturated forms with the carbon nanostructures (CNSs) have been explored by using density functional theory (DFT) calculations. Model systems representing armchair and zigzag carbon nanotubes (CNTs) and graphene have been considered to investigate the effect of chirality and curvature of the CNSs toward these interactions. Results of this study reveal contrasting binding nature of the acyclic and cyclic hydrocarbons toward CNSs. While the saturated molecules show stronger binding affinity in acyclic hydrocarbons; the unsaturated molecules exhibit higher binding affinity in cyclic hydrocarbons. In addition, acyclic hydrocarbons exhibit stronger binding affinity toward the CNSs when compared to their corresponding cyclic counterparts. The computed results excellently corroborate the experimental observations. The interaction of hydrocarbons with graphene is more favorable when compared with CNTs. Bader's theory of atoms in molecules has been invoked to characterize the noncovalent interactions of saturated and unsaturated hydrocarbons. Our results are expected to provide useful insights toward the development of rational strategies for designing complexes with desired noncovalent interaction involving CNSs.

Understanding the fundamental role of π/π, σ/σ, and σ/π dispersion interactions in shaping carbon-based materials

Noncovalent interactions involving aromatic rings, such as π-stacking and CH/π interactions, are central to many areas of modern chemistry. However, recent studies proved that aromaticity is not required for stacking interactions, since similar interaction energies were computed for several aromatic and aliphatic dimers. Herein, the nature and origin of π/π, σ/σ, and σ/π dispersion interactions has been investigated by using dispersion-corrected density functional theory, energy decomposition analysis, and the recently developed noncovalent interaction (NCI) method. Our analysis shows that π/π and σ/σ stacking interactions are equally important for the benzene and cyclohexane dimers, explaining why both compounds have similar boiling points. Also, similar dispersion forces are found in the benzene⋅⋅⋅methane and cyclohexane⋅⋅⋅methane complexes. However, for systems larger than naphthalene, there are enhanced stacking interactions in the aromatic dimers adopting a parallel-displaced configuration compared to the analogous saturated systems. Although dispersion plays a decisive role in stabilizing all the complexes, the origin of the π/π, σ/σ, and σ/π interactions is different. The NCI method reveals that the dispersion interactions between the hydrogen atoms are responsible for the surprisingly strong aliphatic interactions. Moreover, whereas σ/σ and σ/π interactions are local, the π/π stacking are inherently delocalized, which give rise to a non-additive effect. These new types of dispersion interactions between saturated groups can be exploited in the rational design of novel carbon materials.

Functionalization of homodiamantane: oxygen insertion reactions without rearrangement with dimethyldioxirane

Homodiamantane bromination and nitroxylation are accompanied by contraction of the seven-membered ring to give the corresponding substituted 1-diamantylmethyl derivatives. In contrast, CH-bond hydroxylations with dimethyldioxirane retain the cage and give both apically and medially substituted homodiamantanes. The product ratios are in accord with the barriers for the oxygen insertion computed with density functional theory methods only if solvation is included through a polarizable continuum model. B3LYP-D3 and M06-2X computations with a 6-31G(d,p) basis set on the oligomeric van der Waals complexes predict the potential of homodiamantane derivatives for surface modifications with conformationally slightly flexible diamondoid homologues.

Polytwistane

Twistane, C10H16, is a classic D2-symmetric chiral hydrocarbon that has been studied for decades due to its fascinating stereochemical and thermodynamic properties. Here we propose and analyze in detail the contiguous linear extension of twistane with ethano (ethane-1,2-diyl) bridges to create a new chiral, C2-symmetric hydrocarbon nanotube called polytwistane. Polytwistane, (CH)n, has the same molecular formula as polyacetylene but is composed purely of C(sp(3))-H units, all of which are chemically equivalent. The polytwistane nanotube has the smallest inner diameter (2.6 Å) of hydrocarbons considered to date. A rigorous topological analysis of idealized polytwistane and a C236H242 prototype optimized by B3LYP density functional theory reveals that the polymer has a nonrepeating, alternating σ-helix, with an irrational periodicity parameter and an instantaneous rise (or lead) angle near 15 °. A theoretical analysis utilizing homodesmotic equations and explicit computations as high as CCSD(T)/cc-pVQZ yields the enthalpies of formation Delta(f)H(0)°(twistane) = -1.7 kcal mol(-1) and Delta(f)H(0)°(polytwistane) = +1.28 kcal (mol CH)(-1), demonstrating that the hypothetical formation of polytwistane from acetylene is highly exothermic. Hence, polytwistane is synthetically viable both on thermodynamic grounds and also because no obvious pathways exist for its rearrangement to lower-lying isomers. The present analysis should facilitate the preparation and characterization of this new chiral hydrocarbon nanotube.

Poly[[di-aqua-(μ4-benzene-1,2,4,5-tetra-carboxyl-ato)tetrakis-(1H-imidazole-κN (3))dicopper(II)] N,N-di-methyl-formamide monosolvate]

The asymmetric unit of the polymeric title compound, {[Cu2(C10H2O8)(C3H4N2)4(H2O)2]·C3H7NO} n , contains two independent Cu(II) ions, each coordinated by one water mol-ecule, two imidazole N atoms and two carboxyl-ate O atoms from benzene-1,2,4,5-tetra-carboxyl-ate anions in a distorted square-pyramidal geometry. The benzene-1,2,4,5-tetra-carboxyl-ate anion bridges four Cu(II) ions, forming a polymeric sheet parallel to (010). In the crystal, extensive N-H⋯O and O-H⋯O hydrogen bonds link the polymeric sheets and di-methyl-formamide solvent mol-ecules into a three-dimensional supra-molecular structure.

Understanding the Nature of the CH···HC Interactions in Alkanes

To understand the dispersion stabilization of hydrocarbons in solids and of encumbered molecules, wherein CH···HC interactions act as sticky fingers, we developed here a valence bond (VB) model and applied it to analyze the H···H interactions in dimers of H2 and alkanes. The VB analysis revealed two distinct mechanisms of "dispersion." In the dimers of small molecules like H-H···H-H and H3CH···HCH3, the stabilization arises primarily due to the increased importance of the VB structures which possess charge alternation, e.g., C(+)H(-)···H(+)C(-) and C(-)H(+)···H(-)C(+), and hence bring about electrostatic stabilization that holds the dimer. This is consistent with the classical mechanism of oscillating dipoles as the source of dispersion interactions. However, in larger alkanes, this mechanism is insufficient to glue the two molecules together. Here, the "dispersion" interaction comes about through perturbational mixing of VB structures, which reorganize the bonding electrons of the two interacting CH bonds via recoupling of these electrons to H···H, C···C, and C···H "bonds." Finally, an attempt is made to create a bridge from VB to molecular orbital (MO) and local pair natural-orbital coupled electron pair approximation (LPNO-CEPA/1) analyses of the interactions, which bring about CH···HC binding.

XH/π (X = C, Si) Interactions in Graphene and Silicene: Weak in Strength, Strong in Tuning Band Structures

The lack of a band gap has greatly hindered the applications of graphene in electronic devices. By means of dispersion-corrected density functional theory computations, we demonstrated that considerable CH/π interactions exist between graphene and its fully (graphane) or patterned partially (C4H) hydrogenated derivatives. Due to the equivalence breaking of two sublattices of graphene, a 90 meV band gap is opened in the graphene/C4H bilayer. The band gap can be further increased to 270 meV by sandwiching graphene between two C4H layers. By taking advantage of the similar SiH/π interactions, a 120 meV band gap also can be opened for silicene. Interestingly, the high carrier mobility of graphene/silicene can be well-preserved. Our theoretical results suggest a rather practical solution for gap opening of graphene and silicene, which would allow them to serve as field effect transistors and other nanodevices.