Researcher profile

Peter Blaha

Peter Blaha contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

Bandgap of two-dimensional materials: Thorough assessment of modern exchange-correlation functionals

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.

preprint2022arXiv

Implementation of self-consistent MGGA functionals in augmented plane wave based methods

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.

preprint2022arXiv

Length-gauge optical matrix elements in WIEN2k

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.

preprint2021arXiv

Density analysis for estimating the degree of on-site correlation on transition-metal atoms in extended systems

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.

preprint2020arXiv

Perturbation approach to ab initio effective mass calculations

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.

preprint2020arXiv

Shortcomings of meta-GGA functionals when describing magnetism

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.

preprint2013arXiv

High-pressure cupric oxide: a room-temperature multiferroic

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.

preprint2013arXiv

Hybrid functionals for solids with an optimized Hartree-Fock mixing parameter

(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.

preprint2013arXiv

Oxide Heterostructures for Efficient Solar Cells

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.

preprint2012arXiv

Application of screened hybrid functionals to the bulk transition metals Rh, Pd, and Pt

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.

preprint2012arXiv

Ferromagnetic insulating state in tensile-strained LaCoO$_3$ thin films

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.

preprint2012arXiv

Three-Dimensional Fermiology by Soft-X-Ray ARPES: Origin of Charge Density Waves in VSe2

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.

preprint2011arXiv

Implementation of screened hybrid functionals based on the Yukawa potential within the LAPW basis set

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.

preprint2010arXiv

Systematic investigation of a family of gradient-dependent functionals for solids

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.