Density-functional approaches to noncovalent interactions: a comparison of dispersion corrections (DFT-D), exchange-hole dipole moment (XDM) theory, and specialized functionals

A systematic study of techniques for treating noncovalent interactions within the computationally efficient density functional theory (DFT) framework is presented through comparison to benchmark-quality evaluations of binding strength compiled for molecular complexes of diverse size and nature. In particular, the efficacy of functionals deliberately crafted to encompass long-range forces, a posteriori DFT+dispersion corrections (DFT-D2 and DFT-D3), and exchange-hole dipole moment (XDM) theory is assessed against a large collection (469 energy points) of reference interaction energies at the CCSD(T) level of theory extrapolated to the estimated complete basis set limit. The established S22 [revised in J. Chem. Phys. 132, 144104 (2010)] and JSCH test sets of minimum-energy structures, as well as collections of dispersion-bound (NBC10) and hydrogen-bonded (HBC6) dissociation curves and a pairwise decomposition of a protein-ligand reaction site (HSG), comprise the chemical systems for this work. From evaluations of accuracy, consistency, and efficiency for PBE-D, BP86-D, B97-D, PBE0-D, B3LYP-D, B970-D, M05-2X, M06-2X, ωB97X-D, B2PLYP-D, XYG3, and B3LYP-XDM methodologies, it is concluded that distinct, often contrasting, groups of these elicit the best performance within the accessible double-ζ or robust triple-ζ basis set regimes and among hydrogen-bonded or dispersion-dominated complexes. For overall results, M05-2X, B97-D3, and B970-D2 yield superior values in conjunction with aug-cc-pVDZ, for a mean absolute deviation of 0.41 - 0.49 kcal/mol, and B3LYP-D3, B97-D3, ωB97X-D, and B2PLYP-D3 dominate with aug-cc-pVTZ, affording, together with XYG3/6-311+G(3df,2p), a mean absolute deviation of 0.33 - 0.38 kcal/mol.     

Acid-base thermochemistry of gaseous oxygen and sulfur substituted amino acids (Ser, Thr, Cys, Met)

Acid-base thermochemistry of isolated amino acids containing oxygen or sulfur in their side chain (serine, threonine, cysteine and methionine) have been examined by quantum chemical computations. Density functional theory (DFT) was used, with B3LYP, B97-D and M06-2X functionals using the 6-31+G(d,p) basis set for geometry optimizations and the larger 6-311++G(3df,2p) basis set for energy computations. Composite methods CBS-QB3, G3B3, G4MP2 and G4 were applied to large sets of neutral, protonated and deprotonated conformers. Conformational analysis of these species, based on chemical approach and AMOEBA force field calculations, has been used to identify the lowest energy conformers and to estimate the population of conformers expected to be present at thermal equilibrium at 298 K. It is observed that G4, G4MP2, G3B3, CBS-QB3 composite methods and M06-2X DFT lead to similar conformer energies. Thermochemical parameters have been computed using either the most stable conformers or equilibrium populations of conformers. Comparison of experimental and theoretical proton affinities and Δ(acid)H shows that the G4 method provides the better agreement with deviations of less than 1.5 kJ mol(-1). From this point of view, a set of evaluated thermochemical quantities for serine, threonine, cysteine and methionine may be proposed: PA = 912, 919, 903, 938; GB = 878, 886, 870, 899; Δ(acid)H = 1393, 1391, 1396, 1411; Δ(acid)G = 1363, 1362, 1367, 1382 kJ mol(-1). This study also confirms that a non-negligible ΔpS° is associated with protonation of methionine and that the most acidic hydrogen of cysteine in the gas phase is that of the SH group. In several instances new conformers were identified thus suggesting a re-examination of several IRMPD spectra.     

Assessment of density functionals with long‐range and/or empirical dispersion corrections for conformational energy calculations of peptides

