粒子法の概観

分子動力学法から SMC サンプラーまで

Particles
Opinion
Author

司馬 博文

Published

4/07/2024

概要

数値実験と LLM とはいずれもシミュレーションに使えるが,用いる形式が違う(数字と文字).これにより,物理的な用途と社会的な用途とに別れている.この形式の違いを超克するのが機械学習の悲願であるとするならば,計算とはなんだろうか? Monte Carlo 法とはシミュレーションにより計算を実行する手段である.

1 導入

分子動力学法 (MD: molecular dynamics) の文脈で 粒子法 (particle methods) と言ったとき,

  • 古典的な質点には Newton 力学
  • 剛体で体積を持つ場合は Euler 方程式
  • 内部構造を持つ場合は Langrange 方程式
  • さらに一般の場合は Hamilton 方程式

に基づいて粒子の動きをシミュレーションすることを指し,多くの場合は Monte Carlo 法と組み合わせることで,粒子の動きを提案し,Metropolis 法によって補正することで,系の物理量の分布を求めたり,平均値を求めたりする営みを指す (Rapaport, 2004, p. 4), (Griebel et al., 2007, p. 17)

そもそも MCMC の手法 (Metropolis et al., 1953) 自体も,Schrödinger 方程式に基づいた基底状態のシミュレーションを実行するために生まれたものである.

これが,強力なのである.極めて強力であるが故に,計算科学の分野ではすでに,そして将来的には確実に,シミュレーションと数値実験が,実際の実験と似たような重要な役割を科学において担うことになると見られている.

この意味での粒子法と,統計学において粒子法と言った場合に想像されるであろう SMC や IPM (interacting particle methods) とは 実は数理的には同一物である ことを見たい.

If, in some cataclysm, all of scientific knowledge were to be destroyed, and only one sentence passed on to the next generations of creatures, what statement would contain the most information in the fewest words? I believe it is the atomic hypothesis (or the atomic fact, or whatever you wish to call it) that all things are made of atoms—little particles that move around in perpetual motion, attracting each other when they are a little distance apart, but repelling upon being squeezed into one another. In that one sentence, you will see, there is an enormous amount of information about the world, if just a little imagination and thinking are applied. (Feynman et al., 1964)

一言で言えば,一般に粒子法はポテンシャル(ハミルトニアン)に従って粒子を伝播する方法であり,粒子フィルターは特にこれが Feynman-Kac 測度で与えられる場合に用いられる粒子法である.

References

Feynman, R. P., Leighton, R. B., and Sands, M. (1964). The feynman lectures of physics. In,Vol. I. Addison-Wesley.
Griebel, M., Zumbusch, G., and Knapek, S. (2007). Numerical simulation in molecular dynamics: Numerics, algorithms, parallelization, applications,Vol. 5. Springer Belin, Heidelberg.
Metropolis, N., Rosenbluth, A. W., Rosenbluth, M. N., Teller, A. H., and Teller, E. (1953). Equation of state calculations by fast computing machines. The Journal of Chemical Physics, 21(6), 1087–1092.
Rapaport, D. C. (2004). The art of molecular dynamics simulation. Cambridge University Press.