IIR with asymptotic impulse response
As per the exact solution of a discrete LTI system: 

$$y(t)=C A^k x_0+C\sum_{k=0}^{t-1}A^{t-k-1} B u(k)+D u(t)$$

and considering that $u(t)=\delta(t)$ and $x_0=0$ for the impulse response, we have:

$$h(t)=C A^{t-1} B+D \delta(t)$$

where $A$,$B$,$C$ and $D$ are the state space discrete matrices:

$$x(t+1)=Ax(t)+Bu(t)\\y(t)=Cx(t)+Du(t)$$

Which in this case are obtained as:

b=[1/6 -1/3 0];
a=[1 1/6 -1/3];
[A,B,C,D]=tf2ss(b,a);


Hence, aside numerical issues, both methods gave the same results:

[h,t] = impz([1 -2],[6 1 -2]);
stem(t,h);

t=(0:23)';
h(1)=D*1;
for k=2:length(t)
    h(k,1)=C*A^(k-2)*B;
end
stem(t,h);


So the $t$=23th term on the $h(t)$ signal is:
$$h(23)=C A^{t-1} B|_{t=23}=C A^{22} B=-3.39702260269318e^{-05} $$

Though small, the $t$th term never vanishes, hence, transfer functions with denominator different than the identity vector, will always be an IIR.