PACS 31.15.ec -Hohenberg-Kohn theorem and formal mathematical properties, completeness theorems PACS 71.10.Fd -Lattice fermion models (Hubbard model, etc.) PACS 71.15.Mb -Density functional theory, local density approximation, gradient and other correctionsAbstract -We demonstrate that, for a fermionic lattice system, the ground-state particle density uniquely determines the external potential except for the sites corresponding to nodes of the wave function, and the limiting case where the Pauli exclusion principle completely determines the occupation of all sites. Our fundamental finding completes, for this general class of systems, the one-to-one correspondence between ground states, their densities, and the external potential at the base of the Hohenberg-Kohn theorem. Moreover we demonstrate that the mapping from wave function to potential is unique not just for the ground state, but also for excited states. To illustrate our findings, we develop a practical inversion scheme to determine the external potential from a given density. Our results hold for a general class of lattice models, which includes the Hubbard model.Mapping many-body systems to a lattice while neglecting interactions beyond a given order has proven a powerful aid to devise tractable models -such as the Hubbard [1] and Heisenberg (see, e.g., Ref.[2]) models -and further the understanding of their properties. Lattice systems are of recent interest for quantum chemistry calculations when the one-particle wave functions can be appropriately mapped to sites in a one-dimensional lattice and DMRG techniques are used to model the system (see, e.g., Refs. 3-5). In addition, lattice systems, both bosonic and fermionic, may be simulated experimentally using an optical lattice (see, e.g., Refs. 6-8). The latter is of great interest to the quantum technology community for developing quantum simulators. Fermionic lattice systems have been used to model the entanglement of nanostructures [9] and, using time-dependent density-functional theory, have been employed to investigate the time evolution of the out-ofequilibrium Mott insulator [10]. Time-dependent densityfunctional theory has also been used to create a method to simulate driven lattice gases with interactions [11] and used with other classical lattice Hamiltonians to model particle transport against a bias [12]. Lattice Hamiltonians have been used to describe the process of excitationenergy transfer and the role of intramolecular vibrations in the photosynthetic process [13]. Lattice systems are also used to simulate molecular systems, as for example the Pariser-Parr-Pople (PPP) model [14,15], similar to the inhomogeneous Hubbard model, can describe molecular systems with delocalised electrons in π-orbitals and has recently been used to model armchair polyacenes in Ref. 16. In the last decades density-functional theory (DFT) [17] in all its flavors has been a very important tool to predict properties of materials and structures. DFT is in principle an exact theory and allows to map many-body p...