Georg Kresse

Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set

We present a detailed description and comparison of algorithms for performing ab-initio quantum-mechanical calculations using pseudopotentials and a plane-wave basis set. We will discuss: (a) partial occupancies within the framework of the linear tetrahedron method and the finite temperature density-functional theory, (b) iterative methods for the diagonalization of the Kohn-Sham Hamiltonian and a discussion of an efficient iterative method based on the ideas of Pulays residual minimization, which is close to an order Natoms2 scaling even for relatively large systems, (c) efficient Broyden-like and Pulay-like mixing methods for the charge density including a new special ‘preconditioning’ optimized for a plane-wave basis set, (d) conjugate gradient methods for minimizing the electronic free energy with respect to all degrees of freedom simultaneously. We have implemented these algorithms within a powerful package called VAMP (Vienna ab-initio molecular-dynamics package). The program and the techniques have been used successfully for a large number of different systems (liquid and amorphous semiconductors, liquid simple and transition metals, metallic and semi-conducting surfaces, phonons in simple metals, transition metals and semiconductors) and turned out to be very reliable.

Norm-conserving and ultrasoft pseudopotentials for first-row and transition elements

The construction of accurate pseudopotentials with good convergence properties for the first-row and transition elements is discussed. We show that by combining an improved description of the pseudowavefunction inside the cut-off radius with the concept of ultrasoft pseudopotentials introduced by Vanderbilt optimal compromise between transferability and plane-wave convergence can be achieved. With the new pseudopotentials, basis sets with no more than 75-100 plane waves per atom are sufficient to reproduce the results obtained with the most accurate norm-conserving pseudopotentials.

Screened hybrid density functionals applied to solids

Hybrid Fock exchange/density functional theory functionals have shown to be very successful in describing a wide range of molecular properties. For periodic systems, however, the long-range nature of the Fock exchange interaction and the resultant large computational requirements present a major drawback. This is especially true for metallic systems, which require a dense Brillouin zone sampling. Recently, a new hybrid functional [HSE03, J. Heyd, G. E. Scuseria, and M. Ernzerhof, J. Chem. Phys. 118, 8207 (2003)] that addresses this problem within the context of methods that evaluate the Fock exchange in real space was introduced. We discuss the advantages the HSE03 functional brings to methods that rely on a reciprocal space description of the Fock exchange interaction, e.g., all methods that use plane wave basis sets. Furthermore, we present a detailed comparison of the performance of the HSE03 and PBE0 functionals for a set of archetypical solid state systems by calculating lattice parameters, bulk moduli, heats of formation, and band gaps. The results indicate that the hybrid functionals indeed often improve the description of these properties, but in several cases the results are not yet on par with standard gradient corrected functionals. This concerns in particular metallic systems for which the bandwidth and exchange splitting are seriously overestimated.

The Perdew–Burke–Ernzerhof exchange-correlation functional applied to the G2-1 test set using a plane-wave basis set

Present local and semilocal functionals show significant errors, for instance, in the energetics of small molecules and in the description of band gaps. One possible solution to these problems is the introduction of exact exchange and hybrid functionals. A plane-wave-based algorithm was implemented in VASP (Vienna ab-initio simulation package) to allow for the calculation of the exact exchange. To systematically assess the precision of the present implementation, calculations for the 55 molecules of the G2-1 quantum chemical test set were performed applying the PBE and PBE0 functionals. Excellent agreement for both atomization energies and geometries compared with the results obtained by GAUSSIAN 03 calculations using large basis sets (augmented correlation consistent polarized valence quadruple zeta for the geometry optimization and augmented correlation-consistent polarized valence quintuple zeta for the energy calculations) was found. The mean absolute error for atomization energies between VASP and the experiment is 8.6 and 3.7 kcalmol, as calculated with the PBE and PBE0 functionals, respectively. The mean deviations between VASP and GAUSSIAN are 0.46 and 0.49 kcalmol for the PBE and PBE0 functionals, respectively.

Ab initio Force Constant Approach to Phonon Dispersion Relations of Diamond and Graphite

The phonon dispersion relations of diamond and graphite are calculated using an ab initio force constant method. The force constants are calculated via a self-consistent supercell approach in the local-density approximation in terms of the Hellmann-Feynman forces induced by the displacement of a single atom in the supercell. For diamond our ab initio results are in very good agreement with the neutron inelastic scattering and Raman data. For graphite we find good agreement with the neutron data for the low-energy modes as well as with the reflections electron energy loss spectroscopy (REELS) and optical data at higher energies. Significant differences to the predictions of semi-empirical models appear.

Hybrid functionals applied to extended systems

We present an overview of the description of structural, thermochemical, and electronic properties of extended systems using several well known hybrid Hartree-Fock/density-functional-theory functionals (PBE0, HSE03, and B3LYP). In addition we address a few aspects of the evaluation of the Hartree-Fock exchange interactions in reciprocal space, relevant to all methods that employ a plane wave basis set and periodic boundary conditions.

