Physical Review B
covering condensed matter and materials physics
- Highlights
- Recent
- Accepted
- Collections
- Authors
- Referees
- Search
- Press
- About
- Editorial Team
Hubbard parameters for transition metals from first principles
Rebekka Tesch and Piotr M. Kowalski
Phys. Rev. B 105, 195153 – Published 31 May 2022
- Article
- References
- Citing Articles (27)
- Supplemental Material
PDFHTMLExport Citation
Abstract
Using the linear response-based constrained local density approximation (cLDA) approach we systematically computed the Hubbard parameters for series of and transition metals. We compare the results with estimations by the constrained random phase approximation (cRPA) method and discuss the performance of the self-consistent density functional theory + () method for prediction of lattice parameters, work functions, -bandwidths and -band centers. Interestingly, we found that blindly applied the standard, fully localized limit (FLL) version of the approach heavily overestimates the positions of -band centers with respect to the Fermi level, but much better agreement with experiment is obtained when applying a more realistic, Wannier-type representation of orbitals for projection of states occupancies. We present another, independent estimate of the Hubbard parameter based on the comparison of Hartree-Fock and DFT eigenvalues, and positions of -band centers. The so-derived estimates are surprisingly well consistent with the ones derived from the above-mentioned first principles approaches, and allow for validation of cRPA or cLDA results for the disputed cases, including Cu, Ag, and Au for which large parameters are obtained from the cLDA method.
- Received 20 July 2021
- Revised 6 May 2022
- Accepted 9 May 2022
DOI:https://doi.org/10.1103/PhysRevB.105.195153
©2022 American Physical Society
Physics Subject Headings (PhySH)
- Physical Systems
Strongly correlated systemsTransition metals
- Techniques
Ab initio calculationsElectron-correlation calculationsRandom phase approximation
Condensed Matter, Materials & Applied Physics
Authors & Affiliations
- Institute of Energy and Climate Research, Theory and Computation of Energy Materials (IEK-13), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany; Chair of Theory and Computation of Energy Materials, Faculty of Georesources and Materials Engineering, RWTH Aachen University, 52062 Aachen, Germany; and Jülich Aachen Research Alliance, JARA-CSD and JARA-ENERGY, 52425 Jülich, Germany
- Institute of Energy and Climate Research, Theory and Computation of Energy Materials (IEK-13), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany and Jülich Aachen Research Alliance, JARA-CSD and JARA-ENERGY, 52425 Jülich, Germany
- *p.kowalski@fz-juelich.de
Article Text (Subscription Required)
Click to Expand
Supplemental Material (Subscription Required)
Click to Expand
References (Subscription Required)
Click to Expand
Issue
Vol. 105, Iss. 19 — 15 May 2022
Access Options
- Buy Article »
- Log in with individual APS Journal Account »
- Log in with a username/password provided by your institution »
- Get access through a U.S. public or high school library »
![Hubbard $U$ parameters for transition metals from first principles (11) Hubbard $U$ parameters for transition metals from first principles (11)](https://i0.wp.com/cdn.journals.aps.org/development/journals/images/author-services-placard.png)
Images
Figure 1
The Hubbard parameters for , and transition metals. The different symbols denote the results obtained with: the cLDA method (red filled and open squares for nonmagnetic and magnetic metals, respectively), the -band center shifts between HF and DFT methods (green circles for nonmagnetic metals), the cRPA method by Şaşıoğlu etal. [15] (blue filled and open diamonds for nonmagnetic and magnetic metals, respectively), the cLDA method by Nakamura etal. [26] (open circles). The black symbols represent the experimental values (the measured correlation energy deduced from Auger and XPS spectroscopy) of de Boer etal. [53], Sawatzky and Post [54], Antonides etal. [55] (stars) and Kaurila etal. [56] (crosses, with uncertainty at the level of ). We note, however, that an exact correspondence between the measured and computed values is not expected and we provide these values only for qualitative comparison. All the reported data are also provided in the ESI, Table S2 [57].
Figure 2
The projected orbitals density of states for fcc Cu metal: (a)computed with the standard DFT (), hybrid functional (PBE0) and Hartree-Fock methods, as well as measured with XPS by Ref.[46], and (b)computed with and different projections of the states occupancy: atomic orbitals (AO), Wannier functions (WF). Gaussian smearing of has been used to match the experimental band broadening. The XPS data are scaled vertically to match the intensity of the computed -band.
Figure 3
The -bandwidths (derived considering: only the occupied part of the -band (left panels) and the entire -bandwidth (right panels)), derived as described in the text. Magnetic states of Cr, Mn, Fe, Co, and Ni metals are considered. Computed reference data (ref. calc.) are taken from Şaşıoğlu etal. [15]. Experimental reference data are taken from Hüfner etal. [66] for Ni, Cu and Ag and from Smith etal. [82] for Rh, Pd, Ir, Pt, and Au. AO and WF indicate calculations performed with atomic orbitals and Wannier functions as projectors, respectively. Uncertainty of the measured values is in the order of , and arises mainly due to unclear definitions of the band limits. All the reported data are also provided in the ESI, Table S4 [57].
Figure 4
The ratio computed with the cLDA and cRPA [15] Hubbard parameter values. The bandwidths are those from standard DFT calculations (), to allow a straightforward comparison to the cRPA reference data [15] computed in such way. These were derived considering: only the occupied part of the -band (left panels) and the entire -bandwidth (right panels). Magnetic states of Cr, Mn, Fe, Co, and Ni metals are considered.
Figure 5
The position of the -band center (dbc) with respect to the Fermi level, derived considering: only the occupied part of the -band (left panels) and the entire -bandwidth (right panels). Magnetic states of Cr, Mn, Fe, Co, and Ni metals are considered. Experimental XPS data (with maximal uncertainty of ), calculated reference data (standard (ref. calc.) and refined (ref. calc. refi., see the text for explanation)) are taken from Hofmann etal. [46] for Fe, Co, Ni, Cu, Pd, Ag, Pt, and Au and from Smith etal. [82] for Rh and Ir metals. Refined DFT calculations of Hofmann etal. [46] are given by blue stars. More details are provided in the ESI, Table S5 [57].
Figure 6
The projected density of states for fcc Cu computed with the (cLDA) method (atomic orbitals as projectors, ), with Gaussian smearing of . The results show a clear hybrid state between the main -band and the Fermi level. XPS data are taken from [46] and scaled vertically to match the intensity of the computed -band.
Figure 7
The projected orbitals density of states for selected metals, computed with different projections of states occupancy: atomic orbitals (AO), Wannier functions (WF). Gaussian smearing of has been applied to match the experimental band broadening. Magnetic and non-magnetic states of Fe, Cr, and Ni metals are considered. XPS data are taken from [46] and are rescaled vertically to match the intensity of the computed -bands. The reference experimental spectra agree well with earlier measurements by Höchst etal. [101], Hüfner etal. [66], Smith etal. [82]. Our computed DOS for (DFT approach) agree well with the calculations by Hofmann etal. [46].