Let $a_k$ a an eigenvalue of the adjacency matrix $A$ of a planar cubic graph with $n$ vertices. For the returning paths without backtracking we get the generating function of
$$G(x,a_k)=\frac{1-x^2}{1-a_kx+2x^2} \tag{1} $$
Now multiply each of the $n$ expressions $G(x,a_k)$ and compare it with Ihara's $\zeta$ function for cubic graphs. Using $\chi-1=F-2$$\;=E-n$, there seems to be a mismatch $$ \zeta _{G}(x)={\frac {(1-x^{2})^{n-E}}{\det(I-Ax+2x^{2}I)}} \neq \frac{(1-x^2)^n}{\prod_{k=1}^n (1-a_kx+2x^2)}, $$ where you have to think of the determinant on the LHS written in a diagonal base.
$\hskip1.3in$Why do I have to divide by $(1-x^2)^{E}$ to get Ihara's formula?
$\hskip2.0in$Is each edge interpreted like a $2$-cycle?
Another thought: The Ihara $\zeta$ function written in another form (rightmost) $$ \zeta _{G}(x) = {\frac {(1-x^{2})^{n-E}}{\det(I-Ax+2x^{2}I)}} = \exp \left( \sum \limits_{k=1}^{\infty} N_k \frac{x^k}{k} \right). $$ The non-negative numbers $N_k$ count non-backtracking, tailless closed paths. Since we divide by $(1-x^2)^{E}$ we kind of remove the edges.
Is this edge-removement what makes the Ihara $\zeta$ function count tailless walks, since an edge is the beginning of a tail?