A density-functional study on pi-aromatic interaction: benzene dimer and naphthalene dimer

The long-range correction (LC) scheme of density-functional theory (DFT) was applied to the calculation of the pi-aromatic interaction of the benzene dimer and naphthalene dimer. In previous calculations, it was confirmed that the LC scheme [Iikura et al., J. Chem. Phys. 115, 3540 (2001)] gives very accurate potential- energy surfaces (PESs) of small van der Waals (vdW) complexes by combining with the Anderson-Langreth-Lundqvist (ALL) vdW correlation functional [Andersson et al., Phys. Rev. Lett. 76, 102 (1996)] (LC-DFT + ALL). In this study, LC-DFT+ALL method was examined by calculating a wide range of PES of the benzene dimer including parallel, T-shaped, and parallel-displaced configurations. As a result, we succeeded in reproducing very accurate PES within the energy deviance of less than 1 kcalmol in comparison with the results of high-level ab initio molecular-orbital methods at all reference points on the PES. It was also found that LC-DFT + ALL gave accurate results independent of exchange-correlation functional used, in contrast with the strong functional dependencies of conventional pure functionals. This indicates that both exchange repulsion and van der Waals attractive interactions should be correctly incorporated in conventional pure functionals in order to calculate accurate pi-aromatic interactions. We also found that LC-DFT + ALL method has a low basis-set dependency in the calculations of pi-aromatic interactions. The present scheme was also successfully applied to the pi,[ellipsis (horizontal)],pi stacking interactions of naphthalene dimer. This may suggest that LC-DFT + ALL method would be a powerful tool in the calculations of large molecules such as biomolecules.     

Modeling chemical reactions on surfaces: The roles of chemical bonding and van der Waals interactions

Chemical reactions on surfaces play central roles in heterogeneous catalysis, and most reactions involve the formation and/or the cleavage of bonds. At present, density functional theory (DFT) has become the workhorse for computational investigation of reaction mechanisms, but its predictive power has been severely limited by the lack of appropriate exchange-correlation functionals. Here, we show that there are many cases where the chemical bonding and van der Waals (vdW) interactions both play a key role in chemical reactions on surfaces. After briefly introducing some DFT methods and basic theory in chemical reactions, we first demonstrate that DFT can help to understand the mechanisms of “classic” reactions that mainly dominated by covalent bonding and vdW forces, as exemplified in electrocatalytic reduction of CO₂ and the fabrication of 2D materials on metal substrates. We next show that DFT calculations can help to uncover the tautomerization reactions of molecules on metal surfaces, wherein the hydrogen bonding and vdW forces would largely affect the reaction process. More importantly, we show that in some cases, the vdW interactions can become the decisive effect that determines the adsorption configuration, energy hierarchy, and the potential-energy surface of chemical reactions, yielding distinct pathways and products. Additionally, we highlight the importance of more realistic conditions, such as surface defects, finite coverage, and temperature effects, in accurate modeling of chemical reactions. Finally, we summarize some challenges in modeling catalysis, which include many-body dispersive correction, strong correlation effect, and non-adiabatic approximations.     

An application of the van der Waals density functional: Hydrogen bonding and stacking interactions between nucleobases

We apply the van der Waals density functional (vdW-DF) to study hydrogen bonding and stacking interactions between nucleobases. The excellent agreement of our results with high level quantum chemical calculations highlights the value of the vdW-DF for first-principles investigations of biologically important molecules. Our results suggest that, in the case of hydrogen-bonded nucleobase pairs, dispersion interactions reduce the cost of propeller twists while having a negligible effect on buckling. Furthermore, the efficient scaling of DFT methods allowed for the easy optimization of separation distance between nucleobase stacks, indicating enhancements in the interaction energy of up to 3 kcalmol over previous fixed distance calculations. We anticipate that these results are significant for extending the vdW-DF method to model larger vdW complexes and biological molecules.     

Free energy decomposition analysis of bonding and nonbonding interactions in solution

A free energy decomposition analysis algorithm for bonding and nonbonding interactions in various solvated environments, named energy decomposition analysis-polarizable continuum model (EDA-PCM), is implemented based on the localized molecular orbital-energy decomposition analysis (LMO-EDA) method, which is recently developed for interaction analysis in gas phase [P. F. Su and H. Li, J. Chem. Phys. 130, 074109 (2009)]. For single determinant wave functions, the EDA-PCM method divides the interaction energy into electrostatic, exchange, repulsion, polarization, desolvation, and dispersion terms. In the EDA-PCM scheme, the homogeneous solvated environment can be treated by the integral equation formulation of PCM (IEFPCM) or conductor-like polarizable continuum model (CPCM) method, while the heterogeneous solvated environment is handled by the Het-CPCM method. The EDA-PCM is able to obtain physically meaningful interaction analysis in different dielectric environments along the whole potential energy surfaces. Test calculations by MP2 and DFT functionals with homogeneous and heterogeneous solvation, involving hydrogen bonding, vdW interaction, metal-ligand binding, cation-π, and ionic interaction, show the robustness and adaptability of the EDA-PCM method. The computational results stress the importance of solvation effects to the intermolecular interactions in solvated environments.     

