化学ポテンシャル (かがくポテンシャル、英語 : chemical potential )は、熱力学 で用いられる示強性 状態量 の一つで、浸透圧 や相平衡 、化学反応 のようなマクロな物質量の移動が伴う現象で重要となる物理量 である。
推奨される量記号は、μ (ミュー)である。
化学ポテンシャルの概念および記号と用語は、ウィラード・ギブズ の1876年の論文『不均一な物質系の平衡に就いて 』で導入された。
化学ポテンシャルは、物質の多寡により系が潜在的に持つエネルギー の大きさの尺度となる量である。
例えば、半透膜 で隔てられた二つの系の間に濃度 差が有った場合、浸透圧 が生じ仕事 を為す事が出来る。
また、物質が増減する化学反応 では熱の出入り(発熱反応 、吸熱反応 )を伴う。
このように、物質が存在することにより系は潜在的にエネルギーを持つ。
その系 に含まれるある成分の単位物質量 あたりのギブスエネルギー がその成分の化学ポテンシャルに相当する。
定義
熱力学的な系 の内部エネルギー U の微小変化dU を考えると、熱力学第一法則 より、
d
U
=
d
′
Q
−
d
′
W
{\displaystyle dU=d'Q-d'W}
が成り立つ。ここで、d'Q は外部から系に流れる熱量、d'W は系が外部にする仕事である。系のエントロピー をS 、熱力学的温度 をT とすると、d'Q=TdS であるから、上式は、
d
U
=
T
d
S
−
d
′
W
{\displaystyle dU=TdS-d'W}
となる。系が外部と粒子の出入りがない場合では、圧力 をP 、体積をV とすると、d'W=PdV であるから、
d
U
=
T
d
S
−
P
d
V
{\displaystyle dU=TdS-PdV}
となるが[ 2] 、粒子の出入りがある場合では、追加の項が必要であり、dU は
d
U
=
T
d
S
−
P
d
V
+
∑
i
μ
i
d
N
i
{\displaystyle dU=TdS-PdV+\sum _{i}\mu _{i}dN_{i}}
(A)
となる[ 3] 。(すなわち、
d
′
W
=
P
d
V
−
∑
i
μ
i
d
N
i
{\displaystyle d'W=PdV-\sum _{i}\mu _{i}dN_{i}}
である)
ここで、
d
N
i
{\displaystyle dN_{i}}
は成分i の物質量
N
i
{\displaystyle N_{i}}
の微小変化を表しており、
μ
i
{\displaystyle \mu _{i}}
は
μ
i
(
S
,
V
,
N
)
=
(
∂
U
∂
N
i
)
S
,
V
,
N
j
≠
i
{\displaystyle \mu _{i}(S,V,{\boldsymbol {N}})=\left({\frac {\partial U}{\partial N_{i}}}\right)_{S,V,N_{j\neq i}}}
で定義される、化学ポテンシャルと呼ばれる量[ 4] である。(全微分 も参照。)括弧に付く添え字はその変数を一定にして偏微分することを意味する。また、
N
j
≠
i
{\displaystyle N_{j\neq i}}
は
i
{\displaystyle i}
以外の全ての成分の物質量を表す。
その他の表現
化学ポテンシャルは様々な変数の組の関数として、また様々な熱力学ポテンシャル の偏微分の関数として表現される。例えば、F を系のヘルムホルツエネルギー とすると、成分
i
{\displaystyle i}
の化学ポテンシャルは
μ
i
(
T
,
V
,
N
)
=
(
∂
F
∂
N
i
)
T
,
V
,
N
j
≠
i
{\displaystyle \mu _{i}(T,V,{\boldsymbol {N}})=\left({\frac {\partial F}{\partial N_{i}}}\right)_{T,V,N_{j\neq i}}}
と表される[ 5] 。これは次のようにして示される。まず、
F
=
U
−
T
S
{\displaystyle F=U-TS}
なので、その微小変化は、
d
F
=
d
U
−
S
d
T
−
T
d
S
{\displaystyle dF=dU-SdT-TdS}
である。ここで(A)のdU を代入すると
d
F
(
T
,
V
,
N
)
=
−
S
d
T
−
P
d
V
+
∑
i
μ
i
d
N
i
{\displaystyle dF\left(T,V,{\boldsymbol {N}}\right)=-SdT-PdV+\sum _{i}\mu _{i}dN_{i}}
であるから、上式が成り立つことが言える。ここで、U からF への変換はルジャンドル変換 となっている。同様に、系のギブズエネルギー G=U-TS+PV とエンタルピー H=U+PV に対して、
μ
i
(
T
,
P
,
N
)
=
(
∂
G
∂
N
i
)
T
,
P
,
N
j
≠
i
{\displaystyle \mu _{i}(T,P,{\boldsymbol {N}})=\left({\frac {\partial G}{\partial N_{i}}}\right)_{T,P,N_{j\neq i}}}
μ
i
(
S
,
P
,
N
)
=
(
∂
H
∂
N
i
)
S
,
P
,
N
j
≠
i
{\displaystyle \mu _{i}(S,P,{\boldsymbol {N}})=\left({\frac {\partial H}{\partial N_{i}}}\right)_{S,P,N_{j\neq i}}}
も示される。また、(A)より、エントロピー S の微小変化は
d
S
=
1
T
d
U
+
P
T
d
V
−
1
T
∑
i
μ
i
d
N
i
{\displaystyle dS={\frac {1}{T}}dU+{\frac {P}{T}}dV-{\frac {1}{T}}\sum _{i}\mu _{i}dN_{i}}
なので、
μ
i
(
U
,
V
,
N
)
=
−
T
(
∂
S
∂
N
i
)
U
,
V
,
N
j
≠
i
{\displaystyle \mu _{i}(U,V,{\boldsymbol {N}})=-T\left({\frac {\partial S}{\partial N_{i}}}\right)_{U,V,N_{j\neq i}}}
も成り立つ。
物理的な意味
化学ポテンシャルの物理的な意味は、同じ示強性である圧力との対応を考えるとわかりやすい。 [独自研究? ] たとえば圧力 (示強性)は、熱力学的な系の体積 (示量性)を少し変えたときに外界が感じる『手ごたえ』 [要追加記述 ] である。この関係性を化学ポテンシャルに当てはめてみると、化学ポテンシャル(示強性)とは、熱力学的な系の物質量(示量性)を少し変えたときの『手ごたえ』と考えることができる[ 6] 。よって [要追加記述 ] 平衡状態に向かうときは、化学ポテンシャルが等しくなるように物質量は移動する[ 7] 。
また電磁気学 において電荷 q とその移動を司る静電ポテンシャル φ との積がポテンシャルエネルギー qφ である。この関係性を化学ポテンシャルに当てはめてみると[ 8] 、マクロな物質量N の移動を司るポテンシャルが化学ポテンシャルμ であり[ 9] 、それらの積であるギブズエネルギーNμ はポテンシャルエネルギーのような量だと考えることもできる。 [誰? ] ただし実際にはミクロな粒子間にある複雑な相互作用などの結果としてマクロな化学ポテンシャルは決まると考えられ、力学におけるポテンシャル [要追加記述 ] と熱力学における化学ポテンシャルはかなり異なり同一視することはできない[ 9] 。
ギブズエネルギーとの関係
部分モルギブズエネルギー
温度 T と圧力 p 、および物質量 の組 N = (N 1 , N 2 ,...) により平衡状態が指定される場合での化学ポテンシャルは
μ
i
(
N
;
T
,
p
)
=
(
∂
G
∂
N
i
)
T
,
p
,
N
j
≠
i
{\displaystyle \mu _{i}({\boldsymbol {N}};T,p)=\left({\frac {\partial G}{\partial N_{i}}}\right)_{T,p,N_{j\neq i}}}
で与えられる[ 10] 。
このように温度と圧力と、1成分を除いた物質量を固定した示量変数の偏微分は部分モル量(partial molal quantity )と呼ばれ[ 12] [ 13] 、この意味で化学ポテンシャルは部分モルギブズエネルギーに等しい。
オイラーの関係式
系のスケール変換を考えれば、ギブスエネルギー と物質量の示量性 、及び温度と圧力の示強性 から、スケール・パラメータ λ に対して
λ
G
(
N
;
T
,
p
)
=
G
(
λ
N
;
T
,
p
)
{\displaystyle \lambda G({\boldsymbol {N}};T,p)=G(\lambda {\boldsymbol {N}};T,p)}
が成り立つ。これを λ について微分すれば
G
(
N
;
T
,
p
)
=
∑
i
N
i
∂
G
∂
N
i
|
N
=
λ
N
=
∑
i
N
i
μ
i
(
λ
N
;
T
,
p
)
{\displaystyle G({\boldsymbol {N}};T,p)=\sum _{i}N_{i}\left.{\frac {\partial G}{\partial N_{i}}}\right|_{N=\lambda N}=\sum _{i}N_{i}\,\mu _{i}(\lambda {\boldsymbol {N}};T,p)}
であり、λ = 1 と置けば
G
(
N
;
T
,
p
)
=
∑
i
N
i
μ
i
(
N
;
T
,
p
)
{\displaystyle G({\boldsymbol {N}};T,p)=\sum _{i}N_{i}\,\mu _{i}({\boldsymbol {N}};T,p)}
の関係が得られる。各成分の化学ポテンシャルとその成分の物質量の積の総和がギブズエネルギーとなる。
特に単一成分系では
G
(
N
;
T
,
p
)
=
N
μ
(
T
,
p
)
{\displaystyle G(N;T,p)=N\mu (T,p)}
であり、ギブズエネルギーは物質量に比例 し、化学ポテンシャルは物質量に依らない。
つまり1成分系では温度と圧力が等しければ化学ポテンシャルは等しい。これは自由に熱を通し自由に動くことができる壁に穴を開けても、平衡状態は変化しない(壁の両側でマクロな物質量は変化しない)ことを意味する[ 7] 。
化学ポテンシャルの偏微分
温度T 、圧力p を変数としたときの化学ポテンシャルの偏微分は
(
∂
μ
∂
T
)
p
=
1
N
(
∂
G
∂
T
)
p
,
N
=
−
S
N
{\displaystyle \left({\frac {\partial \mu }{\partial T}}\right)_{p}={\frac {1}{N}}\left({\frac {\partial G}{\partial T}}\right)_{p,N}=-{\frac {S}{N}}}
(
∂
μ
∂
p
)
T
=
1
N
(
∂
G
∂
p
)
T
,
N
=
V
N
{\displaystyle \left({\frac {\partial \mu }{\partial p}}\right)_{T}={\frac {1}{N}}\left({\frac {\partial G}{\partial p}}\right)_{T,N}={\frac {V}{N}}}
となる。
「ギブズ・デュエムの式 」も参照
具体的な表示
気体の化学ポテンシャル
理想気体 のモル体積は V m = RT /p であり、これを積分すると
μ
(
T
,
p
)
=
μ
∘
(
T
)
+
R
T
ln
p
p
∘
{\displaystyle \mu (T,p)=\mu ^{\circ }(T)+RT\ln {\frac {p}{p^{\circ }}}}
となる。ここで p ° は標準状態圧力 、μ ° は標準化学ポテンシャルである。
実在気体の場合はモル体積をビリアル展開 で
V
m
(
T
,
p
)
=
R
T
p
+
B
(
T
)
+
C
(
T
)
p
+
O
(
p
2
)
{\displaystyle V_{\text{m}}(T,p)={\frac {RT}{p}}+B(T)+C(T)\,p+O(p^{2})}
と表わしたものを積分すれば
μ
(
T
,
p
)
=
μ
∘
(
T
)
+
R
T
ln
p
p
∘
+
B
(
T
)
p
+
C
(
T
)
p
2
2
+
O
(
p
3
)
{\displaystyle \mu (T,p)=\mu ^{\circ }(T)+RT\ln {\frac {p}{p^{\circ }}}+B(T)\,p+C(T)\,{\frac {p^{2}}{2}}+O(p^{3})}
となる。標準化学ポテンシャルは
μ
∘
(
T
)
=
lim
p
→
0
{
μ
(
T
,
p
)
−
R
T
ln
p
p
∘
}
{\displaystyle \mu ^{\circ }(T)=\lim _{p\to 0}\left\{\mu (T,p)-RT\ln {\frac {p}{p^{\circ }}}\right\}}
で定義される[ 14] 。
また、フガシティー f を用いることで実在気体の化学ポテンシャルを
μ
(
T
,
p
)
=
μ
∘
(
T
)
+
R
T
ln
f
(
T
,
p
)
p
∘
{\displaystyle \mu (T,p)=\mu ^{\circ }(T)+RT\ln {\frac {f(T,p)}{p^{\circ }}}}
と表わすこともできる。ビリアル展開の形と比較すれば、フガシティーは
f
(
T
,
p
)
=
p
e
B
p
/
R
T
+
O
(
p
2
)
{\displaystyle f(T,p)=p\,e^{Bp/RT+O(p^{2})}}
である[ 14] 。
混合のポテンシャル
混合物の組成をモル分率 xi の組で表したとき、理想混合系の化学ポテンシャルは
μ
i
id
(
x
;
T
,
p
)
=
μ
i
∗
(
T
,
p
)
+
R
T
ln
x
i
{\displaystyle \mu _{i}^{\text{id}}({\boldsymbol {x}};T,p)=\mu _{i}^{*}(T,p)+RT\ln x_{i}}
で表される[ 15] 。ここで * は純物質 における量を表している。
特に理想混合気体では
μ
i
id
(
x
;
T
,
p
)
=
μ
i
∘
(
T
)
+
R
T
ln
x
i
p
p
∘
{\displaystyle \mu _{i}^{\text{id}}({\boldsymbol {x}};T,p)=\mu _{i}^{\circ }(T)+RT\ln {\frac {x_{i}\,p}{p^{\circ }}}}
となり、純粋系での圧力をその組成の分圧 xi p で置き換えた形となる。
実在気体の混合系では、分圧をフガシティーへ置き換えて表される。
実在の固溶体 や実在溶液 では、モル分率を活量 へ置き換えて
μ
i
(
x
;
T
,
p
)
=
μ
i
∗
(
T
,
p
)
+
R
T
ln
a
i
(
x
;
T
,
p
)
{\displaystyle \mu _{i}({\boldsymbol {x}};T,p)=\mu _{i}^{*}(T,p)+RT\ln a_{i}({\boldsymbol {x}};T,p)}
で表わされる[ 16] 。
溶液の化学ポテンシャル
理想溶液 において、溶質 i の濃度が質量モル濃度 bi で表されるときの化学ポテンシャルは
μ
i
(
b
;
T
,
p
)
=
μ
b
,
i
∘
(
T
,
p
)
+
R
T
ln
b
i
b
∘
{\displaystyle \mu _{i}({\boldsymbol {b}};T,p)=\mu _{b,i}^{\circ }(T,p)+RT\ln {\frac {b_{i}}{b^{\circ }}}}
で表される。ここで b ° は標準質量モル濃度であり、通常 b ° = 1 mol /kg に選ばれる。
溶質の濃度がモル濃度 ci で表されるときの化学ポテンシャルは
μ
i
(
c
;
T
,
p
)
=
μ
c
,
i
∘
(
T
,
p
)
+
R
T
ln
c
i
c
∘
{\displaystyle \mu _{i}({\boldsymbol {c}};T,p)=\mu _{c,i}^{\circ }(T,p)+RT\ln {\frac {c_{i}}{c^{\circ }}}}
で表される[ 17] 。ここで c ° は標準モル濃度であり、通常 c ° = 1 mol/L に選ばれる。
実在溶液の場合は活量 を用いることで、それぞれに
μ
i
(
b
;
T
,
p
)
=
μ
b
,
i
∘
(
T
,
p
)
+
R
T
ln
a
b
,
i
(
b
;
T
,
p
)
{\displaystyle \mu _{i}({\boldsymbol {b}};T,p)=\mu _{b,i}^{\circ }(T,p)+RT\ln a_{b,i}({\boldsymbol {b}};T,p)}
μ
i
(
c
;
T
,
p
)
=
μ
c
,
i
∘
(
T
,
p
)
+
R
T
ln
a
c
,
i
(
c
;
T
,
p
)
{\displaystyle \mu _{i}({\boldsymbol {c}};T,p)=\mu _{c,i}^{\circ }(T,p)+RT\ln a_{c,i}({\boldsymbol {c}};T,p)}
と表わすことができる[ 18] [ 17] 。
無限希釈の極限 b → 0 あるいは c → 0 で理想溶液に漸近するので、標準化学ポテンシャルは
μ
b
,
i
∘
(
T
,
p
)
=
lim
b
→
0
{
μ
i
(
b
;
T
,
p
)
−
R
T
ln
b
i
b
∘
}
{\displaystyle \mu _{b,i}^{\circ }(T,p)=\lim _{b\to 0}\left\{\mu _{i}({\boldsymbol {b}};T,p)-RT\ln {\frac {b_{i}}{b^{\circ }}}\right\}}
μ
c
,
i
∘
(
T
,
p
)
=
lim
c
→
0
{
μ
i
(
c
;
T
,
p
)
−
R
T
ln
c
i
c
∘
}
{\displaystyle \mu _{c,i}^{\circ }(T,p)=\lim _{c\to 0}\left\{\mu _{i}({\boldsymbol {c}};T,p)-RT\ln {\frac {c_{i}}{c^{\circ }}}\right\}}
で定義される[ 18] 。
見かけの化学ポテンシャル
溶質の濃度が質量モル濃度 b で表されるときの溶媒のモル分率は
x
slv
=
1
1
+
M
slv
∑
slt
b
{\displaystyle x_{\text{slv}}={\frac {1}{1+M_{\text{slv}}\sum _{\text{slt}}b}}}
なので、理想溶液における溶媒の化学ポテンシャルは
μ
slv
id
(
b
;
T
,
p
)
=
μ
slv
∗
(
T
,
p
)
−
R
T
ln
(
1
+
M
slv
∑
slt
b
i
)
{\displaystyle \mu _{\text{slv}}^{\text{id}}({\boldsymbol {b}};T,p)=\mu _{\text{slv}}^{*}(T,p)-RT\ln \left(1+M_{\text{slv}}\sum _{\text{slt}}b_{i}\right)}
である。実在溶液においては活量で置き換えて
μ
slv
(
b
;
T
,
p
)
=
μ
slv
∗
(
T
,
p
)
+
R
T
ln
a
slv
(
b
;
T
,
p
)
=
μ
slv
∗
(
T
,
p
)
−
ϕ
R
T
M
slv
∑
slt
b
i
{\displaystyle {\begin{aligned}\mu _{\text{slv}}({\boldsymbol {b}};T,p)&=\mu _{\text{slv}}^{*}(T,p)+RT\ln a_{\text{slv}}({\boldsymbol {b}};T,p)\\&=\mu _{\text{slv}}^{*}(T,p)-\phi RTM_{\text{slv}}\sum _{\text{slt}}b_{i}\end{aligned}}}
となる。
ここで
ϕ
(
b
;
T
,
p
)
=
−
ln
a
slv
(
b
;
T
,
p
)
M
slv
∑
slt
b
{\displaystyle \phi ({\boldsymbol {b}};T,p)=-{\frac {\ln a_{\text{slv}}({\boldsymbol {b}};T,p)}{M_{\text{slv}}\sum _{\text{slt}}b}}}
は浸透係数である。
このときギブズエネルギーは
G
=
N
slt
μ
slv
+
∑
slt
N
i
μ
i
=
N
slt
μ
slv
∗
+
∑
slt
N
i
{
μ
i
−
ϕ
R
T
}
=
N
slt
μ
slv
∗
+
∑
slt
N
i
μ
i
app
{\displaystyle {\begin{aligned}G&=N_{\text{slt}}\,\mu _{\text{slv}}+\sum _{\text{slt}}N_{i}\,\mu _{i}\\&=N_{\text{slt}}\,\mu _{\text{slv}}^{*}+\sum _{\text{slt}}N_{i}\left\{\mu _{i}-\phi RT\right\}\\&=N_{\text{slt}}\,\mu _{\text{slv}}^{*}+\sum _{\text{slt}}N_{i}\,\mu _{i}^{\text{app}}\end{aligned}}}
となる。ここで app は見かけの量 (apparent molar property ) を表している。
化学平衡
化学量論数 νi で表される化学反応 において、反応進行度 を ξ とすれば、物質量は
N
i
(
ξ
)
=
N
0
,
i
+
ν
i
ξ
{\displaystyle N_{i}(\xi )=N_{0,i}+\nu _{i}\,\xi }
と表わされる。等温等圧条件下ではギブズエネルギー が減少する方向に変化が進行し、平衡状態においてギブズエネルギーが極小となる。従って
Δ
r
G
≡
d
d
ξ
G
(
N
(
ξ
)
;
T
,
p
)
=
0
{\displaystyle \Delta _{\text{r}}G\equiv {\frac {d}{d\xi }}G({\boldsymbol {N}}(\xi );T,p)=0}
を満たす ξ において化学平衡となる。
反応のギブズエネルギーは化学ポテンシャルを用いて
Δ
r
G
=
∑
i
d
N
i
d
ξ
(
∂
G
∂
N
i
)
T
,
p
,
N
j
≠
i
=
∑
i
ν
i
μ
i
(
N
(
ξ
)
;
T
,
p
)
{\displaystyle \Delta _{\text{r}}G=\sum _{i}{\frac {dN_{i}}{d\xi }}\,\left({\frac {\partial G}{\partial N_{i}}}\right)_{T,p,N_{j\neq i}}=\sum _{i}\nu _{i}\,\mu _{i}({\boldsymbol {N}}(\xi );T,p)}
と書くことができて、理想混合気体においては
Δ
r
G
=
∑
i
ν
i
μ
i
∘
(
T
)
+
R
T
∑
i
ln
(
p
x
i
(
ξ
)
p
∘
)
ν
i
{\displaystyle \Delta _{\text{r}}G=\sum _{i}\nu _{i}\,\mu _{i}^{\circ }(T)+RT\sum _{i}\ln \left({\frac {px_{i}(\xi )}{p^{\circ }}}\right)^{\nu _{i}}}
となる。
標準平衡定数 を
R
T
ln
K
∘
(
T
;
ν
)
=
−
∑
i
ν
i
μ
i
∘
(
T
)
{\displaystyle RT\ln K^{\circ }(T;{\boldsymbol {\nu }})=-\sum _{i}\nu _{i}\,\mu _{i}^{\circ }(T)}
で定義すれば、平衡の条件は
K
∘
(
T
;
ν
)
=
∏
i
(
p
x
i
(
ξ
)
p
∘
)
ν
i
{\displaystyle K^{\circ }(T;{\boldsymbol {\nu }})=\prod _{i}\left({\frac {px_{i}(\xi )}{p^{\circ }}}\right)^{\nu _{i}}}
となる。
物性物理学への応用
モル数でなく、粒子 数としての化学ポテンシャルμ も考えることができる。固体電子論 における電子系(例:電子ガス )でも化学ポテンシャルを定義することができ、特に温度T = 0 Kにおける化学ポテンシャル μ のことを、フェルミエネルギー ε F と呼ぶ場合がある。
μ
T
=
0
=
ϵ
F
{\displaystyle \mu _{T=0}=\epsilon _{F}}
脚注
^ 北孝文 p.42
^ 北孝文 p.61
^ 北孝文 p.66
^ 田崎 p.174
^ 田崎『熱力学』 p.120
^ a b 佐々真一 『熱力学入門』共立出版 、2000年。ISBN 978-4320033474 。
^ 田崎『熱力学』
^ a b 清水『熱力学の基礎』
^ バーロー『物理化学(上)』 pp.233-235, §8.10
^ バーロー『物理化学(上)』 p.163
^ Kirkwood & Oppenheim p.9
^ a b Kirkwood & Oppenheim p.89-92, §7-6.
^ バーロー『物理化学(下)』 p.612
^ バーロー『物理化学(下)』 pp.619-621, §19.4
^ a b バーロー『物理化学(下)』 pp.621-624, §19.5
^ a b Kirkwood & Oppenheim p.160-164, §11-1.
参考文献
Gibbs, J. Willard (1876-05-01). “On the Equilibrium of Heterogeneous Substances. First Part.”. the Transactions of the Connecticut Academy of Arts and Sciences (Connecticut Academy of Arts and Sciences) 3 : 108–248.
J.G. Kirkwood, I. Oppenheim (1961). Chemical Thermodynamics . McGraw-Hill
G.M.バーロー『物理化学』 上巻、藤代亮一 訳(第1版)、東京化学同人、1968年。
G.M.バーロー『物理化学』 下巻、藤代亮一 訳(第1版)、東京化学同人、1968年。
田崎晴明『熱力学 - 現代的な視点から』培風館〈新物理学シリーズ〉、2000年。ISBN 4-563-02432-5 。
清水明『熱力学の基礎』東京大学出版会、2007年。ISBN 978-4-13-062609-5 。
北孝文、鈴木久男(監修)『演習しよう熱・統計力学』数理工学社、2018年。ISBN 978-4-86481-053-1 。
関連項目
外部リンク