Joint Approximate Diagonalization approach to Quasiparticle Self-Consistent $GW$ calculations

We introduce an alternative route to quasiparticle self-consistent $GW$ calculations ($\mathrm{q}\mathrm{s}GW$) on the basis of a Joint Approximate Diagonalization of the one-body $GW$ Green’s functions $G({\epsilon }_n^{QP})$ taken at the input quasiparticle energies. Such an approach allows working with the full dynamical self-energy, without approximating the latter by a symmetrized static form as in the standard $\mathrm{q}\mathrm{s}GW$ scheme. Calculations on the $GW$100 molecular test set lead nevertheless to a good agreement, at the 65 meV mean-absolute-error accuracy on the ionization potential, with respect to the conventional $\mathrm{q}\mathrm{s}GW$ approach. We show further that constructing the density matrix from the full Green’s function as in the fully self-consistent $\mathrm{s}\mathrm{c}GW$ scheme, and not from the occupied quasiparticle one-body orbitals, allows obtaining a scheme intermediate between $\mathrm{q}\mathrm{s}GW$ and $\mathrm{s}\mathrm{c}GW$ approaches, closer to CCSD(T) reference values.