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

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

Testing the TPSS meta-generalized-gradient-approximation exchange-correlation functional in calculations of transition states and reaction barriers

We have studied the performance of local and semilocal exchange-correlation functionals [meta-generalized-gradient-approximation (GGA)-TPSS, GGA-Perdew-Burke-Ernzerhof (PBE), and local density approximation (LDA)] in the calculation of transition states, reaction energies, and barriers for several molecular and one surface reaction, using the plane-wave pseudopotential approach. For molecular reactions, these results have been compared to all-electron Gaussian calculations using the B3LYP hybrid functional, as well as to experiment and high level quantum chemistry calculations, when available. We have found that the transition state structures are accurately identified irrespective of the level of the exchange-correlation functional, with the exception of a qualitatively incorrect LDA prediction for the H-transfer reaction in the hydrogen bonded complex between a water molecule and a OH radical. Both the meta-GGA-TPSS and the GGA-PBE functionals improve significantly the calculated LDA barrier heights. The meta-GGA-TPSS further improves systematically, albeit not always sufficiently, the GGA-PBE barriers. We have also found that, on the Si(001) surface, the meta-GGA-TPSS barriers for hydrogen adsorption agree significantly better than the corresponding GGA-PBE barriers with quantum Monte Carlo cluster results and experimental estimates.     

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.     

Dynamic Simulation on Surface Hydration and Dehydration of Monoclinic Zirconia

The commonly used oxide-supported metal catalysts are usually prepared in aqueous phase, which then often need to undergo calcination before usage. Therefore, the surface hydration and dehydration of oxide supports are critical for the realistic modeling of supported metal catalysts. In this work, by ab initio molecular dynamics (AIMD) simulations, the initial anhydrous monoclinic ZrO\begin{document}$_2$\end{document}(111) surfaces are evaluated within explicit solvents in aqueous phase at mild temperatures. During the simulations, all the two-fold-coordinated O sites will soon be protonated to form the acidic hydroxyls (HO\begin{document}$_{\rm{L}}$\end{document}), remaining the basic hydroxyls (HO^*) on Zr. The basic hydroxyls (HO^*) can easily diffuse on surfaces via the active proton exchange with the undissociated adsorption water (H\begin{document}$_2$\end{document}O^*). Within the temperatures ranging from 273 K to 373 K, in aqueous phase a certain representative equilibrium hydrated m-ZrO\begin{document}$_2$\end{document}(111) surface is obtained with the coverage (\begin{document}$\theta$\end{document}) of 0.75 on surface Zr atoms. Later, free energies on the stepwise surface water desorption are calculated by density functional theory to mimic the surface dehydration under the mild calcination temperatures lower than 800 K. By obtaining the phase diagrams of surface dehydration, the representative partially hydrated m-ZrO\begin{document}$_2$\end{document}(111) surfaces (0.25\begin{document}$\leq$\end{document}\begin{document}$\theta$\end{document} < 0.75) at various calcination temperatures are illustrated. These hydrated m-ZrO\begin{document}$_2$\end{document}(111) surfaces can be crucial and readily applied for more realistic modeling of ZrO\begin{document}$_2$\end{document} catalysts and ZrO\begin{document}$_2$\end{document}-supported metal catalysts.     

Influence of Van der Waals Interactions on the Solvation Energies of Adsorbates at Pt-Based Electrocatalysts

Solvation can significantly modify the adsorption energy of species at surfaces, thereby influencing the performance of electrocatalysts and liquid-phase catalysts. Thus, it is important to understand adsorbate solvation at the nanoscale. Here we evaluate the effect of van der Waals (vdW) interactions described by different approaches on the solvation energy of *OH adsorbed on near-surface alloys (NSAs) of Pt. Our results show that the studied functionals can be divided into two groups, each with rather similar average *OH solvation energies: (1) PBE and PW91; and (2) vdW functionals, RPBE, PBE-D3 and RPBE-D3. On average, *OH solvation energies are less negative by ∼0.14 eV in group (2) compared to (1), and the values for a given alloy can be extrapolated from one functional to another within the same group. Depending on the desired level of accuracy, these concrete observations and our tabulated values can be used to rapidly incorporate solvation into models for electrocatalysis and liquid-phase catalysis.     

