BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Georg Kresse

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

ResearcherUnclaimed source record
cond-mat.mtrl-scicond-mat.str-elphysics.chem-phphysics.comp-phcond-mat.stat-mechcond-mat.supr-concond-mat.softmath.FAmath.NANumerical Analysisphysics.geo-ph

Catalog footprint

What is connected

31works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

31 published item(s)

preprint2025arXiv

Magnetism-Enhanced Strong Electron-Phonon Coupling in Infinite-Layer Nickelate

Intriguing analogies between the nickelates and the cuprates provide a promising avenue for unraveling the microscopic mechanisms underlying high-$T_c$ superconductivity. While electron correlation effects in the nickelates have been extensively studied, the role of electron-phonon coupling (EPC) remains highly controversial. Here, by taking pristine LaNiO$_2$ as an exemplar nickelate, we present an in-depth study of EPC for both the non-magnetic (NM) and the $C$-type antiferromagnetic ($C$-AFM) phase using advanced density functional theory methods without invoking $U$ or other free parameters. The weak EPC strength $λ$ in the NM phase is found to be greatly enhanced ($\sim$4$\times$) due to the presence of magnetism in the $C$-AFM phase. This enhancement arises from strong interactions between the flat bands associated with the Ni-3$d_{z^2}$ orbitals and the low-frequency phonon modes driven by the vibrations of Ni and La atoms. The resulting phonon softening is shown to yield a distinctive kink in the electronic structure around 15 meV, which would provide an experimentally testable signature of our predictions. Our study highlights the critical role of local magnetic moments and interply EPC in the nickelate.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Accurate Electron-phonon Interactions from Advanced Density Functional Theory

Electron-phonon coupling (EPC) is key for understanding many properties of materials such as superconductivity and electric resistivity. Although first principles density-functional-theory (DFT) based EPC calculations are used widely, their efficacy is limited by the accuracy and efficiency of the underlying exchange-correlation functionals. These limitations become exacerbated in complex $d$- and $f$-electron materials, where beyond-DFT approaches and empirical corrections, such as the Hubbard $U$, are commonly invoked. Here, using the examples of CoO and NiO, we show how the efficient r2scan density functional correctly captures strong EPC effects in transition-metal oxides without requiring the introduction of empirical parameters. We also demonstrate the ability of r2scan to accurately model phonon-mediated superconducting properties of the main group compounds (e.g., MgB$_2$), with improved electronic bands and phonon dispersions over those of traditional density functionals. Our study provides a pathway for extending the scope of accurate first principles modeling of electron-phonon interactions to encompass complex $d$-electron materials.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Core-hole excitations using the projector augmented-wave method and the Bethe-Salpeter equation

We present an implementation of the Bethe-Salpeter equation (BSE) for core-conduction band pairs within the framework of the projector augmented-wave method. For validation, the method is applied to the $K$-edges of diamond, graphite, hexagonal boron-nitride, as well as four lithium-halides (LiF, LiCl, LiI, LiBr). We compare our results with experiment, previous theoretical BSE results, and the density functional theory-based supercell core-hole method. In all considered cases, the agreement with experiment is excellent, in particular for the position of the peaks as well as the fine structure. Comparing BSE to supercell core-hole spectra we find that the latter often qualitatively reproduces the experimental spectrum, however, it sometimes lacks important details. This is shown for the $K$-edges of diamond and nitrogen in hexagonal boron-nitride, where we are capable to resolve within the BSE experimental features that are lacking in the core-hole method. Additionally, we show that in certain systems the supercell core-hole method performs better if the excited electron is added to the background charge. We attribute this improved performance to a reduced self-interaction.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Phase Transitions of Zirconia: Machine-Learned Force Fields Beyond Density Functional Theory

We present an approach to generate machine-learned force fields (MLFF) with beyond density functional theory (DFT) accuracy. Our approach combines on-the-fly active learning and $Δ$-machine learning in order to generate an MLFF for zirconia based on the random phase approximation (RPA). Specifically, an MLFF trained on-the-fly during DFT based molecular dynamics simulations is corrected by another MLFF that is trained on the differences between RPA and DFT calculated energies, forces and stress tensors. Thanks to the relatively smooth nature of the differences, the expensive RPA calculations are performed only on a small number of representative structures of small unit cells. These structures are determined by a singular value decomposition rank compression of the kernel matrix with low spatial resolution. This dramatically reduces the computational cost and allows us to generate an MLFF fully capable of reproducing high-level quantum-mechanical calculations beyond DFT. We carefully validate our approach and demonstrate its success in studying the phase transitions of zirconia.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Accurate optical spectra through time-dependent density functional theory based on screening-dependent hybrid functionals

