Take a sample $X_1 \ldots X_n$ of $n$ independent observations $X_j \in \mathbb{R}$ with zero mean and finite variances $\sigma_j^2$. For $i = 1, 2, \ldots$, define the sums $$S^n_i = \frac{\pm X_1 \pm X_2 \ldots \pm X_n}{s_n} \;,$$ where different $i$'s have different sets of $\pm$ and $s_n^2 = \sum_{j=1}^{n} \sigma_j^2$ is the sum of variances. Assuming suitable conditions on the $X_i$, we have the CLT $$S^n_i \rightsquigarrow N(0,1)$$ for all individual $i = 1 \ldots n$.

**Q:** *What is known about the asymptotic joint distribution of $S^n_i$? Are they jointly normal?*

In particular, let $n$ be a number for which the Hadamard matrix $H_n$ exists, and $X = (X_1 \ldots X_n)$. If the distribution were jointly normal, $\frac{H_n X}{s_n}$ would tend towards the multivariate normal $N(0, I_n)$. That's why I suspect that the distribution is either not jointly normal, or the rate of convergence breaks down.

**Edit:** We can reduce the question to the most basic case. Let $$S^n_1 = \frac{1}{s_n}\sum_{j=1}^{n} X_j \;, \\ S^n_2 = \frac{1}{s_n} \sum_{j=1}^{n} (-1)^j X_j \;,$$ what is the joint limiting distribution of $S^n_1$ and $S^n_2$?