グリーン–久保公式
グリーン–久保公式(Green–Kubo relations)あるいは中野–久保公式とは、線形応答理論における、輸送係数をカレントの時間相関で表す関係式を一般化して定式化されたものである。メルヴィル・S・グリーン、中野藤生、久保亮五らの名前を冠して名付けられている。
外場H'(t) = - A F(t)による物理量Bの線形応答\operatorname{Tr}(\rho'(t)B) は、応答関数を \phi_{BA} (t)とすると
\operatorname{Tr} (\rho'(t) B) = \int_{- \infty}^t \phi_{BA} (t - t') F(t') dt'
\phi_{BA} (t)  \equiv - { 1 \over { i \hbar} } \operatorname{Tr} ( [A, \rho_0] \bar{B} (t) ) = {1 \over { i \hbar} } \operatorname{Tr} ( [\rho_0, A] \bar{B} (t) ) = {1 \over { i \hbar} } \operatorname{Tr} ( [A, \bar{B} (t) ] \rho_0 )
となる。
外場F_\mathrm{ex}が存在するとき、熱伝導率や粘性率などの輸送係数をL(0)、カレントを\boldsymbol{J}とすると、輸送係数は以下のように時間相関関数で表せる。
L(0)
= \lim_{F_\mathrm{ex} \to 0} L(F_\mathrm{ex})
= \beta V \int_{0}^{\infty} \lim_{F_\mathrm{ex} \to 0} \langle J(0) J(t) \rangle_{F_\mathrm{ex} dt
ここで\langle \quad \rangle_{F_\mathrm{ex}は外場F_\mathrm{ex}があるときのアンサンブル平均である。
ある系が時間 t について t = −∞ で熱平衡状態であるとする。この時点で系に外場は印加されていない。時間 t に依存する外場 H′(t)を考え、これが最初の時点(t = −∞)から十分時間が経った段階で系に働くとして、その時の該当する系全体のハミルトニアンは次のように表される。
H_\text{total} = H_0 + H'(t)
ここでH0 は外場のない時の系のハミルトニアンで、これは時間に依存しないとする。
外場H′(t)が次のように表現できるとする。
H'(t) = - A F(t)
ここでA は時間を含まない演算子で、系におけるある物理量を表す。F(t) はこの演算子を通じて系に作用する外場(の大きさ)であり、これは演算子ではないとする。F(t) は次を満たさなければならない。
\lim_{t \to - \infty} F(t) = 0
とする。ここで、系全体を記述する密度行列(統計演算子)を導入し、これを ρtotal  とすると Htotal の式に対応して系全体の密度行列は、
\rho_\text{total} = \rho_0 + \rho'(t)
と表される。系全体の時間発展は次のフォン・ノイマンの式で表される。
i \hbar {\partial \rho_\text{total} \over {\partial t} }  = [H,\rho_\text{total}]
ここで \hbar = h / 2 \pi 、h はプランク定数、上式右辺の括弧は交換関係を表している。 ρ′(t) は外場に対応する密度行列であり、これは次を満たさなければならない。
\lim_{t \to - \infty}\rho'(t) = 0
外場に関係する H′(t), ρ′(t) は、それぞれH0, ρ0 に対し十分に小さいものと考え、2次の項を無視すると以下が得られる。
i \hbar {\partial \rho' \over {\partial t} } = [H_0, \rho'] + [H', \rho_0]
次にρ′(t) を次のように表現し直す。
\bar{\rho}' \equiv \exp \left( { {i H_0 t \over \hbar} } \right) \rho' (t) \exp \left( { - { i H _0 t \over \hbar} } \right)
この時間発展は次のようになる。
i \hbar { \partial \bar{\rho}'(t) \over {\partial t} } = [ \bar{H}' (t) , \rho_0]
\bar{H}'(t) = \exp \left( { {i H_0 t \over \hbar} } \right) H'(t) \exp \left( { - { i H _0 t \over \hbar} } \right)
よって\bar{\rho}' の時間発展は
\bar{\rho}'(t) = {1 \over {i \hbar} } \int_{- \infty}^t [ \bar{H}'(t'), \rho_0] dt'
となり、外場 H′(t) の1次まで考えると、系全体の密度演算子は次のようになる。
\rho_\mathrm{total} (t) = \rho_0 + {1 \over {i \hbar} \int_{- \infty}^t  \exp \left( { - {i H_0 (t - t') \over \hbar} } \right) [ \bar{H}'(t), \rho_0]  \exp \left( { { i H _0 (t-t') \over \hbar} } \right) dt'
ここで、 状態ρtotal(t) について時間に依らないある物理量 B の統計的期待値を取ると
\langle B \rangle = \operatorname{Tr} (\rho (t) B) = \operatorname{Tr} (\rho_0 B) + \operatorname{Tr} (\rho'(t) B)
となる(Tr はトレース:対角和をとることを意味する)。ここで、上式最右辺の第一項の Tr 内は、時間に依存しないのでその期待値をゼロとみなす。従って問題となるのは第二項の部分で、H′(t) = −AF(t) 及び、Tr 内の演算は交換可能で、かつ演算子も循環的に演算順序を変えることができることから、Tr 及び B を積分内に移動し、B、 \exp(-i H_0 (t-t')/ \hbar) の順にこれらを先頭に移動すると、
\begin{align} \operatorname{Tr} (\rho'(t) B)
&= {i \over {\hbar} } \operatorname{Tr} \left( \int_{- \infty}^t  \exp \left( { - {i H_0 (t - t') \over \hbar} } \right) [ A, \rho_0]  \exp \left( { { i H_0 (t-t') \over \hbar} } \right) F(t') dt' B \right)  \\
&= {i \over {\hbar} } \int_{- \infty}^t \operatorname{Tr} \left( B \exp \left( { - {i H_0 (t - t') \over \hbar} } \right) [ A, \rho_0]  \exp \left( { { i H_0 (t-t') \over \hbar} } \right) F(t') \right) dt' \\
&= {i \over {\hbar} } \int_{- \infty}^t \operatorname{Tr} \left( \exp \left( { { i H_0 (t-t') \over \hbar} } \right) B \exp \left( { - {i H_0 (t - t') \over \hbar} } \right) [ A, \rho_0]   F(t') \right) dt'
\end{align}
となる。そして B を
\bar{B} (t-t') = \exp \left( { {i H_0 (t - t') \over \hbar} } \right) B \exp \left( { -{ i H_0 (t-t') \over \hbar} } \right)
と表現し直し、 \, [A, \rho_0] を前に置くと、
\operatorname{Tr} (\rho'(t) B) = - {1 \over { i \hbar} } \int_{- \infty}^t \operatorname{Tr} \left( [A, \rho_0] \bar{B} (t-t') \right) F(t') dt'
と変形できる。以上で、F(t′) は先の定義により単なる大きさを表す量なのでどこにでも置くことができる。
次に応答関数なるものを、
\phi_{BA} (t)  \equiv - { 1 \over { i \hbar} } \operatorname{Tr} ( [A, \rho_0] \bar{B} (t) ) = {1 \over { i \hbar} } \operatorname{Tr} ( [\rho_0, A] \bar{B} (t) ) = {1 \over { i \hbar} } \operatorname{Tr} ( [A, \bar{B} (t) ] \rho_0 )
と定義すると、Tr(ρ′(t)B) は、
\operatorname{Tr} (\rho'(t) B) = \int_{- \infty}^t \phi_{BA} (t - t') F(t') dt'
と表せる。