We investigate optical absorption spectra obtained through time-dependent density functional theory (TD-DFT) based on nonempirical hybrid functionals that are designed to correctly reproduce the dielectric function. The comparison with state-of-the-art $GW$ calculations followed by the solution of the Bethe-Sapeter equation (BSE-$GW$) shows close agreement for both the transition energies and the main features of the spectra. We confront TD-DFT with BSE-$GW$ by focusing on the model dielectric function and the local exchange-correlation kernel. The present TD-DFT approach achieves the accuracy of BSE-$GW$ at a fraction of the computational cost.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Comparative ab initio study of the structural, electronic, magnetic, and dynamical properties of LiOsO$_3$ and NaOsO$_3$

Despite similar chemical compositions, LiOsO$_3$ and NaOsO$_3$ exhibit remarkably distinct structural, electronic, magnetic, and spectroscopic properties. At low temperature, LiOsO$_3$ is a polar bad metal with a rhombohedral $R3c$ structure without the presence of long-range magnetic order, whereas NaOsO$_3$ is a $G$-type antiferromagnetic insulator with an orthorhombic $Pnma$ structure. By means of comparative first-principles DFT+$U$ calculations with the inclusion of the spin-orbit coupling, we ($i$) identify the origin of the different structural ($R3c$ vs. $Pnma$) properties using a symmetry-adapted soft mode analysis, ($ii$) provide evidence that all considered exchange-correlation functionals (LDA, PBE, PBEsol, SCAN, and HSE06) and the spin disordered polymorphous descriptions are unsatisfactory to accurately describe the electronic and magnetic properties of both systems simultaneously, and ($iii$) clarify that the distinct electronic (metallic vs. insulating) properties originates mainly from a cooperative steric and magnetic effect. Finally, we find that although at ambient pressure LiOsO$_3$ with a $Pnma$ symmetry and NaOsO$_3$ with a $R\bar{3}c$ symmetry are energetically unfavorable, they do not show soft phonons and therefore are dynamically stable. A pressure-induced structural phase transition from $R3c$ to $Pnma$ for LiOsO$_3$ is predicted, whereas for NaOsO$_3$ no symmetry change is discerned in the considered pressure range.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Electron-Phonon Interactions Using the PAW Method and Wannier Functions

We present an ab-initio density-functional-theory approach for calculating electron-phonon interactions within the projector augmented-wave method. The required electron-phonon matrix elements are defined as the second derivative of the one-electron energies in the PAW method. As the PAW method leads to a generalized eigenvalue problem, the resulting electron-phonon matrix elements lack some symmetries that are usually present for simple eigenvalue problems and all-electron formulations. We discuss the relation between our definition of the electron-phonon matrix element and other formulations. To allow for efficient evaluation of physical properties, we introduce a Wannier-interpolation scheme, again adapted to generalized eigenvalue problems. To explore the method's numerical characteristics, the temperature-dependent band-gap renormalization of diamond is calculated and compared with previous publications. Furthermore, we apply the method to selected binary compounds and show that the obtained zero-point renormalizations agree well with other values found in literature and experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Local embedding of Coupled Cluster theory into the Random Phase Approximation using plane-waves

We present an embedding approach to treat local electron correlation effects in periodic environments. In a single, consistent framework, our plane-wave based scheme embeds a local high-level correlation calculation (here Coupled Cluster Theory, CC), employing localized orbitals, into a low-level correlation calculation (here the direct Random Phase Approximation, RPA). This choice allows for an accurate and effcient treatment of long-range dispersion effects. Accelerated convergence with respect to the local fragment size can be observed if the low-level and high-level long-range dispersion are quantitatively similar, as is the case for CC in RPA. To demonstrate the capabilities of the introduced embedding approach, we calculate adsorption energies of molecules on a surface and in a chabazite crystal cage, as well as the formation energy of a lattice impurity in a solid at the level of highly accurate many-electron perturbation theories. The absorption energy of a methane molecule in a zeolite chabazite, for instance, is converged with an error well below 20 meV at the CC level. As our largest periodic benchmark system, we apply our scheme to the adsorption of a water molecule on titania in a supercell containing more than 1000 electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Minimax Isometry Method: A compressive sensing approach for Matsubara summation in many-body perturbation theory

