エバルトの方法

エバルト法：クーロン相互作用の効率的な計算手法

エバルト法は、計算物理学、特に分子動力学シミュレーションや第一原理バンド計算において、周期系における原子核やイオン間のクーロン相互作用を効率的に計算するための重要な手法です。クーロン力は距離の逆二乗に比例して働くため、単純に全ての粒子間の相互作用を計算しようとすると、計算量は粒子の数にほぼ比例して増加し、非常に多くの計算時間が必要になります。エバルト法はこの問題を解決するために開発されました。

エバルト法の原理

エバルト法では、クーロン相互作用を、実空間で短距離相互作用、逆格子空間で長距離相互作用という２つの部分に分割して計算します。それぞれの空間での計算は高速に収束するため、計算時間が大幅に短縮されます。

具体的には、単位胞内のイオンが作る静電ポテンシャルを考えます。このポテンシャルは、実空間でのイオン配置と逆格子空間でのフーリエ級数展開の両方で表すことができます。しかし、これらの表現はどちらも無限和であり、そのまま計算すると発散してしまいます。そこでエバルト法では、短距離相互作用と長距離相互作用を分割するために、適切な関数(例えば相補誤差関数)を用いてポテンシャルを２つの項に分割します。

短距離相互作用: 実空間において、短距離で急速に減衰する関数を使うことで、近傍のイオン間の相互作用を効率的に計算します。

長距離相互作用: フーリエ変換を用いて逆格子空間で計算します。この方法では、長距離の相互作用を効率的に計算することができます。

これらの２つの部分の和が、最終的な静電ポテンシャルとなります。分割されたそれぞれの項は、適切なパラメータを選択することで、高速に収束し、有限の計算で十分な精度が得られるようになります。

粒子メッシュエバルト法 (PME)

エバルト法は、計算機の登場以前から理論物理学で用いられていましたが、1970年代以降、コンピュータシミュレーションにおいて広く利用されるようになりました。特に、粒子メッシュエバルト法(PME)は、エバルト法を改良した手法として知られています。

PME法では、電荷密度をメッシュ上に表現することで高速フーリエ変換(FFT)を利用し、長距離相互作用の計算をさらに効率化します。FFTは、フーリエ変換を高速に実行するアルゴリズムであり、多くの粒子を含む系でも効率的に計算が可能です。また、PME法では、短距離相互作用の計算は、直接計算を行うか、あるいは高速多重極法などの高速計算アルゴリズムを用いることで効率化が図られます。

PME法は、分子動力学シミュレーションだけでなく、プラズマ、銀河、生物高分子のシミュレーションなど、様々な分野で応用されています。 特に、レナード・ジョーンズポテンシャルなど、短距離相互作用を含む系にも適用されるようになっています。

エバルト法の適用条件

エバルト法、そしてPME法は、系の周期性を仮定しています。つまり、無限に周期的に繰り返される系を扱うのに適しています。分子動力学シミュレーションでは、周期境界条件を用いることで、この条件を満たすことができます。しかし、系のサイズが小さい場合や、表面効果が重要な場合は、エバルト法の精度が低下する可能性があります。また、電荷密度の揺らぎが大きい系では、高速多重極法などの他の手法の方が効率的である場合があります。

まとめ

エバルト法は、周期的な系におけるクーロン相互作用を効率的に計算するための強力な手法です。PME法などの改良版も開発され、計算コストを削減しつつ高精度な計算を可能にしています。様々なシミュレーションにおいて、不可欠な手法となっています。

もう一度検索