A partir de la independencia tenemos $\varphi_{X+Y+(X-Y)}=\varphi_{2X}$ y por lo tanto
$$
\varphi(2t) = \varphi(t)^2(\varphi(t)\varphi(-t)) = \varphi(t)^3\varphi(-t).
$$
Set $\gamma(t) = \varphi(t)/\varphi(-t)$, luego
$$
\gamma(2t) = \frac{\varphi(2t)}{\varphi(-2t)} = \frac{\varphi(t)^3\varphi (t)}{\varphi(-t)^3\varphi(t)} =\gamma(t)^2.
$$
De ello se desprende que $\gamma(t)=\gamma(t/2)^2$ y por inducción, $\gamma(t) = \gamma(t/2^n)^{2^n}$ para los números enteros no negativos $n$. Por otra parte, $\varphi'(0)=i\mathbb E[X]=0$$\varphi''(0)=i^2\mathbb E[X^2] = -1$, así como el $h\to0$ tenemos por expansiones de Taylor
$$
\gamma(h) = \frac{\varphi(h)}{\varphi(-h)} = \frac{1-\frac12h^2+o(h^2)}{1-\frac12h^2+o(h^2)} = 1 + o(h)^2.
$$
De ello se desprende que para los grandes $n$, $$\gamma(t) = (1+o(t^2/2^{2n})^{2^n}\stackrel{n\to\infty}\longrightarrow1,$$
y por lo $\gamma\equiv1$, produciendo $\varphi(t)=\varphi(-t)$. Por lo tanto,$\varphi(2t)=\varphi(t)^4$, por lo tanto $\varphi(t)=\varphi(t/2)^4$ y por inducción, $$\varphi(t)=\varphi(t/2^n)^{4^n},n\geqslant 1.$$
El uso de Taylor expansiones vemos que un gran $n$,
$$
\varphi(t) = \left(1-\frac12\left(\frac t{2^n}\right)^2+o\left(\left(\frac t{2^n}\right)^2\right)\right)^{4^n}\stackrel{n\to\infty}\longrightarrow \exp\left(-\frac12 t^2\right)
$$
y, por tanto, $X$ $Y$ tiene distribución normal estándar.