We present a compressive sensing approach for the long standing problem of Matsubara summation in many-body perturbation theory. By constructing low-dimensional, almost isometric subspaces of the Hilbert space we obtain optimum imaginary time and frequency grids that allow for extreme data compression of fermionic and bosonic functions in a broad temperature regime. The method is applied to the random phase and self-consistent $GW$ approximation of the grand potential. Integration and transformation errors are investigated for Si and SrVO$_3$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Assessing model-dielectric-dependent hybrid functionals on the antiferromagnetic transition-metal monoxides MnO, FeO, CoO, and NiO

Recently, two nonempirical hybrid functionals, dielectric-dependent range-separated hybrid functional based on the Coulomb-attenuating method (DD-RSH-CAM) and doubly screened hybrid functional (DSH), have been suggested by [Chen et al, Phys. Rev. Mater. 2, 073803 (2018)] and [Cui et al, J. Phys. Chem. Lett. 9, 2338 (2018)], respectively. These two hybrid functionals are both based on a common model dielectric function approach, but differ in the way how to non-empirically obtain the range-separation parameter. By retaining the full short-range Fock exchange and a fraction of the long-range Fock exchange that equals the inverse of the dielectric constant, both DD-RSH-CAM and DSH turn out to perform very well in predicting the band gaps for a large variety of semiconductors and insulators. Here, we assess how these two hybrid functionals perform on challenging antiferromagnetic transition-metal monoxides MnO, FeO, CoO, and NiO by comparing them to other conventional hybrid functionals and the $GW$ method. We find that single-shot DD0-RSH-CAM and DSH0 improve the band gaps towards experiments as compared to conventional hybrid functionals. The magnetic moments are slightly increased, but the predicted dielectric constants are decreased. The valence band density of states (DOS) predicted by DD0-RSH-CAM and DSH0 are as satisfactory as HSE03 in comparison to experimental spectra, however, the conduction band DOS are shifted to higher energies by about 2 eV compared to HSE03. Self-consistent DD-RSH-CAM and DSH deteriorate the results with a significant overestimation of band gaps.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Room temperature dynamic correlation between methylammonium molecules in lead-iodine based perovskites: An ab-initio molecular dynamics perspective

The high efficiency of lead organo-metal-halide perovskite solar cells has raised many questions about the role of the methylammonium (MA) molecules in the Pb-I framework. Experiments indicate that the MA molecules are able to 'freely' spin around at room temperature even though they carry an intrinsic dipole moment. We have performed large supercell (2592 atoms) finite temperature ab-initio molecular dynamics calculations to study the correlation between the molecules in the framework. An underlying long range anti-ferroelectric ordering of the molecular dipoles is observed. The dynamical correlation between neighboring molecules shows a maximum around room temperature in the mid-temperature phase. In this phase, the rotations are slow enough to (partially) couple to neighbors via the Pb-I cage. This results in a collective motion of neighboring molecules in which the cage acts as the mediator. At lower and higher temperatures the motions are less correlated.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

$GW$100: a plane wave perspective for small molecules

In a recent work, van Setten and coworkers have presented a carefully converged $G_0W_0$ study of 100 closed shell molecules [J. Chem. Theory Comput. 11, 5665 (2015)]. For two different codes they found excellent agreement to within few 10 meV if identical Gaussian basis sets were used. We inspect the same set of molecules using the projector augmented wave method and the Vienna ab initio simulation package (VASP). For the ionization potential, the basis set extrapolated plane wave results agree very well with the Gaussian basis sets, often reaching better than 50 meV agreement. In order to achieve this agreement, we correct for finite basis set errors as well as errors introduced by periodically repeated images. For electron affinities below the vacuum level differences between Gaussian basis sets and VASP are slightly larger. We attribute this to larger basis set extrapolation errors for the Gaussian basis sets. For quasi particle (QP) resonances above the vacuum level, differences between VASP and Gaussian basis sets are, however, found to be substantial. This is tentatively explained by insufficient basis set convergence of the Gaussian type orbital calculations as exemplified for selected test cases.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Correlation energy for the homogeneous electron gas: exact Bethe-Salpeter solution and new approximate evaluation

