  We present an extension of a framework for simulating single quasiparticle or
collective excitations on top of strongly correlated quantum many-body ground
states using infinite projected entangled pair states, a tensor network ansatz
for two-dimensional wave functions in the thermodynamic limit. Our approach
performs a systematic summation of locally perturbed states in order to obtain
excited eigenstates localized in momentum space, using the corner transfer
matrix method, and generalizes the framework to arbitrary unit cell sizes, the
implementation of global Abelian symmetries and fermionic systems. Results for
several test cases are presented, including the transverse Ising model, the
spin-$\frac{1}{2}$ Heisenberg model and a free fermionic model, to demonstrate
the capability of the method to accurately capture dispersions. We also provide
insight into the nature of excitations at the $k=(\pi,0)$ point of the
Heisenberg model.