Reproducibility in density functional theory calculations of solids

A comparison of DFT methods Density functional theory (DFT) is now routinely used for simulating material properties. Many software packages are available, which makes it challenging to know which are the best to use for a specific calculation. Lejaeghere et al. compared the calculated values for the equation of states for 71 elemental crystals from 15 different widely used DFT codes employing 40 different potentials (see the Perspective by Skylaris). Although there were variations in the calculated values, most recent codes and methods converged toward a single value, with errors comparable to those of experiment. Science, this issue p. 10.1126/science.aad3000; see also p. 1394 A survey of recent density functional theory methods shows a convergence to more accurate property calculations. [Also see Perspective by Skylaris] INTRODUCTION The reproducibility of results is one of the underlying principles of science. An observation can only be accepted by the scientific community when it can be confirmed by independent studies. However, reproducibility does not come easily. Recent works have painfully exposed cases where previous conclusions were not upheld. The scrutiny of the scientific community has also turned to research involving computer programs, finding that reproducibility depends more strongly on implementation than commonly thought. These problems are especially relevant for property predictions of crystals and molecules, which hinge on precise computer implementations of the governing equation of quantum physics. RATIONALE This work focuses on density functional theory (DFT), a particularly popular quantum method for both academic and industrial applications. More than 15,000 DFT papers are published each year, and DFT is now increasingly used in an automated fashion to build large databases or apply multiscale techniques with limited human supervision. Therefore, the reproducibility of DFT results underlies the scientific credibility of a substantial fraction of current work in the natural and engineering sciences. A plethora of DFT computer codes are available, many of them differing considerably in their details of implementation, and each yielding a certain “precision” relative to other codes. How is one to decide for more than a few simple cases which code predicts the correct result, and which does not? We devised a procedure to assess the precision of DFT methods and used this to demonstrate reproducibility among many of the most widely used DFT codes. The essential part of this assessment is a pairwise comparison of a wide range of methods with respect to their predictions of the equations of state of the elemental crystals. This effort required the combined expertise of a large group of code developers and expert users. RESULTS We calculated equation-of-state data for four classes of DFT implementations, totaling 40 methods. Most codes agree very well, with pairwise differences that are comparable to those between different high-precision experiments. Even in the case of pseudization approaches, which largely depend on the atomic potentials used, a similar precision can be obtained as when using the full potential. The remaining deviations are due to subtle effects, such as specific numerical implementations or the treatment of relativistic terms. CONCLUSION Our work demonstrates that the precision of DFT implementations can be determined, even in the absence of one absolute reference code. Although this was not the case 5 to 10 years ago, most of the commonly used codes and methods are now found to predict essentially identical results. The established precision of DFT codes not only ensures the reproducibility of DFT predictions but also puts several past and future developments on a firmer footing. Any newly developed methodology can now be tested against the benchmark to verify whether it reaches the same level of precision. New DFT applications can be shown to have used a sufficiently precise method. Moreover, high-precision DFT calculations are essential for developing improvements to DFT methodology, such as new density functionals, which may further increase the predictive power of the simulations. Recent DFT methods yield reproducible results. Whereas older DFT implementations predict different values (red darts), codes have now evolved to mutual agreement (green darts). The scoreboard illustrates the good pairwise agreement of four classes of DFT implementations (horizontal direction) with all-electron results (vertical direction). Each number reflects the average difference between the equations of state for a given pair of methods, with the green-to-red color scheme showing the range from the best to the poorest agreement. The widespread popularity of density functional theory has given rise to an extensive range of dedicated codes for predicting molecular and crystalline properties. However, each code implements the formalism in a different way, raising questions about the reproducibility of such predictions. We report the results of a community-wide effort that compared 15 solid-state codes, using 40 different potentials or basis set types, to assess the quality of the Perdew-Burke-Ernzerhof equations of state for 71 elemental crystals. We conclude that predictions from recent codes and pseudopotentials agree very well, with pairwise differences that are comparable to those between different high-precision experiments. Older methods, however, have less precise agreement. Our benchmark provides a framework for users and developers to document the precision of new applications and methodological improvements.

Quasiparticle band structure based on a generalized Kohn-Sham scheme

We present a comparative full-potential study of generalized Kohn-Sham (gKS) schemes with explicit focus on their suitability as starting point for the solution of the quasiparticle equation. We compare

AB initio molecular dynamics for liquid metals

{G}_{0}{W}_{0}

Towards an exact description of electronic wavefunctions in real solids

quasiparticle band structures calculated upon local-density approximation (LDA), screened-exchange, HSE03, PBE0, and Hartree-Fock functionals for exchange and correlation (XC) for Si, InN, and ZnO. Furthermore, the HSE03 functional is studied and compared to the generalized gradient approximation (GGA) for 15 nonmetallic materials for its use as a starting point in the calculation of quasiparticle excitation energies. For this case, the effects of self-consistency in the