Let $\omega = e^{i\frac{2\pi}{3}}$ be a cubic root of unity. For any $z_1,z_2,z_3 \in \mathbb{C}$, let
$$\begin{align}
a &= \frac13 (z_1 + z_2 + z_3)\\
b &= \frac13 (z_1 + z_2\omega + z_3\omega^2)\\
c &= \frac13 (z_1 + z_2\omega^2 + z_3\omega)
\end{align}$$
It is easy to check
$$\begin{align}
z_1 &= a+b+c\\
z_2 &= a+b\omega^2+c\omega\\
z_3 &= a+b\omega+c\omega^2
\end{align}$$
Let $Q(z_1,z_2,z_3)$ be the quadratic form
$$z_1^2+z_2^2+z_3^2-z_1z_2 - z_2z_3 -z_3z_1$$
The reason that $Q(z_1,z_2,z_3) = 0$ iff $\triangle z_1z_2z_3$ is equilateral comes down to the fact
$$Q(z_1,z_2,z_3) = (z_1 + z_2\omega + z_3\omega^2)(z_1 + z_2\omega^2+z_3\omega) = 9bc$$
So $$Q(z_1,z_2,z_3) = 0 \iff bc = 0 \iff b = 0 \lor c = 0$$
When $b = 0$, $(z_1,z_2,z_3) = (a+c,a+c\omega,a+c\omega^2)$ and we can obtain $z_2,z_3$ by rotating $z_1$ with respect to $a$ counter-clockwisely for angles $120^\circ$ and $240^\circ$ respectively. This implies $\triangle z_1,z_2,z_3$ is equilateral and $z_1,z_2,z_3$ are ordered counter-clockwisely with respect to center $a$.
By a similar argument, if $c = 0$,
$\triangle z_1,z_2,z_3$ is again equilateral but $z_1,z_2,z_3$ are ordered clockwisely with respect to $a$.
To show the quadratic form $Q(z_1,z_2,z_3)$ is essentially unique. Let $Q'(z_1,z_2,z_3)$ be another quadratic form which vanishes on and only on equilateral triangles. Expand $Q'(z_1,z_2,z_3)$ in terms of $a,b,c$, there are coefficients $\alpha,\beta_1,\beta_2, \gamma_1,\gamma_2,\gamma_3$ (not al zero) such that
$$Q'(z_1,z_2,z_3) = \alpha a^2 + a (\beta_1 b + \beta_2 c) + (\gamma_1 b^2 + \gamma_2 bc + \gamma_3 c^2)\tag{*1}$$
Pick any $z_1,z_2,z_3$ such that $\triangle z_1z_2z_3$ is equilateral. If one translate all vertices by same amount, $z_k \mapsto z_k + t$, the resulting triangle
is also equilateral. However, $a \mapsto a + t$ but $b,c$ remains the same. This means for that specific $b,c$.
$(*1) = 0$ for all $a$. This forces
$$\alpha = \beta_1 b + \beta_2 c = \gamma_1 b^2 + \gamma_2 bc + \gamma_3 c^2 = 0$$
It is trivial to find equilateral triangles with $b = 0, c\ne 0$ or triangles with $b \ne 0, c = 0$. This forces
$\beta_1 = \beta_2 = \gamma_1 = \gamma_3 = 0$.
As a result, $\gamma_2 \ne 0$ and
$$Q'(z_1,z_2,z_3) = \gamma_2 bc = \frac{\gamma_2}{9} Q(z_1,z_2,z_3)$$
From this, we can deduce the quadratic form $Q(z_1,z_2,z_3)$ is unique up to non-zero scaling constants.
Best Answer
Let $$z_1=z_0+r(\cos\theta+i\sin\theta),$$ $$z_2=z_0+r(\cos(120^{\circ}+\theta)+i\sin(120^{\circ}+\theta))$$ and $$z_3=z_0+r(\cos(240^{\circ}+\theta)+i\sin(240^{\circ}+\theta)).$$ Id est, it's enough to prove that $$\cos2\theta+\cos(240^{\circ}+2\theta)+\cos(480^{\circ}+2\theta)=0,$$ $$\sin2\theta+\sin(240^{\circ}+2\theta)+\sin(480^{\circ}+2\theta)=0,$$ $$\cos\theta+\cos(120^{\circ}+\theta)+\cos(240^{\circ}+\theta)=0$$ and $$\sin\theta+\sin(120^{\circ}+\theta)+\sin(240^{\circ}+\theta)=0.$$ Can you end it now?