Source author record

Christoph Friedrich

Christoph Friedrich appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

17 published item(s)

preprint2022arXiv

Bulk and surface electronic structure of Bi$_4$Te$_3$ from $GW$ calculations and photoemission experiments

We present a combined theoretical and experimental study of the electronic structure of stoichiometric Bi$_4$Te$_3$, a natural superlattice of alternating Bi$_2$Te$_3$ quintuple layers and Bi bilayers. In contrast to the related semiconducting compounds Bi$_2$Te$_3$ and Bi$_1$Te$_1$, density functional theory predicts Bi$_4$Te$_3$ to be a semimetal. In this work, we compute the quasiparticle electronic structure of Bi$_4$Te$_3$ in the framework of the $GW$ approximation within many-body perturbation theory. The quasiparticle corrections are found to modify the dispersion of the valence and conduction bands in the vicinity of the Fermi energy, leading to the opening of a small indirect band gap. Based on the analysis of the eigenstates, Bi$_4$Te$_3$ is classified as a dual topological insulator with bulk topological invariants $\mathbb{Z}_2$ (1;111) and magnetic mirror Chern number $n_M=1$. The bulk $GW$ results are used to build a Wannier-functions based tight-binding Hamiltonian that is further applied to study the electronic properties of the (111) surface. The comparison with our angle-resolved photoemission measurements shows excellent agreement between the computed and measured surface states and indicates the dual topological nature of Bi$_4$Te$_3$.

preprint2021arXiv

$Z_2$ topology of bismuth

While first-principles calculations with different levels of sophistication predict a topologically trivial $Z_2$ state for bulk bismuth, some photoemission experiments show surface states consistent with the interpretation of bismuth being in a topologically non-trivial $Z_2$ state. We resolve this contradiction between theory and experiment by showing, based on quasiparticle self-consistent $GW$ calculations, that the experimental surface states interpreted as supporting a non-trivial phase are actually consistent with a trivial $Z_2$ invariant. We identify this contradiction as the result of a crosstalk effect arising from the extreme penetration depth of the surface states into the bulk of Bi. A film of Bi can be considered bulk-like only for thicknesses of about 1000 bilayers ($\approx$ 400 nm) and more.

preprint2015arXiv

Electronic phase transitions of bismuth under strain from relativistic self-consistent GW calculations

We present quasiparticle self-consistent GW (QSGW) calculations of semimetallic bulk Bi. We go beyond the conventional QSGW method by including the spin-orbit coupling throughout the self-consistency cycle. This approach improves the description of the electron and the hole pockets considerably with respect to standard density functional theory (DFT), leading to excellent agreement with experiment. We employ this relativistic QSGW approach to conduct a study of the semimetal-to-semiconductor and the trivial-to-topological transitions that Bi experiences under strain. DFT predicts that an unphysically large strain is needed for such transitions. We show, by means of the relativistic QSGW description of the electronic structure, that an in-plane tensile strain of only 0.3% and a compressive strain of 0.4% are sufficient to cause the semimetal-to-semiconductor and the trivial-to-topological phase transitions, respectively. Thus, the required strain moves into a regime that is likely to be realizable in experiment, which opens up the possibility to explore bulklike topological behavior of pure Bi.

preprint2013arXiv

Ab-initio calculation of the effective on-site Coulomb interaction parameters for half-metallic magnets

Correlation effects play an important role in the electronic structure of half-metallic (HM) magnets. In particular, they give rise to non-quasiparticle states above (or below) the Fermi energy at finite temperatures that reduce the spin polarization and, as a consequence, the efficiency of spintronics devices. Employing the constrained random-phase approximation (cRPA) within the full-potential linearized augmented-plane-wave (FLAPW) method using maximally localized Wannier functions, we calculate the strength of the effective on-site Coulomb interaction (Hubbard $U$ and Hund exchange $J$) between localized electrons in different classes of HM magnets considering: (i) \emph{sp}-electron ferromagnets in rock-salt structure, (ii) zincblende 3\emph{d} binary ferromagnets, as well as (iii) ferromagnetic and ferrimagnetic semi- and full-Heusler compounds.

preprint2013arXiv

Elimination of the linearization error and improved basis-set convergence within the FLAPW method

