Electric properties of molecules confined by the spherical harmonic potential

The harmonic oscillator potential is very often used in quantum chemical studies of electric properties to model the effect of spatial confinement. In the vast majority of works, the harmonic potential of cylindrical symmetry was applied. Thus far, its spherical counterpart was used mainly to describe properties of spatially restricted atomic systems. Therefore, our main goal was to study the molecular electric properties in the presence of the spherically symmetric harmonic oscillator potential and to characterize the impact of the relative position of the considered molecules and spherical confinement on these properties. Moreover, we analyzed how the topology of confining environment affects the dipole moment and (hyper)polarizability, by comparing the results obtained in the spherical and cylindrical harmonic potential. Based on the conducted research, it was found that the position of the molecules relative to the spherical confinement strongly influences their electric properties. The observed trends of changes in the electric properties, caused by increasing the confinement strength, vary significantly. Moreover, it was shown that in the vast majority of cases, significant differences in the values of electric properties, obtained in the cylindrical and spherical confinement of a given strength, occur.     

Assessing the performance of density functional theory for the dynamic polarizabilities of amino acids: Treatment of correlation and role of exact exchange

The theoretical determination of electric response properties of the biological systems is a field where the application of density functional theory (DFT) appears to be quite promising. In this work, the performance of 41 density functional methods is evaluated in predicting dynamic polarizabilities of an experimental benchmark set of 20 proteinogenic amino acids. The behavior of a large number of density functionals, including various types of the local spin density approximation (LSDA), generalized gradient approximation (GGA), meta‐GGA (m‐GGA), hybrid‐GGA (h‐GGA), hybrid meta‐GGA (hm‐GGA), and range‐separated hybrid‐GGA (rsh‐GGA), has been assessed for the purpose. Analyzing the results of our DFT benchmarking, we found that these computationally economical methods show very diverse predictive capability and a careful selection of DFT functionals is very important in the polarizability calculations. Considering the role of exchange, correlation, dispersion and long‐range corrections, it turned out that in the LSDA class, SVWN3 gives better results than SPL and SVWN5 toward the reference values. Of the GGA methods, OPBE outperforms all other functionals. The M06‐L is the best method of m‐GGA class. The B3LYP and TPSSh are the best functionals of h‐GGA and hm‐GGA lineages, respectively. Finally, CAM‐B3LYP is the best method of rsh‐GGA functionals that predicts the most accurate polarizability for amino acids by a large margin with respect to others. Overall, the best performing functionals turn out to be hm‐GGAs TPSSh, TPSS1KCIS, M05, tau‐HCTHhyb, and h‐GGA B3LYP. Hopefully, the results of this investigation might provide the useful guidance to propose a new exchange‐correlation functional for calculating the optical properties of biomolecular materials. © 2013 Wiley Periodicals, Inc.     

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

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

CO oxidation catalyzed by a single gold atom: benchmark calculations and the performance of DFT methods

Quantum chemical calculations were carried out on CO oxidation catalyzed by a single gold atom. To investigate the performance of density functional theory (DFT) methods, 42 DFT functionals have been evaluated and compared with high-level wavefunction based methods. It was found that in order to obtain accurate results the functionals used must treat long range interaction well. The double-hybrid mPW2PLYP and B2PLYP functionals are the two functionals with best overall performance. CAM-B3LYP, a long range corrected hybrid GGA functional, also performs well. On the other hand, the popular B3LYP, PW91, and PBE functionals do not show good performance and the performance of the latter two are even at the bottom of the 42 functionals. Our accurate results calculated at the CCSD(T)/aug-cc-pVTZ//mPW2PLYP/aug-cc-pVTZ level of theory indicate that Au atom is a good catalysis for CO oxidation. The reaction follows the following mechanism where CO and O(2) adsorb on Au atom forming an Au(OCOO) intermediate and subsequently O(2) transfer one oxygen atom to CO to form CO(2) and AuO. Then AuO reacts with CO to form another CO(2) to complete the catalytic cycle. The overall energy barrier at 0 K for the first CO oxidation step (Au + CO + O(2)→ AuO + CO(2)) is just 4.8 kcal mol(-1), and that for the second CO oxidation step (AuO + CO → Au + CO(2)) is just 1.6 kcal mol(-1).     

Electric dipole (hyper)polarizabilities of spatially confined LiH molecule