On the Role of Charge Transfer in Many-Body Non-Covalent Interactions**

Charge transfer is one of the mechanisms involved in non-covalent interactions. In molecular dimers, its contribution to pairwise interaction energies has been studied extensively using a variety of interaction energy decomposition schemes. In polar interactions such as hydrogen bonds, it can contribute ten or several tens of percent of the interaction energy. Less is known about its importance in higher-order interactions in many-body systems, mainly because of the lack of methods applicable to this problem. In this work, we extend our method for the quantification of the charge-transfer energy based on constrained DFT to many-body cases and apply it to model trimers extracted from molecular crystals. Our calculations show that charge transfer can account for a large fraction of the total three-body interaction energy. This also has implications for DFT calculations of many-body interactions in general as it is known that many DFT functionals struggle to describe charge-transfer effects correctly.     

A density functional study on dielectric properties of acrylic acid grafted polypropylene

Influence of acrylic acid grafting of isotactic polypropylene on the dielectric properties of the polymer is investigated using density functional theory (DFT) calculations, both in the molecular modeling and three-dimensional (3D) bulk periodic system frameworks. In our molecular modeling calculations, polarizability volume, and polarizability volume per mass which reflects the permittivity of the polymer, as well as the HOMO-LUMO gap, one of the important measures indicating the electrical breakdown voltage strength, were examined for oligomers with various chain lengths and carboxyl mixture ratios. When a polypropylene oligomer is grafted with carboxyl groups (cf. acrylic acid), our calculations show that the increase of the polarizability volume α' of the oligomer is proportional to the increase of its mass m, while the ratio α'/m decreases from the value of a pure polymer when increasing the mixture ratio. The decreasing ratio of α'/m under carboxyl grafting indicates that the material permittivity might also decrease if the mass density of the material remains constant. Furthermore, our calculations show that the HOMO-LUMO gap energy decreases by only about 15% in grafting, but this decrease seems to be independent on the mixture ratio of carboxyl. This indicates that by doping polymers with additives better dielectric properties can be tailored. Finally, using the first-principles molecular DFT results for polarizability volume per mass in connection with the classical Clausius-Mossotti relation, we have estimated static permittivity for acrylic acid grafted polypropylene, assuming the structural density keeping constant under grafting. The computed permittivity values are in a qualitative agreement with the recent experiments, showing increasing tendency of the permittivity as a function of the grafting composition. In order to validate our molecular DFT based approach, we have also carried out extensive three-dimensional bulk periodic first-principles total-energy calculations in the frameworks of the density functional theory and density functional perturbation theory (DFPT) for crystalline acrylic acid grafted polypropylene. Interestingly, the computed electronic and dielectric properties behave very similarly between the simplified molecular DFT modeling and the more realistic 3D bulk periodic DFPT method. In particular, the static permittivity values [ε(r)(0)] from the molecular DFT-Clausius-Mossotti modeling are in excellent agreement with the high-frequency dielectric constant values (ε(∞)) from the DFPT method. This obviously implies that the chain-to-chain interaction to dielectric and electronic properties of acrylic acid polypropylene, to a first approximation, can be neglected, therefore justifying the usage of molecular DFT modeling in our calculations.     

Configurational enthalpy of polystyrene

Recent precision measurements of the heats of combustion of atactic and isotactic polystyrene permit an unequivocal calculation of the enthalpy difference in the bulk amorphous forms of the two isomers of this polymer. Contributions to this quantity arise mainly from the differences in nonbonded interactions in the two configurations but may also contain terms relating to higher energy conformations and to intermolecular interactions. The thermochemical and NMR data of specific polymers and of simple molecules are discussed in a comparison with model compound calculations. The thermochemical method has potentially important applications in studying molecular interactions in stereoregular polymers.     

Assessment of density functional theory to calculate the phase transition pressure of ice