We analyze in detail the error that arises from the linearization in linearized augmented-plane-wave (LAPW) basis functions around predetermined energies $E_l$ and show that it can lead to undesirable dependences of the calculated results on method-inherent parameters such as energy parameters $E_l$ and muffin-tin sphere radii. To overcome these dependences, we evaluate approaches that eliminate the linearization error systematically by adding local orbitals (LOs) to the basis set. We consider two kinds of LOs: (i) constructed from solutions $u_l(r,E)$ to the scalar-relativistic approximation of the radial Dirac equation with $E>E_l$ and (ii) constructed from second energy derivatives $\partial^2 u_l(r,E) / \partial E^2$ at $E=E_l$. We find that the latter eliminates the error most efficiently and yields the density functional answer to many electronic and materials properties with very high precision. Finally, we demonstrate that the so constructed LAPW+LO basis shows a more favorable convergence behavior than the conventional LAPW basis due to a better decoupling of muffin-tin and interstitial regions, similarly to the related APW+lo approach, which requires an extra set of LOs to reach the same total energy, though.

preprint2013arXiv

Strong magnon softening in tetragonal FeCo compounds

Magnons play an important role in fast precessional magnetization reversal processes serving as a heat bath for dissipation of the Zeeman energy and thus being responsible for the relaxation of magnetization. Employing \emph{ab initio} many-body perturbation theory we studied the magnon spectra of the tetragonal FeCo compounds considering three different experimental $c/a$ ratios, $c/a=$1.13, 1.18, and 1.24 corresponding to FeCo grown on Pd, Ir, and Rh, respectively. We find that for all three cases the short-wave-length magnons are strongly damped and tetragonal distortion gives rise to a significant magnon softening. The magnon stiffness constant $D$ decreases almost by a factor of two from FeCo/Pd to FeCo/Rh. The combination of soft magnons together with the giant magnetic anisotropy energy suggests FeCo/Rh to be a promising material for perpendicular magnetic recording applications.

preprint2012arXiv

GW study of the half-metallic Heusler compounds Co2MnSi and Co2FeSi

Quasiparticle spectra of potentially half-metallic Co2MnSi and Co2FeSi Heusler compounds have been calculated within the one-shot GW approximation in an all-electron framework without adjustable parameters. For Co2FeSi the many-body corrections are crucial: a pseudogap opens and good agreement of the magnetic moment with experiment is obtained. Otherwise, however, the changes with respect to the density-functional-theory starting point are moderate. For both cases we find that photoemission and x-ray absorption spectra are well described by the calculations. By comparison with the GW density of states, we conclude that the Kohn-Sham eigenvalue spectrum provides a reasonable approximation for the quasiparticle spectrum of the Heusler compounds considered in this work.

preprint2012arXiv

Precise response functions in all-electron methods: Application to the optimized-effective-potential approach

The optimized-effective-potential (OEP) method is a special technique to construct local Kohn-Sham potentials from general orbital-dependent energy functionals. In a recent publication [M. Betzinger, C. Friedrich, S. Blügel, A. Görling, Phys. Rev. B 83, 045105 (2011)] we showed that uneconomically large basis sets were required to obtain a smooth local potential without spurious oscillations within the full-potential linearized augmented-plane-wave method (FLAPW). This could be attributed to the slow convergence behavior of the density response function. In this paper, we derive an incomplete-basis-set correction for the response, which consists of two terms: (1) a correction that is formally similar to the Pulay correction in atomic-force calculations and (2) a numerically more important basis response term originating from the potential dependence of the basis functions. The basis response term is constructed from the solutions of radial Sternheimer equations in the muffin-tin spheres. With these corrections the local potential converges at much smaller basis sets, at much fewer states, and its construction becomes numerically very stable. We analyze the improvements for rock-salt ScN and report results for BN, AlN, and GaN, as well as the perovskites CaTiO3, SrTiO3, and BaTiO3. The incomplete-basis-set correction can be applied to other electronic-structure methods with potential-dependent basis sets and opens the perspective to investigate a broad spectrum of problems in theoretical solid-state physics that involve response functions.

preprint2012arXiv

Strength of the Effective Coulomb Interaction at Metal and Insulator Surfaces

The effective on-site Coulomb interaction (Hubbard $U$) between localized electrons at crystal surfaces is expected to be enhanced due to the reduced coordination number and reduced subsequent screening. By means of first principles calculations employing the constrained random-phase approximation (cRPA) we show that this is indeed the case for simple metals and insulators but not necessarily for transition metals and insulators that exhibit pronounced surface states. In the latter case, the screening contribution from surface states as well as the influence of the band narrowing increases the electron polarization to such an extent as to overcompensate the decrease resulting from the reduced effective screening volume. The Hubbard $U$ parameter is thus substantially reduced in some cases, e.g., by around 30% for the (100) surface of bcc Cr.

preprint2011arXiv

Effective Coulomb interaction in transition metals from constrained random-phase approximation

