Density-functional calculations for large systems: Can GGA functionals Be competitive with hybrid functionals?

Vincent Tognetti,C. Adamo,P. Cortona

DOI: https://doi.org/10.1007/s12539-010-0073-2

2010-05-01

Interdisciplinary Sciences Computational Life Sciences

Abstract:

What problem does this paper attempt to address?

Thirty years of density functional theory in computational chemistry: an overview and extensive assessment of 200 density functionals

Narbe Mardirossian,Martin Head-Gordon

DOI: https://doi.org/10.1080/00268976.2017.1333644

2017-06-21

Molecular Physics

Abstract:In the past 30 years, Kohn–Sham density functional theory has emerged as the most popular electronic structure method in computational chemistry. To assess the ever-increasing number of approximate exchange-correlation functionals, this review benchmarks a total of 200 density functionals on a molecular database (MGCDB84) of nearly 5000 data points. The database employed, provided as Supplemental Data, is comprised of 84 data-sets and contains non-covalent interactions, isomerisation energies, thermochemistry, and barrier heights. In addition, the evolution of non-empirical and semi-empirical density functional design is reviewed, and guidelines are provided for the proper and effective use of density functionals. The most promising functional considered is ωB97M-V, a range-separated hybrid meta-GGA with VV10 nonlocal correlation, designed using a combinatorial approach. From the local GGAs, B97-D3, revPBE-D3, and BLYP-D3 are recommended, while from the local meta-GGAs, B97M-rV is the leading choice, followed by MS1-D3 and M06-L-D3. The best hybrid GGAs are ωB97X-V, ωB97X-D3, and ωB97X-D, while useful hybrid meta-GGAs (besides ωB97M-V) include ωM05-D, M06-2X-D3, and MN15. Ultimately, today's state-of-the-art functionals are close to achieving the level of accuracy desired for a broad range of chemical applications, and the principal remaining limitations are associated with systems that exhibit significant self-interaction/delocalisation errors and/or strong correlation effects.

chemistry, physical,physics, atomic, molecular & chemical
Towards Very Large-Scale Electronic-Structure Calculations

S. Baroni,P. Giannozzi,S Baroni,P Giannozzi

DOI: https://doi.org/10.1209/0295-5075/17/6/012

1992-02-01

Europhysics Letters (EPL)

Abstract:We present a new approach to density functional theory, which does not require the calculation of Kohn-Sham orbitals. The computational workload required by our method—which is based on the calculation of selected elements of the Green's function—scales linearly with the volume of the system, thus opening the way to first-principles calculations for very large systems. Some of the problems which still hinder the achievement of this goal are discussed, and possible solutions are outlined. As an application, we calculate the charge density of a model silicon supercell containing 64 atoms slightly displaced at random from equilibrium.
A New Generation Density Functional XYG3

Zhang Igor Ying,Xu Xin

2012-01-01

Progress in chemistry

Abstract:There is growing evidence, showing that the widely-used approximate functionals, such as B3LYP, degrade as the system becomes large, underestimate reaction barrier heights and fail to bind van der Waals systems, etc. The success of the Kohn-Sham implementation of density functional theory (DFT) depends on the quality of the exchange-correlation functional. This paper provides an overview of the recent progress on the construction of a new generation functional, namely doubly hybrid density functionals (DHDFs). We pointed out that the theoretical basis of DHDFs lies in Coding-Levy (GL) coupling-constant perturbation theory and adiabatic connection method, and we proposed that the currently available DHDFs can be classified into three groups by their different references used to construct the second order perturbation energy. We systematically examined the performance of various DHDFs. Finally, possible directions for future development of DHDFs are forecasted.
A General and Transferable Local Hybrid Functional for Electronic Structure Theory and Many-Fermion Approaches

Christof Holzer,Yannick J. Franzke

2024-10-01