To assess the accuracy of density functional theory (DFT) methods in describing hydrogen bonding in condensed phases, we benchmarked their performance in describing phase transitions among different phases of ice. We performed DFT calculations of ice for phases Ih, II, III, VI and VII using BLYP, PW91, PBE, PBE-D, PBEsol, B3LYP, PBE0, and PBE0-D, and compared the calculated phase transition pressures between Ih-II, Ih-III, II-VI, and VI-VII with the 0 K experimental values of Whalley [J. Chem. Phys., 1984, 81, 4087]. From the geometry optimization of many different candidates, we found that the most stable proton orientation as well as the phase transition pressure does not show much functional dependence for the generalized gradient approximation and hybrid functionals. Although all these methods overestimated the phase transition pressure, the addition of van der Waals (vdW) correction using PBE-D and PBE0-D reduced the transition pressure and improved the agreement for Ih-II. On the other hand, energy ordering between VI and VII reversed and gave an unphysical negative transition pressure. Binding energy profiles of a few conformations of water dimers were calculated to understand the improvement for certain transitions and failures for others with the vdW correction. We conclude that vdW dispersion forces must be considered to accurately describe the hydrogen bond in many different phases of ice, but the simple addition of the R(-6) term with a small basis set tends to over stabilize certain geometries giving unphysical ordering in the high density phases.     

Water monomer interaction with gold nanoclusters from van der Waals density functional theory

We investigate the interaction between water molecules and gold nanoclusters Au(n) through a systematic density functional theory study within both the generalized gradient approximation and the nonlocal van der Waals (vdW) density functional theory. Both planar (n = 6-12) and three-dimensional (3D) clusters (n = 17-20) are studied. We find that applying vdW density functional theory leads to an increase in the Au-Au bond length and a decrease in the cohesive energy for all clusters studied. We classify water adsorption on nanoclusters according to the corner, edge, and surface adsorption geometries. In both corner and edge adsorptions, water molecule approaches the cluster through the O atom. For planar clusters, surface adsorption occurs in a O-up/H-down geometry with water plane oriented nearly perpendicular to the cluster. For 3D clusters, water instead favors a near-flat surface adsorption geometry with the water O atom sitting nearly atop a surface Au atom, in agreement with previous study on bulk surfaces. Including vdW interaction increases the adsorption energy for the weak surface adsorption but reduces the adsorption energy for the strong corner adsorption due to increased water-cluster bond length. By analyzing the adsorption induced charge rearrangement through Bader's charge partitioning and electron density difference and the orbital interaction through the projected density of states, we conclude that the bonding between water and gold nanocluster is determined by an interplay between electrostatic interaction and covalent interaction involving both the water lone-pair and in-plane orbitals and the gold 5d and 6s orbitals. Including vdW interaction does not change qualitatively the physical picture but does change quantitatively the adsorption structure due to the fluxionality of gold nanoclusters.     

Ordering of anisotropic polarizable polymer chains on the full many-body level

We study the effect of dielectric anisotropy of polymers on their equilibrium ordering within mean-field theory, but with a formalism that takes into account the full n-body nature of van der Waals (vdW) forces. Dielectric anisotropy within polymers is to be expected as the electronic properties of the polymer will typically be different along the polymer than across its cross section. It is therefore physically intuitive that larger charge fluctuations can be induced along the chain than perpendicular to it. We show that this dielectric anisotropy leads to n-body interactions which can induce an isotropic-nematic transition. The two body and three body components of the full vdW interaction are extracted and it is shown how the two body term behaves like the phenomenological self-aligning-pairwise nematic interaction. At the three body interaction level we see that the nematic phase that is energetically favorable is discotic, however, on the full n-body interaction level we find that the normal axial nematic phase is always the stable ordered phase. The n-body nature of our approach also shows that the key parameter driving the nematic-isotropic transition is the bare persistence length of the polymer chain.     

Predicting adsorption enthalpies on silicalite and HZSM‐5: A benchmark study on DFT strategies addressing dispersion interactions

We evaluated the accuracy of periodic density functional calculations for adsorption enthalpies of water, alkanes, and alcohols in silicalite and HZSM‐5 zeolites using a gradient‐corrected density functional with empirical dispersion corrections (PBE‐D) as well as a nonlocal correlation functional (vdW‐DF2). Results of both approaches agree in acceptable fashion with experimental adsorption energies of alcohols in silicalite, but the adsorption energies for n‐alkanes in both zeolite models are overestimated, by 21?46 kJ mol^?1. For PBE‐D calculations, the adsorption of alkanes is exclusively determined by the empirical dispersion term, while the generalized gradient approximation‐DFT part is purely repulsive, preventing the molecule to come too close to the zeolite walls. The vdW‐DF2 results are comparable to those of PBE‐D calculations, but the latter values are slightly closer to the experiment in most cases. Thus, both computational approaches are unable to reproduce available experimental adsorption energies of alkanes in silicalite and HZSM‐5 zeolite with chemical accuracy. © 2014 Wiley Periodicals, Inc.     

