【物理シミュレーションに挑戦！】古典力学
様々な力による物理シミュレーション３：クーロン力による運動

本項は物理学における様々な系の物理現象を、解析的に扱うのではなく数値計算による物理シミュレーションを実行するために必要な計算アルゴリズムを示すことを目的とします。 様々な初期条件に対する物理シミュレーションを実現するために、最も汎用的な直交座標系を用います。
・様々な力による物理シミュレーション１：弾性力による運動
・様々な力による物理シミュレーション２：万有引力による運動
・様々な力による物理シミュレーション３：クーロン力による運動
・様々な力による物理シミュレーション４：レナード・ジョーンズ相互作用による運動
・様々な力による物理シミュレーション５：引力＋衝突力

電荷をもつ２つの物体にはクーロン力と呼ばれる引力ないし斥力が働きます。本項では、クーロン力による運動シミュレーションするために必要な計算アルゴリズムを示します。

図は２つの電荷に働くクーロン力の模式図です。クーロン力は２つの電荷の大きさに比例し、２電荷の距離の２乗に反比例し、力の向きは２点間を結ぶ直線上に働きます。 なお、電荷の符号が同じ場合には斥力、異なる場合には斥力となります。以上の結果から、m番目の質点からn番目の質点に与えられるクーロン力は

と与えられます。ただし、 は真空の誘電率と呼ばれる定数です。また、力の向きを与える単位ベクトル は

と定義します。もし、複数の電荷が互いにクーロン力を及ぼし合う場合、n番目の質点に加わる力は各万有引力を足し合わせた

で与えられます。そのため、n番目の質点の時間発展は運動方程式から

となります。

運動エネルギーとポテンシャルエネルギー

系全体の運動エネルギーとポテンシャルエネルギーはそれぞれ次のとおり与えられます。

両者の和系全体の運動エネルギー

は時間に依存しません。

クーロン力にによる運動シミュレーション

赤球：陽子（電荷：（電気素量）、質量：）
青球：電子（電荷：、質量：）
なお、エネルギーの単位は電子ボルト（）です。

古典的系における原子のモデル

電子間相互作用を考慮した古典的原子モデルです。クーロン力だけでは電子の電荷と原子の電荷が釣り合っている場合でも安定しません。

物理シミュレーションについては「HTML5による物理シミュレーション」を参照ください。 数式の表示は「Tex表記によるHTML文書への式の埋め込み」をご覧ください。