Abstract:Density functional theory has become the workhorse of quantum physics, chemistry, and materials science. Within these fields, a broad range of applications needs to be covered. These applications range from solids to molecular systems, from organic to inorganic chemistry, or even from electrons to other fermions such as protons or muons. This is emphasized by the plethora of density functional approximations that have been developed for various cases. In this work, a new local hybrid exchange-correlation density functional is constructed from first principles, promoting generality and transferability. We show that constraint satisfaction can be achieved even for admixtures with full exact exchange, without sacrificing accuracy. The performance of the new functional for electronic structure theory is assessed for thermochemical properties, excitation energies, Mössbauer isomer shifts, NMR spin--spin coupling constants, NMR shieldings and shifts, magnetizabilities, as well as EPR hyperfine coupling constants. Here, the new density functional shows excellent performance throughout all tests and is numerically robust only requiring small grids for converged results. Additionally, the functional can be easily generalized to arbitrary fermions as shown for electron-proton correlation energies. Therefore, we outline that density functionals generated in this way are general purpose tools for quantum mechanical studies.

Chemical Physics
Benchmarking Density Functionals on Structural Parameters of Small-/Medium-Sized Organic Molecules.

Eric Bremond,Marika Savarese,Neil Qiang Su,Angel Jose Perez-Jimenez,Xin Xu,Juan Carlos Sancho-Garcia,Carlo Adamo

DOI: https://doi.org/10.1021/acs.jctc.5b01144

2016-01-01

Journal of Chemical Theory and Computation

Abstract:In this Letter we report the error analysis of 59 exchange-correlation functionals in evaluating the structural parameters of small- and medium-sized organic molecules. From this analysis, recently developed double hybrids, such as xDH-PBE0, emerge as the most reliable methods, while global hybrids confirm their robustness in reproducing molecular structures. Notably the M06-L density functional is the only semilocal method reaching an accuracy comparable to hybrids'. A comparison with errors obtained on energetic databases (including thermochemistry, reaction barriers, and interaction energies) indicate that most of the functionals have a coherent behavior, showing low (or high) deviations on both energy and structure data sets. Only a few of them are more prone toward one of these two properties.
Density Functional Theory for Molecule–Metal Surface Reactions: When Does the Generalized Gradient Approximation Get It Right, and What to Do If It Does Not

Nick Gerrits,Egidius W. F. Smeets,Stefan Vuckovic,Andrew D. Powell,Katharina Doblhoff-Dier,Geert-Jan Kroes

DOI: https://doi.org/10.1021/acs.jpclett.0c02452

IF: 6.888

2020-12-09

The Journal of Physical Chemistry Letters

Abstract:While density functional theory (DFT) is perhaps the most used electronic structure theory in chemistry, many of its practical aspects remain poorly understood. For instance, DFT at the generalized gradient approximation (GGA) tends to fail miserably at describing gas-phase reaction barriers, while it performs surprisingly well for many molecule–metal surface reactions. GGA-DFT also fails for many systems in the latter category, and up to now it has not been clear when one may expect it to work. We show that GGA-DFT tends to work if the difference between the work function of the metal and the molecule's electron affinity is greater than ∼7 eV and to fail if this difference is smaller, with sticking of O<sub>2</sub> on Al(111) being a spectacular example. Using dynamics calculations we show that, for this system, the DFT problem may be solved as done for gas-phase reactions, i.e., by resorting to hybrid functionals, but using screening at long-range to obtain a correct description of the metal. Our results suggest the GGA error in the O<sub>2</sub> + Al(111) barrier height to be functional driven. Our results also suggest the possibility to compute potential energy surfaces for the difficult-to-treat systems with computationally cheap nonself-consistent calculations in which a hybrid functional is applied to a GGA density.This article has not yet been cited by other publications.
A Time-Dependent Random State Approach for Large-Scale Density Functional Calculations

Weiqing Zhou,Shengjun Yuan

DOI: https://doi.org/10.1088/0256-307x/40/2/027101

2023-01-01

Chinese Physics Letters

Abstract:We develop a self-consistent first-principle method based on the density functional theory. Physical quantities such as the density of states, Fermi energy and electron density are obtained using a time-dependent random state method without diagonalization. The numerical error for calculating either global or local variables always scales as 1 / S N e , where N e is the number of electrons and S is the number of random states, leading to a sublinear computational cost with the system size. In the limit of large systems, one random state could be enough to achieve reasonable accuracy. The accuracy and scaling properties of using the method are derived analytically and verified numerically in different condensed matter systems. Our time-dependent random state approach provides a powerful strategy for large-scale density functional calculations.
Workhorse semilocal density functional for condensed matter physics and quantum chemistry

