We use the composite operator method (COM) to analyze the strongly correlated repulsive Hubbard model, investigating the effect of nearest-neighbor hoppings up to fourth order on a square lattice. We consider two sets of self-consistent equations, one enforcing the Pauli principle and the other imposing charge-charge, spin-spin, and pair-pair correlations using a decoupling scheme developed by L. Roth [Phys. Rev. 184, 451–459, (1969)]. We extract three distinct solutions from these equations: COM1 and COM2 by imposing the Pauli principle and one from Roth decoupling. An overview of the method studying the validity of particle-hole symmetry and the Luttinger theorem for each solution is presented. Additionally, we extend the initial basis to study superconductivity, concluding that it is induced by the Van Hove singularity. Finally, we include higher-order hoppings using realistic estimates for tight binding parameters and compare our results with ARPES measurements on cuprates.