Consider the triangular array $X\_{n,k}$ such that, for each $n>0$, the variables $(X\_{n,1},\cdots,X\_{n,n})$ have the following properties:


1. For any given $1 \le L \le n$, all
subsets of
$(X\_{n,1},\cdots,X\_{n,n})$ of size
$L$ have the same joint distribution
(even after applying an arbitrary permutation).
2. Each $X\_{n,k}$ has
zero mean, variance
$0<\sigma^2<\infty$, and third
absolute moment $0<\rho<\infty$
3. Each $(X\_{n,k\_1},X\_{n,k\_2})$ pair
($k\_1 \ne k\_2$) has covariance
$\frac{-C}{n}$ for some $0 < C < \infty$ (hence the variables are all negatively correlated, and the correlation tends to zero)


Let $S\_n = \sum\_{k=1}^{n} X\_{n,k}$. Is $\frac{S\_n}{\sqrt{\mathrm{Var}[S\_n]}}$ asymptotically distributed as $N(0,1)$? If so, is a convergence rate of $O(\frac{1}{\sqrt{n}})$ achieved (cf. Berry-Esseen theorem)? What about in the multivariate setting where each $X\_{n,k}$ is a random vector (and the relevant quantities above are replaced by vectors/matrices)?


If it helps, we can also assume that the $X\_{n,k}$ are uniformly bounded in $n$ and $k$ with probability one. Answers with further assumptions than the ones listed will also be appreciated.



Central Limit Theorem (and Berry-Esseen theorem) for non-independent variables
These types of conditions (exchangeability and $O(1/n)$ covariance) do not even ensure that $S\_n$ is approximately normal.


Ignore odd $n$. (I believe similar examples can be constructed for odd $n$.) Let $n=2m$. Let each $X\_{n,k} \in \lbrace -1,1 \rbrace$. Let each set of $m$ positive signs have probability ${2m \choose m}^{-1}$. Then the sum is always $0$. This satisfies conditions 1 and 2. $P(X\_{n,1} = X\_{n,2}) = {2m-2 \choose m}/{2m-1 \choose m} = \frac{m-1}{2m-1}$ so the covariance is $\frac{m-1}{2m-1} - \frac{m}{2m-1} = \frac{-1}{2m-1}$.


The point mass at $0$ is still normal, but degenerate. However, this example can be modified slightly so that the sum has positive standard deviation. With probability $\frac{1}{4m^2}$ let all $X\_{n,i}$ be equal ($\frac{1}{8m^2}$ chance of all $+1$, $\frac{1}{8m^2}$ all $-1$), and with probability $\frac{4m^2-1}{4m^2}$ let the positive indices be a uniformly random subset of size $m$. Then the covariance $\text{Cov}(X\_{n,1},X\_{n,2})$ is $\frac{-1}{n}$, and $P(S\_n =0) = \frac{n^2-1}{n^2}, P(S\_n = \pm n) = \frac{1}{2n^2}$ so $\text{Var}(S\_n) = 1$, and $\frac{S\_n}{\sqrt{\text{Var}(S\_n)}}$ is far from normal. 