John P Perdew,Adrienn Ruzsinszky,Gábor I Csonka,Lucian A Constantin,Jianwei Sun

DOI: https://doi.org/10.1103/PhysRevLett.103.026403

2009-07-10

Abstract:Semilocal density functionals for the exchange-correlation energy are needed for large electronic systems. The Tao-Perdew-Staroverov-Scuseria (TPSS) meta-generalized gradient approximation (meta-GGA) is semilocal and usefully accurate, but predicts too-long lattice constants. Recent "GGA's for solids" yield good lattice constants but poor atomization energies of molecules. We show that the construction principle for one of them (restoring the density gradient expansion for exchange over a wide range of densities) can be used to construct a "revised TPSS" meta-GGA with accurate lattice constants, surface energies, and atomization energies for ordinary matter.
How Well Can Modern Density Functionals Predict Internuclear Distances at Transition States?

Xuefei Xu,I. M. Alecu,Donald G. Truhlar

DOI: https://doi.org/10.1021/ct2001057

2011-01-01

Journal of Chemical Theory and Computation

Abstract:We introduce a new database called TSG48 containing 48 transition state geometrical data (in particular, internuclear distances in transition state structures) for 16 main group reactions. The 16 reactions are the 12 reactions in the previously published DBH24 database (which includes hydrogen transfer reactions, heavy-atom transfer reactions, nucleophilic substitution reactions, and association reactions plus one unimolecular isomerization) plus four H-transfer reactions in which a hydrogen atom is abstracted by the methyl or hydroperoxyl radical from the two different positions in methanol. The data in TSG48 include data for four reactions that have previously been treated at a very high level in the literature. These data are used to test and validate methods that are affordable for the entire test suite, and the most accurate of these methods is found to be the multilevel BMC-CCSD method. The data that constitute the TSG48 database are therefore taken to consist of these very high level calculations for the four reactions where they are available and BMC-CCSD calculations for the other 12 reactions. The TSG48 database is used to assess the performance of the eight Minnesota density functionals from the M05-M08 families and 26 other high-performance and popular density functionals for locating transition state geometries. For comparison, the MP2 and QCISD wave function methods have also been tested for transition state geometries. The MC3BB and MC3MPW doubly hybrid functionals and the M08-HX and M06-2X hybrid meta-GGAs are found to have the best performance of all of the density functionals tested. M08-HX is the most highly recommended functional due to the excellent performance for all five subsets of TSG48, as well as having a lower cost when compared to doubly hybrid functionals. The mean absolute errors in transition state internuclear distances associated with breaking and forming bonds as calculated by the B2PLYP, MP2, and B3LYP methods are respectively about 2, 3, and 5 times larger than those calculated by MC3BB and M08-HX.
Density functional calculation of many-electron systems in cartesian coordinate grid

Amlan K. Roy

DOI: https://doi.org/10.48550/arXiv.1001.2370

2010-01-14

Abstract:A recently developed density functional method, within Hohenberg-Kohn-Sham framework, is used for faithful description of atoms, molecules in Cartesian coordinate grid, by using an LCAO-MO ansatz. Classical Coulomb potential is obtained by means of a Fourier convolution technique. All two-body potentials (including exchange-correlation (XC)) are constructed directly on real grid, while their corresponding matrix elements are computed from numerical integration. Detailed systematic investigation is made for a representative set of atoms/molecules through a number of properties like total energies, component energies, ionization energies, orbital energies, etc. Two nonlocal XC functionals (FT97 and PBE) are considered for pseudopotential calculation of 35 species while preliminary all-electron results are reported for 6 atoms using the LDA XC density functional. Comparison with literature results, wherever possible, exhibits near-complete agreement. This offers a simple efficient route towards accurate reliable calculation of many-electron systems in the Cartesian grid. Future prospect of this method is also discussed.

Chemical Physics
Amyloidosis presenting as urticaria.

D. D. Kennedy,F. Rosenthal,I. Sneddon

DOI: https://doi.org/10.1136/bmj.1.5478.31

1966-01-01

British Medical Journal

Abstract:Attention has recently been drawn to the syndrome of inherited amyloidosis by Gafni et at. (1964), who suggested that some of the isolated cases of idiopathic primary amyloidosis might be sporadic examples of the inherited form. We record below the facts about a patient who suffered from amyloidosis, apparently of the hereditary type. The early clinical features were so unusual and misleading that we hope
A double-hybrid density functional based on good local physics with outstanding performance on the GMTKN55 database

