Notice that b=0 or d=0 as a2+ac+c2+1=(a+c/2)2+(3/4)c2+1>0. Suppose that c=0. By assumption, we then have a=0 and the second equation becomes 2ab+ad=0, hence b=−d/2. Substituting this into the first equation gives d2=(4/3)(a2+1). So b2=(a2+1)/3 and
(3a2−b2+1)(3c2−d2+1)=91(8a2+2)(−1−4a2)<0,
(2a+c)b=−(a+2c)d.
Since b=0 or d=0, and (2a+c)(a+2c)=2(a+c)2+ac>0, we deduce that bd<0. Consider the complex numbers z=a+ib and w=c+id. Then
z2+zw+w2=(a2−b2+ac−bd+c2−d2)+i(2ab+ad+bc+2cd)=−1.
So
(z3+z)−(w3+w)=(z−w)(z2+zw+w2+1)=0.
Now
z3+zw3+w=a(a2−3b2+1)+ib(3a2−b2+1),=c(c2−3d2+1)+id(3c2−d2+1).
In particular d(3c2−d2+1)=b(3a2−b2+1). But bd<0. So
(3a2−b2+1)(3c2−d2+1)≤0.