<s>[INST] Sum of Squares for $a^2+b^2+c^2+d^2+abcd+1\ge ab+bc+cd+da + ac+bd$ $\color{red}{\textrm{Update}}$
To clarify: I hope to see a (simple) SOS solution similar to the one of the following 3-variable problem, no other assumptions e.g. $a\ge b\ge c\ge d$.
Let us see an example. In the link below, by letting $y = \mathrm{mid}(x, y, z)$, 
two SOS solutions are given. They are not what I want. It is more difficult to get a SOS solution without the assumption $y = \mathrm{mid}(x, y, z)$ for that problem.
Prove $2\left(x^2+y^2+z^2+1)(x^3y+y^3z+z^3x+xyz\right) \le \left(x^2+y^2+z^2+3xyz\right)^2.$
Another example: Mongolian TST 2008, day 2 problem 3, Ji Chen gave a SOS solution with the assumption $x = \min(x, y, z)$.
Actually, I gave a SOS solution without any assumption.
https://artofproblemsolving.com/community/c6h205316p11219067
$\phantom{2}$
Problem 1. Let $a, b, c, d \ge 0$. Prove that $a^2+b^2+c^2+d^2+abcd+1\ge ab+bc+cd+da + ac+bd$.
There are quite a few solutions (including my solution). Here, I am particularly interested in $\color{blue}{\textrm{(simple) Sum of Squares (SOS) solutions}}$.
Denote $f(a, b, c, d) = \mathrm{LHS} - \mathrm{RHS}$.
Then, $f(x^2, y^2, z^2, w^2) \ge 0$ for all real numbers $x, y, z, w$.
However, $f(x^2, y^2, z^2, w^2)$ may not be expressed as SOS (polynomial).
I found that $(x^2+y^2+z^2+w^2)^2f(x^2, y^2, z^2, w^2)$ may be expressed as SOS,
since numerically $(x^2+y^2+z^2+w^2)^2f(x^2, y^2, z^2, w^2) \approx u^\mathsf{T}Q u$ where $Q$ is $185\times 185$ matrix
and $u$ is a vector containing monomials in $x, y, z, w$. But $Q$ is quite large, I have not yet proceeded.
Any comments and solutions are welcome.
Relevant information
For 3-variable problem below, there are quite a few SOS solutions.
Problem 2. Let $a, b, c\ge 0$. Prove that $a^2+b^2+c^2+2abc+1 \ge 2(ab+bc+ca)$.
My SOS solution is
\begin{align}
&a^2+b^2+c^2+2abc+1 - 2(ab+bc+ca) \\
=\ & \frac{1}{2(a+b)^2}\Big[(a^2-ac-b^2+bc-a+b)^2 +(a^2-2ab-ac+b^2-bc+a+b)^2\\
&\qquad\qquad\quad + 4ab(a-b)^2 + 4ab(c-1)^2 + 4abc(a+b-2)^2\Big].
\end{align}
FYI, my non-SOS solution for Problem 1: 
WLOG, assume that $a\ge b\ge c\ge d$.
If $cd \ge 1$, then $ab \ge 1$ and $(ab-1)(cd-1) \ge 0$, i.e., $ab+cd \le abcd+1$.
Also, $bc + da + ac + bd \le \frac{b^2+c^2}{2} + \frac{d^2+a^2}{2}
+ \frac{a^2+c^2}{2} + \frac{b^2+d^2}{2} = a^2+b^2+c^2+d^2$. Add them up to get the desired result.
If $cd < 1$, let $u = a - \frac{c+d}{cd+1}, \ v = b - \frac{c+d}{cd+1}$. We have
\begin{align}
\mathrm{LHS} - \mathrm{RHS} &= u^2 + v^2 - (1-cd)uv +\frac{c^3d+cd^3-c^2d^2-2cd+1}{cd+1}\\
&\ge 2|uv| - (1-cd) |uv| + \frac{c^3d+cd^3-c^2d^2-2cd+1}{cd+1}\\
&\ge \frac{c^3d+cd^3-c^2d^2-2cd+1}{cd+1}\\
&\ge \frac{2c^2d^2-c^2d^2-2cd+1}{cd+1}\\
&= \frac{(cd-1)^2}{cd+1}\\
&\ge 0.
\end{align}
(Q.E.D.)
 [/INST] A very nice SOS solution by KaiRain@AoPS:
We have
\begin{align*}
&a^2+b^2+c^2+d^2+abcd+ 1 - ( ab+bc+cd+da + ac+bd)\\
=\,& \frac{\left(2a^2-2ab+2b^2-ac-ad-bc-bd+2 \right)^2+3 \left(a-b \right)^2 \left(c-d \right)^2 }{4 \left(a^2-ab+b^2+1 \right)}\\
&\quad + \frac{abcd \left(a-b \right)^2+cd \left(ab-1 \right)^2 +\left(c-d \right)^2}{a^2-ab+b^2+1 }.
\end{align*}
</s>