統計物理学 において、ブロッホ=ドミニシスの定理 (英 : Bloch-De Dominicis theorem )とは、量子多体系 における熱平均で定義された多点相関関数 を、2点相関関数の組み合わせ和に分解する定理。場の量子論 の真空期待値 に関するウィックの定理 に対し、有限温度の系での類似版に相当しており、ウィックの定理 とも呼ばれる。物理学者 松原武生 によって、温度グリーン関数 の理論展開ともに導入された[ 1] 。定理の名は、最初に完全な証明を与えた物理学者C. ブロッホとC. T. ドミニシスに因む[ 2] 。
A を生成演算子 aα † 、または消滅演算子 aα 、もしくはそれらを虚時間 で相互作用表示 したものとする。ここで、相互作用表示において、非摂動系の自由ハミルトニアン は
H
0
=
∑
α
ε
α
a
α
†
a
α
{\displaystyle H_{0}=\sum _{\alpha }\varepsilon _{\alpha }a_{\alpha }^{\,\dagger }a_{\alpha }}
のように2次形式で表されているとする。また平均値
⟨
⋅
⟩
0
{\displaystyle \langle \cdot \rangle _{0}}
は
e
−
β
(
H
0
−
μ
N
)
{\displaystyle e^{-\beta (H_{0}-\mu N)}}
による非摂動系でのグランドカノニカル分布 の熱平均を表すものとする。
このとき、この熱平均で定義されるn 点相関関数は、n が偶数である場合のみゼロにならず、
⟨
A
1
A
2
⋯
A
n
⟩
0
=
{
0
(
n
:
odd
)
∑
m
=
2
n
(
±
1
)
m
⟨
A
1
A
m
⟩
0
⟨
A
2
⋯
A
m
+
1
A
m
−
1
⋯
A
n
⟩
0
(
n
:
even
)
{\displaystyle \langle A_{1}A_{2}\cdots A_{n}\rangle _{0}=\left\{{\begin{matrix}0&\quad (n:{\mbox{odd}})\\\sum _{m=2}^{n}(\pm 1)^{m}\langle A_{1}A_{m}\rangle _{0}\langle A_{2}\cdots A_{m+1}A_{m-1}\cdots A_{n}\rangle _{0}&\quad (n:{\mbox{even}})\end{matrix}}\right.}
が成り立つ。ここで現れる2点相関関数は縮約 (contraction)と呼ばれる。また、(±1)m の項はフェルミ粒子 での演算子の順番の並べ替えにおいて、隣合う演算子同士を置き換える際に生じる符号の反転を表しており、符号は正がボーズ粒子 、負がフェルミ粒子 に対応するものとする。以降、本項に現れる複合の符号は全て、上がボーズ粒子、下がフェルミ粒子に対応するものとする。
さらに、n が偶数であるときは、この結果を繰り返し適用することで、
⟨
A
1
A
2
⋯
A
n
⟩
0
=
∑
P
′
(
±
1
)
P
⟨
A
i
1
A
i
2
⟩
0
⟨
A
i
3
A
i
4
⟩
0
⋯
⟨
A
i
n
−
1
A
i
n
⟩
0
{\displaystyle \langle A_{1}A_{2}\cdots A_{n}\rangle _{0}=\sum _{P}'(\pm 1)^{P}\langle A_{i_{1}}A_{i_{2}}\rangle _{0}\langle A_{i_{3}}A_{i_{4}}\rangle _{0}\cdots \langle A_{i_{n-1}}A_{i_{n}}\rangle _{0}}
と全ての演算子の縮約の組み合わせ和に分解できる。ここで、P は(1, 2, …, n )→(i 1 , i 2 , …, in )なる置換 を表し、和Σ'において、ここの縮約で対となる演算子は、ik -1 < ik を満たし、全体としてはi 1 < i 3 < …< in の順序が満たされる項について和をとるものする。また、(±1)P は、フェルミ粒子の隣合う演算子同士の並べ替えの際に付与する正負の符号の変化を表し、ボーズ粒子については+1、フェルミ粒子については演算子の並べ替えの回数に応じた符号(置換の符号)を与えるものとする。これをブロッホ=ドミニシスの定理 またはウィックの定理 と呼ぶ[ 3] [ 4] 。
縮約について、記法
A
∙
B
∙
:=
⟨
A
B
⟩
0
{\displaystyle A^{\bullet }B^{\bullet }:=\langle AB\rangle _{0}}
を導入すれば、全ての可能な縮約の組み合わせをとるというブロッホ=ドミニシスの定理は
⟨
A
1
A
2
⋯
A
n
⟩
0
=
A
1
∙
A
2
∙
A
3
∙
∙
⋯
A
n
∙
∙
∙
+
A
1
∙
A
2
∙
∙
A
3
∙
⋯
A
n
∙
∙
∙
+
A
1
∙
A
2
∙
∙
A
3
∙
∙
⋯
A
n
∙
∙
∙
+
⋯
{\displaystyle \langle A_{1}A_{2}\cdots A_{n}\rangle _{0}=A_{1}^{\bullet }A_{2}^{\bullet }A_{3}^{\bullet \bullet }\cdots A_{n}^{\bullet \bullet \bullet }+A_{1}^{\bullet }A_{2}^{\bullet \bullet }A_{3}^{\bullet }\cdots A_{n}^{\bullet \bullet \bullet }+A_{1}^{\bullet }A_{2}^{\bullet \bullet }A_{3}^{\bullet \bullet }\cdots A_{n}^{\bullet \bullet \bullet }+\cdots }
とも表すことができる。但し、複数個の演算子の縮約については、
A
∙
B
∙
∙
C
∙
D
∙
∙
=
±
(
A
∙
C
∙
)
(
B
∙
∙
D
∙
∙
)
=
±
⟨
A
C
⟩
0
⟨
B
D
⟩
0
A
∙
B
∙
∙
C
∙
∙
D
∙
=
(
A
∙
D
∙
)
(
B
∙
∙
C
∙
∙
)
=
⟨
A
D
⟩
0
⟨
B
C
⟩
0
{\displaystyle {\begin{aligned}A^{\bullet }B^{\bullet \bullet }C^{\bullet }D^{\bullet \bullet }&=\pm (A^{\bullet }C^{\bullet })(B^{\bullet \bullet }D^{\bullet \bullet })=\pm \langle AC\rangle _{0}\langle BD\rangle _{0}\\A^{\bullet }B^{\bullet \bullet }C^{\bullet \bullet }D^{\bullet }&=(A^{\bullet }D^{\bullet })(B^{\bullet \bullet }C^{\bullet \bullet })=\langle AD\rangle _{0}\langle BC\rangle _{0}\end{aligned}}}
のように、同じ右付き添え字の演算子同士の縮約を行い、さらにフェルミ粒子の場合には演算子の入れ替えの回数に応じた符号を与えるものとする。
A =A (τ)が虚時間での相互作用表示の演算子であるとし、虚時間に対する時間順序積 をとる場合にも、同様にブロッホ=ドミニシスの定理
⟨
T
τ
[
A
1
(
τ
1
)
A
2
(
τ
2
)
⋯
A
n
(
τ
n
)
]
⟩
0
=
∑
P
′
(
±
1
)
P
⟨
T
τ
[
A
i
1
(
τ
i
1
)
A
i
2
(
τ
i
2
)
]
⟩
0
⟨
T
τ
[
A
i
3
(
τ
i
3
)
A
i
4
(
τ
i
4
)
]
⟩
0
⋯
⟨
T
τ
[
A
i
n
−
1
(
τ
i
n
−
1
)
A
i
n
(
τ
i
n
)
]
⟩
0
{\displaystyle \langle T_{\tau }[A_{1}(\tau _{1})A_{2}(\tau _{2})\cdots A_{n}(\tau _{n})]\rangle _{0}=\sum _{P}'(\pm 1)^{P}\langle T_{\tau }[A_{i_{1}}(\tau _{i_{1}})A_{i_{2}}(\tau _{i_{2}})]\rangle _{0}\langle T_{\tau }[A_{i_{3}}(\tau _{i_{3}})A_{i_{4}}(\tau _{i_{4}})]\rangle _{0}\cdots \langle T_{\tau }[A_{i_{n-1}}(\tau _{i_{n-1}})A_{i_{n}}(\tau _{i_{n}})]\rangle _{0}}
が成り立つ。
この場合も縮約について、
A
(
τ
a
)
∘
B
(
τ
b
)
∘
:=
⟨
T
τ
[
A
(
τ
a
)
B
(
τ
b
)
]
⟩
0
{\displaystyle A(\tau _{a})^{\circ }B(\tau _{b})^{\circ }:=\langle T_{\tau }[A(\tau _{a})B(\tau _{b})]\rangle _{0}}
となる記法を導入すれば、
⟨
T
τ
[
A
1
(
τ
1
)
A
2
(
τ
2
)
⋯
A
n
(
τ
n
)
]
⟩
0
=
A
1
(
τ
1
)
∘
A
2
(
τ
2
)
∘
A
3
(
τ
3
)
∘
∘
⋯
A
n
(
τ
n
)
∘
∘
∘
+
A
1
(
τ
1
)
∘
A
2
(
τ
2
)
∘
∘
A
3
(
τ
3
)
∘
⋯
A
n
(
τ
n
)
∘
∘
∘
+
A
1
(
τ
1
)
∘
A
2
(
τ
2
)
∘
∘
A
3
(
τ
3
)
∘
∘
⋯
A
n
(
τ
n
)
∘
∘
∘
+
⋯
{\displaystyle {\begin{aligned}\langle T_{\tau }[A_{1}(\tau _{1})A_{2}(\tau _{2})\cdots A_{n}(\tau _{n})]\rangle _{0}=&A_{1}(\tau _{1})^{\circ }A_{2}(\tau _{2})^{\circ }A_{3}(\tau _{3})^{\circ \circ }\cdots A_{n}(\tau _{n})^{\circ \circ \circ }\\&+A_{1}(\tau _{1})^{\circ }A_{2}(\tau _{2})^{\circ \circ }A_{3}(\tau _{3})^{\circ }\cdots A_{n}(\tau _{n})^{\circ \circ \circ }+A_{1}(\tau _{1})^{\circ }A_{2}(\tau _{2})^{\circ \circ }A_{3}(\tau _{3})^{\circ \circ }\cdots A_{n}(\tau _{n})^{\circ \circ \circ }+\cdots \end{aligned}}}
と表すことができる。
ブロッホ=ドミニシスの定理により、3点相関関数、4点相関関数について
⟨
A
B
C
⟩
0
=
0
{\displaystyle \langle ABC\rangle _{0}=0}
⟨
A
B
C
D
⟩
0
=
A
∙
B
∙
C
∙
∙
D
∙
∙
+
A
∙
B
∙
∙
C
∙
D
∙
∙
+
A
∙
B
∙
∙
C
∙
∙
D
∙
=
⟨
A
B
⟩
0
⟨
C
D
⟩
0
±
⟨
A
C
⟩
0
⟨
B
D
⟩
0
+
⟨
A
D
⟩
0
⟨
B
C
⟩
0
{\displaystyle {\begin{aligned}\langle ABCD\rangle _{0}&=A^{\bullet }B^{\bullet }C^{\bullet \bullet }D^{\bullet \bullet }+A^{\bullet }B^{\bullet \bullet }C^{\bullet }D^{\bullet \bullet }+A^{\bullet }B^{\bullet \bullet }C^{\bullet \bullet }D^{\bullet }\\&=\langle AB\rangle _{0}\langle CD\rangle _{0}\pm \langle AC\rangle _{0}\langle BD\rangle _{0}+\langle AD\rangle _{0}\langle BC\rangle _{0}\end{aligned}}}
が成り立つ。
時間順序積をとる場合にも同様に
⟨
T
τ
[
A
(
τ
a
)
B
(
τ
b
)
C
(
τ
c
)
]
⟩
0
=
0
{\displaystyle \langle T_{\tau }[A(\tau _{a})B(\tau _{b})C(\tau _{c})]\rangle _{0}=0}
⟨
T
τ
[
A
(
τ
a
)
B
(
τ
b
)
C
(
τ
c
)
D
(
τ
d
)
]
⟩
0
=
A
(
τ
a
)
∘
B
(
τ
b
)
∘
C
(
τ
c
)
∘
∘
D
(
τ
d
)
∘
∘
+
A
(
τ
a
)
∘
B
(
τ
b
)
∘
∘
C
(
τ
c
)
∘
D
(
τ
d
)
∘
∘
+
A
(
τ
a
)
∘
B
(
τ
b
)
∘
∘
C
(
τ
c
)
∘
∘
D
(
τ
d
)
∘
=
⟨
T
τ
[
A
(
τ
a
)
B
(
τ
b
)
]
⟩
0
⟨
T
τ
[
C
(
τ
c
)
D
(
τ
d
)
]
⟩
0
±
⟨
T
τ
[
A
(
τ
a
)
C
(
τ
c
)
]
⟩
0
⟨
T
τ
[
B
(
τ
b
)
D
(
τ
d
)
]
⟩
0
+
⟨
T
τ
[
A
(
τ
a
)
D
(
τ
d
)
]
⟩
0
⟨
T
τ
[
B
(
τ
b
)
C
(
τ
c
)
]
⟩
0
{\displaystyle {\begin{aligned}\langle T_{\tau }[A(\tau _{a})B(\tau _{b})C(\tau _{c})D(\tau _{d})]\rangle _{0}=&A(\tau _{a})^{\circ }B(\tau _{b})^{\circ }C(\tau _{c})^{\circ \circ }D(\tau _{d})^{\circ \circ }\\&+A(\tau _{a})^{\circ }B(\tau _{b})^{\circ \circ }C(\tau _{c})^{\circ }D(\tau _{d})^{\circ \circ }+A(\tau _{a})^{\circ }B(\tau _{b})^{\circ \circ }C(\tau _{c})^{\circ \circ }D(\tau _{d})^{\circ }\\=&\langle T_{\tau }[A(\tau _{a})B(\tau _{b})]\rangle _{0}\langle T_{\tau }[C(\tau _{c})D(\tau _{d})]\rangle _{0}\\&\pm \langle T_{\tau }[A(\tau _{a})C(\tau _{c})]\rangle _{0}\langle T_{\tau }[B(\tau _{b})D(\tau _{d})]\rangle _{0}+\langle T_{\tau }[A(\tau _{a})D(\tau _{d})]\rangle _{0}\langle T_{\tau }[B(\tau _{b})C(\tau _{c})]\rangle _{0}\end{aligned}}}
が成り立つ。
量子多体系 における温度グリーン関数 の理論では、温度グリーン関数によって系の様々な物理量を求めることができるともに、摂動 計算を系統的に行うことができる。ここで演算子A , B の温度グリーン関数は
G
A
B
(
τ
,
τ
′
)
:=
−
⟨
T
τ
[
A
(
τ
)
B
(
τ
′
)
]
⟩
{\displaystyle G_{AB}(\tau ,\tau '):=-\langle T_{\tau }[A(\tau )B(\tau ')]\rangle }
で定義される2点相関関数である。但し、記号
⟨
⋯
⟩
{\displaystyle \langle \cdots \rangle }
は
⟨
⋯
⟩
=
T
r
{
e
−
β
(
H
−
μ
N
)
⋯
}
T
r
{
e
−
β
(
H
−
μ
N
)
}
=
T
r
{
e
−
β
K
⋯
}
T
r
{
e
−
β
K
}
(
K
:=
H
−
μ
N
)
{\displaystyle \langle \cdots \rangle ={\frac {\mathrm {Tr} \{e^{-\beta (H-\mu N)}\cdots \}}{\mathrm {Tr} \{e^{-\beta (H-\mu N)}\}}}={\frac {\mathrm {Tr} \{e^{-\beta K}\cdots \}}{\mathrm {Tr} \{e^{-\beta K}\}}}\quad (K:=H-\mu N)}
で定義されるグランドカノニカル分布での熱平均であり、H はハミルトニアン 、N は数演算子 、β は逆温度 、μ は化学ポテンシャル を表す。またA (τ )は
A
(
τ
)
=
e
+
K
τ
A
e
−
K
τ
A
(
τ
)
†
=
e
+
K
τ
A
†
e
−
K
τ
{\displaystyle {\begin{aligned}A(\tau )&=e^{+K\tau }Ae^{-K\tau }\\A(\tau )^{\dagger }&=e^{+K\tau }A^{\dagger }e^{-K\tau }\end{aligned}}}
で定義される虚時間 τ =it についてのハイゼンベルク表示 の演算子である。Tτ は虚時間についての時間順序積であり、
T
τ
{
A
1
(
τ
1
)
A
2
(
τ
2
)
}
=
{
A
1
(
τ
1
)
A
2
(
τ
2
)
(
τ
1
>
τ
2
)
±
A
2
(
τ
2
)
A
1
(
τ
1
)
(
τ
2
>
τ
1
)
{\displaystyle T_{\tau }\{A_{1}(\tau _{1})A_{2}(\tau _{2})\}=\left\{{\begin{matrix}A_{1}(\tau _{1})A_{2}(\tau _{2})&\quad (\tau _{1}>\tau _{2})\\\pm A_{2}(\tau _{2})A_{1}(\tau _{1})&\quad (\tau _{2}>\tau _{1})\end{matrix}}\right.}
を意味する。
一般に温度グリーン関数の計算において、相互作用にある系では、
H
=
H
0
+
V
{\displaystyle H=H_{0}+V\,}
とハミルトニアンを可解な非摂動項と相互作用を含む摂動項に分け、相互作用表示 の演算子
A
I
(
τ
)
=
e
+
K
0
τ
A
e
−
K
0
τ
A
I
(
τ
)
†
=
e
+
K
0
τ
A
†
e
−
K
0
τ
(
K
0
:=
H
0
−
μ
N
)
{\displaystyle {\begin{aligned}A_{I}(\tau )&=e^{+K_{0}\tau }Ae^{-K_{0}\tau }\\A_{I}(\tau )^{\dagger }&=e^{+K_{0}\tau }A^{\dagger }e^{-K_{0}\tau }\quad (K_{0}:=H_{0}-\mu N)\end{aligned}}}
に対して、摂動計算を行うことが必要となる。このとき、摂動計算において、
⟨
A
1
I
(
τ
1
)
⋯
A
n
I
(
τ
n
)
⟩
0
=
T
r
{
e
−
β
K
0
A
1
I
(
τ
1
)
⋯
A
n
I
(
τ
n
)
}
T
r
{
e
−
β
K
0
}
{\displaystyle \langle A_{1I}(\tau _{1})\cdots A_{nI}(\tau _{n})\rangle _{0}={\frac {\mathrm {Tr} \{e^{-\beta K_{0}}A_{1I}(\tau _{1})\cdots A_{nI}(\tau _{n})\}}{\mathrm {Tr} \{e^{-\beta K_{0}}\}}}}
という高次の相関関数が現れる。ブロッホ=ドミニシスの定理は、こうした多点相関関数を縮約
⟨
a
α
†
a
α
′
⟩
0
=
δ
α
α
′
1
∓
e
β
ε
α
⟨
a
α
a
α
′
†
⟩
0
=
δ
α
α
′
1
∓
e
−
β
ε
α
⟨
a
α
†
a
α
′
†
⟩
0
=
0
⟨
a
α
a
α
′
⟩
0
=
0
:
{\displaystyle {\begin{aligned}\langle a_{\alpha }^{\,\dagger }a_{\alpha '}\rangle _{0}&={\frac {\delta _{\alpha \alpha '}}{1\mp e^{\beta \varepsilon _{\alpha }}}}\\\langle a_{\alpha }a_{\alpha '}^{\,\dagger }\rangle _{0}&={\frac {\delta _{\alpha \alpha '}}{1\mp e^{-\beta \varepsilon _{\alpha }}}}\\\langle a_{\alpha }^{\,\dagger }a_{\alpha '}^{\,\dagger }\rangle _{0}&=0\\\langle a_{\alpha }a_{\alpha '}\rangle _{0}&=0\end{aligned}}:}
に分解し、実際の計算を可能にする。
ブロッホ=ドミニシスの定理は、古典系におけるガウス過程 の持つ性質を量子系に拡張したものに相当する
[ 5] 。実際、分布が
P
n
(
x
1
,
⋯
,
x
n
)
=
1
det
A
exp
(
−
∑
i
,
j
A
i
j
(
x
i
−
⟨
X
i
⟩
)
(
x
j
−
⟨
X
j
⟩
)
)
{\displaystyle P_{n}(x_{1},\cdots ,x_{n})={\frac {1}{\sqrt {\det {A}}}}\exp {{\bigl (}-\sum _{i,j}A_{ij}(x_{i}-\langle X_{i}\rangle )(x_{j}-\langle X_{j}\rangle ){\bigr )}}}
で与えられるガウス過程[ 6] を考えると
⟨
X
1
X
2
⋯
X
n
⟩
=
∑
a
l
l
p
a
r
t
i
t
i
o
n
∏
⟨
X
i
1
⋯
X
i
m
⟩
c
{\displaystyle \langle X_{1}X_{2}\cdots X_{n}\rangle =\sum _{\operatorname {all\,\,partition} }\prod \langle X_{i_{1}}\cdots X_{i_{m}}\rangle _{c}}
が成り立つ。ここで<…>はこの分布に対する期待値 、<…>c はキュムラント を表すものとする。
また、右辺の和は、X 1 ,…,X n をいくつかの集まりに分割する全ての組み合わせにわたってとるものである。例えば、3点相関関数、4点相関関数については、
⟨
X
1
X
2
⟩
=
⟨
X
1
X
2
⟩
c
+
⟨
X
1
⟩
c
⟨
X
2
⟩
c
⟨
X
1
X
2
X
3
⟩
=
⟨
X
1
X
2
X
3
⟩
c
+
⟨
X
1
⟩
c
⟨
X
2
X
3
⟩
c
+
⟨
X
2
⟩
c
⟨
X
1
X
3
⟩
c
+
⟨
X
3
⟩
c
⟨
X
1
X
2
⟩
c
+
⟨
X
1
⟩
c
⟨
X
2
⟩
c
⟨
X
3
⟩
c
{\displaystyle {\begin{aligned}\langle X_{1}X_{2}\rangle =&\langle X_{1}X_{2}\rangle _{c}+\langle X_{1}\rangle _{c}\langle X_{2}\rangle _{c}\\\langle X_{1}X_{2}X_{3}\rangle =&\langle X_{1}X_{2}X_{3}\rangle _{c}+\langle X_{1}\rangle _{c}\langle X_{2}X_{3}\rangle _{c}+\langle X_{2}\rangle _{c}\langle X_{1}X_{3}\rangle _{c}+\langle X_{3}\rangle _{c}\langle X_{1}X_{2}\rangle _{c}+\langle X_{1}\rangle _{c}\langle X_{2}\rangle _{c}\langle X_{3}\rangle _{c}\end{aligned}}}
である。
さらに全てのXi について、⟨Xi ⟩=0であるとするならば、ブロッホ=ドミニシスの定理と同様にn 点相関関数は、n が偶数である場合のみゼロにならず、
⟨
X
1
X
2
⋯
X
n
⟩
=
∑
a
l
l
p
a
i
r
∏
⟨
X
i
X
j
⟩
=
∑
P
′
⟨
X
i
1
X
i
2
⟩
⟨
X
i
3
X
i
4
⟩
⋯
⟨
X
i
n
−
1
X
i
n
⟩
{\displaystyle {\begin{aligned}\langle X_{1}X_{2}\cdots X_{n}\rangle &=\sum _{\operatorname {all\,\,pair} }\prod \langle X_{i}X_{j}\rangle \\&=\sum _{P}'\langle X_{i_{1}}X_{i_{2}}\rangle \langle X_{i_{3}}X_{i_{4}}\rangle \cdots \langle X_{i_{n-1}}X_{i_{n}}\rangle \end{aligned}}}
と2点相関関数の組み合わせ和に分解される。
^ T. Matsubara, Prog. Theor. Phys. , 14 , p.351 (1955)
^ C. Bloch and C. T. de Dominicis, Nucl. Phys. , 7 , p.459 (1958)
^ A. L. Fetter and J. D. Walecka (2003)
^ 阿部龍蔵 (1992)
^ 今田正俊 (2004)
^ A =(A ij )は正定値行列
論文
G. C. Wick (1950). “The Evaluation of the Collision Matrix”. Phys. Rev. 80 : 268. doi :10.1103/PhysRev.80.268 .
T. Matsubara (1955). “A New Approach to Quantum-Statistical Mechanics”. Prog. Theor. Phys. 14 : 351. doi :10.1143/PTP.14.351 .
C. Bloch; C. T. de Dominicis (1958). “Un développement du potentiel de gibbs d'un système quantique composé d'un grand nombre de particules”. Nucl. Phys. 7 : 459. doi :10.1016/0029-5582(58)90285-2 .
M. Gaudin (1960). “Une démonstration simplifiée du théorème de wick en mécanique statistique”. Nucl. Phys. 15 : 89. doi :10.1016/0029-5582(60)90285-6 .
書籍