The method of dispersion correction as an add-on to standard Kohn-Sham density functional theory (DFT-D) has been refined regarding higher accuracy, broader range of applicability, and less empiricism. The main new ingredients are atom-pairwise specific dispersion coefficients and cutoff radii that are both computed from first principles. The coefficients for new eighth-order dispersion terms are computed using established recursion relations. System (geometry) dependent information is used for the first time in a DFT-D type approach by employing the new concept of fractional coordination numbers (CN). They are used to interpolate between dispersion coefficients of atoms in different chemical environments. The method only requires adjustment of two global parameters for each density functional, is asymptotically exact for a gas of weakly interacting neutral atoms, and easily allows the computation of atomic forces. Three-body nonadditivity terms are considered. The method has been assessed on standard benchmark sets for inter- and intramolecular noncovalent interactions with a particular emphasis on a consistent description of light and heavy element systems. The mean absolute deviations for the S22 benchmark set of noncovalent interactions for 11 standard density functionals decrease by 15%-40% compared to the previous (already accurate) DFT-D version. Spectacular improvements are found for a tripeptide-folding model and all tested metallic systems. The rectification of the long-range behavior and the use of more accurate C(6) coefficients also lead to a much better description of large (infinite) systems as shown for graphene sheets and the adsorption of benzene on an Ag(111) surface. For graphene it is found that the inclusion of three-body terms substantially (by about 10%) weakens the interlayer binding. We propose the revised DFT-D method as a general tool for the computation of the dispersion energy in molecules and solids of any kind with DFT and related (low-cost) electronic structure methods for large systems.

http://dns2.asia.edu.tw/~ysho/YSHO-English/1000%20WC/PDF/J%20Che%20Phy132,%20154104.pdf

A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu

We present two new hybrid meta exchange- correlation functionals, called M06 and M06-2X. The M06 functional is parametrized including both transition metals and nonmetals, whereas the M06-2X functional is a high-nonlocality functional with double the amount of nonlocal exchange (2X), and it is parametrized only for nonmetals.The functionals, along with the previously published M06-L local functional and the M06-HF full-Hartree–Fock functionals, constitute the M06 suite of complementary functionals. We assess these four functionals by comparing their performance to that of 12 other functionals and Hartree–Fock theory for 403 energetic data in 29 diverse databases, including ten databases for thermochemistry, four databases for kinetics, eight databases for noncovalent interactions, three databases for transition metal bonding, one database for metal atom excitation energies, and three databases for molecular excitation energies. We also illustrate the performance of these 17 methods for three databases containing 40 bond lengths and for databases containing 38 vibrational frequencies and 15 vibrational zero point energies. We recommend the M06-2X functional for applications involving main-group thermochemistry, kinetics, noncovalent interactions, and electronic excitation energies to valence and Rydberg states. We recommend the M06 functional for application in organometallic and inorganometallic chemistry and for noncovalent interactions.

/pdf/the-m06-suite-of-density-functionals-for-main-group-1748swzp1l.pdf

The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals

Deposits of clastic carbonate-dominated (calciclastic) sedimentary slope systems in the rock record have been identified mostly as linearly-consistent carbonate apron deposits, even though most ancient clastic carbonate slope deposits fit the submarine fan systems better. Calciclastic submarine fans are consequently rarely described and are poorly understood. Subsequently, very little is known especially in mud-dominated calciclastic submarine fan systems. Presented in this study are a sedimentological core and petrographic characterisation of samples from eleven boreholes from the Lower Carboniferous of Bowland Basin (Northwest England) that reveals a >250 m thick calciturbidite complex deposited in a calciclastic submarine fan setting. Seven facies are recognised from core and thin section characterisation and are grouped into three carbonate turbidite sequences. They include: 1) Calciturbidites, comprising mostly of highto low-density, wavy-laminated bioclast-rich facies; 2) low-density densite mudstones which are characterised by planar laminated and unlaminated muddominated facies; and 3) Calcidebrites which are muddy or hyper-concentrated debrisflow deposits occurring as poorly-sorted, chaotic, mud-supported floatstones. These

Phd by thesis

We report re-optimization of a recently proposed long-range corrected (LC) hybrid density functional [J.-D. Chai and M. Head-Gordon, J. Chem. Phys., 2008, 128, 084106] to include empirical atom–atom dispersion corrections. The resulting functional, ωB97X-D yields satisfactory accuracy for thermochemistry, kinetics, and non-covalent interactions. Tests show that for non-covalent systems, ωB97X-D shows slight improvement over other empirical dispersion-corrected density functionals, while for covalent systems and kinetics it performs noticeably better. Relative to our previous functionals, such as ωB97X, the new functional is significantly superior for non-bonded interactions, and very similar in performance for bonded interactions.