In this study we report on the electronic contributions to the linear and nonlinear static electronic electric dipole properties, namely the dipole moment (μ), the polarizability (α), and the first-hyperpolarizability (β), of spatially confined LiH molecule in its ground X (1)Σ(+) state. The finite-field technique is applied to estimate the corresponding energy and dipole moment derivatives with respect to external electric field. Various forms of confining potential, of either spherical or cylindrical symmetry, are included in the Hamiltonian in the form of one-electron operator. The computations are performed at several levels of approximation including the coupled-cluster methods as well as multi-configurational (full configuration interaction) and explicitly correlated Gaussian wavefunctions. The performance of Kohn-Sham density functional theory for the selected exchange-correlation functionals is also discussed. In general, the orbital compression effects lead to a substantial reduction in all the studied properties regardless of the symmetry of confining potential, however, the rate of this reduction varies depending on the type of applied potential. Only in the case of dipole moment under a cylindrical confinement a gradual increase of its magnitude is observed.     

Parallel implementation of γ‐point pseudopotential plane‐wave DFT with exact exchange

Semi‐local functionals commonly used in density functional theory (DFT) studies of solids usually fail to reproduce localized states such as trapped holes, polarons, excitons, and solitons. This failure is ascribed to self‐interaction which creates a Coulomb barrier to localization. Pragmatic approaches in which the exchange correlation functionals are augmented with small amount of exact exchange (hybrid‐DFT, e.g., B3LYP and PBE0) have shown to promise in rectifying this type of failure, as well as producing more accurate band gaps and reaction barriers. The evaluation of exact exchange is challenging for large, solid state systems with periodic boundary conditions, especially when plane‐wave basis sets are used. We have developed parallel algorithms for implementing exact exchange into pseudopotential plane‐wave DFT program and we have implemented them in the NWChem program package. The technique developed can readily be employed in Γ‐point plane‐wave DFT programs. Furthermore, atomic forces and stresses are straightforward to implement, making it applicable to both confined and extended systems, as well as to Car‐Parrinello ab initio molecular dynamic simulations. This method has been applied to several systems for which conventional DFT methods do not work well, including calculations for band gaps in oxides and the electronic structure of a charge trapped state in the Fe(II) containing mica, annite. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Energy landscapes of nucleophilic substitution reactions: a comparison of density functional theory and coupled cluster methods

We have carried out a detailed evaluation of the performance of all classes of density functional theory (DFT) for describing the potential energy surface (PES) of a wide range of nucleophilic substitution (SN2) reactions involving, amongst others, nucleophilic attack at carbon, nitrogen, silicon, and sulfur. In particular, we investigate the ability of the local density approximation (LDA), generalized gradient approximation (GGA), meta-GGA as well as hybrid DFT to reproduce high-level coupled cluster (CCSD(T)) benchmarks that are close to the basis set limit. The most accurate GGA, meta-GGA, and hybrid functionals yield mean absolute deviations of about 2 kcal/mol relative to the coupled cluster data, for reactant complexation, central barriers, overall barriers as well as reaction energies. For the three nonlocal DFT classes, the best functionals are found to be OPBE (GGA), OLAP3 (meta-GGA), and mPBE0KCIS (hybrid DFT). The popular B3LYP functional is not bad but performs significantly worse than the best GGA functionals. Furthermore, we have compared the geometries from several density functionals with the reference CCSD(T) data. The same GGA functionals that perform best for the energies (OPBE, OLYP), also perform best for the geometries with average absolute deviations in bond lengths of 0.06 A and 0.6 degrees, even better than the best meta-GGA and hybrid functionals. In view of the reduced computational effort of GGAs with respect to meta-GGAs and hybrid functionals, let alone coupled cluster, we recommend the use of accurate GGAs such as OPBE or OLYP for the study of SN2 reactions.     

General performance of density functionals

The density functional theory (DFT) foundations date from the 1920s with the work of Thomas and Fermi, but it was after the work of Hohenberg, Kohn, and Sham in the 1960s, and particularly with the appearance of the B3LYP functional in the early 1990s, that the widespread application of DFT has become a reality. DFT is less computationally demanding than other computational methods with a similar accuracy, being able to include electron correlation in the calculations at a fraction of time of post-Hartree-Fock methodologies. In this review we provide a brief outline of the density functional theory and of the historic development of the field, focusing later on the several types of density functionals currently available, and finishing with a detailed analysis of the performance of DFT across a wide range of chemical properties and system types, reviewed from the most recent benchmarking studies, which encompass several well-established density functionals together with the most recent efforts in the field. Globally, an overall picture of the level of performance of the plethora of currently available density functionals for each chemical property is drawn, with particular attention being dedicated to the relative performance of the popular B3LYP density functional.     

DFT performance for the hole transfer parameters in DNA π stacks