van der Waals corrected density functional calculations of the adsorption of benzene on the Cu (111) surface

We investigate the performance of several van der Waals (vdW) functionals at calculating the interactions between benzene and the copper (111) surface, using the local orbital approach in the SIESTA code. We demonstrate the importance of using surface optimized basis sets to calculate properties of pure surfaces, including surface energies and the work function. We quantify the errors created using (3 × 3) supercells to study adsorbate interactions using much larger supercells, and show non‐negligible errors in the binding energies and separation distances. We examine the eight high‐symmetry orientations of benzene on the Cu (111) surface, reporting the binding energies, separation distance, and change in work function. The optimized vdW‐DF(optB88‐vdW) functional provides superior results to the vdW‐DF(revPBE) and vdW‐DF2(rPW86) functionals, and closely matches the experimental and experimentally deduced values. This work demonstrates that local orbital methods using appropriate basis sets combined with a vdW functional can model adsorption between metal surfaces and organic molecules.     

Assessing Density Functionals for Describing Methane Dissociative Chemisorption on Pt(110)-(2×1) Surface

In this work, we explore the suitability of several density functionals with the generalized gradient approximation (GGA) and beyond for describing the dissociative chemisorption of methane on the reconstructed Pt(110)-(2\begin{document}$ \times $\end{document}1) surface. The bulk and surface structures of the metal, methane adsorption energy, and dissociation barrier are used to assess the functionals. A van der Waals corrected GGA functional (optPBE-vdW) and a meta-GGA functional with van der Waals correction (MS PBEl-rVV10) are selected for ab initio molecular dynamics calculations of the sticking probability. Our results suggest that the use of these two functionals may lead to a better agreement with existing experimental results, thus serving as a good starting point for future development of reliable machine-learned potential energy surfaces for the dissociation of methane on the Pt(110)-(2\begin{document}$ \times $\end{document}1) surface.     

Six-dimensional potential energy surface for H2 at Ru(0001)

The six-dimensional (6D) potential energy surface (PES) for the H(2) molecule interacting with a clean Ru(0001) surface has been computed accurately for the first time. Density functional theory (DFT) and a pseudopotential based periodic plane-wave approach have been used to calculate the electronic interactions between the molecule and the surface. Two different generalized gradient approximation (GGA) exchange-correlation functionals, PW91 and RPBE, have been adopted. Based on the DFT/GGA calculated potential energies, an analytical 6D PES has been constructed using the corrugation reducing procedure. A very accurate representation of the DFT/GGA data has been achieved, with an average error in the interpolation of about 3 meV and a maximum error not larger than about 30 meV. The top site is found to be the most reactive site for both functionals used, but PW91 predicts a higher reactivity than RPBE, with lower-energy and earlier-located dissociation barriers. The energetic corrugation displayed by the RPBE PES is larger than the PW91 PES while the geometric corrugation is smaller. The differences between the two PESs increase as the distance of the molecular center of mass to the surface decreases. A direct comparison with experimental investigations on H(2)/Ru(0001) could shed light on the suitability of these XC potentials often used in DFT calculations.     

First-Principles Evaluation of One-Dimensional Metal-Organic Frameworks for Electrocatalytic C−H Activation of Natural Gas

To replace the oxygen evolution reaction with thermodynamically more favorable and economically more profitable methane and ethane (the major components of natural gas) electrochemical partial oxidation, we employed constant electrode potential density functional theory calculations to screen 20 one-dimensional metal-organic frameworks containing heteroatom-substituted benzene as electrocatalysts. By computing the Pourbaix diagrams, O−H binding energies, and C−H activation barriers, we determined that although none of these catalysts were able to activate methane, one was able to hydroxylate ethane to ethanol with facile kinetics, making it a promising electrocatalyst for natural gas oxidation.     