/pdf/long-range-corrected-hybrid-density-functionals-with-damped-44ti6za9f3.pdf

Long-range corrected hybrid density functionals with damped atom–atom dispersion corrections

Long-Range Corrected Hybrid Density Functionals with Damped Atom-Atom Dispersion Corrections Jeng-Da Chai ∗ and Martin Head-Gordon † Department of Chemistry, University of California and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA (Dated: June 14, 2008) We report re-optimization of a recently proposed long-range corrected (LC) hybrid density func- tionals [J.-D. Chai and M. Head-Gordon, J. Chem. Phys. 128, 084106 (2008)] to include empirical atom-atom dispersion corrections. The resulting functional, ωB97X-D yields satisfactory accuracy for thermochemistry, kinetics, and non-covalent interactions. Tests show that for non-covalent sys- tems, ωB97X-D shows slight improvement over other empirical dispersion-corrected density func- tionals, while for covalent systems and kinetics, it performs noticeably better. Relative to our previous functionals, such as ωB97X, the new functional is signiﬁcantly superior for non-bonded interactions, and very similar in performance for bonded interactions. I. INTRODUCTION Due to its favorable cost-to-performance ratio, Kohn- Sham density-functional theory (KS-DFT) [1, 2] has be- come the most popular electronic structure theory for large-scale ground-state systems [3–5]. Its extension for treating excited-state systems [6, 7], time-dependent den- sity functional theory (TDDFT), has also been developed to the stage where it is now very widely used. The essential ingredient of KS-DFT, the exchange- correlation energy functional E xc , remains unknown and needs to be approximated. Semi-local gradient-corrected density functionals, though successful in many applica- tions, lead to qualitative failures in some circumstances, where the accurate treatment of non-locality of exchange- correlation hole becomes crucial. These situations occur mostly in the asymptotic regions of molecular systems, such as spurious self-interaction eﬀects upon dissociation [8, 9] and dramatic failures for long-range charge-transfer excitations [10–12]. Widely used hybrid density function- als, like B3LYP [13, 14], do not qualitatively resolve these problems. These self-interaction errors can be qualitatively re- solved using the long-range corrected (LC) hybrid density functionals [15, 16, 18], which employ 100% Hartree-Fock (HF) exchange for long-range electron-electron interac- tions. This is accomplished by a partition of unity, using erf(ωr)/r for long-range (treated by HF exchange) and erfc(ωr)/r for short-range (treated by an exchange func- tional), with the parameter ω controlling the partition- ing. Over the past ﬁve years, the LC hybrid scheme has been attracting increasing attention [15] since its compu- tational cost is comparable with standard hybrid func- tionals [13]. However, LC functionals have tended to be inferior to the best hybrids for properties such as ther- mochemistry. ∗ Electronic † Author address: jdchai@berkeley.edu to whom correspondence should be addressed. Electronic address: mhg@cchem.berkeley.edu Recently we have improved the overall accuracy at- tainable with the LC functionals by using a systematic optimization procedure [18]. One important conclusion is that optimizing LC and hybrid functionals with identical numbers of parameters in their GGA exchange and cor- relation terms leads to noticeably better results for all properties using the LC form. The resulting LC func- tional is called ωB97. Further statistically signiﬁcant improvement results from re-optimizing the entire func- tional with one extra parameter corresponding to an ad- justable fraction of short-range exact exchange, deﬁning the ωB97X functional. Independent test sets covering thermochemistry and non-covalent interactions support these conclusions. However, problems associated with the lack of non-locality of the correlation hole, such as the lack of dispersion interactions (London forces), still remain, as the semi-local correlation functionals cannot capture long-range correlation eﬀects [19, 20]. There have been signiﬁcant eﬀorts to develop a frame- work that can account for long-range dispersion eﬀects within DFT. Zaremba and Kohn (ZK) [21] derived an exact expression for the second-order dispersion energy in terms of the exact density-density response functions of the two separate systems. To obtain a tractable non- local dispersion functional, Dobson and Dinite (DD) [22] made local density approximations to the ZK response functions. DD’s non-local correlation functional was ob- tained independently [23] by modifying the eﬀective den- sity deﬁned in the earlier work of Rapcewicz and Ashcroft Starting from the formally exact expression of KS- DFT, the adiabatic connection ﬂuctuation-dissipation theorem (ACFDT), for the ground-state exchange- correlation energy, Langreth and co-workers [25] devel- oped a so-called van der Waals density functional (vdW- DF) by making a series of reasonable approximations to yield a computationally tractable scheme. Recently, Becke and Johnson (BJ) developed a series of post-HF correlation models with a novel treatment for dispersion interactions based on the exchange-hole dipole moment [26]. The origin of dispersion claimed in the BJ models was recently questioned by Alonso, and A.