Axel D. Becke,Golokesh Santra,Jan M. L. Martin

DOI: https://doi.org/10.1063/5.0141238

IF: 4.304

2023-04-21

The Journal of Chemical Physics

Abstract:In two recent papers [A. D. Becke, J. Chem. Phys. 156, 214101 (2022) and A. D. Becke, J. Chem. Phys. 157, 234102 (2022)], we compared two Kohn–Sham density functionals based on physical modeling and theory with the best density-functional power-series fits in the literature. The best error statistics reported to date for a hybrid functional on the general main-group thermochemistry, kinetics, and noncovalent interactions (GMTKN55) chemical database of Goerigk et al. [Phys. Chem. Chem. Phys. 19, 32184 (2017)] were obtained. In the present work, additional second-order perturbation-theory terms are considered. The result is a 12-parameter double-hybrid density functional with the lowest GMTKN55 WTMAD2 “weighted total mean absolute deviation” error (1.76 kcal/mol) yet seen for any hybrid or double-hybrid density-functional approximation. We call it “DH23.”

chemistry, physical,physics, atomic, molecular & chemical
How Evenly Can Approximate Density Functionals Treat the Different Multiplicities and Ionization States of 4d Transition Metal Atoms?

Sijie Luo,Donald G Truhlar,Donald G. Truhlar

DOI: https://doi.org/10.1021/ct300737t

2012-09-27

Journal of Chemical Theory and Computation

Abstract:The ability of density functional theory to predict the relative energies of different spin states, especially for systems containing transition metal atoms, is of great importance for many applications. Here, in order to sort out the key factors determining accuracy, we compare the predictions of 60 density functional approximations of 10 different types [local spin density approximation, generalized gradient approximation (GGA), nonseparable gradient approximation (NGA), global-hybrid GGA, range-separated hybrid GGA, range-separated hybrid GGA plus molecular mechanics, meta-GGA, meta-NGA, global-hybrid meta-GGA, and range-separated hybrid meta-GGA] for their ability to represent the spin-flip transitions of all 4d transition metal atoms of groups 3-10 (Y through Pd) and their singly positive cations. We consider all 16 excitation energies connecting the ground states (of the neutral atoms and the cations) to their first excited states of different multiplicities, and we also consider all eight ionization potentials. We also test the Hartree-Fock method. All density functional and Hartree-Fock calculations are converged to a stable solution, in which the spatial symmetry is allowed to be completely broken to achieve the lowest possible energy solution. By analyzing the fractional subshell occupancies and spin contaminations, we are able to sort out the effects of s orbital vs d orbital bias and high-spin vs low-spin bias. A reliable functional should have little or no bias of either type rather than succeeding for a limited subset of cases by cancellation of errors. We find that the widely used correlations of spin splittings to percentage of Hartree-Fock exchange are not borne out by the data, and the correlation functionals also play a significant role. We eventually conclude that SOGGA11-X, B1LYP, B3V5LYP, and MPW3LYP are the most consistently reliable functionals for balanced treatments of 4d transition metal atoms and their cations.

chemistry, physical,physics, atomic, molecular & chemical
Extending Density-Corrected Density Functional Theory to Large Molecular Systems

Youngsam Kim,Mingyu Sim,Minhyeok Lee,Sehun Kim,Suhwan Song,Kieron Burke,Eunji Sim

DOI: https://doi.org/10.26434/chemrxiv-2024-tn0bc

2024-10-04

Abstract:Practical density-corrected density functional theory (DC-DFT) calculations rely on Hartree-Fock (HF) densities, which can be computationally expensive for systems with over a hundred atoms. We extend the applicability of HF-DFT using the dual-basis method, where the density matrix from a smaller basis set is used to estimate the HF solution on a larger basis set. Benchmarks on many systems, including the GMTKN55 database for main-group chemistry, and the L7 and S6L datasets for large molecular systems demonstrate the efficacy of our approach. We apply the dual-basis method to both DNA and HIV systems, and compare with the literature. A careful reparameterisation of HF-r2SCAN-DC4 eliminates the negative s8 coefficient, with no loss of performance.