Density functionals with long‐range and/or empirical dispersion corrections, including LC‐ωPBE, B97‐D, ωB97X‐D, M06‐2X, B2PLYP‐D, and mPW2PLYP‐D functionals, are assessed for their ability to describe the conformational preferences of Ac‐Ala‐NHMe (the alanine dipeptide) and Ac‐Pro‐NHMe (the proline dipeptide) in the gas phase and in water, which have been used as prototypes for amino acid residues of peptides. For both dipeptides, the mean absolute deviation (MAD) is estimated to be 0.22–0.40 kcal/mol in conformational energy and 2.0–3.2° in torsion angles ? and ψ using these functionals with the 6‐311++G(d,p) basis set against the reference values calculated at the MP2/aug‐cc‐pVTZ//MP2/aug‐cc‐pVDZ level of theory in the gas phase. The overall performance is obtained in the order B2PLYP‐D ≈ mPW2PLYP‐D > ωB97X‐D ≈ M06‐2X > MP2 > LC‐ωPBE > B3LYP with the 6–311++G(d,p) basis set. The SMD model at the M06‐2X/6‐31+G(d) level of theory well reproduced experimental hydration free energies of the model compounds for backbone and side chains of peptides with MADs of 0.47 and 4.3 kcal/mol for 20 neutral and 5 charged molecules, respectively. The B2PLYP‐D/6‐311++G(d,p)//SMD M06‐2X/6‐31+G(d) level of theory provides the populations of backbone and/or prolyl peptide bond for the alanine and proline dipeptides in water that are consistent with the observed values. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Electronic structure and absorption spectra of supramolecular complexes of a fullerene crown ether with a π-extended TTF derivative

The present work is a theoretical investigation on supramolecular complexes of a fullerene crown ether (A and B isomers) with a derivative of π-extended tetrathiafulvalene (T). The geometry and the electronic structure of seven different conformers of the complex of dibenzo-18-crown-6 ether of fullero-N-methylpyrrolidine with a N-benzyl-N-(4-{[9,10-bis(1,3-dithiol-2-ylidene)-9,10-dihydroanthracen-2-yl]ethynyl}benzyl)ammonium cation were determined. We calculated the complexation energies and the absorption spectra, i.e., the lowest 50 excited electronic states of the complexes have been determined at the ground state optimum geometry. All calculations were carried out employing the density functional theory (DFT) and the time-dependent DFT, using the B3LYP, CAM-B3LYP, ωB97X-D, and M06-2X functionals in conjunction with the 6-31G(d,p) basis set. Various types of van der Waals interactions are observed in the complexes. Conformer complexation energies (CE) range from 2.54 to 2.14 eV in the gas phase and from 1.75 to 1.34 eV in CHCl(3) solvent at the ωB97X-D/6-31G(d,p)//M06-2X/6-31G(d,p) level of theory. There are three major features at about 390, 330, and 290 nm in the calculated absorption spectra of all the conformers. The major peaks correspond to T→T, T→T/F (electron density in both T and the fullerene F of B) and to T→F transitions, depending on the particular conformer. Other charge transfer T→F transitions are observed close to the T→T transition, indicating the possibility of photoinduced electron transfer in all these complexes.     

A comparison of the behavior of functional/basis set combinations for hydrogen-bonding in the water dimer with emphasis on basis set superposition error

We evaluate the performance of ten functionals (B3LYP, M05, M05-2X, M06, M06-2X, B2PLYP, B2PLYPD, X3LYP, B97D, and MPWB1K) in combination with 16 basis sets ranging in complexity from 6-31G(d) to aug-cc-pV5Z for the calculation of the H-bonded water dimer with the goal of defining which combinations of functionals and basis sets provide a combination of economy and accuracy for H-bonded systems. We have compared the results to the best non-density functional theory (non-DFT) molecular orbital (MO) calculations and to experimental results. Several of the smaller basis sets lead to qualitatively incorrect geometries when optimized on a normal potential energy surface (PES). This problem disappears when the optimization is performed on a counterpoise (CP) corrected PES. The calculated interaction energies (ΔEs) with the largest basis sets vary from -4.42 (B97D) to -5.19 (B2PLYPD) kcal/mol for the different functionals. Small basis sets generally predict stronger interactions than the large ones. We found that, because of error compensation, the smaller basis sets gave the best results (in comparison to experimental and high-level non-DFT MO calculations) when combined with a functional that predicts a weak interaction with the largest basis set. As many applications are complex systems and require economical calculations, we suggest the following functional/basis set combinations in order of increasing complexity and cost: (1) D95(d,p) with B3LYP, B97D, M06, or MPWB1k; (2) 6-311G(d,p) with B3LYP; (3) D95++(d,p) with B3LYP, B97D, or MPWB1K; (4) 6-311++G(d,p) with B3LYP or B97D; and (5) aug-cc-pVDZ with M05-2X, M06-2X, or X3LYP.     

