一些密度泛函方法预测电子亲和势

选取了杂化泛函B3LYP, B3PW91, O3LYP, PBE0, 以及与之相对应的GGA泛函BLYP, BPW91, OLYP和PBE, 还选取了能更好地兼顾强相互作用和弱相互作用的X3LYP泛函和在预测NMR的化学位移有较好表现的OPBE泛函, 以及两种meta-GGA泛函VSXC和TPSS, 共12种泛函, 详细地考察了这些泛函在预测EA方面的准确性.     

Systematic investigation on the geometric dependence of the calculated nuclear magnetic shielding constants

We perform a systematic examination on the dependence of the calculated nuclear magnetic shielding constants on the chosen geometry for a selective set of density functional methods of B3LYP, PBE0, and OPBE. We find that the OPBE exchange-correlation functional performs remarkably well when either the optimized geometries or the experimental geometries are used. The popular B3LYP and PBE0 functionals have a clear tendency of deshielding, giving shieldings that are usually too low and shifts that are usually too high, at the experimental geometries. Combined with the Hartree-Fock geometries, however, much improved magnetic constants are obtained for B3LYP and PBE0, due to the compensation effect from the systematic underestimation of bond lengths by the Hartree-Fock method.     

Performance of six functionals (LDA, PBE, PBESOL, B3LYP, PBE0, and WC1LYP) in the simulation of vibrational and dielectric properties of crystalline compounds. The case of forsterite Mg2SiO4

The performance of six different density functionals (LDA, PBE, PBESOL, B3LYP, PBE0, and WC1LYP) in describing the infrared spectrum of forsterite, a crystalline periodic system with orthorhombic unit cell (28 atoms in the primitive cell, Pbmn space group), is investigated by using the periodic ab initio CRYSTAL09 code and an all‐electron Gaussian‐type basis set. The transverse optical (TO) branches of the 35 IR active modes are evaluated at the equilibrium geometry together with the oscillator strengths and the high‐frequency dielectric tensor ?_∞. These quantities are essential to compute the dielectric function ?(ν), and then the reflectance spectrum R(ν), which is compared with experiment. It turns out that hybrid functionals perform better than LDA and GGA, in general; that B3LYP overperforms WC1LYP and, in turn, PBE0; that PBESOL is better than PBE; that LDA is the worst performing functional among the six under study. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Benchmarking density functionals in conjunction with Grimme's dispersion correction for noble gas dimers (Ne2, Ar2, Kr2, Xe2, Rn2)

Eleven exchange‐correlational functionals of different types corrected for dispersion by Grimme's D3 correction in conjunction with the aug‐cc‐pVTZ basis set were tested on the following noble gas (Ng) dimers: Ne₂, Ar₂, Kr₂, Xe₂, and Rn₂. For comparison, the D2 and D3BJ corrections were probed with the B3LYP functional. From post‐HF wavefunction methods, CCSD(T) theory was also included. The investigated properties involved potential energy curves, equilibrium bond distances, and interaction energies. The B3LYP‐D3, B3LYP‐D3BJ, and PBE0‐D3 functionals performed overall best for bond distances, while B3LYP‐D3 and B97‐D3 performed best for interaction energies. The importance of fortunate error cancellations was seen in the often reduced agreement with reference data upon correction for BSSE. As several functionals performed well selectively for some noble gases (and poorly for others), we also analysed the performance on the Ng₂ dimers individually and recommended DFT‐D3 functionals for the calculation of large clusters of each Ng.     

Computational study of static first hyperpolarizability of donor-acceptor substituted (E)-benzaldehyde phenylhydrazone

Various hybrid functionals (B3LYP, B97-2, PBE0, BMK, BH&HLYP, CAM-B3LYP, and LC-ωPBE) implemented in density functional theory were applied to give estimate of static first hyperpolarizabilty (β(0)) of (E)-benzaldehyde phenylhydrazone designated as (E)-BPH. Against those of MP2 computations as a function of the underlying density functional, good agreement was obtained with the BH&HLYP and CAM-B3LYP functionals. The LC-ωPBE functional and the B3LYP, PBE0, B97-2, and BMK functionals underestimated and overestimated β(0), respectively. The basis set effect on the calculated β(0) was also investigated. It turned out that the 6-311+G(2d,p) basis set provided excellent converged value of β(0). On the basis of the calculated results, we investigated the substituent effect on β(0) of donor-acceptor (D-A) substituted (E)-BPH systematically by using the BH&HLYP and CAM-B3LYP computations with the 6-311+G(2d,p) basis set. We proposed a Zwitterion structure to explain the calculated trend in the substituent effect and the enhanced hyperpolarizability of type II compounds (A-(E)-BPH-D) than type I compounds (D-(E)-BPH-A). Natural bonding orbital analysis carried out at BH&HLYP/6-311+G(2d,p)//B3LYP/6-31G(2df,p) level of theory substantiated the claim.     