The correlation energy of the homogeneous electron gas is evaluated by solving the Bethe-Salpeter equation (BSE) beyond the Tamm-Dancoff approximation for the electronic polarisation propagator. The BSE is expected to improve upon the random phase approximation, owing to the inclusion of exchange diagrams. For instance, since the BSE reduces in second order to Møller-Plesset perturbation theory, it is self-interaction free in second order. Results for the correlatione energy are compared with Quantum Monte Carlo benchmarks and excellent agreement is observed. For low densities, however, we find imaginary eigenmodes in the polarisation propagator. To avoid the occurence of imaginary eigenmodes, an approximation to the BSE kernel is proposed, which allows to completely remove this issue in the low electron density region. We refer to this approximation as the random phase approximation with screened exchange (RPAsX). We show that this approximation even slightly improves upon the standard BSE kernel.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cubic scaling $GW$: towards fast quasiparticle calculations

Within the framework of the full potential projector-augmented wave methodology, we present a promising low-scaling $GW$ implementation. It allows for quasiparticle calculations with a scaling that is cubic in the system size and linear in the number of $k$ points used to sample the Brillouin zone. This is achieved by calculating the polarizability and self-energy in the real space and imaginary time domain. The transformation from the imaginary time to the frequency domain is done by an efficient discrete Fourier transformation with only a few nonuniform grid points. Fast Fourier transformations are used to go from real space to reciprocal space and vice versa. The analytic continuation from the imaginary to the real frequency axis is performed by exploiting Thiele's reciprocal difference approach. Finally, the method is applied successfully to predict the quasiparticle energies and spectral functions of typical semiconductors (Si, GaAs, SiC, and ZnO), insulators (C, BN, MgO, and LiF), and metals (Cu and SrVO$_3$). The results are compared with conventional $GW$ calculations. Good agreement is achieved, highlighting the strength of the present method.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Electron and hole doping in the relativistic Mott insulator Sr$_2$IrO$_4$: a first-principles study using band unfolding technique

We study the effects of dilute La and Rh doping on the electronic structure of the relativistic Mott insulator Sr$_2$IrO$_4$ using fully relativistic and magnetically non-collinear density functional theory with the inclusion of an on-site Hubbard $U$. To model doping effects, we have adopted the supercell approach that allows for a realistic treatment of structural relaxations and electronic effects beyond a purely rigid band approach. By means of the band unfolding technique we have computed the spectral function and constructed the effective band structure and Fermi surface (FS), which are readily comparable with available experimental data. Our calculations clearly indicate that La and Rh doping can be interpreted as effective electron and (fractional) hole doping, respectively. In Sr$_{2-x}$La$_x$IrO$_4$ the IMT is accompanied by a moderate renormalization of the electronic correlation substantiated by a reduction of the effective on-site Coulomb repulsion $U-J$ from 1.6 eV ($x=0$) to 1.4 eV (metallic regime $x=12.5\%$). The progressive closing of the relativistic Mott gap leads to the emergence of connected elliptical electron pockets at ($π$/2,$π$/2) and less intense features at $X$ in the FS. The substitution of Ir with the nominally isovalent Rh is accompanied by a substantial hole transfer from the Rh site to the nearest neighbor Ir sites. This shifts down the chemical potential, creates almost circular disconnected hole pockets in the FS and establishes the emergence of a two-dimensional metallic state formed by conducting Rh-planes intercalated by insulating Ir-planes. Finally, our data indicate that hole doping causes a flipping of the in-plane net ferromagnetic moment in the Rh plane and induces a magnetic transition from the AF-I to the AF-II ordering.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Role of Polar Phonons in the Photo Excited State of Metal Halide Perovskites

The development of high efficiency perovskite solar cells has sparked a multitude of measurements on the optical properties of these materials. For the most studied methylammonium(MA)PbI$_3$ perovskite, a large range (6-55 meV) of exciton binding energies has been reported by various experiments. The existence of excitons at room temperature is unclear. For the MAPb$X_3$ perovskites we report on relativistic $GW$-BSE calculations. This method is capable to directly calculate excitonic properties from first-principles. At low temperatures it predicts exciton binding energies in agreement with the reported 'large' values. For MAPbI$_3$, phonon modes present in this frequency range have a negligible contribution to the ionic screening. By calculating the polarisation in time from finite temperature molecular dynamics, we show that at room temperature this does not change. We therefore exclude ionic screening as an explanation for the experimentally observed reduction of the exciton binding energy at room temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Screened moments and absence of ferromagnetism in FeAl