CF3CH2CF2CH3（HFC-365mfc）与Cl原子反应的微观机理及动力学性质

采用密度泛函理论方法 M06-2X结合6-31+G(d,p)基组研究了CF3CH2CF2CH3与Cl原子反应的反应机理.计算获得了CF3CH2CF2CH3的两种可区分的稳定几何构象RC1和RC2以及与它们相对应的8条氢提取反应通道和2条取代反应通道.运用改进的正则变分过渡态理论(ICVT)并结合小曲率隧道效应校正(SCT),在M06-2X/6-31+G(d,p)水平上计算了各氢提取通道的速率常数,并由Boltzmann配分函数得到总包反应的速率常数kT(cm3.molecule-1.s-1).计算结果表明,体系的总反应速率常数与已有实验值相吻合,进而给出了该反应在200～1000 K温度区间内反应速率常数kT的三参数表达式kT=1.88×10-22T3.76.exp(-1780.69/T),并讨论了两种构象RC1和RC2对总反应的贡献及各构象中氢提取发生在—CH3或—CH2—基团上的位置选择性.此外,由于缺少相关反应物及产物自由基标准生成焓ΔHf,298 K的数据,利用等化学键法估算了在上述物种的标准生成焓.     

用GGA密度泛函及其长程、色散校正方法计算各类氢键的结合能

应用广义梯度近似(GGA) (PW91和PBE)、含动能密度的广义梯度近似(meta-GGA) (M06-L)、杂化泛函(hyper-GGA)(M06-2X、X3LYP和B3LYP)及其长程校正泛函LC-DFT(CAM-B3LYP、LC-ωPBE和ωB97X)和色散校正密度泛函(DFT-D)(ωB97X-D和B97-D),用多种基函数对15种不同强度的传统氢键和非传统氢键体系的结合能进行了系统的计算与分析.并与高精度的CCSD(T)/aug-cc-pVQZ结果比较发现:在上述各类泛函中,对于氢键结合能的计算M06-2X和ωB97X-D泛函较为精确与可靠,且没有必要使用过大的基函数,6-311++G(2d,2p)或aug-cc-pVDZ水平的基组就已足够,各类泛函所计算结合能的基组重叠误差(BSSE)均较小,除ωB97X和ωB97X-D外,其它9种泛函不经BSSE校正也能得到同样甚至更准确的结果.     

Potential reaction initiation points of polycyclic aromatic hydrocarbons

Reaction initiation points of the 16 priority polycyclic aromatic hydrocarbons (PAHs) have been determined by calculating all the different C-H bond dissociation enthalpy (BDE) values. Six density functional theory methods (B3LYP, B3LYP-D3, B97D3, M06-LD3, M06-2X-D3, and ωB97X-D) in combination with 4 basis sets (6-31G(d), 6-31+G(d,p), 6-311++G(d,p), def2-TZVP) have been applied and the most feasible combination has been selected. The BDE values and the corresponding bond lengths have been used to determine potential attack points on the structures. The studied molecules have been categorized structurally as well, within which the position of the hydrogen atoms is considered. Results show that most of the hydrogens are in zig-zag positions, and the BDE and bond length values for the 16 priority PAHs are in a range between 342.0 and 485.6 kJ/mol and 1.0817–1.952 Å, respectively. Most of the initiation points are represented by armchair and peak hydrogens. The lowest and highest BDE and shortest and longest bond length values belong to fluorene and acenaphthylene where the hydrogens were aliphatic and in peak position, respectively.     

Conformations of the glycine tripeptide analog Ac-Gly-Gly-NHMe: a computational study including aqueous solvation effects

A computational study of the conformational preferences of the glycine tripeptide analog, Ac-Gly-Gly-NHMe, has been carried out. The molecule is considered in isolation as well as with a continuum model of aqueous solvation. In the absence of solvent, several low-energy conformers are found that exhibit turnlike structures including type I and type II β turns. Upon consideration of aqueous solvation, two conformers, corresponding to the type I and II turn structures are found to be significantly lower in energy than all others. Results from ab initio molecular orbital theory calculations at MP2/aug-cc-pVTZ//MP2/6-311+G(d,p) are compared with those from density functional theory with B3LYP, ωB97X-D, B97-D, and M06-2X as well as several empirical force fields.     