Optical rotation calculated with time-dependent density functional theory: the OR45 benchmark

Time-dependent density functional theory (TDDFT) computations are performed for 42 organic molecules and three transition metal complexes, with experimental molar optical rotations ranging from 2 to 2 × 10(4) deg cm(2) dmol(-1). The performances of the global hybrid functionals B3LYP, PBE0, and BHLYP, and of the range-separated functionals CAM-B3LYP and LC-PBE0 (the latter being fully long-range corrected), are investigated. The performance of different basis sets is studied. When compared to liquid-phase experimental data, the range-separated functionals do, on average, not perform better than B3LYP and PBE0. Median relative deviations between calculations and experiment range from 25 to 29%. A basis set recently proposed for optical rotation calculations (LPol-ds) on average does not give improved results compared to aug-cc-pVDZ in TDDFT calculations with B3LYP. Individual cases are discussed in some detail, among them norbornenone for which the LC-PBE0 functional produced an optical rotation that is close to available data from coupled-cluster calculations, but significantly smaller in magnitude than the liquid-phase experimental value. Range-separated functionals and BHLYP perform well for helicenes and helicene derivatives. Metal complexes pose a challenge to first-principles calculations of optical rotation.     

On the performance of quantum chemical methods to predict solvatochromic effects: the case of acrolein in aqueous solution

The performance of the Hartree-Fock method and the three density functionals B3LYP, PBE0, and CAM-B3LYP is compared to results based on the coupled cluster singles and doubles model in predictions of the solvatochromic effects on the vertical n-->pi* and pi-->pi* electronic excitation energies of acrolein. All electronic structure methods employed the same solvent model, which is based on the combined quantum mechanics/molecular mechanics approach together with a dynamical averaging scheme. In addition to the predicted solvatochromic effects, we have also performed spectroscopic UV measurements of acrolein in vapor phase and aqueous solution. The gas-to-aqueous solution shift of the n-->pi* excitation energy is well reproduced by using all density functional methods considered. However, the B3LYP and PBE0 functionals completely fail to describe the pi-->pi* electronic transition in solution, whereas the recent CAM-B3LYP functional performs well also in this case. The pi-->pi* excitation energy of acrolein in water solution is found to be very dependent on intermolecular induction and nonelectrostatic interactions. The computed excitation energies of acrolein in vacuum and solution compare well to experimental data.     

On the applicability of density functional theory to manganese‐based complexes with catalytic activity toward water oxidation

The present contribution assesses the performance of several popular and accurate density functionals, namely B3LYP, BP86, M06, MN12L, mPWPW91, PBE0, and TPSSh toward manganese‐based coordination complexes. These compounds show promising properties toward application to catalytic water oxidation. Although manganese with N‐ and O‐biding ligands tends to give rise to high spin complexes, the results show that BP86, mPWPW91, and specially MN12L, tend to yield low‐spin complexes. The usage of these functionals for such compounds is, thus, discouraged. All the functionals considered deliver accurate geometries. The present results show, however, that B3LYP delivers geometries deviating from experimental values when compared to the other functionals of the set. M06, PBE0, and TPSSh deliver geometries of similar accuracy, PBE0 outstanding slightly with respect to the other two. © 2017 Wiley Periodicals, Inc.     

Performance of density functional theory in computing nonresonant vibrational (hyper)polarizabilities

A set of exchange‐correlation functionals, including BLYP, PBE0, B3LYP, BHandHLYP, CAM‐B3LYP, LC‐BLYP, and HSE, has been used to determine static and dynamic nonresonant (nuclear relaxation) vibrational (hyper)polarizabilities for a series of all‐trans polymethineimine (PMI) oligomers containing up to eight monomer units. These functionals are assessed against reference values obtained using the Møller–Plesset second‐order perturbation theory (MP2) and CCSD methods. For the smallest oligomer, CCSD(T) calculations confirm the choice of MP2 and CCSD as appropriate for assessing the density functionals. By and large, CAM‐B3LYP is the most successful, because it is best for the nuclear relaxation contribution to the static linear polarizability, intensity‐dependent refractive index second hyperpolarizability, static second hyperpolarizability, and is close to the best for the electro‐optical Pockels effect first hyperpolarizability. However, none of the functionals perform satisfactorily for all the vibrational (hyper)polarizabilities studied. In fact, in the case of electric field‐induced second harmonic generation all of them, as well as the Hartree–Fock approximation, yield the wrong sign. We have also found that the Pople 6–31+G(d) basis set is unreliable for computing nuclear relaxation (hyper)polarizabilities of PMI oligomers due to the spurious prediction of a nonplanar equilibrium geometry. © 2013 Wiley Periodicals, Inc.     

