I'm assuming you're restricting your definition of an infinite impulse response filter to one that can be expressed as
$$
y_n = \sum_{k=0}^N b_k x_{n - k} - \sum_{k=1}^N a_k y_{n - k}
$$
where $y_n$ is the current output of the filter and $x_n$ is the current input to the filter.
In this case, there are $N$ possible modes which are a mix of $n^{p_m}d_m^n$ and $\delta(n - p_m)$. $\delta(n - p_m)$ is just a delayed Kroneker delta, and takes into account the fact that you can express a FIR filter or mixed IIR and FIR filter in the above form. Also, $d_m$ can be any complex number, with the constraint that if the coefficients $a_k$ are real, then any $d_m$ with a non-zero imaginary part must be accompanied by its complex conjugate.
This is basically just direct solutions to linear shift-invariant difference equations which are analogous to linear time-invariant differential equations in continuous time.