I think a key is to think in terms of fixed points and cycles. You want an $F$ that either treats $\{1,-1,\sqrt{2},-\sqrt{2}\}$ as fixed points or as parts of two-cycles. Your implication $\implies$ is actually a two-way implication $\iff$. So if you can find a cubic polynomial that either fixes these four points or takes them through a two-cycle, you have such an $F$.
There are $10$ permutations of $\{1,-1,\sqrt{2},-\sqrt{2}\}$ that only involve fixed points and two-cycles. For each one, we can write down the corresponding system of four linear equations in $F$s coefficients. Each system has the same coefficient matrix, but different constant vector. For instance the system corresponding to $(1\leftrightarrow\sqrt{2},-1\leftrightarrow-\sqrt{2})$ is
$$\begin{aligned}
A+B+C+D&=\sqrt{2}\\
-A+B-C+D&=-\sqrt{2}\\
2\sqrt{2}A+2B+\sqrt{2}C+D&=1\\
-2\sqrt{2}A+2B-\sqrt{2}C+D&=-1\\
\end{aligned}$$
Or in matrix form, is
$$\begin{aligned}
\begin{bmatrix}
1&1&1&1\\
-1&1&-1&1\\
2\sqrt{2}&2&\sqrt{2}&1\\
-2\sqrt{2}&2&-\sqrt{2}&1\\
\end{bmatrix}
\begin{bmatrix}A\\B\\C\\D
\end{bmatrix}
&=\begin{bmatrix}\sqrt{2}\\-\sqrt{2}\\1\\-1\end{bmatrix}\\
M
\begin{bmatrix}A\\B\\C\\D
\end{bmatrix}
&=\begin{bmatrix}\sqrt{2}\\-\sqrt{2}\\1\\-1\end{bmatrix}\\
\end{aligned}$$
and $M$ is the same for each permutation. It is worth our effort then to compute $$M^{-1}=\begin{bmatrix}
-1/2&1/2 &\sqrt{2}/4&-\sqrt{2}/4\\
-1/2&-1/2&1/2&1/2\\
1&-1&-\sqrt{2}/4&\sqrt{2}/4\\
1&1&-1/2&-1/2
\end{bmatrix}$$
And so
$$\begin{aligned}
\begin{bmatrix}A\\B\\C\\D
\end{bmatrix}
&=\begin{bmatrix}
-1/2&1/2 &\sqrt{2}/4&-\sqrt{2}/4\\
-1/2&-1/2&1/2&1/2\\
1&-1&-\sqrt{2}/4&\sqrt{2}/4\\
1&1&-1/2&-1/2
\end{bmatrix}
\begin{bmatrix}\sqrt{2}\\-\sqrt{2}\\1\\-1\end{bmatrix}\\
&=\begin{bmatrix}-\sqrt{2}/2\\0\\3\sqrt{2}/2\\0
\end{bmatrix}
\end{aligned}$$
If this is right, then $F(x)={-\frac{\sqrt{2}}{2}}x^3+\frac{3\sqrt{2}}{2}x$ works. Well, a CAS confirms that $F(F(x))-x$ factors as $1/4 (x-1) x (x+1) (x^2-2) (x^4-6 x^2+7)$.
You can get more solutions for $F$ using other permutations on $\{1,-1,\sqrt{2},-\sqrt{2}\}$ that only involve fixed points and two-cycles. However some of these will lead to $0$ for $A$, and not quite give a cubic $F$. For example, the identity permutation leads to $F(x)=x$ and $(1\leftrightarrow-1,\sqrt{2}\leftrightarrow-\sqrt{2})$ leads to $F(x)=-x$.
There are 10 permutations in $S_4$ using only fixed points and two-cycles; two of them certainly lead to degenerate "cubics" as cited above. At least one has now been confirmed to lead to a truly cubic $F$, and my guess would be that the same would happen for the remaining $7$. Thus this method provides 8 $F$s that work, plus two more that might count if degenerate cubics are allowed.
I think to find all solutions to this problem, you could invoke a system of eight linear equations and then find under what conditions that system is consistent. We would allow each of $\{1,-1,\sqrt{2},-\sqrt{2}\}$ to go to arbitrary real numbers $\{m,n,p,q\}$ and then demand that they come back. (These equations will be linear in the coefficients of $F$, but not in $\{m,n,p,q\}$.)
$$\begin{aligned}
\begin{bmatrix}
1&1&1&1\\
-1&1&-1&1\\
2\sqrt{2}&2&\sqrt{2}&1\\
-2\sqrt{2}&2&-\sqrt{2}&1\\
m^3&m^2&m&1\\
n^3&n^2&n&1\\
p^3&p^2&p&1\\
q^3&q^2&q&1\\
\end{bmatrix}
\begin{bmatrix}A\\B\\C\\D
\end{bmatrix}
&=\begin{bmatrix}m\\n\\p\\q\\1\\-1\\\sqrt{2}\\-\sqrt{2}\end{bmatrix}\\
\end{aligned}$$
You could continue to investigate this system and see what conditions on $\{m,n,p,q\}$ make the system consistent, and then see what solution sets for $\{A,B,C,D\}$ there are. That could get messy though since the equations are not linear in $\{m,n,p,q\}$.
The above system simplifies to
$$\begin{aligned}
\begin{bmatrix}
m^3&m^2&m&1\\
n^3&n^2&n&1\\
p^3&p^2&p&1\\
q^3&q^2&q&1\\
\end{bmatrix}
\begin{bmatrix}
-1/2&1/2 &\sqrt{2}/4&-\sqrt{2}/4\\
-1/2&-1/2&1/2&1/2\\
1&-1&-\sqrt{2}/4&\sqrt{2}/4\\
1&1&-1/2&-1/2
\end{bmatrix}
\begin{bmatrix}m\\n\\p\\q
\end{bmatrix}
&=\begin{bmatrix}1\\-1\\\sqrt{2}\\-\sqrt{2}\end{bmatrix}\\
\end{aligned}$$
$$\begin{aligned}
\begin{bmatrix}A\\B\\C\\D
\end{bmatrix}
&=\begin{bmatrix}
-1/2&1/2 &\sqrt{2}/4&-\sqrt{2}/4\\
-1/2&-1/2&1/2&1/2\\
1&-1&-\sqrt{2}/4&\sqrt{2}/4\\
1&1&-1/2&-1/2
\end{bmatrix}\begin{bmatrix}m\\n\\p\\q
\end{bmatrix}
\end{aligned}$$
So there is a system of quartic equations in $m,n,p,q$ that once solved, provides solutions for $A,B,C,D$.
Having four polynomial equations in four variables $m,n,p,q$, one would expect only a finite number of solutions. We've already accounted for $10$ solutions; it's possible that there are no more. Assuming these four equations have the appropriate algebraic independence, Bézout's Theorem says there are at most 64 complex solutions.