Long-Range Corrected Hybrid Density Functionals with Damped Atom-Atom Dispersion Corrections

Popular modern generalized gradient approximations are biased toward the description of free-atom energies. Restoration of the first-principles gradient expansion for exchange over a wide range of density gradients eliminates this bias. We introduce a revised Perdew-Burke-Ernzerhof generalized gradient approximation that improves equilibrium properties of densely packed solids and their surfaces.

/pdf/restoring-the-density-gradient-expansion-for-exchange-in-4uz2ddboii.pdf

Restoring the Density-Gradient Expansion for Exchange in Solids and Surfaces

The ground-state energy, electron density, and related properties of ordinary matter can be computed efficiently when the exchange-correlation energy as a functional of the density is approximated semilocally. We propose the first meta-generalized-gradient approximation (meta-GGA) that is fully constrained, obeying all 17 known exact constraints that a meta-GGA can. It is also exact or nearly exact for a set of "appropriate norms," including rare-gas atoms and nonbonded interactions. This strongly constrained and appropriately normed meta-GGA achieves remarkable accuracy for systems where the exact exchange-correlation hole is localized near its electron, and especially for lattice constants and weak interactions.

/pdf/strongly-constrained-and-appropriately-normed-semilocal-1f2s80n9n2.pdf

Strongly Constrained and Appropriately Normed Semilocal Density Functional.

The ground-state energy, electron density, and related properties of ordinary matter can be computed efficiently when the exchange-correlation energy as a functional of the density is approximated semilocally. We propose the first meta-GGA (meta-generalized gradient approximation) that is fully constrained, obeying all 17 known exact constraints that a meta-GGA can. It is also exact or nearly exact for a set of appropriate norms, including rare-gas atoms and nonbonded interactions. This SCAN (strongly constrained and appropriately normed) meta-GGA achieves remarkable accuracy for systems where the exact exchange-correlation hole is localized near its electron, and especially for lattice constants and weak interactions.

Strongly Constrained and Appropriately Normed Semilocal Density Functional

We present the case for the nonempirical construction of density functional approximations for the exchange-correlation energy by the traditional method of “constraint satisfaction” without fitting to data sets, and present evidence that this approach has been successful on the first three rungs of “Jacob’s ladder” of density functional approximations [local spin-density approximation (LSD), generalized gradient approximation (GGA), and meta-GGA]. We expect that this approach will also prove successful on the fourth and fifth rungs (hyper-GGA or hybrid and generalized random-phase approximation). In particular, we argue for the theoretical and practical importance of recovering the correct uniform density limit, which many semiempirical functionals fail to do. Among the beyond-LSD functionals now available to users, we recommend the nonempirical Perdew–Burke–Ernzerhof (PBE) GGA and the nonempirical Tao–Perdew–Staroverov–Scuseria (TPSS) meta-GGA, and their one-parameter hybrids with exact exchange. TPSS improvement over PBE is dramatic for atomization energies of molecules and surface energies of solids, and small or moderate for other properties. TPSS is now or soon will be available in standard codes such as GAUSSIAN, TURBOMOLE, NWCHEM, ADF, WIEN, VASP, etc. We also discuss old and new ideas to eliminate the self-interaction error that plagues the functionals on the first three rungs of the ladder, bring up other related issues, and close with a list of “do’s and don’t’s” for software developers and users.

/pdf/prescription-for-the-design-and-selection-of-density-535wtfz6b0.pdf

Prescription for the design and selection of density functional approximations: more constraint satisfaction with fewer fits.

One atom or molecule binds to another through various types of bond, the strengths of which range from several meV to several eV. Although some computational methods can provide accurate descriptions of all bond types, those methods are not efficient enough for many studies (for example, large systems, ab initio molecular dynamics and high-throughput searches for functional materials). Here, we show that the recently developed non-empirical strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation (meta-GGA) within the density functional theory framework predicts accurate geometries and energies of diversely bonded molecules and materials (including covalent, metallic, ionic, hydrogen and van der Waals bonds). This represents a significant improvement at comparable efficiency over its predecessors, the GGAs that currently dominate materials computation. Often, SCAN matches or improves on the accuracy of a computationally expensive hybrid functional, at almost-GGA cost. SCAN is therefore expected to have a broad impact on chemistry and materials science.

Adrienn Ruzsinszky

Papers

Restoring the Density-Gradient Expansion for Exchange in Solids and Surfaces

Strongly Constrained and Appropriately Normed Semilocal Density Functional.

Strongly Constrained and Appropriately Normed Semilocal Density Functional

Prescription for the design and selection of density functional approximations: more constraint satisfaction with fewer fits.

Accurate first-principles structures and energies of diversely bonded systems from an efficient density functional