The effective on-site Coulomb interaction (Hubbard $U$) between localized \textit{d} electrons in 3\textit{d}, 4\textit{d}, and 5\textit{d} transition metals is calculated employing a new parameter-free realization of the constrained random-phase approximation using Wannier functions within the full-potential linearized augmented-plane-wave method. The $U$ values lie between 1.5 and 5.7 eV and depend on the crystal structure, spin polarization, \textit{d} electron number, and \textit{d} orbital filling. On the basis of the calculated $U$ parameters, we discuss the strength of the electronic correlations and the instability of the paramagnetic state towards the ferromagnetic one for 3\textit{d} metals.

preprint2011arXiv

First-principles calculation of electronic excitations in solids with SPEX

We describe the software package SPEX, which allows first-principles calculations of quasiparticle and collective electronic excitations in solids using techniques from many-body perturbation theory. The implementation is based on the full-potential linearized augmented-plane-wave (FLAPW) method, which treats core and valence electrons on an equal footing and can be applied to a wide range of materials, including transition metals and rare earths. After a discussion of essential features that contribute to the high numerical efficiency of the code, we present illustrative results for quasiparticle band structures calculated within the GW approximation for the electronic self-energy, electron-energy-loss spectra with inter- and intraband transitions as well as local-field effects, and spin-wave spectra of itinerant ferromagnets. In all cases the inclusion of many-body correlation terms leads to very good quantitative agreement with experimental spectroscopies.

preprint2011arXiv

Local exact exchange potentials within the all-electron FLAPW method and a comparison with pseudopotential results

We present a general numerical approach to construct local Kohn-Sham potentials from orbital-dependent functionals within the all-electron full-potential linearized augmented-plane-wave (FLAPW) method, in which core and valence electrons are treated on an equal footing. As a practical example, we present a treatment of the orbital-dependent exact-exchange (EXX) energy and potential. A formulation in terms of a mixed product basis, which is constructed from products of LAPW basis functions, enables a solution of the optimized-effective-potential (OEP) equation with standard numerical algebraic tools and without shape approximations for the resulting potential. We find that the mixed product and LAPW basis sets must be properly balanced to obtain smooth and converged EXX potentials without spurious oscillations. The construction and convergence of the exchange potential is analyzed in detail for diamond. Our all-electron results for C, Si, SiC, Ge, GaAs semiconductors as well as Ne and Ar noble-gas solids are in very favorable agreement with plane-wave pseudopotential calculations. This confirms the adequacy of the pseudopotential approximation in the context of the EXX-OEP formalism and clarifies a previous contradiction between FLAPW and pseudopotential results.

preprint2011arXiv

The HSE hybrid functional within the FLAPW method and its application to GdN

We present an implementation of the Heyd-Scuseria-Ernzerhof (HSE) hybrid functional within the full-potential linearized augmented-plane-wave (FLAPW) method. Pivotal to the HSE functional is the screened electron-electron interaction, which we separate into the bare Coulomb interaction and the remainder, a slowly varying function in real space. Both terms give rise to exchange potentials, which sum up to the screened nonlocal exchange potential of HSE. We evaluate the former with the help of an auxiliary basis, defined in such a way that the bare Coulomb matrix becomes sparse. The latter is computed in reciprocal space, exploiting its fast convergence behavior in reciprocal space. This approach is general and can be applied to a whole class of screened hybrid functionals. We obtain excellent agreement of band gaps and lattice constants for prototypical semiconductors and insulators with electronic-structure calculations using plane-wave or Gaussian basis sets. We apply the HSE hybrid functional to examine the ground-state properties of rocksalt GdN, which have been controversially discussed in literature. Our results indicate that there is a half-metal to insulator transition occurring between the theoretically optimized lattice constant at 0 K and the experimental lattice constant at room temperature. Overall, we attain good agreement with experimental data for band transitions, magnetic moments, and the Curie temperature.

preprint2010arXiv

Efficient implementation of the GW approximation within the all-electron FLAPW method

We present an implementation of the GW approximation for the electronic self-energy within the full-potential linearized augmented-plane-wave (FLAPW) method. The algorithm uses an all-electron mixed product basis for the representation of response matrices and related quantities. This basis is derived from the FLAPW basis and is exact for wave-function products. The singularity of the bare and screened interaction potentials gives rise to a numerically important self-energy contribution, which we treat analytically to achieve good convergence with respect to the k-point sampling. As numerical realizations of the GW approximation typically suffer from the high computational expense required for the evaluation of the nonlocal and frequency-dependent self-energy, we demonstrate how the algorithm can be made very efficient by exploiting spatial and time-reversal symmetry as well as by applying an optimization of the mixed product basis that retains only the numerically important contributions of the electron-electron interaction. Furthermore, we demonstrate that one can employ an extrapolar approximation for high-lying states to reduce the number of empty states that must be taken into account explicitly in the construction of the polarization function and the self-energy. We show convergence tests, CPU timings, and results for prototype semiconductors and insulators as well as ferromagnetic nickel.

