-
How to verify the precision of density-functional-theory implementations via reproducible and universal workflows
Authors:
Emanuele Bosoni,
Louis Beal,
Marnik Bercx,
Peter Blaha,
Stefan Blügel,
Jens Bröder,
Martin Callsen,
Stefaan Cottenier,
Augustin Degomme,
Vladimir Dikan,
Kristjan Eimre,
Espen Flage-Larsen,
Marco Fornari,
Alberto Garcia,
Luigi Genovese,
Matteo Giantomassi,
Sebastiaan P. Huber,
Henning Janssen,
Georg Kastlunger,
Matthias Krack,
Georg Kresse,
Thomas D. Kühne,
Kurt Lejaeghere,
Georg K. H. Madsen,
Martijn Marsman
, et al. (20 additional authors not shown)
Abstract:
In the past decades many density-functional theory methods and codes adopting periodic boundary conditions have been developed and are now extensively used in condensed matter physics and materials science research. Only in 2016, however, their precision (i.e., to which extent properties computed with different codes agree among each other) was systematically assessed on elemental crystals: a firs…
▽ More
In the past decades many density-functional theory methods and codes adopting periodic boundary conditions have been developed and are now extensively used in condensed matter physics and materials science research. Only in 2016, however, their precision (i.e., to which extent properties computed with different codes agree among each other) was systematically assessed on elemental crystals: a first crucial step to evaluate the reliability of such computations. We discuss here general recommendations for verification studies aiming at further testing precision and transferability of density-functional-theory computational approaches and codes. We illustrate such recommendations using a greatly expanded protocol covering the whole periodic table from Z=1 to 96 and characterizing 10 prototypical cubic compounds for each element: 4 unaries and 6 oxides, spanning a wide range of coordination numbers and oxidation states. The primary outcome is a reference dataset of 960 equations of state cross-checked between two all-electron codes, then used to verify and improve nine pseudopotential-based approaches. Such effort is facilitated by deploying AiiDA common workflows that perform automatic input parameter selection, provide identical input/output interfaces across codes, and ensure full reproducibility. Finally, we discuss the extent to which the current results for total energies can be reused for different goals (e.g., obtaining formation energies).
△ Less
Submitted 26 May, 2023;
originally announced May 2023.
-
Topological spiral magnetism in the Weyl semimetal SmAlSi
Authors:
Xiaohan Yao Jonathan Gaudet,
Rahul Verma,
David E. Graf,
Hung-Yu Yang,
Faranak Bahrami,
Ruiqi Zhang,
Adam A. Aczel,
Sujan Subedi,
Darius H. Torchinsky,
Jianwei Sun,
Arun Bansil,
Shin-Ming Huang,
Bahadur Singh,
Predrag Nikolic,
Peter Blaha,
Fazel Tafti
Abstract:
Weyl electrons are intensely studied due to novel charge transport phenomena such as chiral anomaly, Fermi arcs, and photogalvanic effect. Recent theoretical works suggest that Weyl electrons can also participate in magnetic interactions, and the Weyl-mediated indirect exchange coupling between local moments is proposed as a new mechanism of spiral magnetism that involves chiral electrons. Despite…
▽ More
Weyl electrons are intensely studied due to novel charge transport phenomena such as chiral anomaly, Fermi arcs, and photogalvanic effect. Recent theoretical works suggest that Weyl electrons can also participate in magnetic interactions, and the Weyl-mediated indirect exchange coupling between local moments is proposed as a new mechanism of spiral magnetism that involves chiral electrons. Despite reports of incommensurate and non-collinear magnetic ordering in Weyl semimetals, an actual spiral order has remained hitherto undetected. Here, we present evidence of Weyl-mediated spiral magnetism in SmAlSi from neutron diffraction, transport, and thermodynamic data. We show that the spiral order in SmAlSi results from the nesting between topologically non-trivial Fermi pockets and weak magnetocrystalline anisotropy, unlike related materials (Ce,Pr,Nd)AlSi, where a strong anisotropy prevents the spins from freely rotating. We map the magnetic phase diagram of SmAlSi and reveal an A-phase where topological magnetic excitations may exist. This is corroborated by the observation of a topological Hall effect within the A-phase.
△ Less
Submitted 10 June, 2022;
originally announced June 2022.
-
Length-gauge optical matrix elements in WIEN2k
Authors:
Oleg Rubel,
Peter Blaha
Abstract:
Hybrid exchange-correlation functionals provide superior electronic structure and optical properties of semiconductors or insulators as compared to semilocal exchange-correlation potentials due to admixing a portion of the non-local exact exchange potential from a Hartree-Fock theory. Since the non-local potential does not commute with the position operator, the momentum matrix elements do not ful…
▽ More
Hybrid exchange-correlation functionals provide superior electronic structure and optical properties of semiconductors or insulators as compared to semilocal exchange-correlation potentials due to admixing a portion of the non-local exact exchange potential from a Hartree-Fock theory. Since the non-local potential does not commute with the position operator, the momentum matrix elements do not fully capture the oscillator strength, while the length-gauge velocity matrix elements do. So far, length-gauge velocity matrix elements were not accessible in the all-electron full-potential WIEN2k package. We demonstrate the feasibility of computing length-gauge matrix elements in WIEN2k for a hybrid exchange-correlation functional based on a finite difference approach. To illustrate the implementation we determined matrix elements for optical transitions between the conduction and valence bands in GaAs, GaN, (CH$_3$NH$_3$)PbI$_3$ and a monolayer MoS$_2$. The non-locality of the Hartree-Fock exact exchange potential leads to a strong enhancement of the oscillator strength as noticed recently in calculations employing pseudopotentials [Laurien and Rubel: arXiv:2111.14772 (2021)]. We obtained an analytical expression for the enhancement factor in terms of the difference in eigenvalues not captured by the kinetic energy. It is expected that these results can also be extended to other non-local potentials, e.g., a many-body $GW$ approximation.
△ Less
Submitted 29 January, 2022; v1 submitted 30 December, 2021;
originally announced December 2021.
-
Implementation of self-consistent MGGA functionals in augmented plane wave based methods
Authors:
Jan Doumont,
Fabien Tran,
Peter Blaha
Abstract:
Functionals of the meta-generalized gradient approximation (MGGA) are nowadays widely used in chemistry and solid-state physics for the simulation of electronic systems like molecules, solids, or surfaces. Due to their dependency on the kinetic energy density, they are in principle more accurate than GGA functionals for various properties (geometry, binding energy, electronic structure, etc.), whi…
▽ More
Functionals of the meta-generalized gradient approximation (MGGA) are nowadays widely used in chemistry and solid-state physics for the simulation of electronic systems like molecules, solids, or surfaces. Due to their dependency on the kinetic energy density, they are in principle more accurate than GGA functionals for various properties (geometry, binding energy, electronic structure, etc.), while being nearly as fast since they are still of the semilocal form. Thus, when an accuracy better than GGA is required, one may consider using a MGGA instead of the much more costly hybrid functionals or methods like the random-phase approximation or $GW$. In this work, the self-consistent implementation of MGGA functionals in APW based methods is presented. Technical aspects of the implementation are discussed, and calculations of band gaps, lattice constants, and magnetic moments are presented in order to validate our implementation. To test the changes of the electron density due to a MGGA, the electric field gradient on transition-metal atoms is calculated.
△ Less
Submitted 27 May, 2022; v1 submitted 24 December, 2021;
originally announced December 2021.
-
Correct and accurate polymorphic energy ordering of transition-metal monoxides obtained from semilocal and onsite-hybrid exchange-correlation approximations
Authors:
Arghya Ghosh,
Subrata Jana,
Manish K Niranjan,
Fabien Tran,
David Wimberger,
Peter Blaha,
Lucian A. Constantin,
Prasanjit Samal
Abstract:
The relative energetic stability of the structural phases of common antiferromagnetic transition-metal oxides (MnO, FeO, CoO, and NiO) within the semilocal and hybrid density functionals are fraught with difficulties. In particular, MnO is known to be the most difficult case for almost all common semilocal and hybrid density approximations. Here, we show that the meta-generalized gradient approxim…
▽ More
The relative energetic stability of the structural phases of common antiferromagnetic transition-metal oxides (MnO, FeO, CoO, and NiO) within the semilocal and hybrid density functionals are fraught with difficulties. In particular, MnO is known to be the most difficult case for almost all common semilocal and hybrid density approximations. Here, we show that the meta-generalized gradient approximation (meta-GGA) constructed from the cuspless hydrogen model and Pauli kinetic energy density (MGGAC) can lead to the correct ground state of MnO. The relative energy differences of zinc-blende (zb) and rock-salt (rs) structures as computed using MGGAC are found to be in nice agreement with those obtained from high-level correlation methods like the random phase approximation or quantum Monte Carlo techniques. Besides, we have also applied the onsite hybrid functionals (closely related to DFT+U ) based on GGA and meta-GGA functionals, and it is shown that a relatively high amount of Hartree-Fock exchange is necessary to obtain the correct ground-state structure. Our present investigation suggests that the semilocal MGGAC and onsite hybrids, both being computationally cheap, as methods of choice for the calculation of the relative stability of antiferromagnetic transition-metal oxides having potential applications in solid-state physics and structural chemistry.
△ Less
Submitted 16 December, 2021;
originally announced December 2021.
-
Density analysis for estimating the degree of on-site correlation on transition-metal atoms in extended systems
Authors:
Leila Kalantari,
Fabien Tran,
Peter Blaha
Abstract:
In the context of the modified Becke-Johnson (mBJ) potential, we recently underlined that $\bar{g}$, the average of $\left\vert\nablaρ\right\vert/ρ$ in the unit cell, has markedly different values in transition-metal oxides and pure transition metals [Tran et al., J. Appl. Phys. 126, 110902 (2019)]. However, since $\bar{g}$ is a constant it is not able to provide local information about a particul…
▽ More
In the context of the modified Becke-Johnson (mBJ) potential, we recently underlined that $\bar{g}$, the average of $\left\vert\nablaρ\right\vert/ρ$ in the unit cell, has markedly different values in transition-metal oxides and pure transition metals [Tran et al., J. Appl. Phys. 126, 110902 (2019)]. However, since $\bar{g}$ is a constant it is not able to provide local information about a particular atom in the system. Furthermore, while $\overline{g}$ can be used only for periodic bulk solids, a local (i.e., position-dependent) version would allow us to consider also low-dimensional systems and interfaces. Such a local function has been proposed by Rauch et al. [J. Chem. Theory Comput. 16, 2654 (2020)] for the local mBJ potential. Actually, a local version of $\overline{g}$, or of another similar quantity like the reduced density gradient $\overline{s}$, could also be used in the framework of other methods. Here, we explored the idea to use such a local function $\tilde{g}$ (or $\tilde{s}$), defined as the average of $g$ (or $s$) over a certain region around a transition-metal atom, to estimate the degree of on-site correlation on this atom. We found a large difference in our correlation estimators between non-correlated and correlated materials, proving its usefulness and reliability. Our estimators can subsequently be used to determine whether or not a Hubbard $U$ on-site correction in the DFT+$U$ method should be applied to a particular atom. This is particularly interesting in cases where the degree of correlation of the transition-metal atoms is not clear, like interfaces between correlated and non-correlated materials or oxygen-covered metal surfaces. In such cases, our estimators could also be used for an interpolation of $U$ between correlated and non-correlated atoms.
△ Less
Submitted 17 August, 2021;
originally announced August 2021.
-
Bandgap of two-dimensional materials: Thorough assessment of modern exchange-correlation functionals
Authors:
Fabien Tran,
Jan Doumont,
Leila Kalantari,
Peter Blaha,
Tomáš Rauch,
Pedro Borlido,
Silvana Botti,
Miguel A. L. Marques,
Abhilash Patra,
Subrata Jana,
Prasanjit Samal
Abstract:
The density functional theory (DFT) approximations that are the most accurate for the calculation of band gap of bulk materials are hybrid functionals like HSE06, the MBJ potential, and the GLLB-SC potential. More recently, generalized gradient approximations (GGA), like HLE16, or meta-GGAs, like (m)TASK, have proven to be also quite accurate for the band gap. Here, the focus is on 2D materials an…
▽ More
The density functional theory (DFT) approximations that are the most accurate for the calculation of band gap of bulk materials are hybrid functionals like HSE06, the MBJ potential, and the GLLB-SC potential. More recently, generalized gradient approximations (GGA), like HLE16, or meta-GGAs, like (m)TASK, have proven to be also quite accurate for the band gap. Here, the focus is on 2D materials and the goal is to provide a broad overview of the performance of DFT functionals by considering a large test set of 298 2D systems. The present work is an extension of our recent studies [Rauch et al., Phys. Rev. B 101, 245163 (2020) and Patra et al., J. Phys. Chem. C 125, 11206 (2021)]. Due to the lack of experimental results for the band gap of 2D systems, $G_{0}W_{0}$ results were taken as reference. It is shown that the GLLB-SC potential and mTASK functional provide the band gaps that are the closest to $G_{0}W_{0}$. Following closely, the local MBJ potential has a pretty good accuracy that is similar to the accuracy of the more expensive hybrid functional HSE06.
△ Less
Submitted 3 January, 2022; v1 submitted 25 May, 2021;
originally announced May 2021.
-
WloopPHI: A tool for ab initio characterization of Weyl semimetals
Authors:
Himanshu Saini,
Magdalena Laurien,
Peter Blaha,
Oleg Rubel
Abstract:
WloopPHI is a Python code that expands the features of WIEN2k, a full-potential all-electron density functional theory package, by the characterization of Weyl semimetals. It enables the calculation of the chirality (or "monopole charge") associated with Weyl nodes and nodal lines. The theoretical methodology for the calculation of the chirality is based on an extended Wilson loop method and a Ber…
▽ More
WloopPHI is a Python code that expands the features of WIEN2k, a full-potential all-electron density functional theory package, by the characterization of Weyl semimetals. It enables the calculation of the chirality (or "monopole charge") associated with Weyl nodes and nodal lines. The theoretical methodology for the calculation of the chirality is based on an extended Wilson loop method and a Berry phase approach. We validate the code using TaAs, which is a well-characterized Weyl semimetal, both theoretically and experimentally. Afterwards, we applied the method to the characterization of YRh$_6$Ge$_4$ and found two sets of Weyl points (ca. 0.2 eV below the Fermi energy) together with a topological nodal line (protected by mirror symmetry) crossing the Fermi energy and mapped their chiralities.
△ Less
Submitted 8 September, 2021; v1 submitted 18 August, 2020;
originally announced August 2020.
-
Pristine quantum criticality in a Kondo semimetal
Authors:
W. T. Fuhrman,
A. Sidorenko,
J. Hänel,
H. Winkler,
A. Prokofiev,
J. A. Rodriguez-Rivera,
Y. Qiu,
P. Blaha,
Q. Si,
C. L. Broholm,
S. Paschen
Abstract:
The observation of quantum criticality in diverse classes of strongly correlated electron systems has been instrumental in establishing ordering principles, discovering new phases, and identifying the relevant degrees of freedom and interactions. At focus so far have been insulators and metals. Semimetals, which are of great current interest as candidate phases with nontrivial topology, are much l…
▽ More
The observation of quantum criticality in diverse classes of strongly correlated electron systems has been instrumental in establishing ordering principles, discovering new phases, and identifying the relevant degrees of freedom and interactions. At focus so far have been insulators and metals. Semimetals, which are of great current interest as candidate phases with nontrivial topology, are much less explored in experiments. Here we study the Kondo semimetal CeRu$_4$Sn$_6$ by magnetic susceptibility, specific heat, and inelastic neutron scattering experiments. The power-law divergence of the magnetic Grünesien ratio reveals that, surprisingly, this compound is quantum critical without tuning. The dynamical energy over temperature scaling in the neutron response, seen throughout the Brillouin zone, as well as the temperature dependence of the static uniform susceptibility indicate that temperature is the only energy scale in the criticality. Such behavior, which has been associated with Kondo destruction quantum criticality in metallic systems, may well be generic in the semimetal setting.
△ Less
Submitted 28 June, 2021; v1 submitted 18 July, 2020;
originally announced July 2020.
-
Perturbation approach to ab initio effective mass calculations
Authors:
Oleg Rubel,
Fabien Tran,
Xavier Rocquefelte,
Peter Blaha
Abstract:
A degenerate perturbation $k\cdot p$ approach for effective mass calculations is implemented in the all-electron density functional theory (DFT) package WIEN2k. The accuracy is tested on major group IVA, IIIA-VA, and IIB-VIA semiconductor materials. Then, the effective mass in graphene and CuI with defects is presented as illustrative applications. For states with significant Cu-d character additi…
▽ More
A degenerate perturbation $k\cdot p$ approach for effective mass calculations is implemented in the all-electron density functional theory (DFT) package WIEN2k. The accuracy is tested on major group IVA, IIIA-VA, and IIB-VIA semiconductor materials. Then, the effective mass in graphene and CuI with defects is presented as illustrative applications. For states with significant Cu-d character additional local orbitals with higher principal quantum numbers (more radial nodes) have to be added to the basis set in order to converge the results of the perturbation theory. Caveats related to a difference between velocity and momentum matrix elements are discussed in the context of application of the method to non-local potentials, such as Hartree-Fock/DFT hybrid functionals and DFT+U.
△ Less
Submitted 7 July, 2020;
originally announced July 2020.
-
Shortcomings of meta-GGA functionals when describing magnetism
Authors:
Fabien Tran,
Guillaume Baudesson,
Jesús Carrete,
Georg K. H. Madsen,
Peter Blaha,
Karlheinz Schwarz,
David J. Singh
Abstract:
Several recent studies have shown that SCAN, a functional belonging to the meta-generalized gradient approximation (MGGA) family, leads to significantly overestimated magnetic moments in itinerant ferromagnetic metals. However, this behavior is not inherent to the MGGA level of approximation since TPSS, for instance, does not lead to such severe overestimations. In order to provide a broader view…
▽ More
Several recent studies have shown that SCAN, a functional belonging to the meta-generalized gradient approximation (MGGA) family, leads to significantly overestimated magnetic moments in itinerant ferromagnetic metals. However, this behavior is not inherent to the MGGA level of approximation since TPSS, for instance, does not lead to such severe overestimations. In order to provide a broader view of the accuracy of MGGA functionals for magnetism, we extend the assessment to more functionals, but also to antiferromagnetic solids. The results show that to describe magnetism there is overall no real advantage in using a MGGA functional compared to GGAs. For both types of approximation, an improvement in ferromagnetic metals is necessarily accompanied by a deterioration (underestimation) in antiferromagnetic insulators, and vice-versa. We also provide some analysis in order to understand in more detail the relation between the mathematical form of the functionals and the results.
△ Less
Submitted 2 June, 2020; v1 submitted 9 April, 2020;
originally announced April 2020.
-
On the calculation of the bandgap of periodic solids with MGGA functionals using the total energy
Authors:
Fabien Tran,
Jan Doumont,
Peter Blaha,
Miguel A. L. Marques,
Silvana Botti,
Albert P. Bartók
Abstract:
During the last few years, it has become more and more clear that functionals of the meta generalized gradient approximation (MGGA) are more accurate than GGA functionals for the geometry and energetics of electronic systems. However, MGGA functionals are also potentially more interesting for the electronic structure, in particular when the potential is non-multiplicative (i.e., when MGGAs are imp…
▽ More
During the last few years, it has become more and more clear that functionals of the meta generalized gradient approximation (MGGA) are more accurate than GGA functionals for the geometry and energetics of electronic systems. However, MGGA functionals are also potentially more interesting for the electronic structure, in particular when the potential is non-multiplicative (i.e., when MGGAs are implemented in the generalized Kohn-Sham framework), which may help to get more accurate bandgaps. Here, we show that the calculation of bandgap of solids with MGGA functionals can be done very accurately also in a non-self-consistent manner. This scheme uses only the total energy and can, therefore, be very useful when the self-consistent implementation of a particular MGGA functional is not available. Since self-consistent MGGA calculations may be difficult to converge, the non-self-consistent scheme may also help to speed-up the calculations. Furthermore, it can be applied to any other types of functionals, for which the implementation of the corresponding potential is not trivial.
△ Less
Submitted 1 November, 2019; v1 submitted 2 September, 2019;
originally announced September 2019.
-
Semilocal exchange-correlation potentials for solid-state calculations: Current status and future directions
Authors:
Fabien Tran,
Jan Doumont,
Leila Kalantari,
Ahmad W. Huran,
Miguel A. L. Marques,
Peter Blaha
Abstract:
Kohn-Sham (KS) density functional theory (DFT) is a very efficient method for calculating various properties of solids as, for instance, the total energy, the electron density, or the electronic band structure. The KS-DFT method leads to rather fast calculations, however the accuracy depends crucially on the chosen approximation for the exchange and correlation (xc) functional $E_{\text{xc}}$ and/…
▽ More
Kohn-Sham (KS) density functional theory (DFT) is a very efficient method for calculating various properties of solids as, for instance, the total energy, the electron density, or the electronic band structure. The KS-DFT method leads to rather fast calculations, however the accuracy depends crucially on the chosen approximation for the exchange and correlation (xc) functional $E_{\text{xc}}$ and/or potential $v_{\text{xc}}$. Here, an overview of xc methods to calculate the electronic band structure is given, with the focus on the so-called semilocal methods that are the fastest in KS-DFT and allow to treat systems containing up to thousands of atoms. Among them, there is the modified Becke-Johnson potential that is widely used to calculate the fundamental band gap of semiconductors and insulators. The accuracy for other properties like the magnetic moment or the electron density, that are also determined directly by $v_{\text{xc}}$, is also discussed.
△ Less
Submitted 18 September, 2019; v1 submitted 5 August, 2019;
originally announced August 2019.
-
Nonlocal van der Waals functionals for solids: Choosing an appropriate one
Authors:
Fabien Tran,
Leila Kalantari,
Boubacar Traoré,
Xavier Rocquefelte,
Peter Blaha
Abstract:
The nonlocal van der Waals (NL-vdW) functionals [Dion et al., Phys. Rev. Lett. 92, 246401 (2004)] are being applied more and more frequently in solid-state physics, since they have shown to be much more reliable than the traditional semilocal functionals for systems where weak interactions play a major role. However, a certain number of NL-vdW functionals have been proposed during the last few yea…
▽ More
The nonlocal van der Waals (NL-vdW) functionals [Dion et al., Phys. Rev. Lett. 92, 246401 (2004)] are being applied more and more frequently in solid-state physics, since they have shown to be much more reliable than the traditional semilocal functionals for systems where weak interactions play a major role. However, a certain number of NL-vdW functionals have been proposed during the last few years, such that it is not always clear which one should be used. In this work, an assessment of NL-vdW functionals is presented. Our test set consists of weakly bound solids, namely rare gases, layered systems like graphite, and molecular solids, but also strongly bound solids in order to provide a more general conclusion about the accuracy of NL-vdW functionals for extended systems. We found that among the tested functionals, rev-vdW-DF2 [Hamada, Phys. Rev. B 89, 121103(R) (2014)] is very accurate for weakly bound solids, but also quite reliable for strongly bound solids.
△ Less
Submitted 5 June, 2019; v1 submitted 1 March, 2019;
originally announced March 2019.
-
Local geometry around B atoms in B/Si(111) from polarized x-ray absorption spectroscopy
Authors:
Saleem Ayaz Khan,
Martin Vondráček,
Peter Blaha,
Kateřina Horáková,
Jan Minár,
Ondřej Šipr,
Vladimír Cháb
Abstract:
The arrangement of B~atoms in a doped Si(111)-$(\sqrt{3}\times\sqrt{3})R30^{\circ}$:B system was studied using near-edge x-ray absorption fine structure (NEXAFS). Boron atoms were deposited via segregation from the bulk by flashing the sample repeatedly. The positions of B~atoms are determined by comparing measured polarized (angle-dependent) NEXAFS spectra with spectra calculated for various stru…
▽ More
The arrangement of B~atoms in a doped Si(111)-$(\sqrt{3}\times\sqrt{3})R30^{\circ}$:B system was studied using near-edge x-ray absorption fine structure (NEXAFS). Boron atoms were deposited via segregation from the bulk by flashing the sample repeatedly. The positions of B~atoms are determined by comparing measured polarized (angle-dependent) NEXAFS spectra with spectra calculated for various structural models based on ab-initio total energy calculations. It is found that most of boron atoms are located in sub-surface L$_{1}^{c}$ positions, beneath a Si atom. However, depending on the preparation method a significant portion of B~atoms may be located elsewhere. A possible location of these non-L$_{1}^{c}$-atoms is at the surface, next to those Si atoms which form the $(\sqrt{3}\times\sqrt{3})R30^{\circ}$ reconstruction.
△ Less
Submitted 17 February, 2019;
originally announced February 2019.
-
Limitations of the DFT-1/2 method for covalent semiconductors and transition-metal oxides
Authors:
Jan Doumont,
Fabien Tran,
Peter Blaha
Abstract:
The DFT-1/2 method in density functional theory [L. G. Ferreira et al., Phys. Rev. B 78, 125116 (2008)] aims to provide accurate band gaps at the computational cost of semilocal calculations. The method has shown promise in a large number of cases, however some of its limitations or ambiguities on how to apply it to covalent semiconductors have been pointed out recently [K.-H. Xue et al., Comput.…
▽ More
The DFT-1/2 method in density functional theory [L. G. Ferreira et al., Phys. Rev. B 78, 125116 (2008)] aims to provide accurate band gaps at the computational cost of semilocal calculations. The method has shown promise in a large number of cases, however some of its limitations or ambiguities on how to apply it to covalent semiconductors have been pointed out recently [K.-H. Xue et al., Comput. Mater. Science 153, 493 (2018)]. In this work, we investigate in detail some of the problems of the DFT-1/2 method with a focus on two classes of materials: covalently bonded semiconductors and transition-metal oxides. We argue for caution in the application of DFT-1/2 to these materials, and the condition to get an improved band gap is a spatial separation of the orbitals at the valence band maximum and conduction band minimum.
△ Less
Submitted 15 January, 2019;
originally announced January 2019.
-
Giant spontaneous Hall effect in a nonmagnetic Weyl-Kondo semimetal
Authors:
Sami Dzsaber,
Xinlin Yan,
Mathieu Taupin,
Gaku Eguchi,
Andrey Prokofiev,
Toni Shiroka,
Peter Blaha,
Oleg Rubel,
Sarah E. Grefe,
Hsin-Hua Lai,
Qimiao Si,
Silke Paschen
Abstract:
Nontrivial topology in condensed matter systems enriches quantum states of matter, to go beyond either the classification into metals and insulators in terms of conventional band theory or that of symmetry broken phases by Landau's order parameter framework. So far, focus has been on weakly interacting systems, and little is known about the limit of strong electron correlations. Heavy fermion syst…
▽ More
Nontrivial topology in condensed matter systems enriches quantum states of matter, to go beyond either the classification into metals and insulators in terms of conventional band theory or that of symmetry broken phases by Landau's order parameter framework. So far, focus has been on weakly interacting systems, and little is known about the limit of strong electron correlations. Heavy fermion systems are a highly versatile platform to explore this regime. Here we report the discovery of a giant spontaneous Hall effect in the Kondo semimetal Ce3Bi4Pd3 that is noncentrosymmetric but preserves time reversal symmetry. We attribute this finding to Weyl nodes - singularities of the Berry curvature - that emerge in the immediate vicinity of the Fermi level due to the Kondo interaction. We stress that this phenomenon is distinct from the previously detected anomalous Hall effect in materials with broken time reversal symmetry; instead, it manifests an extreme topological response that requires a beyond-perturbation-theory description of the previously proposed nonlinear Hall effect. The large magnitude of the effect in even tiny electric and zero magnetic fields, as well as its robust bulk nature may aid the exploitation in topological quantum devices.
△ Less
Submitted 1 February, 2021; v1 submitted 7 November, 2018;
originally announced November 2018.
-
Density-gradient-free variable in exchange-correlation functionals for detecting inhomogeneities in the electron density
Authors:
Fabien Tran,
Peter Blaha
Abstract:
An alternative type of approximation for the exchange and correlation functional in density functional theory is proposed. This approximation depends on a variable $u$ that is able to detect inhomogeneities in the electron density $ρ$ without using derivatives of $ρ$. Instead, $u$ depends on the orbital energies which can also be used to measure how a system differs from the homogeneous electron g…
▽ More
An alternative type of approximation for the exchange and correlation functional in density functional theory is proposed. This approximation depends on a variable $u$ that is able to detect inhomogeneities in the electron density $ρ$ without using derivatives of $ρ$. Instead, $u$ depends on the orbital energies which can also be used to measure how a system differs from the homogeneous electron gas. Starting from the functional of Perdew, Burke, and Ernzerhof (PBE) [Phys. Rev. Lett. 77, 3865 (1996)], a functional depending on $u$ is constructed. Tests on the lattice constant, bulk modulus, and cohesive energy of solids show that this $u$-dependent PBE-like functional is on average as accurate as the original PBE or its solid-state version PBEsol. Since $u$ carries more nonlocality than the reduced density gradient $s$ used in functionals of the generalized gradient approximation (GGA) like PBE and $α$ used in meta-GGAs, it will be certainly useful for the future development of more accurate exchange-correlation functionals.
△ Less
Submitted 23 December, 2018; v1 submitted 27 August, 2018;
originally announced August 2018.
-
Orbital-free approximations to the kinetic-energy density in exchange-correlation MGGA functionals: tests on solids
Authors:
Fabien Tran,
Péter Kovács,
Leila Kalantari,
Georg K. H. Madsen,
Peter Blaha
Abstract:
A recent study of Mejia-Rodriguez and Trickey [Phys. Rev. A 96, 052512 (2017)] showed that the deorbitalization procedure (replacing the exact Kohn-Sham kinetic-energy density by an approximate orbital-free expression) applied to exchange-correlation functionals of the meta-generalized gradient approximation (MGGA) can lead to important changes in the results for molecular properties. For the pres…
▽ More
A recent study of Mejia-Rodriguez and Trickey [Phys. Rev. A 96, 052512 (2017)] showed that the deorbitalization procedure (replacing the exact Kohn-Sham kinetic-energy density by an approximate orbital-free expression) applied to exchange-correlation functionals of the meta-generalized gradient approximation (MGGA) can lead to important changes in the results for molecular properties. For the present work, the deorbitalization of MGGA functionals is further investigated by considering various properties of solids. It is shown that depending on the MGGA, common orbital-free approximations to the kinetic-energy density can be sufficiently accurate for the lattice constant, bulk modulus, and cohesive energy. For the band gap, calculated with the modified Becke-Johnson MGGA potential, the deorbitalization has a larger impact on the results.
△ Less
Submitted 19 July, 2018;
originally announced July 2018.
-
Interplay of Magnetism and Transport in HoBi
Authors:
H. -Y. Yang,
J. Gaudet,
A. A. Aczel,
D. E. Graf,
P. Blaha,
B. D. Gaulin,
Fazel Tafti
Abstract:
We report the observation of an extreme magnetoresistance (XMR) in HoBi with a large magnetic moment from Ho f-electrons. Neutron scattering is used to determine the magnetic wave vectors across several metamagnetic (MM) transitions on the phase diagram of HoBi. Unlike other magnetic rare-earth monopnictides, the field dependence of resistivity in HoBi is non-monotonic and reveals clear signatures…
▽ More
We report the observation of an extreme magnetoresistance (XMR) in HoBi with a large magnetic moment from Ho f-electrons. Neutron scattering is used to determine the magnetic wave vectors across several metamagnetic (MM) transitions on the phase diagram of HoBi. Unlike other magnetic rare-earth monopnictides, the field dependence of resistivity in HoBi is non-monotonic and reveals clear signatures of every metamagnetic transition in the low-temperature and low-field regime, at T < 2 K and H < 2.3 T. The XMR appears at H > 2.3 T after all the metamagnetic transitions are complete and the system is spin-polarized by the external magnetic field. The existence of an onset field for XMR and the intimate connection between magnetism and transport in HoBi are unprecedented among the magnetic rare-earth monopnictides. Therefore, HoBi provides a unique opportunity to understand the electrical transport in magnetic XMR semimetals.
△ Less
Submitted 30 May, 2018;
originally announced May 2018.
-
Subsurface Cation Vacancy Stabilization of the Magnetite (001) Surface
Authors:
R. Bliem,
E. McDermott,
P. Ferstl,
M. Setvin,
O. Gamba,
J. Pavelec,
M. A. Schneider,
M. Schmid,
U. Diebold,
P. Blaha,
L. Hammer,
G. S. Parkinson
Abstract:
Iron oxides play an increasingly prominent role in heterogeneous catalysis, hydrogen production, spintronics and drug delivery. The surface or material interface can be performance limiting in these applications, so it is vital to determine accurate atomic-scale structures for iron oxides and understand why they form. Using a combination of quantitative low-energy electron diffraction, scanning tu…
▽ More
Iron oxides play an increasingly prominent role in heterogeneous catalysis, hydrogen production, spintronics and drug delivery. The surface or material interface can be performance limiting in these applications, so it is vital to determine accurate atomic-scale structures for iron oxides and understand why they form. Using a combination of quantitative low-energy electron diffraction, scanning tunneling microscopy, and density functional theory calculations, we show that an ordered array of subsurface iron vacancies and interstitials underlies the well-known (rt2xrt2)R45° reconstruction of Fe3O4(001). This hitherto unobserved stabilization mechanism occurs because the iron oxides prefer to redistribute cations in the lattice in response to oxidizing or reducing environments. Many other metal oxides also achieve stoichiometric variation in this way, so such surface structures are likely commonplace.
△ Less
Submitted 13 April, 2018;
originally announced April 2018.
-
Assessment of the GLLB-SC potential for solid-state properties and attempts for improvement
Authors:
Fabien Tran,
Sohaib Ehsan,
Peter Blaha
Abstract:
Based on the work of Gritsenko et al. (GLLB) [Phys. Rev. A 51, 1944 (1995)], the method of Kuisma et al. [Phys. Rev. B 82, 115106 (2010)] to calculate the band gap in solids was shown to be much more accurate than the common local density approximation (LDA) and generalized gradient approximation (GGA). The main feature of the GLLB-SC potential (SC stands for solid and correlation) is to lead to a…
▽ More
Based on the work of Gritsenko et al. (GLLB) [Phys. Rev. A 51, 1944 (1995)], the method of Kuisma et al. [Phys. Rev. B 82, 115106 (2010)] to calculate the band gap in solids was shown to be much more accurate than the common local density approximation (LDA) and generalized gradient approximation (GGA). The main feature of the GLLB-SC potential (SC stands for solid and correlation) is to lead to a nonzero derivative discontinuity that can be conveniently calculated and then added to the Kohn-Sham band gap for a comparison with the experimental band gap. In this work, a thorough comparison of GLLB-SC with other methods, e.g., the modified Becke-Johnson (mBJ) potential [F. Tran and P. Blaha, Phys. Rev. Lett. 102, 226401 (2009)], for electronic, magnetic, and density-related properties is presented. It is shown that for the band gap, GLLB-SC does not perform as well as mBJ for systems with a small band gap and strongly correlated systems, but is on average of similar accuracy as hybrid functionals. The results on itinerant metals indicate that GLLB-SC overestimates significantly the magnetic moment (much more than mBJ does), but leads to excellent results for the electric field gradient, for which mBJ is in general not recommended. In the aim of improving the results, variants of the GLLB-SC potential are also tested.
△ Less
Submitted 18 December, 2017;
originally announced December 2017.
-
Simple way to apply nonlocal van der Waals functionals within all-electron methods
Authors:
Fabien Tran,
Julia Stelzl,
David Koller,
Thomas Ruh,
Peter Blaha
Abstract:
The method based on fast Fourier transforms proposed by G. Román-Pérez and J. M. Soler [Phys. Rev. Lett. 103, 096102 (2009)], which allows for a computationally fast implementation of the nonlocal van der Waals (vdW) functionals, has significantly contributed to making the vdW functionals popular in solid-state physics. However, the Román-Pérez-Soler method relies on a plane-wave expansion of the…
▽ More
The method based on fast Fourier transforms proposed by G. Román-Pérez and J. M. Soler [Phys. Rev. Lett. 103, 096102 (2009)], which allows for a computationally fast implementation of the nonlocal van der Waals (vdW) functionals, has significantly contributed to making the vdW functionals popular in solid-state physics. However, the Román-Pérez-Soler method relies on a plane-wave expansion of the electron density; therefore it can not be applied readily to all-electron densities for which an unaffordable number of plane waves would be required for an accurate expansion. In this work, we present the results for the lattice constant and binding energy of solids that were obtained by applying a smoothing procedure to the all-electron density calculated with the linearized augmented plane-wave method. The smoothing procedure has the advantages of being very simple to implement, basis-set independent, and allowing the calculation of the potential. It is also shown that the results agree very well with those from the literature that were obtained with the projector augmented wave method.
△ Less
Submitted 27 July, 2017; v1 submitted 22 May, 2017;
originally announced May 2017.
-
Local environment effects in the magnetic properties and electronic structure of disordered FePt
Authors:
Saleem Ayaz Khan,
Jan Minár,
Hubert Ebert,
Peter Blaha,
Ondřej Šipr
Abstract:
Local aspects of magnetism of disordered FePt are investigated by ab initio fully relativistic full potential calculations, employing the supercell approach and the coherent potential approximation (CPA). The focus is on trends of the spin and orbital magnetic moments with chemical composition and with bond lengths around the Fe and Pt atoms. A small but distinct difference between average magneti…
▽ More
Local aspects of magnetism of disordered FePt are investigated by ab initio fully relativistic full potential calculations, employing the supercell approach and the coherent potential approximation (CPA). The focus is on trends of the spin and orbital magnetic moments with chemical composition and with bond lengths around the Fe and Pt atoms. A small but distinct difference between average magnetic moments obtained when using the supercells and when relying on the CPA is identified and linked to the neglect of the Madelung potential in the CPA.
△ Less
Submitted 14 October, 2016;
originally announced October 2016.
-
NMR shieldings from density functional perturbation theory: GIPAW versus all-electron calculations
Authors:
G. A. de Wijs,
R. Laskowski,
P. Blaha,
R. W. A. Havenith,
G. Kresse,
M. Marsman
Abstract:
We present a benchmark of the density functional linear response calculation of NMR shieldings within the Gauge-Including Projector-Augmented-Wave method against all-electron Augmented-Plane-Wave$+$local-orbital and uncontracted Gaussian basis set results for NMR shieldings in molecular and solid state systems. In general, excellent agreement between the aforementioned methods is obtained. Scalar…
▽ More
We present a benchmark of the density functional linear response calculation of NMR shieldings within the Gauge-Including Projector-Augmented-Wave method against all-electron Augmented-Plane-Wave$+$local-orbital and uncontracted Gaussian basis set results for NMR shieldings in molecular and solid state systems. In general, excellent agreement between the aforementioned methods is obtained. Scalar relativistic effects are shown to be quite large for nuclei in molecules in the deshielded limit. The small component makes up a substantial part of the relativistic corrections.
△ Less
Submitted 13 September, 2016;
originally announced September 2016.
-
Approximations to the exact exchange potential: KLI versus semilocal
Authors:
Fabien Tran,
Peter Blaha,
Markus Betzinger,
Stefan Blügel
Abstract:
In the search for an accurate and computationally efficient approximation to the exact exchange potential of Kohn-Sham density functional theory, we recently compared various semilocal exchange potentials to the exact one [F. Tran et al., Phys. Rev. B 91, 165121 (2015)]. It was concluded that the Becke-Johnson (BJ) potential is a very good starting point, but requires the use of empirical paramete…
▽ More
In the search for an accurate and computationally efficient approximation to the exact exchange potential of Kohn-Sham density functional theory, we recently compared various semilocal exchange potentials to the exact one [F. Tran et al., Phys. Rev. B 91, 165121 (2015)]. It was concluded that the Becke-Johnson (BJ) potential is a very good starting point, but requires the use of empirical parameters to obtain good agreement with the exact exchange potential. In the present work, we extend the comparison by considering the Krieger-Li-Iafrate (KLI) approximation, which is a beyond-semilocal approximation. It is shown that overall the KLI and BJ-based potentials are the most reliable approximations to the exact exchange potential, however, sizeable differences, especially for the antiferromagnetic transition-metal oxides, can be obtained.
△ Less
Submitted 30 August, 2016;
originally announced August 2016.
-
Magnetocrystalline anisotropy of FePt: a detailed view
Authors:
Saleem Ayaz Khan,
Peter Blaha,
Hubert Ebert,
Jan Minár,
Ondřej Šipr
Abstract:
To get a reliable ab-initio value for the magneto-crystalline anisotropy (MCA) energy of FePt, we employ the full-potential linearized augmented plane wave (FLAPW) method and the full-potential Korringa-Kohn-Rostoker (KKR) Green function method. The MCA energies calculated by both methods are in a good agreement with each other. As the calculated MCA energy significantly differs from experiment, i…
▽ More
To get a reliable ab-initio value for the magneto-crystalline anisotropy (MCA) energy of FePt, we employ the full-potential linearized augmented plane wave (FLAPW) method and the full-potential Korringa-Kohn-Rostoker (KKR) Green function method. The MCA energies calculated by both methods are in a good agreement with each other. As the calculated MCA energy significantly differs from experiment, it is clear that many-body effects beyond the local density approximation are essential. It is not really important whether relativistic effects for FePt are accounted for by solving the full Dirac equation or whether the spin-orbit coupling (SOC) is treated as a correction to the scalar-relativistic Hamiltonian. From the analysis of the dependence of the MCA energy on the magnetization angle and on the SOC strength it follows that the main mechanism of MCA in FePt can be described within second order perturbation theory. However, a distinct contribution not accountable for by second order perturbation theory is present as well.
△ Less
Submitted 14 October, 2016; v1 submitted 1 April, 2016;
originally announced April 2016.
-
Rungs 1 to 4 of DFT Jacob's ladder: extensive test on the lattice constant, bulk modulus, and cohesive energy of solids
Authors:
Fabien Tran,
Julia Stelzl,
Peter Blaha
Abstract:
A large panel of old and recently proposed exchange-correlation functionals belonging to rungs 1 to 4 of Jacob's ladder of density functional theory are tested (with and without a dispersion correction term) for the calculation of the lattice constant, bulk modulus, and cohesive energy of solids. Particular attention will be paid to the functionals MGGA_MS2 [J. Sun et al., J. Chem. Phys. 138, 0441…
▽ More
A large panel of old and recently proposed exchange-correlation functionals belonging to rungs 1 to 4 of Jacob's ladder of density functional theory are tested (with and without a dispersion correction term) for the calculation of the lattice constant, bulk modulus, and cohesive energy of solids. Particular attention will be paid to the functionals MGGA_MS2 [J. Sun et al., J. Chem. Phys. 138, 044113 (2013)], mBEEF [J. Wellendorff et al., J. Chem. Phys. 140, 144107 (2014)], and SCAN [J. Sun et al., Phys. Rev. Lett. 115, 036402 (2015)] that are approximations of the meta-generalized gradient type and were developed with the goal to be universally good. Another goal is also to determine for which semilocal functionals and groups of solids it is beneficial (or not necessary) to use the Hartree-Fock exchange or a dispersion correction term.
△ Less
Submitted 4 March, 2016;
originally announced March 2016.
-
Transition from reconstruction towards thin film on the (110) surface of strontium titanate
Authors:
Z. Wang,
A. Loon,
A. Subramanian,
S. Gerhold,
E. McDermott,
J. A. Enterkin,
M. Hieckel,
B. C. Russell,
R. J. Green,
A. Moewes,
J. Guo,
P. Blaha,
M. R. Castell,
U. Diebold,
L. D. Marks
Abstract:
The surfaces of metal oxides often are reconstructed with a geometry and composition that is considerably different from a simple termination of the bulk. Such structures can also be viewed as ultrathin films, epitaxed on a substrate. Here, the reconstructions of the SrTiO3 (110) surface are studied combining scanning tunneling microscopy, transmission electron diffraction, and X-ray absorption sp…
▽ More
The surfaces of metal oxides often are reconstructed with a geometry and composition that is considerably different from a simple termination of the bulk. Such structures can also be viewed as ultrathin films, epitaxed on a substrate. Here, the reconstructions of the SrTiO3 (110) surface are studied combining scanning tunneling microscopy, transmission electron diffraction, and X-ray absorption spectroscopy, and analyzed with density functional theory calculations. While SrTiO3 (110) invariably terminates with an overlayer of titania, with increasing density its structure switches from nx1 and 2xn. At the same time the coordination of the Ti atoms changes from a network of corner-sharing tetrahedra to a double layer of edge-shared octahedra with bridging units of octahedrally coordinated strontium. This transition from the nx1 to 2xn reconstructions is a transition from a pseudomorphically stabilized tetrahedral network towards an octahedral titania thin film with stress-relief from octahedral strontia units at the surface.
△ Less
Submitted 10 February, 2016;
originally announced February 2016.
-
Woptic: optical conductivity with Wannier functions and adaptive k-mesh refinement
Authors:
E. Assmann,
P. Wissgott,
J. Kuneš,
A. Toschi,
P. Blaha,
K. Held
Abstract:
We present an algorithm for the adaptive tetrahedral integration over the Brillouin zone of crystalline materials, and apply it to compute the optical conductivity, dc conductivity, and thermopower. For these quantities, whose contributions are often localized in small portions of the Brillouin zone, adaptive integration is especially relevant. Our implementation, the woptic package, is tied into…
▽ More
We present an algorithm for the adaptive tetrahedral integration over the Brillouin zone of crystalline materials, and apply it to compute the optical conductivity, dc conductivity, and thermopower. For these quantities, whose contributions are often localized in small portions of the Brillouin zone, adaptive integration is especially relevant. Our implementation, the woptic package, is tied into the wien2wannier framework and allows including a many-body self energy, e.g. from dynamical mean-field theory (DMFT). Wannier functions and dipole matrix elements are computed with the DFT package Wien2k and Wannier90. For illustration, we show DFT results for fcc-Al and DMFT results for the correlated metal SrVO$_3$.
△ Less
Submitted 19 January, 2016; v1 submitted 17 July, 2015;
originally announced July 2015.
-
Ab-initio perspective on the Mollwo-Ivey relation for F-centers in alkali halides
Authors:
P. Tiwald,
F. Karsai,
R. Laskowski,
S. Gräfe,
P. Blaha,
J. Burgdörfer,
L. Wirtz
Abstract:
We revisit the well-known Mollwo-Ivey relation that describes the "universal" dependence of the absorption energies of F-type color centers on the lattice constant $a$ of the alkali-halide crystals, $E_{\mbox{abs}}\propto a^{-n}.$ We perform both state-of-the-art ab-initio Quantum Chemistry and post-DFT calculations of F-center absorption spectra. By "tuning" independently the lattice constant and…
▽ More
We revisit the well-known Mollwo-Ivey relation that describes the "universal" dependence of the absorption energies of F-type color centers on the lattice constant $a$ of the alkali-halide crystals, $E_{\mbox{abs}}\propto a^{-n}.$ We perform both state-of-the-art ab-initio Quantum Chemistry and post-DFT calculations of F-center absorption spectra. By "tuning" independently the lattice constant and the atomic species we show that the scaling of the lattice constant alone (keeping the elements fixed) would yield $n=2$ in agreement with the "particle-in-the-box" model. Keeping the lattice constant fixed and changing the atomic species enables us to quantify the ion-size effects which are shown to be responsible for the exponent $n \approx 1.8$.
△ Less
Submitted 8 July, 2015;
originally announced July 2015.
-
Fermi surface of three-dimensional La1-xSrxMnO3 explored by soft-X-ray ARPES: Rhombohedral lattice distortion and its effect on magnetoresistance
Authors:
L. L. Lev,
J. Krempaský,
U. Staub,
V. A. Rogalev,
T. Schmitt,
M. Shi,
P. Blaha,
A. S. Mishchenko,
A. A. Veligzhanin,
Y. V. Zubavichus,
M. B. Tsetlin,
H. Volfová,
J. Braun,
J. Minár,
V. N. Strocov
Abstract:
Electronic structure of the three-dimensional colossal magnetoresistive perovskite La1-xSrxMnO3 has been established using soft-X-ray ARPES with its intrinsically sharp definition of three-dimensional electron momentum. The experimental results show much weaker polaronic coupling compared to the bilayer manganites and are consistent with the GGA+U band structure. The experimental Fermi surface unv…
▽ More
Electronic structure of the three-dimensional colossal magnetoresistive perovskite La1-xSrxMnO3 has been established using soft-X-ray ARPES with its intrinsically sharp definition of three-dimensional electron momentum. The experimental results show much weaker polaronic coupling compared to the bilayer manganites and are consistent with the GGA+U band structure. The experimental Fermi surface unveils the canonical topology of alternating three-dimensional electron spheres and hole cubes, with their shadow contours manifesting the rhombohedral lattice distortion. This picture has been confirmed by one-step photoemission calculations including displacement of the apical oxygen atoms. The rhombohedral distortion is neutral to the Jahn-Teller effect and thus polaronic coupling, but affects the double-exchange electron hopping and thus the colossal magnetoresistance effect.
△ Less
Submitted 1 June, 2015;
originally announced June 2015.
-
Quantized electronic fine structure with large anisotropy in ferromagnetic Fe films
Authors:
C. Seibel,
A. Nuber,
H. Bentmann,
M. Mulazzi,
P. Blaha,
G. Sangiovanni,
F. Reinert
Abstract:
We report on the spectroscopic observation of a quantized electronic fine structure near the Fermi energy in thin Fe films grown on W(110). The quantum well states are detected down to binding energies of $\sim$10 meV by angle-resolved photoelectron spectroscopy. The band dispersion of these states is found to feature a pronounced anisotropy within the surface plane: It is free-electron like along…
▽ More
We report on the spectroscopic observation of a quantized electronic fine structure near the Fermi energy in thin Fe films grown on W(110). The quantum well states are detected down to binding energies of $\sim$10 meV by angle-resolved photoelectron spectroscopy. The band dispersion of these states is found to feature a pronounced anisotropy within the surface plane: It is free-electron like along the $\bar{Γ\rm{H}}$-direction while it becomes heavy along $\bar{Γ\rm{N}}$. Density functional theory calculations identify the observed states to have both majority and minority spin character and indicate that the large anisotropy can be dependent on the number of Fe-layers and coupling to the substrate.
△ Less
Submitted 9 March, 2015;
originally announced March 2015.
-
Comparison between exact and semilocal exchange potentials: An all-electron study for solids
Authors:
Fabien Tran,
Markus Betzinger,
Peter Blaha,
Stefan Blügel
Abstract:
The exact-exchange (EXX) potential, which is obtained by solving the optimized-effective potential (OEP) equation, is compared to various approximate semilocal exchange potentials for a set of selected solids (C, Si, BN, MgO, Cu$_{2}$O, and NiO). This is done in the framework of the linearized augmented plane-wave method, which allows for a very accurate all-electron solution of electronic structu…
▽ More
The exact-exchange (EXX) potential, which is obtained by solving the optimized-effective potential (OEP) equation, is compared to various approximate semilocal exchange potentials for a set of selected solids (C, Si, BN, MgO, Cu$_{2}$O, and NiO). This is done in the framework of the linearized augmented plane-wave method, which allows for a very accurate all-electron solution of electronic structure problems in solids. In order to assess the ability of the semilocal potentials to approximate the EXX-OEP, we considered the EXX total energy, electronic structure, electric-field gradient, and magnetic moment. An attempt to parameterize a semilocal exchange potential is also reported.
△ Less
Submitted 28 January, 2015;
originally announced January 2015.
-
Quantum oscillations of the superconductor LaRu$_2$P$_2$ : comparable mass enhancement $λ\approx 1$ in Ru and Fe phosphides
Authors:
Philip J. W. Moll,
Jakob Kanter,
Ross D. McDonald,
Fedor F. Balakirev,
Peter Blaha,
Karlheinz Schwarz,
Zbigniew Bukowski,
Nikolai D. Zhigadlo,
Sergiy Katrych,
Kurt Mattenberger,
Janusz Karpinski,
Bertram Batlogg
Abstract:
We have studied the angular dependent de Haas-van Alphen oscillations of LaRu$_2$P$_2$ using magnetic torque in pulsed magnetic fields up to 60T. The observed oscillation frequencies are in excellent agreement with the geometry of the calculated Fermi surface. The temperature dependence of the oscillation amplitudes reveals effective masses m*($α$)=0.71 and m*($β$)=0.99 m$_e$, which are enhanced o…
▽ More
We have studied the angular dependent de Haas-van Alphen oscillations of LaRu$_2$P$_2$ using magnetic torque in pulsed magnetic fields up to 60T. The observed oscillation frequencies are in excellent agreement with the geometry of the calculated Fermi surface. The temperature dependence of the oscillation amplitudes reveals effective masses m*($α$)=0.71 and m*($β$)=0.99 m$_e$, which are enhanced over the calculated band mass by $λ^{cyc}$ of 0.8. We find a similar enhancement $λ^γ \approx 1$ in comparing the measured electronic specific heat ($γ= 11.5$ mJ/mol K$^2$) with the total DOS from band structure calculations. Remarkably, very similar mass enhancements have been reported in other pnictides LaFe$_2$P$_2$, LaFePO ($T_c \approx 4K$), and LaRuPO, independent of whether they are superconducting or not. This is contrary to the common perceptions that the normal state quasi-particle renormalizations reflect the strength of the superconducting paring mechanism and leads to new questions about pairing in isostructural and isoelectronic Ru- and Fe-pnictide superconductors.
△ Less
Submitted 28 October, 2014;
originally announced October 2014.
-
Nonmagnetic and ferromagnetic fcc cerium studied with one-electron methods
Authors:
Fabien Tran,
Ferenc Karsai,
Peter Blaha
Abstract:
Density functional theory was used to study the nonmagnetic (NM) and ferromagnetic (FM) phases of face-centered cubic cerium. Functionals of four levels of approximations for the exchange-correlation energy were used: LDA, PBE, LDA/PBE+$U$, and YS-PBEh. The latter two contain an adjustable parameter, the onsite Coulomb repulsion parameter $U$ for LDA/PBE+$U$ and the fraction $α_x$ of Hartree-Fock…
▽ More
Density functional theory was used to study the nonmagnetic (NM) and ferromagnetic (FM) phases of face-centered cubic cerium. Functionals of four levels of approximations for the exchange-correlation energy were used: LDA, PBE, LDA/PBE+$U$, and YS-PBEh. The latter two contain an adjustable parameter, the onsite Coulomb repulsion parameter $U$ for LDA/PBE+$U$ and the fraction $α_x$ of Hartree-Fock exchange for YS-PBEh, which were varied in order to study their influence on the results. By supposing that, as a first approximation, the NM and FM solutions can be identified to the observed $α$ and $γ$ phases, respectively, it is concluded that while a small value of $U$ or $α_x$ leads to the correct trend for the stability ordering of the two phases, larger values are necessary for a more appropriate (but still not satisfying) description of the electronic structure.
△ Less
Submitted 7 April, 2014; v1 submitted 19 December, 2013;
originally announced December 2013.
-
Application of Finite Strain Landau Theory To High Pressure Phase Transitions
Authors:
A. Tröster,
W. Schranz,
F. Karsai,
P. Blaha
Abstract:
In this paper we explain how to set up what is in fact the only possible consistent construction scheme for a Landau theory of high pressure phase transitions that systematically allows to take into account elastic nonlinearities. We also show how to incorporate available information on the pressure dependence of elastic constants taken from experiment or simulation. We apply our new theory to the…
▽ More
In this paper we explain how to set up what is in fact the only possible consistent construction scheme for a Landau theory of high pressure phase transitions that systematically allows to take into account elastic nonlinearities. We also show how to incorporate available information on the pressure dependence of elastic constants taken from experiment or simulation. We apply our new theory to the example of the high pressure cubic-tetragonal phase transition in Strontium Titanate, a model perovskite that has played a central role in the development of the theory of structural phase transitions. Armed with pressure dependent elastic constants calculated by density functional theory, we give a both qualitatively as well as quantitatively satisfying description of recent high precision experimental data. Our nonlinear theory also allows to predict a number of additional elastic transition anomalies that are accessible to experiment.
△ Less
Submitted 20 September, 2013;
originally announced September 2013.
-
High-pressure cupric oxide: a room-temperature multiferroic
Authors:
Xavier Rocquefelte,
Karlheinz Schwarz,
Peter Blaha,
Sanjeev Kumar,
Jeroen van den Brink
Abstract:
Multiferroic materials, in which ferroelectric and magnetic ordering coexist, are of fundamental interest for the development of multi-state memory devices that allow for electrical writing and non-destructive magnetic read-out operation. The great challenge is to create multiferroic materials that operate at room-temperature and have a large ferroelectric polarization P. Cupric oxide, CuO, is pro…
▽ More
Multiferroic materials, in which ferroelectric and magnetic ordering coexist, are of fundamental interest for the development of multi-state memory devices that allow for electrical writing and non-destructive magnetic read-out operation. The great challenge is to create multiferroic materials that operate at room-temperature and have a large ferroelectric polarization P. Cupric oxide, CuO, is promising because of its large P ~ 10^{2} μC.m^{-2}, but is unfortunately only multiferroic in a temperature range of 20 K, from 210 to 230 K. Here, using a combination of density functional theory and Monte Carlo calculations, we establish that pressure-driven phase competition induces a giant stabilization of the multiferroic phase of CuO, which at 20-40 GPa becomes stable in a domain larger than 300 K, from 0 to T > 300 K. Thus, under high-pressure, CuO is predicted to be a room-temperature multiferroic with large polarization.
△ Less
Submitted 1 September, 2013;
originally announced September 2013.
-
Hybrid functionals for solids with an optimized Hartree-Fock mixing parameter
Authors:
David Koller,
Peter Blaha,
Fabien Tran
Abstract:
(Screened) hybrid functionals are being used more and more for solid-state calculations. Usually the fraction alpha of Hartree-Fock exchange is kept fixed during the calculation, however there is no single (universal) value for alpha which systematically leads to satisfying accuracy. Instead, one could use a property of the system under consideration to determine alpha and in this way the function…
▽ More
(Screened) hybrid functionals are being used more and more for solid-state calculations. Usually the fraction alpha of Hartree-Fock exchange is kept fixed during the calculation, however there is no single (universal) value for alpha which systematically leads to satisfying accuracy. Instead, one could use a property of the system under consideration to determine alpha and in this way the functional would be more flexible and potentially more accurate. Recently, it was proposed to use the static dielectric constant epsilon for the calculation of alpha [Shimazaki and Asai, Chem. Phys. Lett. 466, 91 (2008) and Marques et al., Phys. Rev. B 83, 035119 (2011)]. We explore this idea further and propose a scheme where the connection between epsilon and alpha is optimized based on experimental band gaps. epsilon, and thus alpha, is recalculated at each iteration of the self-consistent procedure. We present results for the band gap and lattice constant of various semiconductors and insulators with this procedure. In addition, we show that this approach can also be combined with a non-self-consistent hybrid approximation to speed up the calculations considerably, while retaining an excellent accuracy in most cases.
△ Less
Submitted 16 April, 2013;
originally announced April 2013.
-
CO Induced Adatom Sintering in a Model Catalyst: Pd/Fe3O4
Authors:
Gareth S. Parkinson,
Zbynek Novotny,
Giacomo Argentero,
Michael Schmid,
Jiří Pavelec,
Rukan Kosak,
Peter Blaha,
Ulrike Diebold
Abstract:
The coarsening of catalytically-active metal clusters is often accelerated by the presence of gases through the formation of mobile intermediates, though the exact mechanism through which this happens is often subject to debate. We use scanning tunneling microscopy (STM) to follow the CO induced coalescence of Pd adatoms supported on the Fe3O4(001) surface at room temperature. We show that highly-…
▽ More
The coarsening of catalytically-active metal clusters is often accelerated by the presence of gases through the formation of mobile intermediates, though the exact mechanism through which this happens is often subject to debate. We use scanning tunneling microscopy (STM) to follow the CO induced coalescence of Pd adatoms supported on the Fe3O4(001) surface at room temperature. We show that highly-mobile Pd-carbonyl species, formed via the so-called skyhook effect, are temporarily trapped at other Pd adatoms. Once these reach a critical density, clusters nucleate; subsequent coarsening occurs through cluster diffusion and coalescence. While CO increases the mobility in the Pd/Fe3O4 system, surface hydroxyls have the opposite effect. Pd atoms transported to surface OH groups are no longer susceptible to the skyhook effect and remain isolated. Following the evolution from well-dispersed metal adatoms into clusters, atom-by-atom, allows identification of the key processes that underlie gas-induced mass transport.
△ Less
Submitted 4 March, 2013;
originally announced March 2013.
-
Strain-induced topological insulator phase transition in HgSe
Authors:
Lars Winterfeld,
Luis A. Agapito,
Jin Li,
Nicholas Kioussis,
Peter Blaha,
Yong P. Chen
Abstract:
Using ab initio electronic structure calculations we investigate the change of the band structure and the nu_0 topological invariant in HgSe (non-centrosymmetric system) under two different type of uniaxial strain along the [001] and [110] directions, respectively. Both compressive [001] and [110] strain leads to the opening of a (crystal field) band gap (with a maximum value of about 37 meV) in t…
▽ More
Using ab initio electronic structure calculations we investigate the change of the band structure and the nu_0 topological invariant in HgSe (non-centrosymmetric system) under two different type of uniaxial strain along the [001] and [110] directions, respectively. Both compressive [001] and [110] strain leads to the opening of a (crystal field) band gap (with a maximum value of about 37 meV) in the vicinity of Gamma, and the concomitant formation of a camel-back- (inverse camel-back-) shape valence (conduction) band along the direction perpendicular to the strain with a minimum (maximum) at Gamma. We find that the Z_2 invariant nu_0=1, which demonstrates conclusively that HgSe is a strong topological insulator (TI). With further increase of the strain the band gap decreases vanishing at a critical strain value (which depends on the strain type) where HgSe undergoes a transition from a strong TI to a trivial (normal) insulator. HgSe exhibits a similar behavior under a tensile [110] uniaxial strain. On the other hand, HgSe remains a normal insulator by applying a [001] tensile uniaxial strain. Complementary electronic structure calculations of the non-polar (110) surface under compressive [110] tensile strain show two Dirac cones at the Gamma point whose spin chiral states are associated with the top and bottom slab surfaces.
△ Less
Submitted 17 February, 2013;
originally announced February 2013.
-
Oxide Heterostructures for Efficient Solar Cells
Authors:
Elias Assmann,
Peter Blaha,
Robert Laskowski,
Karsten Held,
Satoshi Okamoto,
Giorgio Sangiovanni
Abstract:
We propose an unexplored class of absorbing materials for high-efficiency solar cells: heterostructures of transition-metal oxides. In particular, LaVO_3 grown on SrTiO_3 has a direct band gap ~1.1 eV in the optimal range as well as an internal potential gradient, which can greatly help to separate the photo-generated electron-hole pairs. Furthermore, oxide heterostructures afford the flexibility…
▽ More
We propose an unexplored class of absorbing materials for high-efficiency solar cells: heterostructures of transition-metal oxides. In particular, LaVO_3 grown on SrTiO_3 has a direct band gap ~1.1 eV in the optimal range as well as an internal potential gradient, which can greatly help to separate the photo-generated electron-hole pairs. Furthermore, oxide heterostructures afford the flexibility to combine LaVO_3 with other materials such as LaFeO_3 in order to achieve even higher efficiencies with band-gap graded solar cells. We use density-functional theory to demonstrate these features.
△ Less
Submitted 8 April, 2013; v1 submitted 3 January, 2013;
originally announced January 2013.
-
Application of screened hybrid functionals to the bulk transition metals Rh, Pd, and Pt
Authors:
Fabien Tran,
David Koller,
Peter Blaha
Abstract:
We present the results of calculations on bulk transition metals Rh, Pd, and Pt using the screened hybrid functional YS-PBE0 [F. Tran and P. Blaha, Phys. Rev. B \textbf{83}, 235118 (2011)]. The results for the equilibrium geometry are compared with those obtained from (semi)local functionals, namely, the local density approximation and the generalized gradient approximation PBE of Perdew \textit{e…
▽ More
We present the results of calculations on bulk transition metals Rh, Pd, and Pt using the screened hybrid functional YS-PBE0 [F. Tran and P. Blaha, Phys. Rev. B \textbf{83}, 235118 (2011)]. The results for the equilibrium geometry are compared with those obtained from (semi)local functionals, namely, the local density approximation and the generalized gradient approximation PBE of Perdew \textit{et al}. [J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. \textbf{77}, 3865 (1996)]. It is shown that the screened hybrid functional yields more accurate equilibrium geometry than PBE, but, overall, it is not more accurate than LDA. However, in contradiction with experiment, we find that the screened hybrid functional favors a ferromagnetic state as the ground state for all three transition metals. Therefore, the use of hybrid functionals for, e.g., the study of catalytically active systems with correlated oxides on a metal support is questionable.
△ Less
Submitted 27 September, 2012;
originally announced September 2012.
-
Accounting for spin fluctuations beyond LSDA in the density functional theory
Authors:
L. Ortenzi,
I. I. Mazin,
P. Blaha,
L. Boeri
Abstract:
We present a method to correct the magnetic properties of itinerant systems in local spin density approximation (LSDA) and we apply it to the ferromagnetic-paramagnetic transition under pressure in a typical itinerant system, Ni$_{3}$Al. We obtain a scaling of the critical fluctuations as a function of pressure equivalent to the one obtained within Moryia's theory. Moreover we show that in this ma…
▽ More
We present a method to correct the magnetic properties of itinerant systems in local spin density approximation (LSDA) and we apply it to the ferromagnetic-paramagnetic transition under pressure in a typical itinerant system, Ni$_{3}$Al. We obtain a scaling of the critical fluctuations as a function of pressure equivalent to the one obtained within Moryia's theory. Moreover we show that in this material the role of the bandstructure is crucial in driving the transition. Finally we calculate the magnetic moment as a function of pressure, and find that it gives a scaling of the Curie temperature that is in good agreement with the experiment. The method can be easily extended to the antiferromagnetic case and applied, for instance, to the Fe-pnictides in order to correct the LSDA magnetic moment.
△ Less
Submitted 25 May, 2012;
originally announced May 2012.
-
Electronic structure of CrN: A comparison between different exchange correlation potentials
Authors:
A. S. Botana,
F. Tran,
V. Pardo,
D. Baldomir,
P. Blaha
Abstract:
We report a series of electronic structure calculations for CrN using different exchange correlation potentials: PBE, LDA+$U$, the Tran-Blaha modified Becke-Johnson, and hybrid functionals. In every case, our calculations show that the onset of magnetism in CrN should be accompanied by a gap opening. The experimentally found antiferromagnetic order always leads to an insulating behavior. Our resul…
▽ More
We report a series of electronic structure calculations for CrN using different exchange correlation potentials: PBE, LDA+$U$, the Tran-Blaha modified Becke-Johnson, and hybrid functionals. In every case, our calculations show that the onset of magnetism in CrN should be accompanied by a gap opening. The experimentally found antiferromagnetic order always leads to an insulating behavior. Our results give further evidence that the Tran-Blaha functional is very useful for treating the electronic structure of correlated semiconductors allowing a parameter free description of the system. Hybrid functionals are also well capable of describing the electronic structure of CrN. The analysis of the system is complemented with our calculations of the thermopower that are in agreement with the experimental data.
△ Less
Submitted 11 May, 2012;
originally announced May 2012.
-
Three-Dimensional Fermiology by Soft-X-Ray ARPES: Origin of Charge Density Waves in VSe2
Authors:
Vladimir N. Strocov,
Ming Shi,
Masaki Kobayashi,
Claude Monney,
Xiaoqiang Wang,
Juraj Krempasky,
Thorsten Schmitt,
Luc Patthey,
Helmuth Berger,
Peter Blaha
Abstract:
Electronic structure of crystalline materials is their fundamental characteristic which is the basis of almost all their physical and chemical properties. Angle-resolved photoemission spectroscopy (ARPES) is the main experimental tool to study all electronic structure aspects with resolution in k-space. However, its application to three-dimensional (3D) materials suffers from a fundamental problem…
▽ More
Electronic structure of crystalline materials is their fundamental characteristic which is the basis of almost all their physical and chemical properties. Angle-resolved photoemission spectroscopy (ARPES) is the main experimental tool to study all electronic structure aspects with resolution in k-space. However, its application to three-dimensional (3D) materials suffers from a fundamental problem of ill-defined surface-perpendicular wavevector kz. Here, we achieve sharp definition of kz to enable precise navigation in 3D k space by pushing ARPES into the soft-X-ray photon energy range. Essential to break through the notorious problem of small photoexcitation cross-section was an advanced photon flux performance of our instrumentation. We explore the electronic structure of a transition metal dichalcogenide VSe2 which develops charge density waves (CDWs) possessing exotic 3D character. We experimentally identify nesting of its 3D Fermi surface (FS) as the precursor for these CDWs. Our study demonstrates an immense potential of soft-X-ray ARPES (SX-ARPES) to resolve various aspects of 3D electronic structure.
△ Less
Submitted 1 May, 2012;
originally announced May 2012.
-
Ferromagnetic insulating state in tensile-strained LaCoO$_3$ thin films
Authors:
Han Hsu,
Peter Blaha,
Renata M. Wentzcovitch
Abstract:
With local density approximation + Hubbard $U$ (LDA+$U$) calculations, we show that the ferromagnetic (FM) insulating state observed in tensile-strained LaCoO$_3$ epitaxial thin films is most likely a mixture of low-spin (LS) and high-spin (HS) Co, namely, a HS/LS mixture state. Compared with other FM states, including the intermediate-spin (IS) state (\textit{metallic} within LDA+$U$), which cons…
▽ More
With local density approximation + Hubbard $U$ (LDA+$U$) calculations, we show that the ferromagnetic (FM) insulating state observed in tensile-strained LaCoO$_3$ epitaxial thin films is most likely a mixture of low-spin (LS) and high-spin (HS) Co, namely, a HS/LS mixture state. Compared with other FM states, including the intermediate-spin (IS) state (\textit{metallic} within LDA+$U$), which consists of IS Co only, and the insulating IS/LS mixture state, the HS/LS state is the most favorable one. The FM order in HS/LS state is stabilized via the superexchange interactions between adjacent LS and HS Co. We also show that Co spin state can be identified by measuring the electric field gradient (EFG) at Co nucleus via nuclear magnetic resonance (NMR) spectroscopy.
△ Less
Submitted 8 March, 2012;
originally announced March 2012.
-
Implementation of screened hybrid functionals based on the Yukawa potential within the LAPW basis set
Authors:
Fabien Tran,
Peter Blaha
Abstract:
The implementation of screened hybrid functionals into the WIEN2k code, which is based on the LAPW basis set, is reported. The Hartree-Fock exchange energy and potential are screened by means of the Yukawa potential as proposed by Bylander and Kleinman [Phys. Rev. B 41, 7868 (1990)] for the calculation of the electronic structure of solids with the screened-exchange local density approximation. De…
▽ More
The implementation of screened hybrid functionals into the WIEN2k code, which is based on the LAPW basis set, is reported. The Hartree-Fock exchange energy and potential are screened by means of the Yukawa potential as proposed by Bylander and Kleinman [Phys. Rev. B 41, 7868 (1990)] for the calculation of the electronic structure of solids with the screened-exchange local density approximation. Details of the formalism, which is based on the method of Massidda, Posternak, and Baldereschi [Phys. Rev. B 48, 5058 (1993)] for the unscreened Hartree-Fock exchange are given. The results for the transition-energy and structural properties of several test cases are presented. Results of calculations of the Cu electric-field gradient in Cu2O are also presented, and it is shown that the hybrid functionals are much more accurate than the standard local-density or generalized gradient approximations.
△ Less
Submitted 3 May, 2011; v1 submitted 23 March, 2011;
originally announced March 2011.
-
Spin-state crossover and hyperfine interactions of ferric iron in MgSiO$_3$ perovskite
Authors:
Han Hsu,
Peter Blaha,
Matteo Cococcioni,
Renata M. Wentzcovitch
Abstract:
Using density functional theory plus Hubbard $U$ calculations, we show that the ground state of (Mg,Fe)(Si,Fe)O$_3$ perovskite, a major mineral phase in the Earth's lower mantle, has high-spin ferric iron ($S=5/2$) at both the dodecahedral (A) and octahedral (B) site. As the pressure increases, the B-site iron undergoes a spin-state crossover to the low-spin state ($S=1/2$), while the A-site iron…
▽ More
Using density functional theory plus Hubbard $U$ calculations, we show that the ground state of (Mg,Fe)(Si,Fe)O$_3$ perovskite, a major mineral phase in the Earth's lower mantle, has high-spin ferric iron ($S=5/2$) at both the dodecahedral (A) and octahedral (B) site. As the pressure increases, the B-site iron undergoes a spin-state crossover to the low-spin state ($S=1/2$), while the A-site iron remains in the high-spin state. Our calculation shows that the B-site spin-state crossover in the pressure range of 40-70 GPa is accompanied by a noticeable volume reduction and an increase in quadrupole splitting, consistent with recent X-ray diffraction and Mössbauer spectroscopy measurements. The volume reduction leads to a significant softening in the bulk modulus, which suggests a possible source of seismic velocity anomalies in the lower mantle.
△ Less
Submitted 1 February, 2011; v1 submitted 19 January, 2011;
originally announced January 2011.
-
Systematic investigation of a family of gradient-dependent functionals for solids
Authors:
Philipp Haas,
Fabien Tran,
Peter Blaha,
Luana S. Pedroza,
Antonio J. R. da Silva,
Mariana M. Odashima,
Klaus Capelle
Abstract:
Eleven density functionals are compared with regard to their performance for the lattice constants of solids. We consider standard functionals, such as the local-density approximation and the Perdew-Burke-Ernzerhof (PBE) generalized-gradient approximation (GGA), as well as variations of PBE GGA, such as PBEsol and similar functionals, PBE-type functionals employing a tighter Lieb-Oxford bound, and…
▽ More
Eleven density functionals are compared with regard to their performance for the lattice constants of solids. We consider standard functionals, such as the local-density approximation and the Perdew-Burke-Ernzerhof (PBE) generalized-gradient approximation (GGA), as well as variations of PBE GGA, such as PBEsol and similar functionals, PBE-type functionals employing a tighter Lieb-Oxford bound, and combinations thereof. Several of these variations are proposed here for the first time. On a test set of 60 solids we perform a system-by-system analysis for selected functionals and a full statistical analysis for all of them. The impact of restoring the gradient expansion and of tightening the Lieb-Oxford bound is discussed, and confronted with previous results obtained from other codes, functionals or test sets. No functional is uniformly good for all investigated systems, but surprisingly, and pleasingly, the simplest possible modifications to PBE turn out to have the most beneficial effect on its performance. The atomization energy of molecules was also considered and on a testing set of six molecules, we found that the PBE functional is clearly the best, the others leading to strong overbinding.
△ Less
Submitted 12 March, 2010;
originally announced March 2010.