[論文レビュー] Coarse-graining schemes and a posteriori error estimates for stochastic lattice systems

Markos A. Katsoulakis, Petr Plecháč+2 more|ArXiv.org|Aug 1, 2006

Advanced Mathematical Modeling in Engineering参考文献 22被引用数 4

ひとこと要約

本稿では、初期のモンテカルロで計算可能な近似の周囲におけるクラスターエクスパンションを用いて、確率的格子系のための適応的粗粒度化スキームを開発する。相対エントロピーに基づく事後誤差推定が導出され、中程度および長距離相互作用を有する系における臨界的挙動やヒステリシスの予測を向上させ、短距離相互作用においても有効ハミルトニアンの高次補正を組み込むことで、従来の手法に比べて顕著に精度が向上する。

ABSTRACT

The primary objective of this work is to develop coarse-graining schemes for stochastic many-body microscopic models and quantify their effectiveness in terms of a priori and a posteriori error analysis. In this paper we focus on stochastic lattice systems of interacting particles at equilibrium. %such as Ising-type models. The proposed algorithms are derived from an initial coarse-grained approximation that is directly computable by Monte Carlo simulations, and the corresponding numerical error is calculated using the specific relative entropy between the exact and approximate coarse-grained equilibrium measures. Subsequently we carry out a cluster expansion around this first-and often inadequate-approximation and obtain more accurate coarse-graining schemes. The cluster expansions yield also sharp a posteriori error estimates for the coarse-grained approximations that can be used for the construction of adaptive coarse-graining methods. We present a number of numerical examples that demonstrate that the coarse-graining schemes developed here allow for accurate predictions of critical behavior and hysteresis in systems with intermediate and long-range interactions. We also present examples where they substantially improve predictions of earlier coarse-graining schemes for short-range interactions.

研究の動機と目的

• 平衡状態における確率的多体系の計算的に効率的な粗粒度化スキームの開発。
• 相対エントロピーに基づく事後誤差推定を用いて、粗粒度化近似の精度を定量化すること。
• 一次近似を超えるために、クラスターエクスパンションを用いて有効ハミルトニアンの高次補正を適用すること。
• 鋭い、計算可能な誤差境界を提供することにより、粗粒度化の適応的改善を可能にすること。
• 長距離および中距離相互作用を有する系における相転移やヒステリシスといった臨界現象を正確に予測すること。

提案手法

• 最初の一次粗粒度化ハミルトニアン \bar{H}^{(0)}_M(\eta) を、粗格子内の磁化 \eta が与えられたもとでの微視的状態の平均をとることで、モンテカルロシミュレーションにより計算可能とする。
• 条件付き確率測度 PN(dσ|η) を用いて、正確な粗粒度化ハミルトニアン \bar{H}_M(\eta) を \bar{H}^{(0)}_M(\eta) の摂動として再定式化する。
• ポリマー・モデルに基づくクラスターエクスパンション技法を用い、級数展開を導出する： \bar{H}_M(\eta) = \bar{H}^{(0)}_M(\eta) + \bar{H}^{(1)}_M(\eta) + \cdots + \bar{H}^{(p)}_M(\eta) + O(\epsilon^{p+1})。
• 粗粒度化比 q、相互作用範囲 L、逆温度 β に依存する小さなパラメータ \epsilon を収束性および誤差推定の制御に用いる。
• 正確なおよび近似された粗粒度化ギブス測度間の相対エントロピーを用いて、鋭い事後誤差推定を導出する。
• 誤差推定を基に、特に誤差が大きく、臨界付近にある領域での粗粒度化スキームの適応的改善を実現する。

実験結果

リサーチクエスチョン

• RQ1 確率的格子系の一次粗粒度化近似を摂動補正を用いて体系的に改善する方法は何か？
• RQ2 粗粒度化の文脈において、クラスターエクスパンションが有効ハミルトニアンの高次補正を導出する役割を果たす仕組みは何か？
• RQ3 相対エントロピーに基づく事後誤差推定を用いて、保証された誤差許容範囲を持つ適応的粗粒度化手法を構築できるか？
• RQ4 高次補正は、短距離、中距離、長距離相互作用を有する系における臨界的挙動およびヒステリシスの予測にどのように影響を与えるか？
• RQ5 導出されたスキームは、既存の粗粒度化手法に比べて、相転移や準安定状態をどれほどよく捉えられるか？

主な発見

• 最近接相互作用を有するイジング模型において、有効ハミルトニアンに三次補正を組み込むことで、高い粗粒度化比であっても誤ったヒステリシスが完全に除去される。
• 中距離相互作用（例：L = 8, q = 8）の状況では、二次近似では正しく相転移を捉えられないが、三次補正によりヒステリシスおよび臨界的挙動の正確な予測が回復する。
• 長距離相互作用領域（例：L = 32, q = 32）では、二次近似でも平均場解と良好に一致し、三次近似ではほぼ同一の結果が得られる。
• 小さなパラメータ \epsilon = q/L が小さくない場合（例：q = 32, L = 8）でも、高次補正により顕著に精度が向上し、漸近的極限の外でも頑健であることが示された。
• 相対エントロピーから導出された事後誤差推定は、計算可能であり、誤差が大きい領域を的確に特定することができ、適応的改善戦略の実現に有効である。
• 数値シミュレーションでは、正確な解と完全に解像度の高い微視的シミュレーションとの間で視覚的に区別がつかない（q = 1）ことが確認され、全テスト粗粒度化レベルにおいて粗粒度化近似の正確性が裏付けられた。