Patrik Thunström

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.

Electron correlations in MnxGa1-xAs as seen by resonant electron spectroscopy and dynamical mean field theory

After two decades since the discovery of ferromagnetism in manganese-doped gallium arsenide, its origin is still debated, and many doubts are related to the electronic structure. Here we report an experimental and theoretical study of the valence electron spectrum of manganese-doped gallium arsenide. The experimental data are obtained through the differences between off- and on-resonance photo emission data. The theoretical spectrum is calculated by means of a combination of density-functional theory in the local density approximation and dynamical mean field theory, using exact diagonalization as impurity solver. Theory is found to accurately reproduce measured data and illustrates the importance of correlation effects. Our results demonstrate that the manganese states extend over a broad range of energy, including the top of the valence band, and that no impurity band splits-off from the valence band edge, whereas the induced holes seem located primarily around the manganese impurity.

Electronic entanglement in late transition metal oxides.

We present a study of the entanglement in the electronic structure of the late transition metal monoxides--MnO, FeO, CoO, and NiO--obtained by means of density-functional theory in the local density approximation combined with dynamical mean-field theory. The impurity problem is solved through exact diagonalization, which grants full access to the thermally mixed many-body ground state density operator. The quality of the electronic structure is affirmed through a direct comparison between the calculated electronic excitation spectrum and photoemission experiments. Our treatment allows for a quantitative investigation of the entanglement in the electronic structure. Two main sources of entanglement are explicitly resolved through the use of a fidelity based geometrical entanglement measure, and additional information is gained from a complementary entropic entanglement measure. We show that the interplay of crystal field effects and Coulomb interaction causes the entanglement in CoO to take a particularly intricate form.

Charge self-consistent dynamical mean-field theory based on the full-potential linear muffin-tin orbital method: methodology and applications.

Full charge self-consistence (CSC) over the electron density has been implemented into the local density approximation plus dynamical mean-field theory (LDA + DMFT) scheme based on a full-potential ...

Identifying the electronic character and role of the Mn states in the valence band of (Ga,Mn)As.

We report high-resolution hard x-ray photoemission spectroscopy results on (Ga,Mn)As films as a function of Mn doping. Supported by theoretical calculations we identify, for both low (1%) and high (13%) Mn doping values, the electronic character of the states near the top of the valence band. Magnetization and temperature-dependent core-level photoemission spectra reveal how the delocalized character of the Mn states enables the bulk ferromagnetic properties of (Ga,Mn)As.

Chebyshev expansion for impurity models using matrix product states

We improve a recently developed expansion technique for calculating real-frequency spectral functions of any one-dimensional model with short-range interactions, by postprocessing computed Chebyshev moments with linear prediction. This can be achieved at virtually no cost, and in sharp contrast to existing methods based on the dampening of the moments, improves the spectral resolution rather than lowering it. We validate the method for the exactly solvable resonating level model and the single impurity Anderson model. It is capable of resolving sharp Kondo resonances, as well as peaks within the Hubbard bands when employed as an impurity solver for dynamical mean-field theory. Our method works at zero temperature and allows for arbitrary discretization of the bath spectrum. It achieves similar precision as the dynamical density matrix renormalization group, at lower cost. We also propose an alternative expansion, of 1 - exp(-tau H) instead of the usual H, which opens the possibility of using established methods for the time evolution of matrix product states to calculate the spectral functions directly.

Efficient DMFT impurity solver using real-time dynamics with matrix product states

We propose to calculate spectral functions of quantum impurity models using the time evolving block decimation (TEBD) for matrix product states. The resolution of the spectral function is improved by a so-called linear prediction approach. We apply the method as an impurity solver within the dynamical mean-field theory (DMFT) for the single- and two-band Hubbard model on the Bethe lattice. For the single-band model, we observe sharp features at the inner edges of the Hubbard bands. A finite-size scaling shows that they remain present in the thermodynamic limit. We analyze the real time-dependence of the double occupation after adding a single electron and observe oscillations at the same energy as the sharp feature in the Hubbard band, indicating a long-lived coherent superposition of states that correspond to the Kondo peak and the side peaks. For a two-band Hubbard model, we observe an even richer structure in the Hubbard bands, which cannot be related to a multiplet structure of the impurity, in addition to sharp excitations at the band edges of a type similar to the single-band case.

Correlated electronic structure and chemical bonding of cerium pnictides and γ -Ce

We present calculated spectral properties and lattice parameters for cerium pnictides (CeN, CeP, CeAs, CeSb, CeBi) and

Ab initio dynamical vertex approximation

\ensuremath{\gamma}

Treatment of 4 f states of the rare earths: The case study of tbn

-Ce, within the LDA/GGA+DMFT (local density approximation/generalized gradient approximation+dynamical mean field theory) approach. The effective impurity model arising in the DMFT is solved by using the spin-polarized