原子間ポテンシャル（げんしかんポテンシャル、英: Interatomic potential）は、与えられた原子の空間位置から原子系のポテンシャルエネルギーを計算するための関数である[1][2][3][4]。原子間ポテンシャルは化学や分子物理学、材料物理学の分野で、凝集 (en:英語版) 、熱膨張、材料の弾性特性などの効果に関する分子力学法や分子動力学法 (MD) によるシミュレーションの物理的基盤として広く用いられる[5][6][7][8][9][10]。
関数形

原子間ポテンシャルは1つの原子位置の関数、もしくは2つ、3つ…の原子位置の組の関数を項とする級数展開で表すことができる。このとき系のポテンシャルの総和 VTOT は以下のように書ける[3]。 V T O T = ∑ i N V 1 ( r → i ) + ∑ i , j N V 2 ( r → i , r → j ) + ∑ i , j , k N V 3 ( r → i , r → j , r → k ) + ⋯ {\displaystyle V_{\mathrm {TOT} }=\sum _{i}^{N}V_{1}({\vec {r}}_{i})+\sum _{i,j}^{N}V_{2}({\vec {r}}_{i},{\vec {r}}_{j})+\sum _{i,j,k}^{N}V_{3}({\vec {r}}_{i},{\vec {r}}_{j},{\vec {r}}_{k})+\cdots }

ここで V1 は一体項、V2 は二体項、V3 は三体項を表し、N は系に含まれる原子の数、 r → i {\displaystyle {\vec {r}}_{i}} は i 番目の原子の位置である。指標 i, j, k …の総和はすべての原子（の組）について行う。

対ポテンシャルが原子対ごとに与えられているなら、それに1/2をかけたものが級数展開の二体項となることに注意が必要である。さもなければ一つの対のポテンシャルが2回数えられてしまう。同じように三体項には1/6の係数がかかる[3]。あるいは、ポテンシャルの関数形が指標 i, j, k …の交換に対して対称である場合には、二体項の総和を i < j の場合のみに、三体項の総和を i < j < k の場合のみに限定する方法もある（原子の種類が複数ある場合にはこのような対称性が存在しないこともある）。

一体項は原子が外場（電場など）の中にある場合にのみ意味を持つ。外場がなければ、ポテンシャル V は個々の原子の絶対位置ではなく原子間の相対位置にのみ依存するはずである。それはつまり、ポテンシャルの表式を原子間距離 r i j = 。 r → i − r → j 。 {\displaystyle \textstyle r_{ij}=|{\vec {r}}_{i}-{\vec {r}}_{j}|} および結合角（ある原子から隣接する複数の原子に向けて引いたベクトルの間の角度）θijk の関数として書き換えられるということである。よって外力がない場合の一般形は以下のようになる。 V T O T = ∑ i , j N V 2 ( r i j ) + ∑ i , j , k N V 3 ( r i j , r i k , θ i j k ) + ⋯ {\displaystyle V_{\mathrm {TOT} }=\sum _{i,j}^{N}V_{2}(r_{ij})+\sum _{i,j,k}^{N}V_{3}(r_{ij},r_{ik},\theta _{ijk})+\cdots }

3次元空間においてi, j, kの3原子の相対位置を確定するには rij, rjk, θijk の3項のみで十分であるため、上式の三体項 V3 に原子間距離 r j k {\displaystyle \textstyle r_{jk}} は含まれない。3次以上の全ての項は多体ポテンシャルとも呼ばれる。ある種の原子間ポテンシャルでは対アポテンシャルの中に多体相互作用が埋め込まれている（埋め込み原子型ポテンシャルおよび結合次数ポテンシャルについては後に論じる）。

原理的には、これらの数式に含まれる総和は N 個の原子すべてについて取る。しかし、原子間相互作用が及ぶ範囲が有限であるならば（すなわち、あるカットオフ距離 rcut より大きい r に対してポテンシャルが V(r) ≡ 0 なら）、総和はカットオフ距離より近い原子の組だけに限られる。またセル法によって隣接原子を選び出すことで[1]、MDシミュレーションの計算アルゴリズムのオーダーを O(N) にできる。