Consider a $W \sim \text{Pareto}(\beta)$ and let $\mathbf{X}=(X_1,X_2)^T$ be a bivariate vector whose components are given by:
$$X_1 = \sqrt{W}(Z_1+Z_2) \quad \text{and} \quad X_2 = \sqrt{W}(Z_1-Z_2)$$
where $Z_1$ and $Z_2$ iid standard normal, independent of $W$.
Now since it is obvious that the margins of $\mathbf{X}$ are not independent, I would like to show that formally. Natural way of showing this would be to show that $\mathrm{E}\left[X_1X_2\right] \neq \mathrm{E}[X_1]\mathrm{E}[X_2]$, however this approach doesn't produce anything useful.
How does one go about proving random variables are $\textit{not}$ independent in general?
Showing that their joint PDF/CDF do not factorize into a product form seems difficult in general.