Density-functional expansion methods: evaluation of LDA, GGA, and meta-GGA functionals and different integral approximations.
Ontology highlight
ABSTRACT:
SUBMITTER: Giese TJ
PROVIDER: S-EPMC3026014 | biostudies-other | 2010 Dec
REPOSITORIES: biostudies-other
Similar Datasets
Project description:This work examines the suitability of meta-GGA functionals for symmetry-adapted perturbation theory (SAPT) calculations. The assessment is based on the term-by-term comparison with the benchmark SAPT variant based on coupled-cluster singles and doubles description of monomers, SAPT(CCSD). Testing systems include molecular complexes ranging from strong to weak and the He dimer. The following nonempirical meta-GGAs are examined: TPSS, revTPSS, MVS, SCAN, and SCAN0 with and without the asymptotic correction (AC) of the exchange-correlation potential. One range-separated meta-GGA functional, LC-PBETPSS, is also included. The AC-corrected pure meta-GGAs (with the exception of MVS) represent a definite progress in SAPT(DFT) compared to pure GGA, such as PBEAC, with their more consistent predictions of energy components. However, none of the meta-GGAs is better than the hybrid GGA approach SAPT(PBE0AC). The SAPT(DFT) electrostatic energy offers the most sensitive probe of the quality of the underlying DFT density. Both SCAN- and TPSS-based electrostatic energies agree with reference to within 5% or better which is an excellent result. We find that SCAN0 can be used in SAPT without the AC correction. The long-range corrected LC-PBETPSS is a reliable performer both for the components and total interaction energies.
| S-EPMC7256086 | biostudies-literature
Project description:We have used ab initio molecular dynamics (AIMD) to characterize water properties using two meta-generalized gradient approximation (meta-GGA) functionals, M06-L-D3 and B97M-rV, and compared their performance against a standard GGA corrected for dispersion, revPBE-D3, at ambient conditions (298 K, and 1 g cm-3 or 1 atm). Simulations of the equilibrium density, radial distribution functions, self-diffusivity, the infrared spectrum, liquid dipole moments, and characterizations of the hydrogen bond network show that all three functionals have overcome the problem of the early AIMD simulations that erroneously found ambient water to be highly structured, but they differ substantially among themselves in agreement with experiment on this range of water properties. We show directly using water cluster data up through the pentamer that revPBE-D3 benefits from a cancellation of its intrinsic functional error by running classical trajectories, whereas the meta-GGA functionals are demonstrably more accurate and would require the simulation of nuclear quantum effects to realize better agreement with all cluster and condensed phase properties.
| S-EPMC6092720 | biostudies-literature
Project description:The homogeneous electron gas (HEG) is a key ingredient in the construction of most exchange-correlation functionals of density-functional theory. Often, the energy of the HEG is parameterized as a function of its spin density nσ, leading to the local density approximation (LDA) for inhomogeneous systems. However, the connection between the electron density and kinetic energy density of the HEG can be used to generalize the LDA by evaluating it on a geometric average nσavg(r) = nσ1-x(r)ñσx(r) of the local spin density nσ(r) and the spin density ñσ(r) of a HEG that has the local kinetic energy density τσ(r) of the inhomogeneous system. This leads to a new family of functionals that we term meta-local density approximations (meta-LDAs), which are still exact for the HEG, which are derived only from properties of the HEG and which form a new rung of Jacob's ladder of density functionals [ AIP Conf. Proc. 2001, 577, 1]. The first functional of this ladder, the local τ approximation (LTA) of Ernzerhof and Scuseria [ J. Chem. Phys. 1999, 111, 911] that corresponds to x = 1 is unfortunately not stable enough to be used in self-consistent field calculations because it leads to divergent potentials, as we show in this work. However, a geometric averaging of the LDA and LTA densities with smaller values of x not only leads to numerical stability of the resulting functional but also yields more accurate exchange energies in atomic calculations than the LDA, the LTA, or the tLDA functional (x = 1/4) of Eich and Hellgren [ J. Chem. Phys. 2014, 141, 224107]. We choose x = 0.50, as it gives the best total energy in self-consistent exchange-only calculations for the argon atom. Atomization energy benchmarks confirm that the choice x = 0.50 also yields improved energetics in combination with correlation functionals in molecules, almost eliminating the well-known overbinding of the LDA and reducing its error by two thirds.
| S-EPMC8023657 | biostudies-literature
Project description:Creatinine concentrations in blood and urine can be used to detect renal insufficiencies and muscle diseases, but current chemical sensors cannot measure creatinine with sufficient selectivity and robustness because they lack a receptor that binds protonated creatinine (creatininium) selectively enough. As a first step toward identifying potential receptors for creatininium, we examine the accuracy of density functional theory (DFT) and wave function theory (WFT) calculations for creatininium cation geometries, evaluated against reference parameters from experiment. We tested twenty-one local and nonlocal density functionals, Hartree-Fock theory, four semiempirical molecular orbital (SEMO) methods of the neglect of differential overlap (NDO) type, and one self-consistent-field nonorthogonal tight-binding method (SCC-DFTB) as implemented in two closely related software packages. DFT and HF calculations were carried out with the MG3S augmented polarized triple-zeta basis set. We find that DFT significantly outperforms SEMO methods for our dataset, and both SCC-DFTB releases we tested (which gave almost identical results) were less accurate than 81% of the density functionals evaluated. The top five functionals for the creatininium structures we examined were MPW1B95, PBEh, mPW1PW, SVWN5, and B97-2, with MMUEs in bond length of 0.0126 Å, 0.0129 Å, 0.0133 Å, 0.0142 Å, and 0.0144 Å respectively, which indicates that all five functionals are similarly suitable for creatininium. The popular B3LYP functional has a MMUE of 0.0178 Å, which ranks it 16(th) overall. B3LYP also performs less favorably than the best local functional (SVWN5), which is less expensive.
| S-EPMC2504755 | biostudies-literature
Project description:We have assessed integration grid errors arising from the use of popular DFT quadrature schemes for a set of 34 organic reaction energies. The focus is primarily on M05-2X and the M06 suite of functionals (M06-L, M06, M06-2X, and M06-HF). M05-2X, M06, and M06-2X outperform popular older DFT functionals for the reaction energies studied, and offer accuracies comparable to results from perturbative hybrid DFT functionals. However, these new functionals are more sensitive to the choice of quadrature grid than previous generations of DFT functionals. Errors in predicted reaction energies arising from the use of the popular SG-1 grid, which is the default in the Q-Chem package, are significant. In particular, M06-HF reaction energies computed with the SG-1 grid exhibit errors ranging from -6.7 to 3.2 kcal mol(-1) relative to results computed with a very fine integration grid. This grid-sensitivity is not a problem for meta-GGA functionals in general, but is instead due to the specific functional forms used in these functionals. The large grid errors are traced to the kinetic energy density enhancement factor utilized in the exchange component of the M05-2X and M06 functionals. This term contains empirically adjusted parameters that are of large magnitude for all of the M06 functionals and for M06-HF in particular. The product of these large constants with modest integration errors for the kinetic energy density results in very large errors in individual contributions to the exchange energy. This gives rise to the troubling large errors exhibited by these functionals for certain integration grids.
| S-EPMC2840268 | biostudies-other
Project description:We have assessed the accuracy of the magnetic properties of a set of 51 density functional approximations, including both recently published and already established functionals. The accuracy assessment considers a series of 27 small molecules and is based on comparing the predicted magnetizabilities to literature reference values calculated using coupled-cluster theory with full singles and doubles and perturbative triples [CCSD(T)] employing large basis sets. The most accurate magnetizabilities, defined as the smallest mean absolute error, are obtained with the BHandHLYP functional. Three of the six studied Berkeley functionals and the three range-separated Florida functionals also yield accurate magnetizabilities. Also, some older functionals like CAM-B3LYP, KT1, BHLYP (BHandH), B3LYP, and PBE0 perform rather well. In contrast, unsatisfactory performance is generally obtained with Minnesota functionals, which are therefore not recommended for calculations of magnetically induced current density susceptibilities and related magnetic properties such as magnetizabilities and nuclear magnetic shieldings. We also demonstrate that magnetizabilities can be calculated by numerical integration of magnetizability density; we have implemented this approach as a new feature in the gauge-including magnetically induced current (GIMIC) method. Magnetizabilities can be calculated from magnetically induced current density susceptibilities within this approach even when analytical approaches for magnetizabilities as the second derivative of the energy have not been implemented. The magnetizability density can also be visualized, providing additional information that is not otherwise easily accessible on the spatial origin of magnetizabilities.
| S-EPMC8023670 | biostudies-literature
Project description:A new self-consistent procedure for calculating the total energy with an orbital-dependent density functional approximation (DFA), the generalized optimized effective potential (GOEP), is developed in the present work. The GOEP is a nonlocal Hermitian potential that delivers the sets of occupied and virtual orbitals and minimizes the total energy. The GOEP optimization leads to the same minimum as does the orbital optimization. The GOEP method is promising as an effective optimization approach for orbital-dependent functionals, as demonstrated for the self-consistent calculations of the random phase approximation (RPA) to the correlation functionals in the particle-hole (ph) and particle-particle (pp) channels. The results show that the accuracy in describing the weakly interacting van der Waals systems is significantly improved in the self-consistent calculations. In particular, the important single excitations contribution in non-self-consistent RPA calculations can be captured self-consistently through the GOEP optimization, leading to orbital renormalization, without using the single excitations in the energy functional.
| S-EPMC6209318 | biostudies-literature
Project description:A systematic comparative study has been performed to better understand DFT+U (density functional theory?+?U) method. We examine the effect of choosing different double counting and exchange-correlation functionals. The calculated energy distribution and the Hund-J dependence of potential profile for representative configurations clearly show the different behaviors of each DFT+U formalism. In particular, adopting spin-dependent exchange-correlation functionals likely leads to undesirable magnetic solution. Our analyses are further highlighted by real material examples ranging from insulating oxides (MnO and NiO) to metallic magnetic systems (SrRuO3 and BaFe2As2). The current work sheds new light on understanding DFT+U and provides a guideline to use the related methods.
| S-EPMC6015075 | biostudies-literature
Project description:The formulation of density-functional expansion methods is extended to treat the second and higher-order terms involving the response density and spin densities with an arbitrary single-center auxiliary basis. The two-center atomic orbital products are represented by the auxiliary functions centered about those two atoms, and the mapping coefficients are determined from a local constrained variational procedure. This two-center variational procedure allows the mapping coefficients to be pretabulated and splined as a function of internuclear separation for efficient look up. The splines of mapping coefficients have a range no longer than that of the overlap integrals, and the auxiliary density appears as a single point-multipole expansion to all nonoverlapping atoms, thus allowing for the trivial implementation of a linear-scaling algorithm. The method is tested using Gaussian multipole expansions, and the effect of angular and radial completeness is explored. Several auxiliary basis sets are parametrized and compared to an auxiliary basis analogous to that used in the self-consistent-charge density-functional tight-binding model, and the method is demonstrated to greatly improve the representation of the density response with respect to a reference expansion model that does not use an auxiliary basis.
| S-EPMC3113327 | biostudies-other
Project description:We analyze and suggest improvements to a recently developed approximate continuum-electrostatic model for proteins. The model, called BIBEE/I (boundary-integral based electrostatics estimation with interpolation), was able to estimate electrostatic solvation free energies to within a mean unsigned error of 4% on a test set of more than 600 proteins-a significant improvement over previous BIBEE models. In this work, we tested the BIBEE/I model for its capability to predict residue-by-residue interactions in protein-protein binding, using the widely studied model system of trypsin and bovine pancreatic trypsin inhibitor (BPTI). Finding that the BIBEE/I model performs surprisingly less well in this task than simpler BIBEE models, we seek to explain this behavior in terms of the models' differing spectral approximations of the exact boundary-integral operator. Calculations of analytically solvable systems (spheres and tri-axial ellipsoids) suggest two possibilities for improvement. The first is a modified BIBEE/I approach that captures the asymptotic eigenvalue limit correctly, and the second involves the dipole and quadrupole modes for ellipsoidal approximations of protein geometries. Our analysis suggests that fast, rigorous approximate models derived from reduced-basis approximation of boundary-integral equations might reach unprecedented accuracy, if the dipole and quadrupole modes can be captured quickly for general shapes.
| S-EPMC3897805 | biostudies-literature