While the stoichiometric intermetallic compound FeAl is found to be paramagnetic in experiment, standard band-theory approaches predict the material to be ferromagnetic. We show that this discrepancy can be overcome by a better treatment of electronic correlations with density functional plus dynamical mean field theory. Our results show no ferromagnetism down to 100 K and since the susceptibility is decreasing at the lowest temperatures studied we also do not expect ferromagnetism at even lower temperatures. This behavior is found to originate from temporal quantum fluctuations that screen short-lived local magnetic moments of 1.6 $μ_B$ on Fe.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Anisotropic Magnetic Couplings and Structure-Driven Canted to Collinear Transitions in Spin-orbit Coupled Sr2IrO4

We put forward a scheme to study the anisotropic magnetic couplings in Sr2IrO4 by mapping fully relativistic constrained noncollinear density functional theory including an on-site Hubbard U correction onto a general spin model Hamiltonian. This procedure allows for the simultaneous account and direct control of the lattice, spin and orbital interactions within a fully ab initio scheme. We compute the isotropic, single site anisotropy and Dzyaloshinskii-Moriya (DM) coupling parameters, and clarify that the origin of the canted magnetic state in Sr2IrO4 arises from the interplay between structural distortions and the competition between isotropic exchange and DM interactions. A complete magnetic phase diagram with respect to the tetragonal distortion and the rotation of IrO6 octahedra is constructed, revealing the presence of two types of canted to collinear magnetic transitions: a spin-flop transition with increasing tetragonal distortion and a complete quenching of the basal weak ferromagnetic moment below a critical octahedral rotation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hidden scale invariance of metals

Density functional theory (DFT) calculations of 58 liquid elements at their triple point show that most metals exhibit near proportionality between thermal fluctuations between virial and potential-energy in the isochoric ensemble. This demonstrates a general "hidden" scale invariance of metals making the dense part of the thermodynamic phase diagram effectively one dimensional with respect to structure and dynamics. DFT computed density scaling exponents, related to the Gr{ü}neisen parameter, are in good agreement with experimental values for 16 elements where reliable data were available. Hidden scale invariance is demonstrated in detail for magnesium by showing invariance of structure and dynamics. Computed melting curves of period three metals follow curves with invariance (isomorphs). The experimental structure factor of magnesium is predicted by assuming scale invariant inverse power-law (IPL) pair interactions. However, crystal packings of several transition metals (V, Cr, Mn, Fe, Nb, Mo, Ta, W and Hg), most post-transition metals (Ga, In, Sn, and Tl) and the metalloids Si and Ge cannot be explained by the IPL assumption. Thus, hidden scale invariance can be present even when the IPL-approximation is inadequate. The virial-energy correlation coefficient of iron and phosphorous is shown to increase at elevated pressures. Finally, we discuss how scale invariance explains the Gr{ü}neisen equation of state and a number of well-known empirical melting and freezing rules.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Prediction of new thermodynamically stable aluminum oxides

Recently, it has been shown that under pressure, unexpected and counterintuitive chemical compounds become stable. Laser shock experiments (A. Rode, unpublished) on alumina (Al2O3) have shown non-equilibrium decomposition of alumina with the formation of free Al and a mysterious transparent phase. Inspired by these observations, with have explored the possibility of the formation of new chemical compounds in the system Al-O. Using the variable-composition structure prediction algorithm USPEX, in addition to the well-known Al2O3, we have found two extraordinary compounds Al4O7 and AlO2 to be thermodynamically stable in the pressure range 330-443 GPa and above 332 GPa, respectively. Both of these compounds at the same time contain oxide O2- and peroxide O22- ions, and both are insulating. Peroxo-groups are responsible for gap states, which significantly reduce the electronic band gap of both Al4O7 and AlO2.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Singles correlation energy contributions in solids

