Are decompositions of a random variable into a sum of two IID random variables unique?

Let $Z$ be a real-valued random variable, and suppose that $Z = X_1 + X_2$ where $X_1$ and $X_2$ are i.i.d. random variables. Suppose further that $Z = Y_1 + Y_2$ where $Y_1$ and $Y_2$ i.i.d. random variables. Does it follow that $Y_1$ and $X_1$ are identically distributed?

I've looked into indecomposable distributions and infinitely divisible distributions, but could not find a result/example immediately answering the above question.

The answer is no in general.

In the language of characteristic functions, you have $\phi_Z(t)=\phi^2_X(t)=\phi_Y^2(t)$. Taking the square root implies $\phi_Y(t)=\pm\phi_X(t)$ for every $t$, in particular with $\phi_X(t)=\phi_Y(t)$ for $t$ in a neighborhood of 0 (since $\phi(0)=1$ for any characteristic function and by uniform continuity).

Thus, this is true if you require that $\phi_X$ and $\phi_Y$ are analytic, as they would necessarily agree on an interval, and therefore everywhere.

Otherwise it is not true, meaning that $|\phi_X(t)|=|\phi_Y(t)|$, such that there are $t$ where $\phi_X(t)\neq \phi_Y(t)$. A classic example of this is as follows:

Let $P(X=x)=\frac{2}{\pi^2(2k-1)^2}$ whenever $x=\pm (2k-1)\pi$ for $k=1,2\cdots$ , $P(X=0)=1/2$, and $P(X=x)=0$ otherwise.

Let $P(Y=y)=\frac{4}{\pi^2(2k-1)^2}$ whenever $y=\pm (2k-1)\pi/2$ and $k=1,2,\cdots$ and $P(Y=y)=0$ otherwise.

Both these give rise to "triangle-wave" characteristic functions that are periodic outside the interval specified:

$$\phi_X(t)=1-|t|, t\in[-1,1], \mbox{ periodic otherwise}$$ $$\phi_Y(t)=1-|t|, t\in[-2,2], \mbox{ periodic otherwise},$$

which look like: