グリーン–久保公式

グリーン–久保公式()あるいは中野–久保公式とは、線形応答理論における、輸送係数をカレントの時間相関で表す関係式を一般化して定式化されたものである。メルヴィル・S・グリーン、中野藤生、久保亮五らの名前を冠して名付けられている。
外場formula_1による物理量formula_2の線形応答formula_3 は、応答関数をformula_4とすると
となる。
輸送係数におけるグリーン–久保公式.
外場formula_7が存在するとき、熱伝導率や粘性率などの輸送係数をformula_8、カレントをformula_9とすると、輸送係数は以下のように時間相関関数で表せる。
ここでformula_11は外場formula_7があるときのアンサンブル平均である。
導出.
系のハミルトニアン.
ある系が時間 "t" について "t" = −∞ で熱平衡状態であるとする。この時点で系に外場は印加されていない。時間 "t" に依存する外場 "H"′("t")を考え、これが最初の時点("t" = −∞)から十分時間が経った段階で系に働くとして、その時の該当する系全体のハミルトニアンは次のように表される。
ここで"H"0 は外場のない時の系のハミルトニアンで、これは時間に依存しないとする。
外場H′("t")が次のように表現できるとする。
ここで"A" は時間を含まない演算子で、系におけるある物理量を表す。"F"("t") はこの演算子を通じて系に作用する外場（の大きさ）であり、これは演算子ではないとする。"F"("t") は次を満たさなければならない。
密度行列.
とする。ここで、系全体を記述する密度行列（統計演算子）を導入し、これを "ρ"total とすると "H"total の式に対応して系全体の密度行列は、
と表される。系全体の時間発展は次のフォン・ノイマンの式で表される。
ここでformula_18、"h" はプランク定数、上式右辺の括弧は交換関係を表している。 "ρ"′("t") は外場に対応する密度行列であり、これは次を満たさなければならない。
外場に関係する "H"′("t"), "ρ"′("t") は、それぞれ"H"0, "ρ"0 に対し十分に小さいものと考え、2次の項を無視すると以下が得られる。
次にρ′("t") を次のように表現し直す。
この時間発展は次のようになる。
よってformula_24の時間発展は
となり、外場 "H"′("t") の1次まで考えると、系全体の密度演算子は次のようになる。
物理量の期待値.
ここで、 状態"ρ"total("t") について時間に依らないある物理量 "B" の統計的期待値を取ると
となる（Tr はトレース：対角和をとることを意味する）。ここで、上式最右辺の第一項の Tr 内は、時間に依存しないのでその期待値をゼロとみなす。従って問題となるのは第二項の部分で、"H"′("t") = −"AF"("t") 及び、Tr 内の演算は交換可能で、かつ演算子も循環的に演算順序を変えることができることから、Tr 及び "B" を積分内に移動し、"B"、formula_28の順にこれらを先頭に移動すると、
となる。そして "B" を
と表現し直し、formula_31を前に置くと、
と変形できる。以上で、"F"("t"′) は先の定義により単なる大きさを表す量なのでどこにでも置くことができる。
応答関数・Green–Kubo公式.
次に応答関数なるものを、
と定義すると、Tr("ρ"′("t")"B") は、
と表せる。