Structure and properties of metal-exchanged zeolites studied using gradient-corrected and hybrid functionals. III. Energetics and vibrational spectroscopy of adsorbates

The influence of the exchange-correlation functional (semilocal gradient corrected or hybrid functional) on density-functional studies of the adsorption of CO and NO in Cu- and Co-exchanged chabazite has been investigated, extending the studies of the structural and electronic properties of these materials [F. Go?ltl and J. Hafner, J. Chem. Phys. 136, 064501 (2012); 136, 064502 (2012)] and including for comparison carbonyls and nitrosyls of Cu and Co. Hybrid functionals predict much lower adsorption energies than conventional semilocal functionals, in better agreement with experiment as far as data are available for comparison. The calculated adsorption energies show a strong linear correlation with the stability of the cation sites. For Cu(I)-chabazite the calculated adsorption energies span almost the interval between the adsorption energies calculated for pure neutral and positively charged Cu-carbonyls and nitrosyls. For divalent Cu(II) and Co(II) the adsorption energies at cations in chabazite are much lower than the metal-molecule binding energies in the free carbonyls or nitrosyls, especially for the most stable cation location in a six-membered ring of the chabazite structure. For the stretching modes of adsorbed CO only hybrid functionals reproduce the blueshift of the frequency reported for all Cu(I)- and Co(II)-zeolites. For Cu(II)-chabazite both types of functionals predict a blueshift, the larger value calculated with hybrid functionals being in better agreement with observation. For NO adsorbed on Cu(I)-chabazite all functionals produce a redshift, the smaller value derived with hybrid functionals being in better agreement with experiment. For NO adsorbed in Cu(II)- and Co(II)-chabazite gradient-corrected functionals produce the best agreement with experiment for cations located in a six-membered ring. Semilocal functionals tend to underestimate the frequencies, while hybrid functionals tend to overestimate. The decisive factors determining the influence of the functionals are the larger HOMO-LUMO gap and the larger bandgap of the zeolite host, as well as the larger exchange-splitting of the cation eigenstates predicted with hybrid functionals. For Co(II)-chabazite the tendency to overestimate the exchange-splitting and to stabilize a high-spin state lead to better results with semilocal functionals. Finally, a comprehensive discussion of the influence of the exchange-correlation functional on the physico-chemical properties of these complex systems, based all three papers of this series is presented.     

Assessment of a Coulomb-attenuated exchange-correlation energy functional

The recently proposed CAM-B3LYP exchange-correlation energy functional, based on a partitioning of the r operator in the exchange interaction into long- and short-range components, is assessed for the determination of molecular thermochemistry, structures, and second order response properties. Rydberg and charge transfer excitation energies and static electronic polarisabilities are notably improved over the standard B3LYP functional; classical reaction barriers also improve. Ionisation potentials, bond lengths, NMR shielding constants and indirect spin-spin coupling constants are comparable with the two functionals. CAM-B3LYP atomisation energies and diatomic harmonic vibrational wavenumbers are less accurate than those of B3LYP. Future research directions are outlined.     

Regional self-interaction correction of density functional theory

We propose a new simple scheme for self-interaction correction (SIC) of exchange functionals in the density functional theory. In the new scheme, exchange energies are corrected by substituting exchange self-interactions for exchange functionals in regions of self-interaction. To classify the regions of self-interaction, we take advantage of the property of the total kinetic energy density approaching the Weizs?cker density in the case of electrons in isolated orbitals. The scheme differs from conventional SIC methods in that it produces optimized molecular structures. Applying the scheme to the calculation of reaction energy barriers showed that it provides a clear improvement in cases where the barriers are underestimated by conventional "pure" functionals. In particular, we found that this scheme even reproduces a transition state that is not given by pure functionals.     