The random phase approximation to the correlation energy often yields highly accurate results for condensed matter systems. However, ways how to improve its accuracy are being sought and here we explore the relevance of singles contributions for prototypical solid state systems. We set out with a derivation of the random phase approximation using the adiabatic connection and fluctuation dissipation theorem, but contrary to the most commonly used derivation, the density is allowed to vary along the coupling constant integral. This yields results closely paralleling standard perturbation theory. We re-derive the standard singles of Görling-Levy perturbation theory [Görling and Levy, Phys. Rev. A {\bf 50}, 196 (1994)], highlight the analogy of our expression to the renormalized singles introduced by Ren and coworkers [Ren, Tkatchenko, Rinke, and Scheffler, Phys. Rev. Lett. {\bf 106}, 153003 (2011)], and introduce a new approximation for the singles using the density matrix in the random phase approximation. We discuss the physical relevance and importance of singles alongside illustrative examples of simple weakly bonded systems, including rare gas solids (Ne, Ar, Xe), ice, adsorption of water on NaCl, and solid benzene. The effect of singles on covalently and metallically bonded systems is also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Kohn-Sham band gaps and potentials of solids from the optimised effective potential method within the random phase approximation

We present an implementation of the optimised effective potential (OEP) scheme for the exact-exchange (EXX) and random phase approximation (RPA) energy functionals and apply these methods to a range of bulk materials. We calculate the Kohn-Sham (KS) potentials and the corresponding band gaps and compare them to the potentials obtained by standard local density approximation (LDA) calculations. The KS gaps increase upon going from the LDA to the OEP in the RPA and finally to the OEP for EXX. This can be explained by the different depth of the potentials in the bonding and interstitial regions. To obtain the true quasi-particle gaps the derivative discontinuities or $G_0W_0$ corrections need to be added to the RPA-OEP KS gaps. The predicted $G_0W_0$@RPA-OEP quasi-particle gaps are about 5% too large compared to the experimental values. However, compared to $G_0W_0$ calculations based on local or semi-local functionals, where the errors vary between different materials, we obtain a rather consistent description among all the materials.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Predictive GW calculations using plane waves and pseudopotentials

We show that quasiparticle (QP) energies as calculated in the $GW$ approximation converge to the wrong value using the projector augmented wave (PAW) method, since the overlap integrals between occupied orbitals and high energy, plane wave like orbitals, are incorrectly described. The error is shown to be related to the incompleteness of the partial wave basis set inside the atomic spheres. It can be avoided by adopting norm-conserving partial waves, as shown by analytic expressions for the contribution from unoccupied orbitals with high kinetic energy. Furthermore, $G_0W_0$ results based on norm-conserving PAW potentials are presented for a large set of semiconductors and insulators. Accurate extrapolation procedures to the infinite basis set limit and infinite k-point limit are discussed in detail.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The random phase approximation applied to ice

Standard density functionals without van der Waals interactions yield an unsatisfactory description of ice phases, specifically, high density phases occurring under pressure are too unstable compared to the common low density phase I$_h$ observed at ambient conditions. Although the description is improved by using functionals that include van der Waals interactions, the errors in relative volumes remain sizable. Here we assess the random phase approximation (RPA) for the correlation energy and compare our results to experimental data as well as diffusion Monte Carlo data for ice. The RPA yields a very balanced description for all considered phases, approaching the accuracy of diffusion Monte Carlo in relative energies and volumes. This opens a route towards a concise description of molecular water phases on surfaces and in cavities.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Computing Gibbs free energy differences by interface pinning

We propose an approach for computing the Gibbs free energy difference between phases of a material. The method is based on the determination of the average force acting on interfaces that separate the two phases of interest. This force, which depends on the Gibbs free energy difference between the phases, is computed by applying an external harmonic field that couples to a parameter which specifies the two phases. Validated first for the Lennard-Jones model, we demonstrate the flexibility, efficiency and practical applicability of this approach by computing the melting temperatures of sodium, magnesium, aluminum and silicon at ambient pressure using density functional theory. Excellent agreement with experiment is found for all four elements, except for silicon, for which the melting temperature is, in agreement with previous simulations, seriously underestimated.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Assessment of correlation energies based on the random-phase approximation

The random-phase approximation to the ground state correlation energy (RPA) in combination with exact exchange (EX) has brought Kohn-Sham (KS) density functional theory one step closer towards a universal, "general purpose first principles method". In an effort to systematically assess the influence of several correlation energy contributions beyond RPA, this work presents dissociation energies of small molecules and solids, activation energies for hydrogen transfer and non-hydrogen transfer reactions, as well as reaction energies for a number of common test sets. We benchmark EX+RPA and several flavors of energy functionals going beyond it: second-order screened exchange (SOSEX), single excitation (SE) corrections, renormalized single excitation (rSE) corrections, as well as their combinations. Both the single excitation correction as well as the SOSEX contribution to the correlation energy significantly improve upon the notorious tendency of EX+RPA to underbind. Surprisingly, activation energies obtained using EX+RPA based on a KS reference alone are remarkably accurate. RPA+SOSEX+rSE provides an equal level of accuracy for reaction as well as activation energies and overall gives the most balanced performance, which makes it applicable to a wide range of systems and chemical reactions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Comparing $GW$+DMFT and LDA+DMFT for the testbed material SrVO$_3$