Recently, we showed that unoccupied Kohn‐Sham (KS) orbitals stemming from DFT calculations of a neutral system can be used to derive accurate estimates of the free energy and electronic couplings for excess electron transfer in DNA (Félix and Voityuk, J Phys Chem A 2008, 112, 9043). In this article, we consider the propagation of radical cation states (hole transfer) through DNA π‐stacks and compare the performance of different exchange‐correlation functionals to estimate the hole transfer (HT) parameters. Two different approaches are used: (1) calculations that use occupied KS orbitals of neutral π stacks of nucleobases, and (2) the time‐dependent DFT method which is applied to the radical cation states of these stacks. Comparison of the calculated parameters with the reference data suggests that the best results are provided by the KS scheme with hybrid functionals (B3LYP, PBE0, and BH&HLYP). The TD DFT approach gives significantly less accurate values of the HT parameters. In agreement with high‐level ab initio results, the KS scheme predicts that the hole in π stacks is confined to a single nucleobase; in contrast, the spin‐unrestricted DFT method considerably overestimates the hole delocalization in the radical cations. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Rung 3.5 density functionals: Another step on Jacob's ladder

Applications of density functional theory (DFT) to computational chemistry and solid‐state physics rely on a “Jacob's Ladder” of progressively more complicated approximations to the many‐body exchange‐correlation (XC) density functional. Accurate, computationally tractable DFT calculations on large and periodic systems remain challenging for existing XC functionals. Simple XC functionals on the three lowest rungs of Jacob's Ladder are insufficiently accurate for many properties, while fourth‐rung hybrid functionals incorporating nonlocal information can be prohibitively expensive. This perspective presents our work toward a compromise, a new class of “Rung 3.5” functionals that incorporate a linear dependence on the nonlocal one‐particle density matrix. This work reviews these functionals' formal underpinning, numerical performance, and prospects for modeling solids and surfaces. © 2012 Wiley Periodicals, Inc.     

Ab initio studies on the polarizability of lithium clusters: Some unusual results

The polarizability of lithium clusters was studied using ab initio as well as Density Functional Theory (DFT)‐based methods. The electron correlation was found not to play any major role in the evaluation of polarizability of the lithium clusters, contrary to the case of sodium and other types of clusters. The theoretical methods used in the present study are not able to reproduce the unusual features in the size evolution of the polarizability observed experimentally. An interesting feature is that the polarizability of lithium clusters calculated by various methods at 0K temperature is found to be higher than the experimental results measured at higher temperature, which is quite unusual. The significance of these issues is highlighted and the need of highly accurate DFT exchange‐correlation functionals and ab initio methods in the study of the electronic properties of lithium clusters is emphasized. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Is the uniform electron gas limit important for small Ag clusters? Assessment of different density functionals for Ag(n) (n < or = 4)

Twenty-three density functional theory (DFT) methods, including the second- and the third-generation functionals, are tested in conjunction with two basis sets (LANL2DZ and SDD) for studying the properties of neutral and ionic silver clusters. We find that DFT methods incorporating the uniform electron gas limit in the correlation functional, namely, those with Perdew's correlation functionals (PW91, PBE, P86, and TPSS), Becke's B95, and the Van Voorhis-Scuseria functional VSXC, generally perform better than the other group of functionals, e.g., those incorporating the LYP correlation functional and variations of the B97 functional. Strikingly, these two groups of functionals can produce qualitatively different results for the Ag3 and Ag4 clusters. The energetic properties and vibrational frequencies of Ag(n) are also evaluated by the different functionals. The present study shows that the choice of DFT methods for heavy metals may be critical. It is found that the exact-exchange-incorporated PBE functional (PBE1PBE) is among the best for predicting the range of properties.     

On the accuracy of density-functional theory exchange-correlation functionals for H bonds in small water clusters: benchmarks approaching the complete basis set limit

The ability of several density-functional theory (DFT) exchange-correlation functionals to describe hydrogen bonds in small water clusters (dimer to pentamer) in their global minimum energy structures is evaluated with reference to second order Moller-Plesset perturbation theory (MP2). Errors from basis set incompleteness have been minimized in both the MP2 reference data and the DFT calculations, thus enabling a consistent systematic evaluation of the true performance of the tested functionals. Among all the functionals considered, the hybrid X3LYP and PBE0 functionals offer the best performance and among the nonhybrid generalized gradient approximation functionals, mPWLYP and PBE1W perform best. The popular BLYP and B3LYP functionals consistently underbind and PBE and PW91 display rather variable performance with cluster size.     

用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校正也能得到同样甚至更准确的结果.     

A density functional study of chemical reactions

