I'd like to optimize $x_1 x_2 x_3 + x_1 x_2 x_4 + x_1 x_3 x_4 + x_2 x_3 x_4$ on the unit 4-sphere. I'm thinking I should do lagrangian optimization, but I'm having trouble solving the resulting system.
We may assume that the maximum is attained when $x_i \gt 0$ for all $i$.
I happened to "guess" one of the critical points which is at $x_1 = x_2 = x_3 = x_4$, so any argument which shows there is only one critical point will suffice.
Your sum is a sum of geometric means of $x_1^3, x_2^3, x_3^3,x_4^3$ you know each is less than the arithmetic means so the max, $M$, satisfies
$$M\le x_1^3+x_2^3+x_3^3+x_4^3=f(x_1,x_2,x_3,x_4)$$
On the sphere, $f$ has a maximum at $(x_1^2,x_2^2,x_3^2,x_4^2)=\lambda(x_1,x_2,x_3,x_4)$ implying $x_i=x_j$ for all $i,j$. Since the AM is maximized there and GM=AM at this point, you know that there is no means to get a larger value than there.