In this paper, we investigate evolution systems in two components, characterized by higher-order spatial derivatives and the presence of two arbitrary functions. Our study begins with an analysis of a fourth-order system. We perform a detailed group classification and identify specific forms of the constitutive functions that allow the system to exhibit additional symmetries in addition to spatial and temporal translations. We extend these results to nth-order systems. Moreover, we derive invariant solutions for these systems. Finally, for each order n, we are able to find non-negative solutions.