This paper describes an all-electron implementation of the self-consistent GW (sc-GW ) approach-i.e., based on the solution of the Dyson equation-in an all-electron numeric atom-centered orbital basis set. We cast Hedin's equations into a matrix form that is suitable for numerical calculations by means of (i) the resolution-of-identity technique to handle four-center integrals and (ii) a basis representation for the imaginary-frequency dependence of dynamical operators. In contrast to perturbative G 0 W 0 , sc-GW provides a consistent framework for ground-and excited-state properties and facilitates an unbiased assessment of the GW approximation. For excited states, we benchmark sc-GW for five molecules relevant for organic photovoltaic applications: thiophene, benzothiazole, 1,2,5-thiadiazole, naphthalene, and tetrathiafulvalene. At self-consistency, the quasiparticle energies are found to be in good agreement with experiment and, on average, more accurate than G 0 W 0 based on Hartree-Fock or density-functional theory with the Perdew-Burke-Ernzerhof exchange-correlation functional. Based on the Galitskii-Migdal total energy, structural properties are investigated for a set of diatomic molecules. For binding energies, bond lengths, and vibrational frequencies sc-GW and G 0 W 0 achieve a comparable performance, which is, however, not as good as that of exact-exchange plus correlation in the random-phase approximation and its advancement to renormalized second-order perturbation theory. Finally, the improved description of dipole moments for a small set of diatomic molecules demonstrates the quality of the sc-GW ground-state density. Many-body perturbation theory (MBPT) 1 in the GW approach for the electron self-energy 2-4 provides a natural framework for an ab initio, parameter-free description of photo-ionization processes and charged excitations. 5 In recent years, the GW approach has become a popular method for the computation of band gaps and charged excitation energies for extended 6,7 and finite systems. 8,9 In numerical implementations, following Hybertsen and Louie, 10 it is standard practice to treat the GW self-energy as a single-shot perturbation (G 0 W 0 ) acting on a Kohn-Sham (KS) or Hartree-Fock (HF) reference system. Thus, excitation energies are evaluated from first-order Feynman-Dyson perturbation theory as corrections to a set of single-particle eigenvalues.The popularity of the G 0 W 0 approximation stems from the substantial reduction in the complexity of Hedin's equations at first-order perturbation theory: The KS or HF eigenstates from a self-consistent field calculation can be used as basis functions and provide a convenient representation in which the noninteracting Green's function is diagonal. In this basis, only diagonal matrix elements of the self-energy are needed to evaluate quasiparticle corrections at first order. Thus, G 0 W 0 grants a considerable simplification of the linear algebra operations which is decisive for applying the theory to large molecules and solids.Although num...