We have implemented the $GW$+dynamical mean field theory (DMFT) approach in the Vienna ab initio simulation package. Employing the interaction values obtained from the locally unscreened random phase approximation (RPA), we compare $GW$+DMFT and LDA+DMFT against each other and against experiment for SrVO$_3$. We observed a partial compensation of stronger electronic correlations due to the reduced $GW$ bandwidth and weaker correlations due to a larger screening of the RPA interaction, so that the obtained spectra are quite similar and well agree with experiment. Noteworthily, the $GW$+DMFT better reproduces the position of the lower Hubbard side band.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Convergence of many-body wavefunction expansions using a plane wave basis: from the homogeneous electron gas to the solid state

Using the finite simulation-cell homogeneous electron gas (HEG) as a model, we investigate the convergence of the correlation energy to the complete basis set (CBS) limit in methods utilising plane-wave wavefunction expansions. Simple analytic and numerical results from second-order Møller-Plesset theory (MP2) suggest a 1/M decay of the basis-set incompleteness error where M is the number of plane waves used in the calculation, allowing for straightforward extrapolation to the CBS limit. As we shall show, the choice of basis set truncation when constructing many-electron wavefunctions is far from obvious, and here we propose several alternatives based on the momentum transfer vector, which greatly improve the rate of convergence. This is demonstrated for a variety of wavefunction methods, from MP2 to coupled-cluster doubles theory (CCD) and the random-phase approximation plus second-order screened exchange (RPA+SOSEX). Finite basis-set energies are presented for these methods and compared with exact benchmarks. A transformation can map the orbitals of a general solid state system onto the HEG plane wave basis and thereby allow application of these methods to more realistic physical problems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effective Onsite Interaction for Dynamical Mean-Field Theory

A scheme to incorporate non-local polarizations into the dynamical mean-field theory (DMFT) and a tailor-made way to determine the effective interaction for the DMFT are systematically investigated. Applying it to the two-dimensional Hubbard model, we find that non-local polarizations induce a non-trivial filling-dependent anti-screening effect for the effective interaction. The present scheme combined with density functional theory offers an ab initio way to derive effective onsite interactions for the impurity problem in DMFT. We apply it to SrVO3 and find that the anti-screening competes with the screening caused by the off-site interaction.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The multiferroic phase of DyFeO$_{3}$:an ab--initio study

By performing accurate ab-initio density functional theory calculations, we study the role of $4f$ electrons in stabilizing the magnetic-field-induced ferroelectric state of DyFeO$_{3}$. We confirm that the ferroelectric polarization is driven by an exchange-strictive mechanism, working between adjacent spin-polarized Fe and Dy layers, as suggested by Y. Tokunaga [Phys. Rev. Lett, \textbf{101}, 097205 (2008)]. A careful electronic structure analysis suggests that coupling between Dy and Fe spin sublattices is mediated by Dy-$d$ and O-$2p$ hybridization. Our results are robust with respect to the different computational schemes used for $d$ and $f$ localized states, such as the DFT+$U$ method, the Heyd-Scuseria-Ernzerhof (HSE) hybrid functional and the GW approach. Our findings indicate that the interaction between the $f$ and $d$ sublattice might be used to tailor ferroelectric and magnetic properties of multiferroic compounds.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Theory of Spin-Conserving Excitation of the $N-V^-$ Center in Diamond

The negatively charged nitrogen-vacancy defect ($N-V^-$ center) in diamond is an important atomic-scale structure that can be used as a qubit in quantum computing and as a marker in biomedical applications. Its usefulness relies on the ability to optically excite electrons between well-defined gap states, which requires clear and detailed understanding of the relevant states and excitation processes. Here we show that by using hybrid density-functional-theory calculations in a large supercell we can reproduce the zero-phonon line and the Stokes and anti-Stokes shifts, yielding a complete picture of the spin-conserving excitation of this defect.

0 citations0 reviewsNo code linkedOpen access