Acceleration of catalyst discovery with easy,fast, and reproducible computational alchemy

The expense of quantum chemistry calculations significantly hinders the search for novel catalysts. Here, we provide a tutorial for using an easy and highly cost‐efficient calculation scheme, called alchemical perturbation density functional theory (APDFT), for rapid predictions of binding energies of reaction intermediates and reaction barrier heights based on the Kohn‐Sham density functional theory (DFT) reference data. We outline standard procedures used in computational catalysis applications, explain how computational alchemy calculations can be carried out for those applications, and then present benchmarking studies of binding energy and barrier height predictions. Using a single OH binding energy on the Pt(111) surface as a reference case, we use computational alchemy to predict binding energies of 32 variations of this system with a mean unsigned error of less than 0.05 eV relative to single‐point DFT calculations. Using a single nudged elastic band calculation for CH₄ dehydrogenation on Pt(111) as a reference case, we generate 32 new pathways with barrier heights having mean unsigned errors of less than 0.3 eV relative to single‐point DFT calculations. Notably, this easy APDFT scheme brings no appreciable computational cost once reference calculations are performed, and this shows that simple applications of computational alchemy can significantly impact DFT‐driven explorations for catalysts. To accelerate computational catalysis discovery and ensure computational reproducibility, we also include Python modules that allow users to perform their own computational alchemy calculations.     

甲烷晶体的晶格能和弹性性质: 不同方法及泛函的评估

通过对甲烷晶体进行结构、晶格能和弹性特性的研究, 评估了不包含和包含色散能量修正的密度泛函理论的性能. 我们分别利用不包含色散能量修正的密度泛函理论(DFT) (包含不同的标准泛函和杂化泛函)和包含色散能量修正的密度泛函理论(DFT-D)计算了甲烷晶体特性, 并与实验作对比. 尽管DFT-D 与传统密度泛函理论及杂化密度泛函理论相比, 修正了甲烷晶体中的范德华(vdW)相互作用, 但是一些修正方案过分修正了这种相互作用. 因此, 人们在使用DFT-D方法时务必谨慎.     

Large changes of static electric properties induced by hydrogen bonding: an ab initio study of linear HCN oligomers

We report the partitioning of the interaction-induced static electronic dipole (hyper)polarizabilities for linear hydrogen cyanide complexes into contributions arising from various interaction energy terms. We analyzed the nonadditivities of the studied properties and used these data to predict the electric properties of an infinite chain. The interaction-induced static electric dipole properties and their nonadditivities were analyzed using an approach based on numerical differentiation of the interaction energy components estimated in an external electric field. These were obtained using the hybrid variational-perturbational interaction energy decomposition scheme, augmented with coupled-cluster calculations, with singles, doubles, and noniterative triples. Our results indicate that the interaction-induced dipole moments and polarizabilities are primarily electrostatic in nature; however, the composition of the interaction hyperpolarizabilities is much more complex. The overlap effects substantially quench the contributions due to electrostatic interactions, and therefore, the major components are due to the induction and exchange-induction terms, as well as the intramolecular electron-correlation corrections. A particularly intriguing observation is that the interaction first hyperpolarizability in the studied systems not only is much larger than the corresponding sum of monomer properties, but also has the opposite sign. We show that this effect can be viewed as a direct consequence of hydrogen-bonding interactions that lead to a decrease of the hyperpolarizability of the proton acceptor and an increase of the hyperpolarizability of the proton donor. In the case of the first hyperpolarizability, we also observed the largest nonadditivity of interaction properties (nearly 17%) which further enhances the effects of pairwise interactions.     

The Nature of the Binding of Au, Ag, and Pd to Benzene, Coronene, and Graphene: From Benchmark CCSD(T) Calculations to Plane-Wave DFT Calculations

The adsorption of Ag, Au, and Pd atoms on benzene, coronene, and graphene has been studied using post Hartree-Fock wave function theory (CCSD(T), MP2) and density functional theory (M06-2X, DFT-D3, PBE, vdW-DF) methods. The CCSD(T) benchmark binding energies for benzene-M (M = Pd, Au, Ag) complexes are 19.7, 4.2, and 2.3 kcal/mol, respectively. We found that the nature of binding of the three metals is different: While silver binds predominantly through dispersion interactions, the binding of palladium has a covalent character, and the binding of gold involves a subtle combination of charge transfer and dispersion interactions as well as relativistic effects. We demonstrate that the CCSD(T) benchmark binding energies for benzene-M complexes can be reproduced in plane-wave density functional theory calculations by including a fraction of the exact exchange and a nonempirical van der Waals correction (EE+vdW). Applying the EE+vdW method, we obtained binding energies for the graphene-M (M = Pd, Au, Ag) complexes of 17.4, 5.6, and 4.3 kcal/mol, respectively. The trends in binding energies found for the benzene-M complexes correspond to those in coronene and graphene complexes. DFT methods that use empirical corrections to account for the effects of vdW interactions significantly overestimate binding energies in some of the studied systems.     

