Is there a way to construct homogenous polynomials of small degree over a certain finite field having only trivial zero?
For instance, the polynomial $f (a, b, c) = a^3 + b^3 + c^3 - 3abc - 3a^2b - 3b^2c - 3c^2a$ has only a trivial zero over $\mathbb{F}_{13}^3$ (I found this by trial and error). Is there a nice way to show that this polynomial has indeed the forementioned property? Can we find other polynomials over $\mathbb{F}_{13}$ with this property?
If you can find a matrix $M$ of integers such that $$ \det (aI + b M + c M^2) $$ is your polynomial, you have a method. For example, with $$ M \; = \; \left( \begin{array}{rrr} 0 & 1 & 0 \\ 0 & 0 & 1 \\ 1 & 0 & 0 \end{array} \right) , $$ we get $$ aI + b M + c M^2 = \left( \begin{array}{rrr} a & b & c \\ c & a & b \\ b & c & a \end{array} \right) . $$ In this case $$ \det ( aI + b M + c M^2) = a^3 + b^3 + c^3 - 3 a b c. $$ As a result, you can see at a glance that this is singular when $a=b=c,$ but it is also singular when $a+b+c = 0.$ So the determinant is $0$ in those two cases. In your case, if you find an appropriate matrix $M,$ which is the low-budget version of Jyrki's comment, you can do Gauss-Jordan elimination in the field of 13 elements to find out when $ aI + b M + c M^2$ is singular, which will (as you say) be only when $$ a \equiv b \equiv c \equiv 0 \pmod {13}. $$
I will see if I can come up with a matrix $M$ for your problem. I think Jyrki would know an elegant method for doing that, but I will do trial and error.
EEDDIITTTT: if absolutely necessary, one may use a matrix where the determinant does not come out the same but is equivalent to your polynomial mod 13; so a $-3$ coefficient in your polynomial could become $-16,$ for instance. This variant is limited to one prime at a time...
EEDDDIIIIITT, Thursday, 6 March. It did not turn out as cleanly as I had expected, but it did have an identity matrix and (the negation of) a companion matrix, so everything still holds. In this case we find $$f (a, b, c) = a^3 + b^3 + c^3 - 3abc - 3a^2b - 3b^2c - 3c^2a = \det (a X + b Y + c I),$$ with $$ Y \; = \; \left( \begin{array}{rrr} 0 & -1 & 0 \\ 0 & 0 & -1 \\ 1 & -3 & 0 \end{array} \right) , $$ also $X = I + Y - Y^2,$ so $$ X \; = \; \left( \begin{array}{rrr} 1 & -1 & -1 \\ 1 & -2 & -1 \\ 1 & -2 & -2 \end{array} \right) . $$
All we need for the multiplication property to hold is to be able to express $Y^2, XY,YX, X^2$ as sums of $X,Y,I$ with integer coefficients. Note first that $Y^3 = 3 Y + I.$ Then we get $$ Y^2 = Y - X + I, $$ $$ XY = YX = -Y - X, $$ $$ X^2 = - Y - 2 X + I. $$ Done. Everything works. So, there is a middle ground. The identity matrix and a companion matrix, then some polynomial in the companion matrix. Live and learn.