Broken symmetry solutions in one-dimensional lattice models via many-body perturbation theory

In this work we study self-consistent solutions in one-dimensional lattice models obtained via many-body perturbation theory. The Dyson equation is solved in a fully self-consistent manner via the algorithmic inversion method based on the sum-over-pole representation (AIM-SOP) of dynamical operators. We start by validating our self-consistent AIM-SOP implementation by taking as test case the one-dimensional Hubbard model. We then move to the study of antiferromagnetic and charge density wave solutions in one-dimensional lattice models, taking into account a long-range Coulomb interaction between the electrons. Complementary, by solving the Sham-Schl\"uter equation, we can compute the non-interacting Green's function reproducing the same charge density of the interacting system. In turn, this allows for the evaluation of the derivative discontinuity in the Kohn-Sham potential, which gives a measure of how the Kohn-Sham gap approximates the many-body fundamental gap of these systems.