[論文レビュー] Finite Gelfand Pair Approaches for Ehrenfest Diffusion Model

Hiroshi Mizukawa|arXiv (Cornell University)|Sep 7, 2010

Theoretical and Computational Physics参考文献 2被引用数 6

ひとこと要約

本稿は、複素反射群 $G(r,1,n)$ および対称群 $S_n$ を用いた有限ゲルファンド対に基づき、古典的なエーレンフェスト拡散モデルを多 urn システムに一般化する。zonal spherical 関数およびこれらの群上のフーリエ解析を用いて、配置空間 $B(r,n)$ 上のランダムウォークの混合時間と極限分布を分析し、$r > 2$ の場合に一様分布に収束すること、$r=2$ の場合に周期的挙動を示し、偶数または奇数ステップに応じて収束が異なることを示している。

ABSTRACT

A classical diffusion model of Ehrenfest which consists of $2$-urns and $n$-balls is realized by a finite Gelfand pair $(H_{n},S_{n})$, where $H_{n}$ is the hyperoctahedral group and $S_{n}$ is the symmetric group. This fact can be generalized to multi-urn version by using Gelfand pairs of complex reflection groups.

研究の動機と目的

• 2 つの urn に限定された古典的エーレンフェスト拡散モデルを、$r \geq 2$ 個の urn に拡張する。
• $B(r,n)$、すなわち $r$ 個の urn 間の玉の配置空間上での確率過程を、$G(r,1,n)$ の作用に関して不変なランダムウォークとしてモデル化する。
• 有限ゲルファンド対上の調和解析を用いて、過程の定常分布への収束挙動を分析する。
• 特に $N$ が大きい場合に、定常分布への全 Variation 距離の明示的評価を導出する。
• 巡回的および双方向移動などの代替シャッフルルールに対応する遷移カーネルのフーリエ係数を計算する。

提案手法

• 配置空間 $B(r,n)$ を、$G(r,1,n)$ と $S_n$ の商空間 $G(r,1,n)/S_n$ として実現する。ここで $G(r,1,n)$ は複素反射群 $\mathbb{Z}/r\mathbb{Z} \wr S_n$ である。
• ゲルファンド対 $(G(r,1,n), S_n)$ を用いて、置換表現を、$\mathbf{k} \in N(r,n)$ でインデックス付けられた既約表現に分解する。
• zonal spherical 関数 $\omega_{\mathbf{k}}(x)$ を、単項対称多項式を単位根で評価することで表現する：$\omega_{\mathbf{k}}(x) = \frac{m_{\lambda(\mathbf{k})}(\xi^{x_1}, \dots, \xi^{x_n})}{\binom{n}{k_0, \dots, k_{r-1}}}$。
• 遷移カーネル $\tilde{\nu}$ を zonal spherical 関数の線形結合として表現する：$\tilde{\nu} = \sum_{\mathbf{k}} a_{\mathbf{k}} \omega_{\mathbf{k}}$。これによりスペクトル解析が可能になる。
• 正規直交性および球関数の減衰見積もりを用いて、$\tilde{\nu}^{*N}$ の $N$ 重畳みを計算する。
• 二項係数および指数的減衰項を含む和の上限を用いて、全 Variation 距離 $||\nu^{*N} - \pi||_{TV}^2$ を分析する。

実験結果

リサーチクエスチョン

• RQ1 多 urn エーレンフェスト拡散モデルの混合時間は、$n$ および $r$ に対してどのようにスケーリングされるか？
• RQ2 極限 $N \to \infty$ における過程の極限分布は、$r > 2$ および $r = 2$ の場合にそれぞれどうなるか？
• RQ3 巡回的移動や左右隣接への移動といった異なるシャッフルルールは、収束速度にどのように影響するか？
• RQ4 群表現論を用いて、遷移カーネルのフーリエ係数を明示的に計算できるか？
• RQ5 全 Variation 距離の上界を $n$、$r$、$N$ の関数として導出できるか？

主な発見

• $r > 2$ の場合、$N \to \infty$ のとき、過程は $B(r,n)$ 上の一様分布 $\pi$ に収束する。すなわち $\lim_{N \to \infty} \tilde{\nu}^{*N} = \frac{1}{r^n n!}$。
• $r = 2$ の場合、過程は周期 2 の周期的挙動を示す：$\lim_{N \to \infty} \tilde{\nu}^{*2N}$ は偶数パリティの配置上で一様であり、$\lim_{N \to \infty} \tilde{\nu}^{*2N-1}$ は奇数パリティの配置上で一様である。
• 全 Variation 距離に対する上界が導出される：$r \geq 3$ の場合、$||\nu^{*N} - \pi||_{TV}^2 \leq \frac{1}{4} e^{-c}$ であり、$r = 2$ の場合、$||\nu^{*N} - \pi||_{TV}^2 - \frac{1}{4} \leq \frac{1}{4} e^{-c}$ が成り立つ。ここで $N = \frac{n(r-1)}{2r}(\log r^n + c)$ を仮定する。
• 巡回シャッフル (1) のフーリエ係数は $f_k = \sum_{i=0}^{r-1} \frac{k_i}{n} \zeta^{-i}$ であり、ここで $\zeta$ は原始 $r$ 乗根である。
• 双方向巡回シャッフル (2) のフーリエ係数は $f_k = \frac{1}{r-1}\left(\frac{r k_0}{n} - 1\right)$ であり、これは urn 0 に含まれる玉の数に線形に依存する。
• $r \geq 3$ の場合、全 Variation 上界は $N$ に対して指数的に減少し、減少率は $\frac{2r}{n(r-1)}$ である。この上界は、$r^n$ 要因の程度までタイトである。