We present theoretical calculations of quasiparticle energies in closed-shell molecules using the GW method. We compare three different approaches: a full-frequency $G_0W_0$ (FF-$G_0W_0$) method with density functional theory (DFT-PBE) used as a star
ting mean field; a full-frequency $GW_0$ (FF-$GW_0$) method where the interacting Greens function is approximated by replacing the DFT energies with self-consistent quasiparticle energies or Hartree-Fock energies; and a $G_0W_0$ method with a Hybertsen-Louie generalized plasmon-pole model (HL GPP-$G_0W_0$). While the latter two methods lead to good agreement with experimental ionization potentials and electron affinities for methane, ozone, and beryllium oxide molecules, FF-$G_0W_0$ results can differ by more than one electron volt from experiment. We trace this failure of the FF-$G_0W_0$ method to the occurrence of incorrect self-energy poles describing shake-up processes in the vicinity of the quasiparticle energies.
Using many-body perturbation theory within the $G_0W_0$ approximation, we explore routes for computing the ionization potential (IP), electron affinity (EA), and fundamental gap of three gas-phase molecules -- benzene, thiophene, and (1,4) diamino-be
nzene -- and compare with experiments. We examine the dependence of the IP on the number of unoccupied states used to build the dielectric function and the self energy, as well as the dielectric function plane-wave cutoff. We find that with an effective completion strategy for approximating the unoccupied subspace, and a converged dielectric function kinetic energy cutoff, the computed IPs and EAs are in excellent quantitative agreement with available experiment (within 0.2 eV), indicating that a one-shot $G_0W_0$ approach can be very accurate for calculating addition/removal energies of small organic molecules. Our results indicate that a sufficient dielectric function kinetic energy cutoff may be the limiting step for a wide application of $G_0W_0$ to larger organic systems.
