エネルギー保存の法則

エネルギー保存の法則（エネルギーほぞんのほうそく、英: law of the conservation of energy）とは、「孤立系のエネルギーの総量は変化しない」という物理学における保存則の一つである。エネルギー保存則とも呼ばれる。 任意の異なる二つの状態について、それらのエネルギー総量の差がゼロであることをいう。 例えば、取り得る状態が全て分かっているとして、全部で 3 つの状態があったとき、それらの状態のエネルギーを A, B, C と表す。 エネルギー保存の法則が成り立つことは、それらの差について、 A − B = 0, B − C = 0, C − A = 0 が成り立っていることをいう。 時間が導入されている場合には、任意の時刻でエネルギー総量の時間変化量がゼロであることをいい、時間微分を用いて表現される。 エネルギー保存の法則は、物理学の様々な分野で扱われる。特に、熱力学におけるエネルギー保存の法則は熱力学第一法則 (英: first law of thermodynamics) と呼ばれ、熱力学の基本的な法則となっている。 熱力学第一法則は、熱力学において基本的な要請として認められるものであり、あるいは熱力学理論を構築する上で成立すべき定理の一つである。第一法則の成立を前提とする根拠は、一連の実験や観測事実のみに基づいており、この意味で第一法則はいわゆる経験則であるといえる。 一方でニュートン力学や量子力学など一般の力学において、エネルギー保存の法則は必ずしも前提とされない。 ルネ・デカルトやゴットフリート・ライプニッツが、それぞれの仕方でこれを主張し、それぞれの支持者によって議論が長年に渡り行われた。 19世紀の中ごろ、ユリウス・ロベルト・フォン・マイヤー、ジェームズ・プレスコット・ジュール、ヘルマン・フォン・ヘルムホルツらによって、「力学的・熱・化学・電気・光などのエネルギーは、それぞれの形態に移り変わるが、エネルギーの総和は変化しない（保存される）」と主張された。 20世紀にアルベルト・アインシュタインによって、質量とエネルギーの等価性という考え方が提唱され、別の形での保存が主張されたが、その有効性や有効範囲については、疑問視されることも多かった。 現在ではエネルギー保存の法則は、しばしば「最も基本的な物理法則の一つ」と考えられている。多くの物理学者が、自然はこの法則にしたがっているはずだ、と信じているのである。 ルネ・デカルトは、1644年に出版した自身の著作『哲学の原理』(Principia Philosophiæ ) で、宇宙においてquantitas motus（運動の量）の総和が保たれている、と主張した。 Deum esse primariam motus causam: et eandem semper motus quantitatem in niverso conservare. — Principia philosophiae, Pars secunda, 36（デカルト『哲学の原理』第二章 36） デカルトが主張した quantitas motus（quantity of motion, 運動の量）という概念は、現代の運動量とある程度似てはいるが、厳密には異なる概念である。デカルトは「質量」という概念を持っていなかったし、デカルトは速度の大きさだけを重視し、向きが変わることについては考慮していなかった。したがって、デカルトの quantitas motus を現代の運動量に対応する量と見なすことはできない。 ゴットフリート・ライプニッツは、運動の量というのを初めて数式で表現してみようと試みたが、デカルトとは異なって mv2 の総和が保存されている、と主張した。ライプニッツはこの量を vis viva（living force, 活力）と呼んだ。この vis viva という概念は、釣り合いなどの場面で想定される動きとしては見えない vis mortua（dead force, 潜在的な力）と対比しつつ置かれた概念である。 デカルトの考え方を支持する人々と、ライプニッツの考え方を支持する人々で議論が起こるようになった。これを「活力論争 (vis viva controversy)」という。議論は長年に渡って続いた。18世紀半ばになって、ジョゼフ＝ルイ・ラグランジュやジャン・ル・ロン・ダランベールらが、両概念の明確化を試み、それらを区別したことによって、ようやく論争は沈静化した。 1807年にトマス・ヤングは、vis viva という用語で表されていた運動の概念を、"energy" と呼んだ。energy はギリシア語の ἐνέργεια（羅: energeia, エネルゲイア）という語を基にした造語である。ギリシャ語のἐνέργεια (energeia ) というのは語の構成としては εν + εργον (en +ergon ) であり、εργον (ergon ) は「仕事」、εν— (en —) は「～の状態」という意味である。よって「仕事をしている状態」といったような意味である。アリストテレスの哲学において ἐνέργεια は、ものが持つ「可能態」の中から現実化された「現実態」を意味する。つまり、energy という用語を用いている背景には、眼には見えない「活力」が具体的な「仕事」に変化したのだ、という発想がある。 ヤングが energy という用語を用いたからといって、それが人々にすぐに用いられるようになったわけでもなく、人々の間に定着するようになったのは、あくまで後のことである。vis viva 相当の概念は、19世紀半ばでもしばしば、英語圏では "force" と呼ばれていたし、ドイツ語圏では „Kraft” と呼ばれていた。 現代的な意味で energy の語が用いられるようになったのは、ヤングより後のことで、1850年頃にウィリアム・トムソンによって kinetic energy（運動エネルギー）、1853年にウィリアム・ランキンによって potential energy（位置エネルギー）の語が定義された。 19世紀前半のドイツの自然哲学では、「破壊されることもなく、形態が様々に変換する根源的な何か」を „Kraft”（力）と呼んでいた。この自然哲学の概念は、現在の「エネルギー保存の法則」という概念の成立に大きな影響を与えている。 19世紀の中ごろ、ロベルト・マイヤー、ジェームズ・プレスコット・ジュール、ヘルマン・フォン・ヘルムホルツらが、それぞれ独立して「エネルギー保存の法則」という考え方に辿りついた。 マイヤーは、ドイツの医者で、船医としてジャワに行った時に熱量とエネルギーとの関係を考察するようになった。船が熱帯を航海すると水夫らの静脈の血液の赤みが増すことに気付き、気温が上昇したことで体温維持のために酸素が使われる量が減るのだ、と解釈した。 そして1842年、「熱」と「仕事」の関係に関する論文 „Bemerkung über die Kräfte der unbelebten Natur”を発表した。 ジュールは1843年に熱の仕事当量の測定を行い、その後も様々な方法で熱の仕事当量を計測した。 ヘルムホルツは、サディ・カルノーやエミール・クラペイロン、ジュールらの仕事について整理し、1847年に著した „Über die Erhaltung der Kraft”で様々な状況でエネルギー保存の法則が成り立つことを示した。 マイヤーやジュールが熱の仕事当量に関する考察をした頃は、1798年のベンジャミン・トンプソン（ランフォード伯）による指摘などがあったものの、アントワーヌ・ラヴォアジエとピエール＝シモン・ラプラスに始まるカロリック説が有力であり、熱は物質であると見なされ、熱は単独で保存されると考えられていた。そのため、熱が仕事に変わり得ることの発見とその事実の定量的評価をすることは、熱力学第一法則を構成する上で重要な仕事だった。 1850年、ルドルフ・クラウジウスは論文 „Über die bewegende Kraft der Wärme”の中で熱力学第一法則について完全な形で述べた。 1905年にアルベルト・アインシュタインは、Annus Mirabilis papers の一つの „Ist die Trägheit eines Körpers von seinem Energieinhalt abhängig?”において、質量とエネルギーが交換可能なのではないか、という提案を行った。これをきっかけとして、物理学が大きく変容していくことになった。「エネルギー」や「物質」という概念自体が大きく変わっていくことになったのである。 特殊相対性理論において、質量はエネルギーの一形態であり、E=mc² という式の関係が成り立っている。したがって相対論の立場では、エネルギー保存の法則は「質量を含めたエネルギーの総和が保存されている」という主張になる。 他の物理学の様々な主張同様に、このアインシュタインの主張も最初は受け入れられなかったり疑問視されたが、原子核反応や電子対生成などの実験において成立していることが確認されると、アインシュタインの考えが次第に受け入れられるようになっていった。 なおそれに伴って、「質量保存則は（厳密に言えば）成り立っていない」と考えられるようになった。特に、原子核反応を扱う場合においては、質量のエネルギーへの変換は無視できないほど大きく、質量は保存されていない、として計算するようになっている。 ただし、この法則を一応受け入れるとしても、一体どの程度まで受け入れてよいのかということについて見解はバラバラであった。例えばニールス・ボーアは、ベータ崩壊をエネルギー保存の法則が成立していない事例だと考えていた。 ただしそのような状況の中で、1932年にヴォルフガング・パウリとエンリコ・フェルミが、β崩壊の事例でも、仮にエネルギー保存の法則が成立していると仮定して計算したところ、中性の粒子が存在しているだろう、と予想することができた。彼らはその粒子の存在を主張したものの具体的な物証は無く、長らく認められなかったが、1956年になり実験によってその粒子（ニュートリノ）が確認された。この出来事によって、有効範囲については疑問視されることも多かったものの、エネルギー保存の法則が成り立つと仮定してみることが、科学的発見につながるひとつの指針にもなり得ることが知られるようになった。 1918年、エミー・ネーターは論文 „Invariante Variationsprobleme”を出版した。この論文の中で、ネーターが1915年に得た、今日ネーターの第一定理と呼ばれる定理の証明が与えられた。ネーターの定理から、作用積分が不変であるような無限小変換が存在する場合、系はその変換に対して対称であるという。このとき系の対称性に対応した量が保存する。特にエネルギー保存の法則は、時間の並進対称性に対応していることが知られる。 熱力学におけるエネルギー保存の法則は、熱力学第一法則である。熱力学第一法則は次のように表現される。 d U = δ Q − δ W . {\displaystyle dU=\delta Q-\delta W.} ここで dU は系の内部エネルギー U の変化量、δQ は系に与えられた熱量、δW は系から取り出された仕事を表す（d は完全微分を、δ は不完全微分（英語版）を表す）。仕事は熱力学的系に繋がっている力学的系へのエネルギーの移動を表し、熱はそれ以外の熱力学的系へのエネルギーの移動を表している。 熱力学第一法則は、エネルギーがひとりでに消えたり生じたりすることはない、という経験的事実を法則化したものであり、上述の定式化では、エネルギーの変化が熱と仕事の和として与えられることで表現されている。 熱力学において第一法則は、上式を満たす状態量（すなわち系に対する操作の方法や途中経過に依存しない量） U が存在することを主張する法則とみなされている。 古典力学におけるエネルギー保存の法則は、力学的エネルギー保存の法則と呼ばれる。力学的エネルギーは位置エネルギーと運動エネルギーに分類され、それらの和が一定であることをいう。 以下に一粒子系の場合についての力学的エネルギー保存の法則を述べる。 一粒子の運動について、粒子に働く力 F ( r ( t ) , t ) {\displaystyle {\boldsymbol {F}}({\boldsymbol {r}}(t),t)} がポテンシャル V ( r ( t ) ) {\displaystyle V({\boldsymbol {r}}(t))} を用いて、 F ( r ( t ) , t ) = − ∇ V ( r ( t ) ) + f ( t ) {\displaystyle {\boldsymbol {F}}({\boldsymbol {r}}(t),t)=-\nabla V({\boldsymbol {r}}(t))+{\boldsymbol {f}}(t)} と表される場合について、ニュートン力学の運動の第2法則、 m d 2 r d t 2 ( t ) = F ( r ( t ) , t ) {\displaystyle m{\frac {d^{2}{\boldsymbol {r}}}{dt^{2}}}(t)={\boldsymbol {F}}({\boldsymbol {r}}(t),t)} より次の運動方程式が得られる。 m d 2 r d t 2 ( t ) = − ∇ V ( r ( t ) ) + f ( t ) . {\displaystyle m{\frac {d^{2}{\boldsymbol {r}}}{dt^{2}}}(t)=-\nabla V({\boldsymbol {r}}(t))+{\boldsymbol {f}}(t).} ここで、 m {\displaystyle m} は質量、 r {\displaystyle {\boldsymbol {r}}} は粒子の位置、 t {\displaystyle t} は時刻をそれぞれ表し、ナブラ ∇ {\displaystyle \nabla } とポテンシャル V {\displaystyle V} の積 ∇ V {\displaystyle \nabla V} はポテンシャルの勾配を意味する。 ∇ V ( r ( t ) ) = ( ∂ ∂ x V ( r ( t ) ) , ∂ ∂ y V ( r ( t ) ) , ∂ ∂ z V ( r ( t ) ) ) T . {\displaystyle \nabla V({\boldsymbol {r}}(t))=\left({\frac {\partial }{\partial x}}V({\boldsymbol {r}}(t)),{\frac {\partial }{\partial y}}V({\boldsymbol {r}}(t)),{\frac {\partial }{\partial z}}V({\boldsymbol {r}}(t))\right)^{T}.} このとき仕事は、 W ( t 1 , t 2 ) = ∫ t 1 t 2 F ( r ( t ) , t ) ⋅ d r ( t ) = ∫ t 1 t 2 { − ∇ V ( r ( t ) ) + f ( t ) } ⋅ d r ( t ) {\displaystyle {\begin{aligned}W(t_{1},t_{2})&=\int _{t_{1}}^{t_{2}}{\boldsymbol {F}}({\boldsymbol {r}}(t),t)\cdot d{\boldsymbol {r}}(t)\\&=\int _{t_{1}}^{t_{2}}\left\{-\nabla V({\boldsymbol {r}}(t))+{\boldsymbol {f}}(t)\right\}\cdot d{\boldsymbol {r}}(t)\end{aligned}}} と r ( t ) {\displaystyle {\boldsymbol {r}}(t)} についての線積分で表される。ここで中黒 '・' はベクトルの内積を意味する。線積分を時間についての積分に直せば、 W ( t 1 , t 2 ) = − ∫ t 1 t 2 ∇ V ( r ( t ) ) ⋅ d r ( t ) + ∫ t 1 t 2 f ( t ) ⋅ d r ( t ) = − ∫ t 1 t 2 { ∂ V ( r ( t ) ) ∂ x d x ( t ) d t + ∂ V ( r ( t ) ) ∂ y d y ( t ) d t + ∂ V ( r ( t ) ) ∂ z d z ( t ) d t } d t + ∫ t 1 t 2 f ( t ) ⋅ d r ( t ) , {\displaystyle {\begin{aligned}W(t_{1},t_{2})&=-\int _{t_{1}}^{t_{2}}\nabla V({\boldsymbol {r}}(t))\cdot d{\boldsymbol {r}}(t)+\int _{t_{1}}^{t_{2}}{\boldsymbol {f}}(t)\cdot d{\boldsymbol {r}}(t)\\&=-\int _{t_{1}}^{t_{2}}\left\{{\frac {\partial V({\boldsymbol {r}}(t))}{\partial x}}{\frac {dx(t)}{dt}}+{\frac {\partial V({\boldsymbol {r}}(t))}{\partial y}}{\frac {dy(t)}{dt}}+{\frac {\partial V({\boldsymbol {r}}(t))}{\partial z}}{\frac {dz(t)}{dt}}\right\}dt+\int _{t_{1}}^{t_{2}}{\boldsymbol {f}}(t)\cdot d{\boldsymbol {r}}(t),\end{aligned}}} となるので、ポテンシャルの時間についての全微分、 d V d t ( r ( t ) ) = { d x ( t ) d t ∂ ∂ x + d y ( t ) d t ∂ ∂ y + d z ( t ) d t ∂ ∂ z } V ( r ( t ) ) {\displaystyle {\frac {dV}{dt}}({\boldsymbol {r}}(t))=\left\{{\frac {dx(t)}{dt}}{\frac {\partial }{\partial x}}+{\frac {dy(t)}{dt}}{\frac {\partial }{\partial y}}+{\frac {dz(t)}{dt}}{\frac {\partial }{\partial z}}\right\}V({\boldsymbol {r}}(t))} を用いて、 W ( t 1 , t 2 ) = − ∫ t 1 t 2 d V ( r ( t ) ) d t d t + ∫ t 1 t 2 f ( t ) ⋅ d r ( t ) = − { V ( r ( t 2 ) ) − V ( r ( t 1 ) ) } + ∫ t 1 t 2 f ( t ) ⋅ d r ( t ) {\displaystyle {\begin{aligned}W(t_{1},t_{2})&=-\int _{t_{1}}^{t_{2}}{\frac {dV({\boldsymbol {r}}(t))}{dt}}dt+\int _{t_{1}}^{t_{2}}{\boldsymbol {f}}(t)\cdot d{\boldsymbol {r}}(t)\\&=-\left\{V({\boldsymbol {r}}(t_{2}))-V({\boldsymbol {r}}(t_{1}))\right\}+\int _{t_{1}}^{t_{2}}{\boldsymbol {f}}(t)\cdot d{\boldsymbol {r}}(t)\end{aligned}}} と書ける。もし、粒子が受ける力がポテンシャルのみによる場合、 f ( t ) {\displaystyle {\boldsymbol {f}}(t)} は存在しないので、粒子に与えられた仕事 W ( t 1 , t 2 ) {\displaystyle W(t_{1},t_{2})} はポテンシャルの差 − { V ( r ( t 2 ) ) − V ( r ( t 1 ) ) } {\displaystyle -\left\{V({\boldsymbol {r}}(t_{2}))-V({\boldsymbol {r}}(t_{1}))\right\}} に等しい。このときポテンシャル　 V ( r ( t ) ) {\displaystyle V({\boldsymbol {r}}(t))} は位置エネルギーと呼ばれる。 再び仕事の定義に戻ると、粒子の運動方程式より、次のように書き換えられる。 W ( t 1 , t 2 ) = ∫ t 1 t 2 F ( r ( t ) , t ) ⋅ d r ( t ) = ∫ t 1 t 2 m d 2 d t 2 r ( t ) ⋅ d r ( t ) = m ∫ t 1 t 2 d 2 d t 2 r ( t ) ⋅ d r ( t ) d t d t {\displaystyle {\begin{aligned}W(t_{1},t_{2})&=\int _{t_{1}}^{t_{2}}{\boldsymbol {F}}({\boldsymbol {r}}(t),t)\cdot d{\boldsymbol {r}}(t)\\&=\int _{t_{1}}^{t_{2}}m{\frac {d^{2}}{dt^{2}}}{\boldsymbol {r}}(t)\cdot d{\boldsymbol {r}}(t)\\&=m\int _{t_{1}}^{t_{2}}{\frac {d^{2}}{dt^{2}}}{\boldsymbol {r}}(t)\cdot {\frac {d{\boldsymbol {r}}(t)}{dt}}dt\end{aligned}}} ここで、ベクトルの内積の微分について、 d d t ( u ( t ) ⋅ v ( t ) ) = ∑ α = x , y , z d d t ( u α ( t ) v α ( t ) ) = ∑ α = x , y , z ( d u α ( t ) d t v α ( t ) + u α ( t ) d v α ( t ) d t ) = d u ( t ) d t ⋅ v ( t ) + u ( t ) ⋅ d v ( t ) d t {\displaystyle {\begin{aligned}{\frac {d}{dt}}\left({\boldsymbol {u}}(t)\cdot {\boldsymbol {v}}(t)\right)&=\sum _{\alpha =x,y,z}{\frac {d}{dt}}\left(u_{\alpha }(t)v_{\alpha }(t)\right)\\&=\sum _{\alpha =x,y,z}\left({\frac {du_{\alpha }(t)}{dt}}v_{\alpha }(t)+u_{\alpha }(t){\frac {dv_{\alpha }(t)}{dt}}\right)\\&={\frac {d{\boldsymbol {u}}(t)}{dt}}\cdot {\boldsymbol {v}}(t)+{\boldsymbol {u}}(t)\cdot {\frac {d{\boldsymbol {v}}(t)}{dt}}\end{aligned}}} という公式が成り立つので、 W ( t 1 , t 2 ) = m ∫ t 1 t 2 d 2 d t 2 r ( t ) ⋅ d r ( t ) d t d t = m ∫ t 1 t 2 { d d t ( d r ( t ) d t ⋅ d r ( t ) d t ) − d r ( t ) d t ⋅ d 2 r ( t ) d t 2 } d t = 1 2 m ∫ t 1 t 2 d d t ( d r ( t ) d t ⋅ d r ( t ) d t ) d t = 1 2 m | d r ( t 2 ) d t | 2 − 1 2 m | d r ( t 1 ) d t | 2 {\displaystyle {\begin{aligned}W(t_{1},t_{2})&=m\int _{t_{1}}^{t_{2}}{\frac {d^{2}}{dt^{2}}}{\boldsymbol {r}}(t)\cdot {\frac {d{\boldsymbol {r}}(t)}{dt}}dt\\&=m\int _{t_{1}}^{t_{2}}\left\{{\frac {d}{dt}}\left({\frac {d{\boldsymbol {r}}(t)}{dt}}\cdot {\frac {d{\boldsymbol {r}}(t)}{dt}}\right)-{\frac {d{\boldsymbol {r}}(t)}{dt}}\cdot {\frac {d^{2}{\boldsymbol {r}}(t)}{dt^{2}}}\right\}dt\\&={\frac {1}{2}}m\int _{t_{1}}^{t_{2}}{\frac {d}{dt}}\left({\frac {d{\boldsymbol {r}}(t)}{dt}}\cdot {\frac {d{\boldsymbol {r}}(t)}{dt}}\right)dt\\&={\frac {1}{2}}m\left|{\frac {d{\boldsymbol {r}}(t_{2})}{dt}}\right|^{2}-{\frac {1}{2}}m\left|{\frac {d{\boldsymbol {r}}(t_{1})}{dt}}\right|^{2}\end{aligned}}} が得られる。ここで得られた関数 1 2 m | d r ( t ) d t | 2 {\displaystyle {\frac {1}{2}}m\left|{\frac {d{\boldsymbol {r}}(t)}{dt}}\right|^{2}} は粒子の運動エネルギーと呼ばれ、この差分は粒子になされた仕事を表す。 ポテンシャルと仕事、運動エネルギーと仕事の関係をそれぞれ見比べると、 − { V ( r ( t 2 ) ) − V ( r ( t 1 ) ) } + ∫ t 1 t 2 f ( t ) ⋅ d r ( t ) = 1 2 m | d r ( t 2 ) d t | 2 − 1 2 m | d r ( t 1 ) d t | 2 ∫ t 1 t 2 f ( t ) ⋅ d r ( t ) = { 1 2 m | d r ( t 2 ) d t | 2 + V ( r ( t 2 ) ) } − { 1 2 m | d r ( t 1 ) d t | 2 + V ( r ( t 1 ) ) } {\displaystyle {\begin{aligned}-\left\{V({\boldsymbol {r}}(t_{2}))-V({\boldsymbol {r}}(t_{1}))\right\}+\int _{t_{1}}^{t_{2}}{\boldsymbol {f}}(t)\cdot d{\boldsymbol {r}}(t)&={\frac {1}{2}}m\left|{\frac {d{\boldsymbol {r}}(t_{2})}{dt}}\right|^{2}-{\frac {1}{2}}m\left|{\frac {d{\boldsymbol {r}}(t_{1})}{dt}}\right|^{2}\\\int _{t_{1}}^{t_{2}}{\boldsymbol {f}}(t)\cdot d{\boldsymbol {r}}(t)&=\left\{{\frac {1}{2}}m\left|{\frac {d{\boldsymbol {r}}(t_{2})}{dt}}\right|^{2}+V({\boldsymbol {r}}(t_{2}))\right\}-\left\{{\frac {1}{2}}m\left|{\frac {d{\boldsymbol {r}}(t_{1})}{dt}}\right|^{2}+V({\boldsymbol {r}}(t_{1}))\right\}\end{aligned}}} という等式が得られる。 f ( t ) {\displaystyle {\boldsymbol {f}}(t)} を粒子に対する力学的な操作によって生じる力だとすれば、それがなす仕事は操作の前後での粒子の力学的エネルギー、すなわち粒子の位置エネルギーと運動エネルギーの和、の差に等しい。 特に、外部から力学的操作を行わない場合には、粒子にはポテンシャルによる力しか働かないので、系の力学的エネルギーは保存されることになる。 また、操作の前後で粒子の速度を変えないようにすれば、操作の前後では粒子の運動エネルギーが変化しないので、外部から与えられた仕事は粒子のポテンシャルの差に等しくなる。 こうして得られた等式が成り立つことを、力学的エネルギー保存の法則と呼ぶ。保存則が成り立っているかどうかは系の設定により、外界の力学的エネルギーを考慮しない場合には、保存則は成り立たないが、外界の力学的エネルギーを考慮するのであれば、外界への仕事を付け加える形で、保存則が成立する。 外界に及ぼされる力は − f ( t ) {\displaystyle -{\boldsymbol {f}}(t)} で表され、摩擦などによる抗力を考える場合には、粒子の速度の関数になる。 以上のことは多粒子系の場合にも成り立つ。一粒子系の場合との変更点は、各粒子に対して力と運動方程式が与えられることと、ポテンシャルがすべての粒子の位置の関数になることである。以下にN 個の粒子がある場合について示す。 力： F i ( { r ( t ) } , t ) = − ∇ i V ( { r i ( t ) } ) + f i ( t ) ( i = 1 , … , N ) . {\displaystyle {\boldsymbol {F}}_{i}(\{{\boldsymbol {r}}(t)\},t)=-\nabla _{i}V(\{{\boldsymbol {r}}_{i}(t)\})+{\boldsymbol {f}}_{i}(t)\quad (i=1,\dots ,N).} 運動方程式： m i d 2 r i d t 2 ( t ) = F i ( { r ( t ) } , t ) ( i = 1 , … , N ) . {\displaystyle m_{i}{\frac {d^{2}{\boldsymbol {r}}_{i}}{dt^{2}}}(t)={\boldsymbol {F}}_{i}(\{{\boldsymbol {r}}(t)\},t)\quad (i=1,\dots ,N).} ナブラ ∇ i {\displaystyle \nabla _{i}} は、粒子 i {\displaystyle i} の位置に対する偏微分を表し、ポテンシャルの勾配は次のように変更される。 ∇ i V ( { r i ( t ) } ) = ( ∂ ∂ x i V ( { r i ( t ) } ) , ∂ ∂ y i V ( { r i ( t ) } ) , ∂ ∂ z i V ( { r i ( t ) } ) ) T ( i = 1 , … , N ) . {\displaystyle \nabla _{i}V(\{{\boldsymbol {r}}_{i}(t)\})=\left({\frac {\partial }{\partial x_{i}}}V(\{{\boldsymbol {r}}_{i}(t)\}),{\frac {\partial }{\partial y_{i}}}V(\{{\boldsymbol {r}}_{i}(t)\}),{\frac {\partial }{\partial z_{i}}}V(\{{\boldsymbol {r}}_{i}(t)\})\right)^{T}\quad (i=1,\dots ,N).} また、ポテンシャルの時間微分は、それぞれの粒子の速度と粒子が感じるポテンシャルの勾配の内積をすべて足しあわせたものになる。 d V d t ( { r i ( t ) } ) = ∑ i = 1 N d r i ( t ) d t ⋅ ∇ i V ( { r i ( t ) } ) . {\displaystyle {\frac {dV}{dt}}(\{{\boldsymbol {r}}_{i}(t)\})=\sum _{i=1}^{N}{\frac {d{\boldsymbol {r}}_{i}(t)}{dt}}\cdot \nabla _{i}V(\{{\boldsymbol {r}}_{i}(t)\}).} 系になされる仕事は、各粒子に対する仕事の和になる。 W ( t 1 , t 2 ) = ∑ i = 1 N ∫ t 1 t 2 F i ( { r i ( t ) } , t ) ⋅ d r i ( t ) . {\displaystyle W(t_{1},t_{2})=\sum _{i=1}^{N}\int _{t_{1}}^{t_{2}}{\boldsymbol {F}}_{i}(\{{\boldsymbol {r}}_{i}(t)\},t)\cdot d{\boldsymbol {r}}_{i}(t).} 以上のことから、力学的エネルギー保存の法則は次のように表される。 ∑ i = 1 N ∫ t 1 t 2 f i ( t ) ⋅ d r i ( t ) = { V ( { r i ( t 2 ) } ) + ∑ i = 1 N ( 1 2 m i | d r i ( t 2 ) d t | 2 ) } − { V ( { r i ( t 1 ) } ) + ∑ i = 1 N ( 1 2 m i | d r i ( t 1 ) d t | 2 ) } . {\displaystyle \sum _{i=1}^{N}\int _{t_{1}}^{t_{2}}{\boldsymbol {f}}_{i}(t)\cdot d{\boldsymbol {r}}_{i}(t)=\left\{V(\{{\boldsymbol {r}}_{i}(t_{2})\})+\sum _{i=1}^{N}\left({\frac {1}{2}}m_{i}\left|{\frac {d{\boldsymbol {r}}_{i}(t_{2})}{dt}}\right|^{2}\right)\right\}-\left\{V(\{{\boldsymbol {r}}_{i}(t_{1})\})+\sum _{i=1}^{N}\left({\frac {1}{2}}m_{i}\left|{\frac {d{\boldsymbol {r}}_{i}(t_{1})}{dt}}\right|^{2}\right)\right\}.} 一粒子の場合と異なり、各粒子の運動エネルギーの総和と系のポテンシャルの和が系の力学的エネルギーの役割を果たしている。 量子力学においてもエネルギー保存の法則は厳密に成立する。量子力学において、あらゆる物理量はそれに対応する自己共役作用素として定義される。閉じた系のエネルギーを与える作用素は、古典力学のハミルトニアンに対応する作用素 H