  The high-dimensional generalization of the one-dimensional Kuramoto paradigm
has been an essential step in bringing about a more faithful depiction of the
dynamics of real-world systems. Despite the multi-dimensional nature of the
oscillators in these generalized models, the interacting schemes so far have
been dominated by a scalar factor unanimously between any pair of oscillators
that leads eventually to synchronization on all dimensions. As a natural
extension of the scalar coupling befitting for the one-dimensional case, we
take a tentative step in studying numerically and theoretically the coupling
mechanism of $2\times2$ real matrices on two-dimensional Kuramoto oscillators.
One of the features stemmed from this new mechanism is that the matrix coupling
enables the two dimensions of the oscillators to separate their transitions to
either synchronization or desynchronization which has not been seen in other
high-dimensional generalizations. Under various matrix configurations, the
synchronization and desynchronization of the two dimensions combine into four
qualitatively distinct modes of position and motion of the system. We
demonstrate that as one matrix is morphed into another in a specific manner,
the system mode also switches correspondingly either through continuous or
explosive transitions of the order parameters, thus mimicking a range of
behaviors in information science and biology.
