Recently I have come across a ranked version of Dirichlet distribution whose densitydefined on the open ranked simplex $\Delta^{n-1}_+$$\nabla^{n-1}_+ = \{\vec x \in \mathbb{R}^n:\sum_{k=1}^n x_k =1, x_1 \geq x_2 \geq \cdots \geq x_n > 0\}$, whose density is given by $$ p(\vec x) = B^{-1} \prod_{k=1}^n x_{(k)}^{\alpha_k -1},\quad \vec x \in \Delta^{n-1}_+, $$$$ p(\vec x) = B^{-1} \prod_{k=1}^n x_{k}^{\alpha_k -1},\quad \vec x \in \nabla^{n-1}_+, $$ where $B = \int_{\Delta^{n-1}} \prod_{k=1}^n x_{(k)}^{\alpha_k -1} dx$$B = \int_{\nabla^{n-1}} \prod_{k=1}^n x_{k}^{\alpha_k -1} dx$ is the normalizing factor, $x_{(k)}$ is the $k-$th largest among $x_1, x_2,\cdots, x_n$ and $\alpha_k$'s positive. (The only difference from the usual Dirichlet distribution is the usual domain $x_k$$\Delta^{n-1}_+$ being replaced by the rank-based one $x_{(k)}$ in pdf$\nabla^{n-1}_+$, and also a different normalizing factor.)
We are interested in under what conditions on $\alpha_k$'s will make the following statement holds true: Assume $X^{n}$ follows the ranked Dirichlet distribution in dimension $n$ with parameters $\{\alpha^n_k\}_{k=1}^n$, for any fixed small $\epsilon >0$, $$ \sup_{n} \left\{\frac{1}{n^2}\sum_{k=1}^{n-1} \sum_{l>k} \mathbb{P} \left( X^n_{(k)} - X^n_{(l)} < \frac{\epsilon}{n}\right) \right\} \leq f(\epsilon) $$$$ \sup_{n} \left\{\frac{1}{n^2}\sum_{k=1}^{n-1} \sum_{l>k} \mathbb{P} \left( X^n_{k} - X^n_{l} < \frac{\epsilon}{n}\right) \right\} \leq f(\epsilon) $$ for some function $f$ such that $\lim_{\epsilon \rightarrow 0} f(\epsilon) = 0$.
Intuitively, this says the averaging probability of gaps between coordinates being less than a dimension-scaled threshold is small. While it seems a pure calculus problem, the computation of the target probability becomes tricky once involved with this "rank" structure. I greatly appreciate any suggestions on how to bound this probability with some choice of $\alpha$ or any reference where you may have seen similar stuff.