Chemistry
GW with hybrid functionals for large molecular systems

Tucker Allen,Minh Nguyen,Daniel Neuhauser

2024-05-21

Abstract:A low-cost approach for stochastically sampling static exchange during TDHF-type propagation is presented. This enables the use of an excellent hybrid DFT starting point for stochastic GW quasiparticle energy calculations. Generalized Kohn-Sham molecular orbitals and energies, rather than those of a local-DFT calculation, are used for building the Green's function and effective Coulomb interaction. The use of an optimally tuned hybrid diminishes the starting point dependency in one-shot stochastic GW, effectively avoiding the need for self-consistent GW iterations.

Chemical Physics
Non-empirical Double-Hybrid Density Functionals as Reliable Tools for Electronic Structure Calculations

Juan C Sancho-Garcia,Éric Brémond,Angel J Perez-Jimenez,Ilaria Ciofini,Carlo Adamo

DOI: https://doi.org/10.1088/2516-1075/ac9942

2022-10-13

Electronic Structure

Abstract:The development of universal and accurate approximations for electronic structure calculations lies at the central core of (past and modern) research in theoretical and computational chemistry. For that purpose, any reliable method needs to treat in a balanced way exchange and correlation effects arising from the intrincate structure of matter at the nanoscopic level. Following this principle, we have developed a set of non-empirical (double-hybrid) density functional expressions, minimizing the parameterization and also widely applicable even for systems of considerable size, while being accurate enough to compete with wavefunction methods or even matching experimental information. The underlying expressions are now implemented in many available codes worldwide, then allowing the access to the whole set of key properties needed for addressing chemical structure, reactivity, and bonding, at all nanostructured levels and/or states of matter. Additionally, the recent extension to excited states through a time-dependent (linear-response) formalism also allows one to deal with photochemistry, photophysical, and related properties. Therefore, this family of methods can now be successfully applied to organic, inorganic, or biomolecular compounds, or any other complex system, within an affordable computational effort.

English Else
Efficient treatment of local meta-generalized gradient density functionals via auxiliary density expansion: the density fitting (DF) J+X approximation

Alyssa V. Bienvenu,Gerald Knizia

DOI: https://doi.org/10.1021/acs.jctc.7b01083

2018-02-06

Abstract:We report an efficient technique to treat density functionals of the meta-generalized gradient approximation (mGGA) class in conjunction with density fitting of Coulomb terms (DF-J) and exchange-correlation terms (DF-X). While the kinetic energy density $\tau$ cannot be computed in the context of a DF-JX calculation, we show that the Laplacian of the density $\upsilon$ can be computed with almost no extra cost. With this technique, $\upsilon$-form mGGAs become only slightly more expensive (10%--20%) than GGAs in DF-JX treatment---and several times faster than regular $\tau$-based mGGA calculations with DF-J and regular treatment of the density functional. We investigate the translation of $\upsilon$-form mGGAs into $\tau$-form mGGAs by employing a kinetic energy functional, but find this insufficiently reliable at this moment. However, $\upsilon$ and $\tau$ are believed to carry essentially equivalent information beyond $\rho$ and $\Vert\vec\nabla\rho\Vert$ [Phys. Rev. B 2007, 75, 155109], so a reparametrization of accurate mGGAs from the $\tau$-form into the $\upsilon$-form should be possible. Once such functionals become available, we expect the presented technique to become a powerful tool in the computation of reaction paths, intermediates, and transition states of medium sized molecules.

Chemical Physics
Hybrid density functional calculations of the band gap of Ga$_x$In$_{1-x}$N

Xifan Wu,Eric J. Walter,Andrew M. Rappe,Roberto Car,Annabella Selloni

DOI: https://doi.org/10.48550/arXiv.0907.2001

2009-07-12

Materials Science

Abstract:Recent theoretical work has provided evidence that hybrid functionals, which include a fraction of exact (Hartree Fock) exchange in the density functional theory (DFT) exchange and correlation terms, significantly improve the description of band gaps of semiconductors compared with local and semilocal approximations. Based on a recently developed order-$N$ method for calculating the exact exchange in extended insulating systems, we have implemented an efficient scheme to determine the hybrid functional band gap. We use this scheme to study the band gap and other electronic properties of the ternary compound In$_{1-x}$Ga$_{x}$N using a 64-atom supercell model.
Hybrid MPI and OpenMP Parallel Implementation of Large-Scale Linear-Response Time-Dependent Density Functional Theory with Plane-Wave Basis Set

