Look at a point of maximum. If one of the variables equals 0, we deal with $n-1$ instead $n$. If all are positive, by Lagrange multipliers the partial derivatives in all variables are equal. Look at a partial derivative in $x_n$, it is $$2x_n\left(\sum_{i
Assume that they take three distinct values $a,b,c$ with multiplicities $\alpha,\beta,\gamma$. By Vieta theorem for the above cubic we get $a+b+c=3\sigma_1/4=3/4(\alpha a+\beta b+\gamma c)$, and without loss of generality we may suppose that $\gamma=1$, i.e., $c=(3\alpha-1)a+(3\beta-1)b$. Substitute it to $ab+bc+ac=2\sigma_2-p_2=\sigma_1^2-2p_2$ we quickly get a contradiction (please check this).
Now assume that $x_i$'s take two different values $a$, $b$ with multiplicities $\alpha,\beta$. We have $4(a^3-b^3)+(\sigma_1^2-2p_2)(a-b)=3\sigma_1(a^2-b^2)$, dividing by $a-b$ we get $$(4-5\alpha+\alpha^2)a^2+(4-5\beta+\beta^2)b^2=(3\alpha+3\beta-2\alpha\beta-4)ab.$$ It is not hard to see that this is impossible when $\alpha\geqslant 4$, $\beta\geqslant 4$. If, say, $\alpha=1$, we get $b=\beta-4$, so this is possible for $n-1=\beta\geqslant 6$, the value of our sum equals $\beta(\beta-1)(\beta-3)/(32(\beta-2)^3)$, worse than $6/125$.The cases when $\alpha=2$ and $\alpha=3$ are computationally painful, but it looks that again the value is less than $6/125$.