The general recipe for a G_{0}W_{0} calculation is as follows:
Ad 1) Symmetry has not yet been implemented for GW, sensible calculations can only be done in C_{1}. The def2-TZVPP basis seems to be the most useful, it comes for all tested systems within 0.1 eV of the def2-QZVP result with about half the number of basis functions. In the define session the calculation of the response function needs to be defined. In the final define menu select the ex menu and select the calculation of RPA singlet excitations, or urpa in case of open shell. Select "all all" to get all excitations. For systems and basis set having less than 4000 rpas excitations just set for all excitations. For large systems start to run G_{0}W_{0} with ˜4000 rpas excitations. In subsequent runs add more excitations until a converged result is reached. escf will keep the converged roots, so not much time is lost using this restart approach.
Ad 2) If $gw is set in the control file the quasi particle energies will be evaluated according to equation 13.3. Additional options are described in the keyword section, section 21.2.13.
Ad 4) In the escf run the response function is calculated which is needed to determine the screened coulomb interaction. At the end of this run the actual GW calculation is performed.
Possible source of errors: When dscf or ridft is repeated after escf the sing_a file may not be correct anymore, this may happen when degenerate levels are present. escf will however not recognize this and continue using the previously converged data in sing_a leading to nonsense values for Σ_{c}. Before running escf the old sing_a file has to be removed.