Construction of a generalized gradient approximation by restoring the density-gradient expansion and enforcing a tight Lieb-Oxford bound

Recently, a generalized gradient approximation (GGA) to the density functional, called PBEsol, was optimized (one parameter) against the jellium-surface exchange-correlation energies, and this, in conjunction with changing another parameter to restore the first-principles gradient expansion for exchange, was sufficient to yield accurate lattice constants of solids. Here, we construct a new GGA that has no empirical parameters, that satisfies one more exact constraint than PBEsol, and that performs 20% better for the lattice constants of 18 previously studied solids, although it does not improve on PBEsol for molecular atomization energies (a property that neither functional was designed for). The new GGA is exact through second order, and it is called the second-order generalized gradient approximation (SOGGA). The SOGGA functional also differs from other GGAs in that it enforces a tighter Lieb-Oxford bound. SOGGA and other functionals are compared to a diverse set of lattice constants, bond distances, and energetic quantities for solids and molecules (this includes the first test of the M06-L meta-GGA for solid-state properties). We find that classifying density functionals in terms of the magnitude mu of the second-order coefficient of the density gradient expansion of the exchange functional not only correlates their behavior for predicting lattice constants of solids versus their behavior for predicting small-molecule atomization energies, as pointed out by Perdew and co-workers [Phys. Rev. Lett. 100, 134606 (2008); Perdew ibid. 80, 891 (1998)], but also correlates their behavior for cohesive energies of solids, reaction barriers heights, and nonhydrogenic bond distances in small molecules.     

Electronic Structure Characteristics of ESIPT and TICT Fluorescence Emissions and Calculations of Emitting Energies

A scheme of time-dependent density functional theory (TDDFT) combined with the single-excitation configuration interaction (CIS) approach was employed to investigate the first excited singlet state (S₁) for eight salicylanilide derivatives and analogues, which have similar structural formulas. The results showed that fluorescence-emitting mechanisms of these molecules were in two distinct manners (excited-state intramolecular proton transfer (ESIPT) and twisted intramolecular charge transfer (TICT)), which agreed with the well-known experiments. For ESIPT compounds with inconspicuous charge transfer (CT) during electron transition, pure functionals without Hartree-Fock (HF) exchange energy, such as OLYP and BLYP, were suitable to calculate emitting energies. For TICT compounds with large CT during electron transition, hybrid functionals with about 37% HF exchange energy, such as mPW1B95 and MPW1K, performed well. On condition that the exchange-correlation (XC) functionals were chosen properly according to the rules above, reliable emitting energies for salicylanilide derivatives and analogues could be obtained at the TDDFT/6-31G(d)//CIS/3-21G(d) level. The average accuracy reached about 0.2 eV. For the salicylanilides with both proton transfer (PT) and CT reaction channels, only one channel occurred actually according to the principle of energy minimum. This actual reaction decided proper XC functionals, whereas the reaction that did not occur actually was trivial. Eight appendent compounds were calculated to prove that this successful scheme is expected to be suitable for other ESIPT and TICT compounds.     

Density functional theory study of water activation and COads + OHads reaction on pure platinum and bimetallic platinum/ruthenium nanoclusters

A density functional theory study of the elementary steps that lead to the removal of CO(ads(Pt)) over alloyed and sequentially deposited Pt/Ru bimetallic nanoclusters is presented. The reaction energies and activation barriers for the H2O(ads(Ru)) dissociation and CO(ads(Pt)) + OH(ads(Ru)) reaction are estimated in solid-gas interface and in a microsolvated environment to determine which surface morphology is more tolerant to COads poisoning. On the basis of the energetics, the sequentially deposited Pt/Ru nanocluster is predicted to be a much more promising anode catalyst than the alloy cluster surface in fuel cell applications.