Formation of self-assembled chains of tetrathiafulvalene on a Cu(100) surface

Formation of self-assembled chains of tetrathiafulvalene (TTF) on the Cu(100) surface has been investigated by scanning tunneling microscopy and density functional theory calculations that include semiempirical van der Waals (vdW) interaction corrections. The calculations show that the chain structures observed in the experiments can only be explained by including the vdW interactions. The molecules are tilted along the chain in order to achieve maximal intermolecular interaction. The chains are metastable on the surface, which is consistent with the experimental observation that they disappear after annealing. The fact that all TTF chains observed in the experiment are short might be possibly explained by the interplay between the stabilizing vdW molecule-molecule interaction and the destabilizing rearrangement of surface atoms due to the strong molecule-substrate interaction.     

Energy decomposition analysis based on a block-localized wavefunction and multistate density functional theory

An interaction energy decomposition analysis method based on the block-localized wavefunction (BLW-ED) approach is described. The first main feature of the BLW-ED method is that it combines concepts of valence bond and molecular orbital theories such that the intermediate and physically intuitive electron-localized states are variationally optimized by self-consistent field calculations. Furthermore, the block-localization scheme can be used both in wave function theory and in density functional theory, providing a useful tool to gain insights on intermolecular interactions that would otherwise be difficult to obtain using the delocalized Kohn-Sham DFT. These features allow broad applications of the BLW method to energy decomposition (BLW-ED) analysis for intermolecular interactions. In this perspective, we outline theoretical aspects of the BLW-ED method, and illustrate its applications in hydrogen-bonding and π-cation intermolecular interactions as well as metal-carbonyl complexes. Future prospects on the development of a multistate density functional theory (MSDFT) are presented, making use of block-localized electronic states as the basis configurations.     

Bonds or not bonds? Pancake bonding in 1,2,3,5-dithiadiazolyl and 1,2,3,5-diselenadiazolyl radical dimers and their derivatives

Unusually long bonds or short intermolecular contacts occur in the title compounds reminiscent of pancake bonding. Pancake bonding interactions seem analogous to π-stacking interactions, but they display much shorter contact distances than normally seen in van der Waals (vdW) dimers. The interpretation of these SN and SeN containing structures has been an outstanding challenge for some time. The antibonding (π*) singly occupied molecular orbital (SOMO) of the radical is the source of two-electron multicenter bonding (2e/mc). Preferred conformations thus can be traced back to SOMO-SOMO overlap. We used several computational methods to understand the nature of pancake bonding in the title compounds including four wave function methods (WFT) and a dozen density functional theories (DFT) including empirical dispersion corrections. We used experimental data and high level CCSD(T)/6-311++G(d,p) and MRPT2/6-311++G(d,p) calculations for comparison. The analysis provided the interpretation a wealth of experimental data including conformational preferences of these SN and SeN containing radical dimers leading to a better overall understanding of pancake bonding. Analysis of the various components of the inter-radical interactions showed that SOMO-SOMO bonding interaction and dispersion interaction contribute to the binding energy and neither of these interactions alone is sufficient to bind the dimer. The dimer is predicted to show weak diradical character.     

On the effect of image charges and ion-wall dispersion forces on electric double layer interactions

Two effects of interactions between polarizable ions and polarizable walls in electric double layers are investigated: ionic image charge forces and ion-wall dispersion forces. The first must be included for a consistent treatment of the wall-wall van der Waals (vdW) interaction, since it contains the effect of screening of the static part of the vdW interaction. The second has been suggested to give rise to ion specificity in double layer interactions. The strength of the ion-wall dispersion forces are estimated from quantum mechanical calculations of ionic polarizability and from experimental data for the dielectric functions of the media. The ion density profiles and the anisotropic ion-ion distribution functions in the double layer are calculated in the highly accurate anisotropic hypernetted chain approximation, which allows the correct treatment of the image charge forces. The double layer interactions are evaluated from these distribution functions. It is found that it is important to include both kinds of ion-wall forces. Quantitative and sometimes even qualitative differences occur in the double layer interactions depending on the ionic species of the electrolyte due to different strengths of the ion-wall dispersion interactions.