The relative energies of polypeptide conformers predicted by linear scaling second-order Møller-Plesset perturbation theory

We describe an implementation of the cluster-in-molecule(CIM) resolution of the identity(RI) approximation second-order M?ller–Plesset perturbation theory(CIM-RI-MP2), with the purpose of extending RI-MP2 calculations to very large systems. For typical conformers of several large polypeptides, we calculated their conformational energy differences with the CIM-RI-MP2 and the generalized energy-based fragmentation MP2(GEBF-MP2) methods, and compared these results with the density functional theory(DFT) results obtained with several popular functionals. Our calculations show that the conformational energy differences obtained with CIM-RI-MP2 and GEBF-MP2 are very close to each other. In comparison with the GEBF-MP2 and CIM-RI-MP2 relative energies, we found that the DFT functionals(CAM-B3LYP-D3, LC-?PBE-D3, M05-2X, M06-2X and ?B97XD) can give quite accurate conformational energy differences for structurally similar conformers, but provide less-accurate results for structurally very different conformers.     

Computational methods to calculate accurate activation and reaction energies of 1,3-dipolar cycloadditions of 24 1,3-dipoles

Theoretical calculations were performed on the 1,3-dipolar cycloaddition reactions of 24 1,3-dipoles with ethylene and acetylene. The 24 1,3-dipoles are of the formula X≡Y(+)-Z(-) (where X is HC or N, Y is N, and Z is CH(2), NH, or O) or X═Y(+)-Z(-) (where X and Z are CH(2), NH, or O and Y is NH, O, or S). The high-accuracy G3B3 method was employed as the reference. CBS-QB3, CCSD(T)//B3LYP, SCS-MP2//B3LYP, B3LYP, M06-2X, and B97-D methods were benchmarked to assess their accuracies and to determine an accurate method that is practical for large systems. Several basis sets were also evaluated. Compared to the G3B3 method, CBS-QB3 and CCSD(T)/maug-cc-pV(T+d)Z//B3LYP methods give similar results for both activation and reaction enthalpies (mean average deviation, MAD, < 1.5 kcal/mol). SCS-MP2//B3LYP and M06-2X give small errors for the activation enthalpies (MAD < 1.5 kcal/mol), while B3LYP has MAD = 2.3 kcal/mol. SCS-MP2//B3LYP and B3LYP give the reasonable reaction enthalpies (MAD < 5.0 kcal/mol). The B3LYP functional also gives good results for most 1,3-dipoles (MAD = 1.9 kcal/mol for 17 common 1,3-dipoles), but the activation and reaction enthalpies for ozone and sulfur dioxide are difficult to calculate by any of the density functional methods.     

Quantum-chemical DFT calculations of the energies of the O—H bond of phenols

The energies of homolytic cleavage of the O—H bonds in 36 phenols with various structures were calculated within the scope of density functional theory using a series of exchange-correlation functionals and basis sets. The best agreement between the calculated and experimental data is given by the M06-2X/6-311++G(2d,2p)//B3LYP/6-31G(d) method.     

Substituent effects on singlet-triplet gaps and mechanisms of 1,2-rearrangements of 1,3-oxazol-2-ylidenes to 1,3-oxazoles

