Using the inequality of the means for 2 variables, $(\sqrt{ab}\leq\frac{a+b}{2})$ for positive $a$ and $b$, equality only when $a=b$ we have
\begin{equation}
\begin{aligned}
&\sqrt{(x_1+x_2)(x_1+x_3)(x_1+x_4)} + \sqrt{(x_2+x_1)(x_2+x_3)(x_2+x_4)} \\&=\sqrt{(x_1+x_2)}(\sqrt{(x_1+x_3)(x_1+x_4)}+\sqrt{(x_2+x_3)(x_2+x_4)})\\ &\leq \sqrt{(x_1+x_2)}(\frac{(x_1+x_3)+(x_1+x_4)}{2}+\frac{(x_2+x_3)+(x_2+x_4)}{2})\\&=\sqrt{(x_1+x_2)}(x_1+x_2+x_3+x_4)\\&=\sqrt{(x_1+x_2)}
\end{aligned}
\end{equation}
with equality attained only when $x_3=x_4$.
Similarly or by swapping $x_1$ and $x_3$, $x_2$ and $x_4$ we have
\begin{equation}
\begin{aligned}
&\sqrt{(x_3+x_1)(x_3+x_2)(x_3+x_4)} + \sqrt{(x_4+x_1)(x_4+x_2)(x_4+x_3)} \\&\leq\sqrt{(x_3+x_4)}
\end{aligned}
\end{equation}
with equality attained only when $x_1=x_2$.
Adding the two inequalities we obtain
\begin{equation}
\begin{aligned}
f(x_1,~x_2,~x_3,~x_4) ~&\leq \sqrt{(x_1+x_2)}+\sqrt{(x_3+x_4)}
\end{aligned}
\end{equation}
By Jensen's Inequality since $-\sqrt{x}$ is convex for $0\leq x \leq 1$ we have
\begin{equation}
\begin{aligned}
\frac{\sqrt{(x_1+x_2)}+\sqrt{(x_3+x_4)}}{2}\leq \sqrt{\frac{(x_1+x_2)+(x_3+x_4)}{2}}=\sqrt{\frac{1}{2}}
\end{aligned}
\end{equation}
with equality attained only when $(x_1+x_2)=(x_3+x_4)=\frac{1}{2}$.
Hence
\begin{equation}
\begin{aligned}
f(x_1,~x_2,~x_3,~x_4) ~&\leq \sqrt{2}
\end{aligned}
\end{equation}
with equality only attained if $x_1=x_2$, $x_3=x_4$ and $(x_1+x_2)=(x_3+x_4)=\frac{1}{2}$ which is equivalent to $x_1=x_2=x_3=x_4=\frac{1}{4}.$