Lingyun Wan,Xiaofeng Liu,Jie Liu,Xinming Qin,Wei Hu,Jinlong Yang

DOI: https://doi.org/10.1088/2516-1075/abfd1f

2021-01-01

Electronic Structure

Abstract:High performance computing is a powerful tool to accelerate the Kohn-Sham density functional theory calculations on modern heterogeneous supercomputers. Here, we describe a massively parallel implementation of large-scale linear-response time-dependent density functional theory (LR-TDDFT) to calculate the excitation energies and wave functions of solids with plane-wave basis set. We adopt a two-level parallelization strategy that combines the message passing interface with open multi-processing parallel programming to deal with the matrix operations and data communications of constructing and diagonalizing the LR-TDDFT Hamiltonian matrix. Numerical results illustrate that the LR-TDDFT calculations can scale up to 24 576 processing cores on modern heterogeneous supercomputers to study the excited state properties of bulky silicon systems containing thousands of atoms (4,096 atoms). We demonstrate that the LR-TDDFT calculations can be used to investigate the photoinduced charge separation of water molecule adsorption on rutile TiO2(110) surface from an excitonic perspective.
Double-hybrid density-functional theory applied to molecular crystals

Kamal Sharkas,Julien Toulouse,Lorenzo Maschio,Bartolomeo Civalleri

DOI: https://doi.org/10.1063/1.4890439

2014-07-04

Abstract:We test the performance of a number of two- and one-parameter double-hybrid approximations, combining semilocal exchange-correlation density functionals with periodic local second-order Møller-Plesset (LMP2) perturbation theory, for calculating lattice energies of a set of molecular crystals: urea, formamide, ammonia, and carbon dioxide. All double-hybrid methods perform better on average than the corresponding Kohn-Sham calculations with the same functionals, but generally not better than standard LMP2. The one-parameter double-hybrid approximations based on the PBEsol density functional gives lattice energies per molecule with an accuracy of about 6 kJ/mol, which is similar to the accuracy of LMP2. This conclusion is further verified on molecular dimers and on the hydrogen cyanide crystal.

Chemical Physics,Materials Science,Computational Physics

Density-functional calculations for large systems: Can GGA functionals Be competitive with hybrid functionals?

Thirty years of density functional theory in computational chemistry: an overview and extensive assessment of 200 density functionals

Towards Very Large-Scale Electronic-Structure Calculations

A New Generation Density Functional XYG3

A General and Transferable Local Hybrid Functional for Electronic Structure Theory and Many-Fermion Approaches

Benchmarking Density Functionals on Structural Parameters of Small-/Medium-Sized Organic Molecules.

Density Functional Theory for Molecule–Metal Surface Reactions: When Does the Generalized Gradient Approximation Get It Right, and What to Do If It Does Not

A Time-Dependent Random State Approach for Large-Scale Density Functional Calculations

Workhorse semilocal density functional for condensed matter physics and quantum chemistry

How Well Can Modern Density Functionals Predict Internuclear Distances at Transition States?

Density functional calculation of many-electron systems in cartesian coordinate grid

Amyloidosis presenting as urticaria.

A double-hybrid density functional based on good local physics with outstanding performance on the GMTKN55 database

How Evenly Can Approximate Density Functionals Treat the Different Multiplicities and Ionization States of 4d Transition Metal Atoms?

Extending Density-Corrected Density Functional Theory to Large Molecular Systems

GW with hybrid functionals for large molecular systems

Non-empirical Double-Hybrid Density Functionals as Reliable Tools for Electronic Structure Calculations

Efficient treatment of local meta-generalized gradient density functionals via auxiliary density expansion: the density fitting (DF) J+X approximation

Hybrid density functional calculations of the band gap of Ga$_x$In$_{1-x}$N

Hybrid MPI and OpenMP Parallel Implementation of Large-Scale Linear-Response Time-Dependent Density Functional Theory with Plane-Wave Basis Set

Double-hybrid density-functional theory applied to molecular crystals