Electronic structures, partial atomic charges, singlet-triplet gaps (Delta E ST), substituent effects, and mechanisms of 1,2-rearrangements of 1,3-oxazol-2-ylidene ( 5) and 4,5-dimethyl- ( 6), 4,5-difluoro- ( 7), 4,5-dichloro- ( 8), 4,5-dibromo- ( 9), and 3-methyl-1,3-oxazol-2-ylidene ( 10) to the corresponding 1,3-oxazoles have been studied using complete-basis-set methods (CBS-QB3, CBS-Q, CBS-4M), second-order M?ller-Plesset perturbation method (MP2), hybrid density functionals (B3LYP, B3PW91), coupled-cluster theory with single and double excitations (CCSD) and CCSD plus perturbative triple excitations [CCSD(T)], and the quadratic configuration interaction method including single and double excitations (QCISD) and QCISD plus perturbative triple excitations [QCISD(T)]. The 6-311G(d,p), 6-31+G(d,p), 6-311+G(d,p), and correlation-consistent polarized valence double-xi (cc-pVDZ) basis sets were employed. The carbenes have singlet ground states, and the CBS-QB3 and CBS-Q methods predict Delta E ST values for 5- 8 and 10 of 79.9, 79.8, 74.7, 77.0, and 82.0 kcal/mol, respectively. CCSD(T), QCISD(T), B3LYP, and B3PW91 predict smaller Delta E ST values than CBS-QB3 and CBS-Q, with the hybrid density functionals predicting the smallest values. The concerted unimolecular exothermic out-of-plane 1,2-rearrangements of singlet 1,3-oxazol-2-ylidenes to their respective 1,3-oxazoles proceed via cyclic three-center transition states. The CBS-predicted barriers to the 1,2-rearrangements of singlet carbenes 5- 9 to their respective 1,3-oxazoles are 41.4, 40.4, 37.8, 40.4, and 40.5 kcal/mol, respectively. During the 1,2-rearrangements of singlet 1,3-oxazol-2-ylidenes 5- 9, there is a decrease in electron density at oxygen, N3 (the migration origin), and C5 and an increase in electron density at C2 (the migration terminus), C4, and the partially positive migrating hydrogen.     

Comparison of aromatic NH···π, OH···π, and CH···π interactions of alanine using MP2, CCSD, and DFT methods

A comparison of the performance of various density functional methods including long‐range corrected and dispersion corrected methods [MPW1PW91, B3LYP, B3PW91, B97‐D, B1B95, MPWB1K, M06‐2X, SVWN5, ωB97XD, long‐range correction (LC)‐ωPBE, and CAM‐B3LYP using 6‐31+G(d,p) basis set] in the study of CH···π, OH···π, and NH···π interactions were done using weak complexes of neutral (A) and cationic (A⁺) forms of alanine with benzene by taking the Møller–Plesset (MP2)/6‐31+G(d,p) results as the reference. Further, the binding energies of the neutral alanine–benzene complexes were assessed at coupled cluster (CCSD)/6‐31G(d,p) method. Analysis of the molecular geometries and interaction energies at density functional theory (DFT), MP2, CCSD methods and CCSD(T) single point level reveal that MP2 is the best overall performer for noncovalent interactions giving accuracy close to CCSD method. MPWB1K fared better in interaction energy calculations than other DFT methods. In the case of M06‐2X, SVWN5, and the dispersion corrected B97‐D, the interaction energies are significantly overrated for neutral systems compared to other methods. However, for cationic systems, B97‐D yields structures and interaction energies similar to MP2 and MPWB1K methods. Among the long‐range corrected methods, LC‐ωPBE and CAM‐B3LYP methods show close agreement with MP2 values while ωB97XD energies are notably higher than MP2 values. © 2010 Wiley Periodicals, Inc. J Comput Chem 2010     

A DFT study of substituent effects in corannulene dimers

Corannulene dimers made up of corannulene monomers with different curvature and substituents were studied using M06-2X, B97D and ωB97XD functionals and 6-31+G* basis set. Corannulene molecules were substituted with five alternating Br, Cl, CH(3), C(2)H or CN units. Geometric results showed that substituents gave rise to small changes in the curvature of corannulene bowls. So, there was not a clear relationship between the curvature of bowls and the changes on interaction energy generated by addition of substituents in the bowl. Electron withdrawing substituents gave rise to a more positive molecular electrostatic potential (MEP) of the bowl, which was able to get a strong interaction with the negative MEP at the surface of a fullerene. Substitution with CN caused the largest effect, giving rise to the most positive MEP and to a large interaction energy of -24.64 kcal mol(-1), at the ωB97XD/6-31+G* level. Dispersive effects must be taken into account to explain the catching ability of the different substituted corannulenes. For unsubstituted dimers, calculations with DFT-D methods employing ωB97XD and B97D functionals led to similar results to those previously reported at the SCS-MP2/cc-pVTZ level for corannulene dimers (A. Sygula and S. Saeb?, Int. J. Quant. Chem., 2009, 109, 65). In particular, the ωB97XD functional led to a difference of only 0.35 kcal mol(-1), regarding MP2 interaction energy for corannulene dimers. On the other hand, the M06-2X functional showed a general considerable underestimation of interaction energies. This functional worked quite well to study trends, but not to obtain absolute interaction energies.     

Assessment of ab initio MP2 and density functionals for characterizing the potential energy profiles of the SN2 reactions at N center