preprint2010arXiv

Hybrid functionals within the all-electron FLAPW method: implementation and applications of PBE0

We present an efficient implementation of the PBE0 hybrid functional within the full-potential linearized augmented-plane-wave (FLAPW) method. The Hartree-Fock exchange term, which is a central ingredient of hybrid functionals, gives rise to a computationally expensive nonlocal potential in the one-particle Schroedinger equation. The matrix elements of this exchange potential are calculated with the help of an auxiliary basis that is constructed from products of FLAPW basis functions. By representing the Coulomb interaction in this basis the nonlocal exchange term becomes a Brillouin-zone (BZ) sum over vector-matrix-vector products. We show that the Coulomb matrix can be made sparse by a suitable unitary transformation of the auxiliary basis, which accelerates the computation of the vector-matrix-vector products considerably. Additionally, we exploit spatial and time-reversal symmetry to identify the nonvanishing exchange matrix elements in advance and to restrict the k summations for the nonlocal potential to an irreducible set of k points. Favorable convergence of the self-consistent-field cycle is achieved by a nested density-only and density-matrix iteration scheme. We discuss the convergence with respect to the parameters of our numerical scheme and show results for a variety of semiconductors and insulators, including oxide materials, where the PBE0 hybrid functional improves the band gaps and the description of localized states in comparison with the PBE functional. Furthermore, we find that in contrast to conventional local exchange-correlation functionals ferromagnetic EuO is correctly predicted to be a semiconductor.

preprint2010arXiv

Wannier-function approach to spin excitations in solids

We present a computational scheme to study spin excitations in magnetic materials from first principles. The central quantity is the transverse spin susceptibility, from which the complete excitation spectrum, including single-particle spin-flip Stoner excitations and collective spin-wave modes, can be obtained. The susceptibility is derived from many-body perturbation theory and includes dynamic correlation through a summation over ladder diagrams that describe the coupling of electrons and holes with opposite spins. In contrast to earlier studies, we do not use a model potential with adjustable parameters for the electron-hole interaction but employ the random-phase approximation. To reduce the numerical cost for the calculation of the four-point scattering matrix we perform a projection onto maximally localized Wannier functions, which allows us to truncate the matrix efficiently by exploiting the short spatial range of electronic correlation in the partially filled d or f orbitals. Our implementation is based on the FLAPW method. Starting from a ground-state calculation within the LSDA, we first analyze the matrix elements of the screened Coulomb potential in the Wannier basis for the 3d transition-metal series. In particular, we discuss the differences between a constrained nonmagnetic and a proper spin-polarized treatment for the ferromagnets Fe, Co, and Ni. The spectrum of single-particle and collective spin excitations in fcc Ni is then studied in detail. The calculated spin-wave dispersion is in good overall agreement with experimental data and contains both an acoustic and an optical branch for intermediate wave vectors along the [100] direction. In addition, we find evidence for a similar double-peak structure in the spectral function along the [111] direction.

preprint2008arXiv

Efficient calculation of the Coulomb matrix and its expansion around k=0 within the FLAPW method

We derive formulas for the Coulomb matrix within the full-potential linearized augmented-plane-wave (FLAPW) method. The Coulomb matrix is a central ingredient in implementations of many-body perturbation theory, such as the Hartree-Fock and GW approximations for the electronic self-energy or the random-phase approximation for the dielectric function. It is represented in the mixed product basis, which combines numerical muffin-tin functions and interstitial plane waves that are here expanded with the Rayleigh formula. The resulting algorithm is very efficient in terms of both computational cost and accuracy and is superior to an implementation with the Fourier transform of the step function. In order to allow an analytic treatment of the divergence at k=0 in reciprocal space, we expand the Coulomb matrix analytically around this point without resorting to a projection onto plane waves. We then apply a basis transformation that diagonalizes the Coulomb matrix and confines the divergence to a single eigenvalue. At the same time, response matrices like the dielectric function separate into head, wings, and body with the same mathematical properties as in a plane-wave basis. As an illustration we apply the formulas to electron-energy-loss spectra for nickel at different k vectors including k=0. The convergence of the spectra towards the result at k=0 is clearly seen. Our treatment also allows to include transitions from core states that give rise to a shallow peak at high energies and lead to good agreement with experiment.