Of course $X=0$ works, but I am looking for a non-singular solution. I haven't made much progress to solve this problem. However, let $\mu_2 = E(X^2)$ and $\mu_1 = E(X)$. For the equality to hold, we must have $\mu_2 = \mu_1 (2-\mu_1) > 0$. This is clearly impossible if $\mu_1 < 0$, thus we can focus on the case $0< \mu_1 < 1$.
Background
Let's $X_1,X_2,X_3$ and so on be i.i.d. with same distribution as $X$. Let's define the following infinite sums: $$Z = X_1 + X_1 X_2 + X_1 X_2 X_3 +\cdots \\ Y=X_1 + X_2 X_3 + X_4 X_5 X_6 +\cdots$$ We have (see here why):
$$ \mbox{Var}(Z) = \frac{\mbox{Var}(X)}{(1-\mu_1)^2(1-\mu_2)} , \mbox{ Var}(Y)=\frac{\mbox{Var}(X)}{(1-\mu_1^2)(1-\mu_2)}$$
Let's use the symbols $\mu$ and $\sigma^2$ to denote the expectation and variance of the infinite sum, regardless as to whether it comes from model $Z$, or from model $Y$. To test whether a data fits the model $Z$ or $Y$, the statistic of the test is
$$T = \sigma^2\cdot\frac{(1-\mu_2)(1-\mu_1^2)}{\mu_2-\mu_1^2}$$ Here $\sigma^2$ is the empirical variance computed on the observations modeled by the infinite series ($Z$ or $Y$ depending on the model). $T$ is expected to be equal to $1$ if the data matches model $Y$. But both models result in the same $T$ only if $(1-\mu_1^2) = (1-\mu_1)^2$. Note that $\mu_X = \mu_1$ and $\sigma_X^2$ are easy to estimate, using some formulas, for instance $\mu_X = \mu_1 = \mu/(1+\mu)$, valid for both models. Also:
- For model $Z$:
$$\sigma_X^2 = \frac{(1-\mu_1)^2(1-\mu_1^2)\sigma^2}{1+\sigma^2(1-\mu_1)^2} $$
- For model $Y$:
$$\sigma_X^2 = \frac{(1-\mu_1^2)^2\sigma^2}{1+\sigma^2(1-\mu_1^2)}$$
Note: If correct, it would imply that $\mu > -\frac{1}{2}$ in all cases where convergence (for the infinite sum) occurs, whether you use model $Z$ or $Y$. Also, if $\mu_X = 0$ then $\mu =0$ (the converse is also true) and $\sigma_X^2 = \sigma^2/(1+\sigma^2)$ regardless of the model.
I realize that I posted the wrong question, due to a typo when copy/pasting a formula. It should have been "can we have $1-E^2(X) =(1-E(X))^2$ which has the obvious answer "yes only if $E(X) = 0$" (since the case $E(X) =1$ must be excluded.)The issue is still the same, that is, getting a statistical test that can discriminate between model $Y$ and model $Z$, and the answers posted by @knRumsey and @Henry to my question are correct, it's just that I posted the wrong question. Not sure how to best handle this. It definitely makes my problem easier, but I need somehow to update my question.