The potential energy profiles of five selected bimolecular nucleophilic substitution (S_N2) reactions at nitrogen (N) center have been reinvestigated with the CCSD(T), G3[MP2,CCSD(T)], MP2, and some density functional methods. The basis sets of 6‐31+G(d,p) and 6‐311+G(3d,2p) are used for the MP2 and density functional calculations. Taking the relative energies at the CCSD(T)/CBS level of theory as benchmarks, we recommend the MP2, B97‐K, B2K‐PLYP, BMK, ωB97X‐D, M06‐2X, M05‐2X, CAM‐B3LYP, M08‐SO, and ωB97X methods to generally characterize the potential energy profiles for the S_N2 reactions at N center. Furthermore, these recommended methods with the relatively small 6‐31+G(d,p) basis set may also be used to perform direct classical trajectory simulations to uncover the dynamic behaviors of the S_N2 reactions at N center. © 2012 Wiley Periodicals, Inc.     

Benchmark of density functional theory methods on the prediction of bond energies and bond distances of noble-gas containing molecules

We have tested three pure density functional theory (DFT) functionals, BLYP, MPWPW91, MPWB95, and ten hybrid DFT functionals, B3LYP, B3P86, B98, MPW1B95, MPW1PW91, BMK, M05-2X, M06-2X, B2GP-PLYP, and DSD-BLYP with a series of commonly used basis sets on the performance of predicting the bond energies and bond distances of 31 small neutral noble-gas containing molecules. The reference structures were obtained using the CCSD(T)∕aug-cc-pVTZ theory and the reference energies were based on the calculation at the CCSD(T)∕CBS level. While in general the hybrid functionals performed significantly better than the pure functionals, our tests showed a range of performance by these hybrid functionals. For the bond energies, the MPW1B95∕6-311+G(2df,2pd), BMK∕aug-cc-pVTZ, B2GP-PLYP∕aug-cc-pVTZ, and DSD-BLYP∕aug-cc-pVTZ methods stood out with mean unsigned errors of 2.0-2.3 kcal∕mol per molecule. For the bond distances, the MPW1B95∕6-311+G(2df,2pd), MPW1PW91∕6-311+G(2df,2pd), and B3P86∕6-311+G(2df,2pd), DSD-BLYP∕6-311+G(2df,2pd), and DSD-BLYP∕aug-cc-pVTZ methods stood out with mean unsigned errors of 0.008-0.013 A? per bond. The current study showed that a careful selection of DFT functionals is very important in the study of noble-gas chemistry, and the most recommended methods are MPW1B95∕6-311+G(2df,2pd) and DSD-BLYP∕aug-cc-pVTZ.     

σ/σ- And π/π-interactions are equally important: multilayered graphanes

The properties of single-sheet [n]graphanes, their double-layered forms (diamondoids), and their van der Waals (vdW) complexes (multilayered [n]graphanes) were studied for n = 10-97 at the dispersion-corrected density functional theory (DFT) level utilizing B97D with a 6-31G(d,p) basis set; for comparison, we also computed a series of structures at M06-2X/6-31G(d,p) as well as B3LYP-D3/6-31G(d,p) and evaluated SCS-MP2/cc-pVDZ single-point energies. The association energies for the vdW complexes reach 120 kcal mol(-1) already at 2 nm particle size ([97]graphane dimer), and graphanes adopt layered structures similar to that of graphenes. The association energies of multilayered graphanes per carbon atom are rather similar and independent of the number of layers (ca. 1.2 kcal mol(-1)). Graphanes show quantum confinement effects as the HOMO-LUMO gaps decrease from 8.2 eV for [10]graphane to 5.7 eV for [97]graphane, asymptotically approaching 5.4 eV previously obtained for bulk graphane. Similar trends were found for layered graphanes, where the differences in the electronic properties of double-sheet CH/σ vdW and double-layered CC/σ diamondoids vanish at particles sizes of 1 nm. For comparison, we studied the parent CC/π systems, i.e., the single- and double-sheet [n]graphenes (n = 10-130) for which the association energies demonstrate the same trends as in the case of [n]graphanes; in both cases the band gaps decrease with an increase in system size. The [112]graphene dimer (HOMO-LUMO gap = 0.5 eV) already approaches the 2D metallic properties of graphite.