Density functional theory (DFT ) was used to study reactions involving small molecules. Relative energies of isomers and transition structures of diazene, formaldehyde, and methylenimine were determined using various DFT functionals and results were compared with MP 2 and MP 4 calculations. DFT reaction barriers were found to be consistently lower. For some reactions, such as OH + H₂→ H₂O + H, gradient-corrected functionals predict very low or nonexistent barriers. The hybrid Hartree–Fock–DFT adiabatic connection method (ACM ) often provides much better results in such cases. The performance of several density functionals, including ACM , was tested in calculations on over 100 atomization, hydrogenation, bond dissociation, and isodesmic reactions. The ACM functional provides consistently better geometries and reaction energetics than does any other functional studied. In cases where both HF and gradient-corrected DFT methods underestimate bond distances, the ACM geometries may be inferior to those predicted by gradient-corrected DFT methods. © 1995 John Wiley & Sons, Inc.     

Exact density functionals for two-electron systems in an external magnetic field

In principle, the extension of density functional theory (DFT) to Coulombic systems in a nonvanishing magnetic field is via current DFT (CDFT). Though CDFT is long established formally, relatively little is known with respect to any generally applicable, reliable approximate E(XC) and A(XC) functionals analogous with the workhorse approximate functionals (local density approximation and generalized gradient approximation) of ordinary DFT. Progress can be aided by having benchmark studies on a solvable correlated system. At zero field, the best-known finite system for such purposes is Hooke's atom. Recently we extended the exact ground state solutions for this two-electron system to certain combinations of nonzero external magnetic fields and confinement strengths. From those exact solutions, as well as high-accuracy numerical results for other field and confinement combinations, we construct the correlated electron density and paramagnetic current density, the exact Kohn-Sham orbitals, and the exact DFT and CDFT exchange-correlation energies and potentials. We compare with results from several widely used approximate functionals, all of which exhibit major qualitative failures, whether in CDFT or in naive application of ordinary DFT. We also illustrate how the CDFT vorticity variable nu is a computationally difficult quantity which may not be appropriate in practice to describe the external B field effects on E(XC) and A(XC).     

Insight and performance of LC-DFT vs DFT in the NMR shielding and chemical shift calculations: Case of CHClCH?CF3

Fifteen density functional theory (DFT) methods and fifteen long-range corrected density functional theory (LC-DFT) methods were used in the present work to assess nuclear magnetic resonance parameters such as nuclear shielding constant (NSC), nuclear chemical shift (NCS), and nuclear anisotropic shielding constant (NAS). These different methods were associated with the full basis set 6-311++G(3df,3pd). The gauge-independent atomic orbital was used for the calculation of nuclear shielding tensors of the nuclei contained in the stereoisomers cis- and trans-CHClCHCF₃. Thus, the effects of LC are clearly observed for heavy nuclei (¹³C, ¹⁹F, ³⁵Cl). The results of NSC, NCS, and NAS from DFT are better described than LC-DFT with regard to the KT3 method. Moreover, the results from the LC-DFT are better described than the standard DFT with regard to CCSD(T). Based on the latter method used as the benchmark, the NSCs of nuclei are well fitted by the competitive functionals LC-TPSSTPSS and LC-PKZBPKZB. In the particular case of the trans-isomer, mPWPKZB was found to be the best method. For the NCSs, the more accurate methods include the latter two LC functionals and the non-LC functionals TPSSTPSS and mPWPKZB. The accuracy of NAS depends strongly on the nuclei. Thus, CAM-B3LYP describes it well for ¹⁹F and LC-PKZBPKZB for ³⁵Cl. The rest of nuclei are well fitted by all the methods except ¹³C₁ and ¹³C₂, which are better reproduced by the LC-DFT except the LC-PKZBPKZB, LC-TPSSTPSS, and CAM-B3LYP functionals.     

Scaling dynamical correlation energy from density functional theory correlation functionals

The scaling of dynamical correlation energy in molecules obtained by the correlation functionals of density functional theory (DFT) is examined. The approach taken is very similar to the scaled external correlation method of Brown and Truhlar but is based on the observation that DFT correlation functionals, especially the LYP, appear to represent the dynamical portion of the correlation energy in molecules. We examine whether higher accuracy in atomization energies can be gained by scaling without significant deterioration of the structural and spectroscopic properties of the molecules using four DFT functionals (BLYP, OLYP, B3LYP, and O3LYP) on 19 molecules including the six molecule AE6 database, the latter being representative of a much larger, 109 molecule training set. We show that, with molecule specific scale factors, nearly perfect agreement with experiment can be achieved in atomization energies without increasing the average errors in other molecular properties relative to the DFT calculation. We further show that it is possible to find optimal scale factors which reduce the mean unsigned error per bond to levels comparable to those of some multilevel multicoefficient methods.