Structure and stability of the (001) alpha-quartz surface

The structure and surface energies of the cleaved, reconstructed, and fully hydroxylated (001) alpha-quartz surface of various thicknesses are investigated with periodic density functional theory (DFT). The properties of the cleaved and hydroxylated surface are reproduced with a slab thickness of 18 atomic layers, while a thicker 27-layer slab is necessary for the reconstructed surface. The performance of the hybrid DFT functional B3LYP, using an atomic basis set, is compared with the generalised gradient approximation, PBE, employing plane waves. Both methodologies give similar structures and surface energies for the cleaved and reconstructed surfaces, which validates studying these surfaces with hybrid DFT. However, there is a slight difference between the PBE and B3LYP approach for the geometry of the hydrogen bonded network on the hydroxylated surface. The PBE adsorption energy of CO on a surface silanol site is in good agreement with experimental values, suggesting that this method is more accurate for hydrogen bonded structures than B3LYP. New hybrid functionals, however, yield improved weak interactions. Since these functionals also give superior activation energies, we recommend applying the new functionals to contemporary issues involving the silica surface and adsorbates on this surface.     

Performance of recently developed kinetic energy density functionals for the calculation of hydrogen binding strengths and hydrogen-bonded structures

 The accuracy of predicted hydrogen binding energies and equilibrium structures for a benchmark set of molecules is compared for some recently developed density functionals, Becke's three parameter hybrid method with the Lee, Yang, and Parr (LYP) correlation functional (B3LYP), Becke's half and half functional combined with the LYP correlation functional (BHLYP), Perdew, Burke and Ernzerhof functional (PBE), Van Voorhis, Scuseria exchange correlation functional (VSXC), the hybrid Perdew, Burke and Ernzerhof functional (PBE1PBE), and meta-generalized gradient approximation (meta-GGA). Overall, the hybrid functionals which contain a portion of Hartree–Fock exchange (B3LYP, BHLYP, and PBE1PBE) yield the most accurate results. The kinetic-energy-density-dependent functionals, VSXC and meta-GGA, are significantly less accurate. Received: 10 December 1999 / Accepted: 5 March 2000 / Published online: 21 June 2000     

Calculations of ionization energies and electron affinities for atoms and molecules: A comparative study with different methods

In the present work, we examined the performance of 36 density functionals, including the newly developed doubly hybrid density functional XYG3 (Y. Zhang, X. Xu, and W. A. Goddard III, Proc. Natl. Acad. Sci, USA, 2009, 106, 4963), to calculate ionization energies (IEs) and electron affinities (EAs). We used the well-established G2-1 set as reference, which contains 14 atoms and 24 molecules for IE, along with 7 atoms and 18 molecules for EA. XYG3 leads to mean absolute deviations (MADs) of 0.057 and 0.080 eV for IEs and EAs, respectively, using the basis set of 6-311 + G (3df,2p). In comparison with some other functionals, MADs for IEs are 0.109 (B2PLYP), 0.119 (M06-2X), 0.159 (X3LYP), 0.161 (PBE), 0.162 (B3LYP), 0.165 (PBE0), 0.173 (TPSS), 0.200 (BLYP), and 0.215 eV (LC-BLYP). MADs for EAs are 0.090 (X3LYP), 0.090 (B2PLYP), 0.102 (PBE), 0.103 (M06-2X), 0.104 (TPSS), 0.105 (BLYP), 0.106 (B3LYP), 0.126 (LC-BLYP), and 0.128 eV (PBE0).     

Systematic studies on the computation of nuclear magnetic resonance shielding constants and chemical shifts: the density functional models

We present a systematic density functional investigation on the prediction of the 13C, 15N, 17O, and 19F NMR properties of 23 molecules with 21 density functionals. Extensive comparisons are made for both 13C magnetic shieldings and chemical shifts with respect to the gas phase experimental data and the best CCSD(T) results. We find that the OPBE and OPW91 exchange-correlation functionals perform significantly better than some popular functionals such as B3LYP and PBE1PBE, even surpassing, in many cases, the standard wavefunction-based method MP2. Further analysis has been performed to explore the individual role played by various exchange and correlation functionals. We find that the B88 and PBE exchange functionals have a too strong tendency of deshielding, leading to too deshielded magnetic shielding constants; whereas the OPTX exchange functional performs remarkably well. We claim that the main source of error arises from the exchange functional, but correlation functional also makes important contribution. We find that the correlation functionals may be grouped into two classes. class A, such as LYP and B98, leads to deshielded NMR values, deteriorating the overall performance; whereas class B, such as PW91 and PBE, generally increases the absolute shieldings, which complements the exchange functionals, leading to improved results in the calculation of NMR data.     

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.