※本記事は
インスタントンとトンネル効果 (1) : 経路積分とWKB近似
の続きです。
今回は二重井戸型ポテンシャルにおける基底状態のエネルギーの評価法について書きます。この系にはEuclid空間における古典解であるインスタントン解が存在します。この解による分配関数を用い、WKB近似におけるの表式を求めます。
具体的なインスタントン分配関数の計算は次回行います。
本記事はRef.Sakitaに基づきます。またこの話題に関してはRef.Colemanが有名ですので文献として載せておきます。
二重井戸型ポテンシャルにおけるインスタントン
次の量子系を考えます(質量は1にしている):
ポテンシャルは以下のような形をしています(図1左):
二重井戸型ポテンシャル。が小さくなると最小点を隔てるポテンシャル障壁が大きくなる
このポテンシャルはにおいてでのポテンシャル障壁が無限に高くなり、2つの独立した無限に深い井戸の系になります(図1右)。この場合のエネルギー準位は簡単に求めることができます。としてでのポテンシャルはのとき
です。これは調和振動子のポテンシャルであり、エネルギー準位は
となります。すなわちでは、に周期の調和振動子が独立に存在します。よってこの極限における基底状態のエネルギーは
であり、二重縮退しています。
しかしが有限のとき、実際には付近の波動関数はトンネル効果で双方に漏れ出し、これによりエネルギーの縮退が解けます。以下、この基底状態のエネルギーをインスタントン解を用いてWKB近似で求めます。
基底状態のエネルギーの求め方
ここではの計算を以下のように行います。次の量
を導入します。このとき
ですが、のとき、カッコ内のの和はの寄与で支配されます。よって
ゆえに
です。ここでは虚時間における分配関数であり、経路積分で以下のように表されます(Appendix参照):
前回の記事
より、はと同じ役割を担っており、のときWKB近似が使えることがわかります。WKB近似ではではなくで議論することに注意してください。
インスタントン解
におけるWKB近似を用いるために、に寄与する古典解を求めます。
Eq.から得られる古典的な運動方程式は以下です:
ここで微分はによる二階微分を表します。この運動方程式は、通常のポテンシャルを逆さまにした系の運動に対応します(図2)。
ポテンシャルを逆さまにしたグラフ
さらに境界条件より、時刻がからまで経過すると、元の位置に戻る運動を考えます。
この系に存在する古典解を考察します。
まず
という解が存在します。これはずっと2つの山の頂上に留まる解です(図3)。
に留まる解
これは上記方程式の運動としては不安定なので、物理的ではないと思うかもしれませんが、元の系では安定な停留点の寄与を与えます。ポテンシャルがひっくり返っているのはあくまで計算上のテクニックであり、実際に物理量に寄与するか否かはこの系の運動の直感から類推してはいけないことに注意してください。
この系には更に、無限の過去にどちらかの山頂に存在し、無限の未来にもうひとつの山頂にたどり着く解が存在します(図4)。
無限の過去にに存在し、無限の未来でに辿り着く解
この運動に対応する解は以下のように表されます:
グラフにすると下図のようになります:
インスタントン解(左図)と反インスタントン解(右図)
からに辿り着く解をインスタントン、からに辿り着く解を反インスタントンと言います。これらが運動方程式を満たすことはすぐに確かめられます。またこの解にはという任意パラメータが存在し、谷を通過する時刻に対応します。この任意パラメータの存在は時間並進不変性の反映です。
ただし、より、片方の山から転がり落ちてもう一度同じ山の頂上に辿り着く解のみが分配関数に寄与します。インスタントンではそのような解を完全に記述することはできませんが、近似的にはインスタントンと反インスタントンをかけたものがそれに対応します。図6からわかるように、インスタントン・反インスタントンの時間方向の距離が十分離れているときに近似が良くなります。
インスタントンと反インスタントンの積
インスタントンペアの分配関数への寄与
ある山頂からもう一方の山頂に戻る解はすべて積分に寄与します。インスタントン・反インスタントンの組の数をとして、それらをすべて足しあげることを考えます。
インスタントンペアの分配関数の足し上げ
の場合に関し、これがどのような積分の寄与をもたらすか考察します。の解は図7の一番左のグラフです。まず、インスタントンと反インスタントの分配関数は同じなので、以下これらをと表記します。ここでインスタントンと反インスタントンが存在している領域を1および2とラベルします(図8)。経路積分およびの肩の虚時間の積分もこれらの領域に分割します。
領域を1と2にわける
では、の範囲に比べ、(反)インスタントンが隣の山に移る時間間隔は非常に狭くなります。この状況では、全領域にインスタントンが1つ存在する状態の分配関数を、に留まり続ける状態の分配関数で割ると、領域1にインスタントンが存在する分配関数を、領域1でに留まる分配関数でわったものになります(図9)。
インスタントンの分配関数をの停留点に留まる分配関数でわる
式で書けば
となります。ここで積分記号の上についている数字は領域を、下についているカッコつきの記号はインスタントンの数を表しています。同様に、全領域に反インスタントンが1つ存在する状態の分配関数をで割ると、領域2に反インスタントンが存在する分配関数を、領域2においてに留まる分配関数でわったものになります。
以上から、インスタントンペアが1組ある場合の分配関数は
になります(図10)。
インスタントンペアの分配関数の計算
この議論は一般のコのインスタントンペアが存在する場合に容易に拡張できて
となります。
ただしこれらの関係には更に係数がつきます。なぜなら、例えば1組のインスタントンペアの場合、では1の領域にインスタントンおよび2の領域に反インスタントンが存在するという制限があるのに対し、ではそのような制限がつかないからです。この制限は、分配関数を計算する際の(反)インスタントンの中心座標の積分の制限に対応します。制限つきのコの(反)インスタントンの中心座標の積分は
です。これと制限なしの積分
との関係は
です。なぜなら、Eq.の積分は一辺がの次元超立方体の体積ですが、これをを満たす領域に分解したとき、その領域の数はからにからを割り当てる通りの数であり、また各分割領域は同じ体積だからです。
以上より
を得ます。最終的に計算したい分配関数は、ペアの個数に関して和をとり、かつ又はから出発して元に戻る2通りの場合があることから
となります。の制限があるため、ペアで足し上げなければならないことに注意してください。
より、は
はに留まり続ける解の周りの量子ゆらぎをまで取り入れた場合の分配関数であり、この運動は調和振動子であるから
であり、基底状態のエネルギーはになります。よって
を得ます。第2項がトンネル効果による基底状態のエネルギーのずれです。
改めて結果を公式の形で記しておきます。
二重井戸型ポテンシャルの基底状態のエネルギーは、でのWKB近似においてインスタントンペアの効果を足し上げることで計算すると
となる。
今回はここまで。おしまい。
【次の記事】
インスタントンとトンネル効果 (3) : ゼロモードの除去、トンネル効果によるエネルギーシフト
Appendix: 統計系の分配関数の経路積分表示
統計力学的分配関数が経路積分表示を持つことを、以下簡単に紹介します。
分配関数は以下のように定義されます:
は同じ状態ではさんで和を取ることを意味します。これは基底に依存しない量です。
基底として位置の固有状態をとります。すると
となります。ここでとすると、は虚時間における時間間隔の時間発展演算子とみなすことができます。よっては
という遷移振幅とみなすことができます。ただし初期状態と終状態は同じの状態であるという条件がつく(およびその状態で和をとる)ことに注意してください。つまり、虚時間をで定義すると、とで周期境界条件が課されることになります。
こうするとあとは通常の経路積分表示の導出と同じ議論により、を経路積分で表すことができます。虚時間を等分し、時間発展演算子を等分した微小時間間隔の時間発展演算子を導入します。ここではの刻み幅であり虚数です。虚時間(実数!)の刻み幅との関係はです。各虚時刻においての完全系を挟み、行列要素を計算すれば
のように表せます。Hamiltonianをと仮定し、の積分を実行すれば、本文のEq.を得ます。詳しくは例えばKugoをご参照ください。
は統計系では(は温度)であり、それゆえは温度を伴う理論の分配関数を表しています。大学で習うような物性系ではスピン系が扱われることが多く、基底として位置の固有状態をとることはあまりないかと思います。しかし上記の議論を用いると、例えば場の量子論にも温度を導入することができます。このようなformalismは松原武生等により議論され、虚時間形式とか松原形式と呼ばれます。