Calculations of charged systems in periodic boundary conditions (PBC) are problematic because there are spurious interactions between the charges in different periodic images that can affect the physical picture. In addition, the intuitive limit of Coulomb interactions decaying to zero as the interacting charges are placed at infinite separation no longer applies, and for example total energies become undefined. Leveraging subsystem density functional theory (also known as density embedding) we define an impurity model that embeds a finite neutral or charged subsystem within an extended (infinite) surrounding subsystem. The combination of the impurity model and a consistent choice of the Coulomb reference provides us with an algorithm for evaluating the ionization potential (IP) in extended systems. We demonstrate our approach in a pilot calculation of the IP of liquid water, based on a configuration from a prior ab initio molecular dynamics (AIMD) simulation of liquid water (Genova et al., J. Chem. Phys. 2016, 144, 234105). The calculations with the impurity model capture the broadening on the ionization energies introduced by the interactions between the water molecules. Furthermore, the calculated average IP value (10.5 eV) compare favorably to experiments (9.9-10.06 eV) and very recent simulations based on the GW approximation (10.55 eV), while at the same time outperforming density embedding calculations carried out with a naïve handling of the electrostatic interactions (about 7 eV). K E Y W O R D S DFT, embedding, periodic boundary conditions, water 1 | INTRODUCTION Atomistic models of extended systems (such as solids and liquids) typically prescribe the use of periodic boundary conditions (PBC). Once PBC are invoked, one has the choice of evaluating Coulomb interactions in real space or in reciprocal space. If real space is chosen, considering r, r 0 2R 3 then the Coulomb kernel is w r, r 0 ð Þ¼ 1 jr− r 0 j hinting that it decays to zero when two charges are far away. Despite the decay feature of the kernel, Coulomb potentials and energies for extended systems involve conditionally convergent integrals [1] -a complication which is dealt with by switching to reciprocal space. From a purely theoretical standpoint, in PBC, distance is not a well-defined quantity [2] and the physical kernel is only the one represented in reciprocal space, w G ð Þ¼ 4π G j j 2 (G 2 R 3 an element of the reciprocal space). Thus, complications arise because when G = 0, as the Coulomb kernel becomes undefined. In Ewald summations, [3] this is dealt with using a two-prongued approach: the short-range interactions are dealt with in real space (using a cut-off function, typically the error function), and the long-range interactions are dealt with in reciprocal space.In practice, however, the kernel singularity at G = 0 is not problematic because, although w(G) is singular for G = 0, the total charge density ρ(G) is zero for G = 0. We note that neutrality of the charge density in each unit cell is the only physical choice in PBC, beca...