# Multitier self-consistent $GW$+EDMFT

## TL;DR

This paper introduces a parameter-free, self-consistent $GW$+EDMFT method for simulating correlated materials, effectively capturing local and nonlocal interactions and screening effects with improved accuracy over previous approaches.

## Contribution

The paper develops and details a multitier $GW$+EDMFT framework that combines ab initio $G^0W^0$ calculations with a self-consistent $GW$ plus EDMFT treatment for better modeling of correlated materials.

## Key findings

- The method accurately captures correlation effects in intermediate regimes.
- It predicts a Mott transition at a specific lattice spacing.
- Impurity interactions show expected trends with lattice spacing.

## Abstract

We discuss a parameter-free and computationally efficient ab initio simulation approach for moderately and strongly correlated materials, the multitier self-consistent $GW$+EDMFT method. This scheme treats different degrees of freedom, such as high-energy and low-energy bands, or local and nonlocal interactions, within appropriate levels of approximation, and provides a fully self-consistent description of correlation and screening effects in the solid. The ab initio input is provided by a one-shot $G^0W^0$ calculation, while the strong-correlation effects originating from narrow bands near the Fermi level are captured by a combined $GW$ plus extended dynamical mean-field (EDMFT) treatment. We present the formalism and technical details of our implementation and discuss some general properties of the effective EDMFT impurity action. In particular, we show that the retarded impurity interactions can have non-causal features, while the physical observables, such as the screened interactions of the lattice system, remain causal. We then turn to stretched sodium as a model system to explore the performance of the multitier self-consistent $GW$+EDMFT method in situations with different degrees of correlation. While the results for the physical lattice spacing $a_0$ show that the scheme is not very accurate for electron-gas like systems, because nonlocal corrections beyond $GW$ are important, it does provide physically correct results in the intermediate correlation regime, and a Mott transition around a lattice spacing of $1.5a_0$. Remarkably, even though the Wannier functions in the stretched compound are less localized, and hence the bare interaction parameters are reduced, the self-consistently computed impurity interactions show the physically expected trend of an increasing interaction strength with increasing lattice spacing.