クリスマスなので温度グリーン関数の話をする
なかしゅです。
はじめに
みなさん、クリスマスですね。クリスマスを象徴する色といえば、やはりクリスマスツリーに代表される緑でしょう。華やかな装飾に彩られた緑の木々に囲まれながら大切な人[1]とひとときを過ごす、そんな特別な日がクリスマスというものではないでしょうか。ところで、冬は寒いものです。人間という生き物は寒さに決して強くはないので、寒い時には心も体も温まるような何かが欲しくなる訳です。
という訳で今回は、温度グリーン関数と呼ばれるものについて解説していきたいと思います。特に物性物理の文脈でそれがどのようなモチベーションで導入されてどのように理論の中で用いられるのかについて、摂動論との関係を中心に説明したいと思います。
なお、記事の長さの関係上ある程度の前提知識を仮定せざるを得なくなりました(それでも長くはなりましたが)。明らかな説明不足と思われる箇所も多々あります(し、その自覚はあります)がご容赦ください。解析力学と場の量子論と統計力学についてのほんの少しの知識があると読みやすいかと思われます[2]。
グリーン関数って何?
グリーン関数とは、一言で言えば「非斉次な微分方程式の一般解を求める上で有用な、微分方程式の特性解」のことです。わかりやすい例を一つ挙げると、電磁気学において、電荷密度の空間分布
この方程式(これはPoisson方程式と呼ばれる形の式です)の一般解は、
を満たすような関数
と書くことができます。実際、上式の両辺に
グリーン関数には「時空上の異なる点どうしの間の相関を担うもの」という役割があります。すなわち、グリーン関数は、場の中のある点における事象の影響がどのように(どれだけの重みで)異なる点に伝播するのかを教えてくれます。
この事情を上の例で確かめましょう。上の解を見てみると、空間上のある点
まとめると、グリーン関数は時空上のある点と別の点とを結びつける(=2点の間に相関を持たせる)役割を果たしていて、すなわちグリーン関数は場の異なる2点どうしがどのようなに相関を持っているのかという情報を持っています。
さらに、ここでは詳細は省きますが、グリーン関数は重要な観測可能量と直接結びついています。例えば、絶対零度の相対論的な場の量子論では、「
グリーン関数がわかると、場を記述する方程式の一般解がわかり、場の相関の様子がわかり、さらに観測可能量の計算もできる。グリーン関数には場に対するこれだけの豊富な情報が含まれているわけです。したがって、ものすごく雑な言い方をしてしまえば、場の理論における主要な目標は与えられた場に対するグリーン関数を求めることだと言えるでしょう。
(絶対零度の)場の量子論におけるグリーン関数
さて、しばらくは「絶対零度の相対論的な場の量子論[3]」におけるグリーン関数の話をします。もっとも簡単な例として、相互作用のない自由なスカラー粒子の運動を考えます。ここでは4次元Minkowski時空の計量を
のように与えられ、これに最小作用の原理を適用する(作用を
が得られます(これはKlein-Gordon方程式と呼ばれる式です)。引数は4次元Minkowski時空における座標です。フーリエ変換を利用して4元運動量の言葉で書き換えてやると、上を満たす
のように書かれます。ここで、波数に依存する展開係数
さて、作用を出発点とするこの場の量子論においては、場の演算子
に対して、次の同時刻交換関係
を課すことによって場の量子化がなされます。場の演算子を平面波展開して生成・消滅演算子の言葉に書き換えると、上の同時刻交換関係は
と、波数に関する同時刻交換関係に言い換えられます。
さて、上で与えた自由場のオイラー-ラグランジュ(E-L)方程式に対するグリーン関数というものを求めてみましょう。グリーン関数は、E-L方程式に登場する微分作用素に対してディラックのデルタ関数を返すもの
として与えられます(
となります(ただし、積分の特異点を「とりあえず」解消するために、虚軸方向に特異点を微小にずらしています)。ゆえ、自由スカラー場のグリーン関数は
と求められます。計算の途中過程で発散する量が出てくると厄介なので、とりあえず暫くの間は特異点をずらして実軸上での発散を形式的に取り除いておいて、計算の最後に
作用(というよりはラグランジアン密度)に相互作用に起因する項がある場合、最小作用の原理から導かれるE-L方程式は非斉次な微分方程式になります。その時の解
ところで、この自由場スカラー場のグリーン関数は別の方法によって与えることもできます。すなわち、「場の演算子の積(時間順序積)の真空期待値=2点相関関数」という量をこの自由スカラー場に対して計算すると、それが微分方程式の特性解という意味でE-L方程式に対するグリーン関数となっていることを確かめることができます。
ここで、「時間順序積」および「真空期待値」というものについて簡単に説明します。まず、時間順序積とは、引数の時間成分=第0成分の大きい演算子から順に左から並べるという順序で演算子の積を定義するということです(場の量子論では演算子同士が非可換なので、何らかの規則によって積の順序を定めてあげる必要があります)。例えば、2つの演算子
と書けます。「階段関数」
を用いれば、時間順序積をまとめて
と書くことができます。
また、任意の演算子
ただし、この
物理量(観測可能な量)と直接結びつくのは演算子そのものではなくあくまでもその期待値であるので、我々が本当に欲しいのは期待値を取ったものになります。
これより、場の演算子の「2点相関関数」を
と与えます。場の演算子をフーリエ変換で波数表示して、上の2点相関関数を
と書けば、これが先に求めた自由スカラー場のグリーン関数と同じ表示を持つことがわかります。階段関数の複素積分表示
を利用すると、
となることが確かめられます。ただし、途中計算では生成/消滅演算子の交換関係から
となることを用いています(2つ目の等号は、真空に消滅演算子を作用させると消える(
ところで、この節では「自由スカラー場のグリーン関数」という言い方をしました。一方、一番初めに述べた静電ポテンシャルを決定するPoisson方程式の文脈では「Poisson方程式のグリーン関数」という言い方をしています。微分方程式の文脈では、グリーン関数は微分演算子の逆演算子として導入され、この時「〇〇方程式のグリーン関数」と呼ばれます。このグリーン関数が座標空間上の異なる2点の間の相関を担う役割を果たしていたことから、これを拡張して場の量子論ないし量子多体系の文脈ではグリーン関数を場の演算子の2点相関関数(時間順序積の期待値)として与えている、という構図になっています。
一般に相互作用のある場合は、場の演算子に相互作用の影響が現れるため(相互作用のある場の演算子は、E-L方程式の非斉次項と自由場のグリーン関数の畳み込みで与えられます)、2点相関関数は微分演算子の逆演算子にはなりません。すなわち、場の量子論におけるグリーン関数が微分方程式の意味でのグリーン関数と一致するのは自由場のケースのみです。以下、本稿で「グリーン関数」といった場合は場の演算子の2点相関関数の意味で用いることとします。
なお、場の演算子の
素粒子実験等で
相互作用のある場と摂動論
今までの計算結果は全て相互作用のない自由場に対するものでした。自由場の理論ではグリーン関数を厳密に計算することができましたが、一般の場合は自由場ほど綺麗に定式化されていません。以下ではラグランジアンに相互作用項があるような場合(
相互作用のある系の場の理論は言ってしまえば不完全です[5]。例えば、自由場の理論のようにフーリエ変換(平面波展開)によって直ちに粒子描像に移行することはできません(相互作用のある場においては、平面波解が固有関数の直交系をなさないため)。より、一般の相互作用のある場をそのまま取り扱うことは困難どころか現在の場の量子論の枠組みでは不可能です。そこで、相互作用の効果を小さいとして摂動論的に取り入れたいという気持ちが湧いてきます。以下では摂動論による相互作用場の取り扱い方をごくごく簡単に紹介します。
相互作用場の2点相関関数はE-L方程式の微分演算子の逆演算子にはなっていないので、グリーン関数を求めるには2点相関関数を具体的に計算しなければなりません。しかし、上述の事情により、一般に相互作用のある場においては状態空間をうまく構成することができません(特に最低エネルギー状態=基底状態を具体的に求めることができません)。そこで、一般の相互作用のある場における真空期待値を、経路積分量子化による重み付けのもとで与えることとします:
これは、真空の具体的な表式によらずに求められます。
経路積分量子化では、相空間内の取りうる可能なあらゆる状態に対して作用の指数関数で重みをつけた上で、全ての状態の寄与を足し合わせて最後に規格化する(=真空の寄与で割る)ことによって、場の量の真空期待値を計算することができます。やっていることとしては、統計力学においてカノニカル分布のもとで熱平均を計算することと近いです(実現可能なあらゆる状態にエネルギーの指数関数で重みをつけて、全ての状態の寄与を足し合わせて規格化する(=分配関数で割る)ことで熱平均を得る)。ただし、上の表式ではすでに
詳細は全て省きますが、経路積分の積分測度の定め方[6]から、経路積分量子化によって場の演算子の積の期待値を計算するとそれは必然的に時間順序積となります。また、真空の寄与を除することによって計算結果から非本質な発散を取り除いています。なお、上の表式は「無限の過去-未来における場の演算子は、漸近的に自由場の演算子と一致する」という、相互作用場の漸近的完全性(という、場の量子論に対する要請)を用いています。
さて、改めてこの経路積分表示による真空期待値の表式を見てみると、これは指数関数の中に作用積分
すなわち、相互作用のある場に対する摂動展開とは、指数関数の肩の作用積分の相互作用項を小さいとして自由場の項の周りで展開することだとわかります。これは、相互作用場のもとでの真空期待値を自由場の真空期待値まわりで展開している、ということになります。
この摂動展開の各項は
のように「場の演算子の積
【定理】(Wickの定理)
ここで、「可能な全ての組み合わせ」という時、重複を避けるため、常に全ての
ところで、相互作用項(被積分関数)に現れる場の演算子は全て時空上の同一の点におけるものであり、これらは全て同じ寄与をもたらします。すなわち、Wickの定理を用いて
有限温度系のグリーン関数
ここまでは、「絶対零度の相対論的な場の量子論」の枠組みの中での話でした。素粒子物理学の文脈ではこのような枠組みでの議論が基本です。一方で、物性物理学においても、物質の中を動くたくさんの電子を場の描像で捉えようという事をよく行います。物質中の電子のように「遅く、たくさん集まっている」ような多体系を場として扱うときは、有限温度の非相対論的な場の量子論の枠組みで議論するべきです。寒い冬には温かみのある話がほしいですよね。より、以下では物質中の電子を有限温度のもとで記述する非相対論的な場を構成し、そこからその場におけるグリーン関数を与えたいと思います。
非相対論的な電子の多体系の場
まずは非相対論的な運動をする電子の集団を記述する場を構成しましょう。電子は(非相対論的な)量子力学に従って運動します。一体のシュレーディンガー方程式
の一般の固有状態は、位置の固有状態を波動関数の重みによって足し合わせたもの
となっています。いま、この一粒子系の位置の固有状態
のように、
絶対零度の自由スカラー場の演算子がそうであったように、この演算子は「エネルギー固有状態
のように書けます。この時の演算子
および
となっています(今考えているのはフェルミオンの多体系なので、Pauliの排他律より反交換関係によって量子化をする必要があります)。場の演算子にスピンなどの内部自由度がある場合、異なる内部自由度の演算子同士は反交換するものとします。
場の演算子を用いることで、場のハミルトニアンおよび個数演算子は
のように書くことができます。相互作用項の形の詳細はひとまずは大事ではないので気にしないことにします。
グリーン関数の定義
場の演算子が与えられたので、ここからこの場に対するグリーン関数を定義したいと思います。ここでは前節の議論とは逆向きに、グリーン関数を「異なる時空上の場の演算子どうしの積の期待値(=場の演算子の2点相関関数)」として与え、それが実際に場の方程式のグリーン関数になっていることを確かめる、という順番で議論します。ここでいう場の運動方程式はシュレディンガー方程式です。
グリーン関数を定義するにあたって、まず場の演算子のハイゼンベルグ表示
を考えます(ここでは熱平衡の下でグランドカノニカル分布に従う系を取り扱うので、
とします(絶対零度の系とは期待値の計算方法が異なることに注意しましょう)。指数関数の中の
と書かれます。
この時、この場のグリーン関数を、異なる時空上の場の演算子同士の時間順序積の期待値として定義します:
諸々の記号(ハイゼンベルグ表示、熱平均、時間順序積)を省略せずに書くと、有限温度系のグリーン関数は
となります。この表式にトレースの巡回性(
となり(第2項も同様に変形できる)、時空間の異なる2点におけるグリーン関数が時間についてその差
さて、このグリーン関数は「異なる時空上の場の演算子どうしの時間順序積の熱平均」として与えたものでしたが、これがちゃんと微分方程式の意味でシュレディンガー方程式のグリーン関数になっていること、すなわち相互作用のない場の演算子に対して
となることが確かめられます(相互作用のある場の演算子に対してはこのように綺麗な形にならないことに注意しましょう)。この計算は、階段関数の微分(正確には「弱微分」)
および、ハイゼンベルグ表示した場の演算子の時間発展の式
を用いることで求められます。階段関数の微分に由来するデルタ関数は、場のハミルトニアンにおける積分の計算の際にうまく消えてくれます。
自由電子ガス模型のグリーン関数
絶対零度の相対論的なスカラー場の量子論の時と同じように、相互作用のない自由場の場合にこのグリーン関数の具体的な形を求めてみます。自由電子ガス模型の場合、エネルギー固有関数の完全系は波数によって指定される平面波解として与えられるので、場の演算子は「波数
のように展開されます(
と書けます(ここで
と計算できます。(2つ目の等号は、Baker-Campbell-Hausdorffの公式と呼ばれる恒等式から導かれます。このようにハイゼンベルグ表示を簡単な形で書けるのは自由場ならではの事です)。よって、グリーン関数の定義に上の表式を素直に代入して時間順序積の熱平均を計算すると
となります。ここで、自由フェルミオン系において、エネルギー
を用いています(記号
このグリーン関数は、時間および空間について並進対称(その差にしか依存しない)となっています。より、そのフーリエ変換を
のように定義すれば、階段関数の複素積分表示を用いることで
と具体的に求まります。
「温度グリーン関数」の導入
さて、有限温度系のグリーン関数というものを前節で定義し、自由電子ガス系の場合にその具体形も与えました。今、例えば系のハミルトニアンが自由電子ガス模型の項
となり、時間発展の部分だけでなく熱平衡分布の重みにも
そこで、グリーン関数の定義を少し変えて、摂動展開ができるような形に書き換えることを考えましょう。指数の中でハミルトニアンの前に実数
有限温度系のグリーン関数に
は温度グリーン関数、あるいは松原グリーン関数と呼ばれます[7]。
によって定義されます(虚時間形式において両者は互いにエルミート共役ではないので、
自由電子ガス模型の場合の表式は、先と同じように計算できて
となります。
温度グリーン関数の周期性
実は、温度グリーン関数は虚時間方向に特徴的な周期性を持っています。実際、虚時間方向の平心対称性から
となり、虚時間方向に
となります。さらに、
のようにできます。
温度グリーン関数の摂動展開
さて、ハミルトニアンに相互作用項がある場合の温度グリーン関数について見ていきましょう。再度述べますが、温度グリーン関数では指数の肩は「ハミルトニアン
ハミルトニアンが
ここで、「(相互作用表示における)虚時間発展演算子」を
として定義すると、相互作用表示とハイゼンベルグ表示の間の関係を
と書くことができます。これによって、温度グリーン関数を相互作用表示した場の演算子によって書き与えることができます。また、ハミルトニアンの相互作用項についても
と相互作用表示することとします。
実は(といっても計算でわかりますが)、虚時間発展演算子
なる式を満たしており、演算子の積の順序に注意してこれを解けば
と書くことができます。
さて、いよいよ相互作用項のある場に対する温度グリーン関数の表式を与えます。温度グリーン関数の定義
に登場するハイゼンベルグ表示の場の演算子を、先の関係によって相互作用表示のものに書き換えると
となります。上で求めた虚時間発展演算子の表式を代入してまとめると、これは
と書かれます。さらに、大分配関数
のように、一粒子項ハミルトニアンと相互作用項を分離した形で書くことができます。ここで、虚時間発展演算子の中の順序積
というわけで、我々の目的はこのようにして書かれる温度グリーン関数を摂動論的に求める事だということがわかりました。ところで、この温度グリーン関数の表式は絶対零度の場の量子論における
これより、この温度グリーン関数の摂動展開は、絶対零度の系における摂動展開と同じ要領で行えるのではないかと期待できます。さらに、この期待に応えるかのように、有限温度系の虚時間形式においては絶対零度の系に対するWickの定理と同じような主張が成立します(Bloch-De Dominicisの定理):
虚時間形式において、
相互作用のある場において, 温度グリーン関数の摂動展開の各項は
のような形になっています。さらに相互作用ハミルトニアンが場の演算子の多項式で書ける場合, これはたくさんの場の演算子の時間順序積の一粒子項ハミルトニアンのもとでの熱平均になっています。Bloch-De Dominicisの定理はまさにこのようなものを計算する時に威力を発揮します。Wickの定理もBloch-De Dominicisの定理も、
温度グリーン関数と解析接続
さて、有限温度系のグリーン関数に対しては摂動計算が困難でしたが、その定義を少し変更して作った温度グリーン関数は絶対零度の場の量子論の時と同じ要領で摂動展開の計算ができるということを見てきました。ですが、この温度グリーン関数は実時間において定義されたグリーン関数とは違って直接的に物理量とは結びつきません(そもそも「虚時間」って何なんでしょうか)。物理的に意味のある結果を得るためには、解析接続によって温度グリーン関数の計算結果を実時間の世界に「持ち込んで」あげる必要があります。以下でもう少し詳しく見てみましょう。
まず、温度グリーン関数に対して、「レーマン表示」と呼ばれる表示に書き直すことを考えます。これは、エネルギー固有状態の完全系を用意してトレースを具体的に書き下した後に(虚)時間方向についてフーリエ変換したものになっています。温度グリーン関数のトレースを固有状態の完全系を基底にとって計算すると、
となります。ここで、場の演算子が粒子の生成/消滅演算子で書けるため、エネルギー固有状態
となります(ただし、松原振動数について
さてここで、「遅延グリーン関数」と呼ばれる次のような量を定義します。
これは
となります。これと温度グリーン関数のレーマン表示を見比べてみると、両者はよく似ています(違うのは、分母が
というものを定義して、温度グリーン関数および遅延グリーン関数を
と書くことにします。これは、複素関数
を用いると、統一的に
と書けます。つまり、温度グリーン関数と遅延グリーン関数はこの
一般の相互作用のある場については、実時間で定義された遅延グリーン関数を摂動論的に計算することはできません。そこで、相互作用場に対する計算は次のようにして行います。
- 虚時間形式のもとで摂動論的に温度グリーン関数を計算し、それをフーリエ変換して
を求める - 1.の手続きで離散的な松原周波数
に対して得られた の値から、解析接続[8]によって関数 の定義域を複素平面 全体に拡張する - この
に対して、 と上半面から実軸に近づけることで の値を得る
このようにして、摂動論の適用できる虚時間形式における計算結果を、実時間に「持ち込む」ことが可能になります。
ところで、温度グリーン関数の値から複素数平面全体に解析接続する方法は一意ではありません。実際、「離散的な点
なる振る舞いを示すことがわかります(スペクトル関数の
となるようなものを選べば、そのような
まとめ
今回の記事の主題を一言でまとめると、「絶対零度の場の量子論におけるグリーン関数の手法と同じようなことを、有限温度の場の量子論に対して行うにはどうすれば良いか?」というものが適切だと思います。絶対零度の場の量子論においてもそうだったように、物性物理でも可観測量がグリーン関数と結びつくので、与えられた場に対するグリーン関数の計算を行うことはとても重要な意味を持ちます。この記事ではダイヤグラムを用いた摂動展開の系統的な計算手法についての説明を全て省いたため、これを読んだだけでは具体的な計算ができるようにはなりません。ですが、それでも「出来る」という事実を強調しておくことは大事だと思うので、この記事にはそういう意義があると勝手に思っています。
改めて述べますが、紙面の都合と私の能力不足によりかなり説明不足かつ冗長な内容になってしまったことをお詫び申し上げます。この記事を読む際は270度ぐらい斜に構えて読んでください。とりあえず、「なんか物性物理(電子の多体系の場)でも絶対零度の場の量子論と同じようなことができそうだ」という気持ちだけが伝われば、執筆者としては幸いです。
というわけで皆さん、良いクリスマスを。僕の恋人は理学書です。うそです[9]。
参考文献
- 小形正男「物性物理のための場の理論・グリーン関数」(SGCライブラリ-142、サイエンス社)
- 坂本眞人「場の量子論 不変性と自由場を中心にして」(量子力学選書、裳華房)
- 坂本眞人「場の量子論(II) ファインマン・グラフとくりこみを中心にして」(量子力学選書、裳華房)
恋愛に限定された話ではありません。 ↩︎
人によってはこれを見て身構えてしまうかもしれませんが、大丈夫です。 ↩︎
すなわち、Lorentz不変(或いはPoincaré不変)な作用によって記述される場の理論の事です。 ↩︎
気になる方は「漸近的完全性」という語で検索してみるとこの辺の事情(相互作用場の理論的取り扱いの難しさ)がわかるような気がします。 ↩︎
相空間内の各経路に対する重みを足し合わせる際、「ある時刻における点とその微小時間
後における点を結ぶ相空間内の経路」からの寄与を離散的に取り込んでから連続極限を取っているので、この積分測度のもとで計算すると必然的に時間順序積が得られる、という仕組みです。 ↩︎ 松原というのは松原武生先生(京都大学名誉教授)のことです。松原先生は、このようなグリーン関数を考えることによって場の量子論における摂動論の手法を統計力学に適用しました。 ↩︎
「解析接続って物理的な変数の定義域を拡張する操作だから、解析接続で得られるものって物理的に意味あるの?」とお思いの方もいるかもしれません。確かに。でも、関数の性質が良い限りは解析接続の手法は問題ないし、実際この手法で物理的に意味のある計算結果が得られているので、まあ大丈夫なのでしょうと思っています(数学的な部分の詳細には全く詳しくないです、ごめんなさい)。 ↩︎
ほんとだよ。 ↩︎