特集／

非可換ベリー位相

初 貝

安 弘
ことを考えると「位相」こそが量子効果の本質的

1. 量子力学と位相部分であるということができるでしょう。何らか日常の世界での古典的な粒子（質点）の運動の実験的観測を行ったときその時刻、場所に粒子はニュートンの運動方程式により十分正確に記が存在すると観測される確率は波動関数の絶対値

述されていますが、小さな系、極めて外乱を受の２乗として与えられます。その意味で絶対値はけにくい系などの極限的な状況における運動は

観測可能量となりますが、通常位相そのものは古

シュレディンガー方程式 と呼ばれる、より基本的

典的な対応物がある物理量として観測されることな法則で支配されていて、ニュートン方程式は近はありません。しかしこれは何も波動関数の位相

似的に成立するにすぎないことは実験的な事実では観測できないことを意味するのではなく、干渉す。金属中の電子は豆電球を点灯させるときには実験のような真に量子的な実験を行えば観測する

古典的な電荷の運動と考えられるわけですが、最ことができるわけです。例えば荷電粒子がソレノ

近話題となっているナノ構造デバイス中の運動はイド周りを運動する際、磁場に直接触れることは量子力学による記述を必要とするわけです。マクないにも関わらず、粒子の軌道が囲む磁束に対応

ロな世界での古典的運動が物理系の大きさが小さした干渉効果を予言するアハロノフ・ボーム効果

くなるにつれ量子力学による補正（量子補正）をの実験的検証等がこの種の実験に対応します。1, 2)

うけ、より小さなナノワールドでは補正では済ま

また近年その有効性が確認され、実現の現実的可ず、量子効果が本質的となります。今後の２１世能性が議論されるようになった量子計算機の実現

紀においては量子効果を本質的にもちいた、いわのためにはこの量子干渉効果が観測にかかること

ゆる量子デバイスが実現するでしょう。そこで実が必須の条件です。

現する高度情報化社会においては量子効果は机上の空論ではなく、本質的に現象を決定する重要な因子となることに疑いはありません。

2. 定常状態と断熱定理

この量子的な世界観にたてば、粒子の運動はニュートン方程式が粒子の座標の時間変化を記

波動関数 と呼ばれる「複素数」で表現されること

述するように、シュレディンガー方程式 (1) は波動となります。古典的な粒子の運動は、ある時刻の関数 Ψ(t, r) の時刻 t における空間の各点 r での粒子の場所を指定する座標である実数で表現され値を規定します。この波動関数を Dirac の記法でることと対比的です。複素数 z = reΨ(t, r) = !r|Ψ" と書いたときシュレディンガー

iθ

はその絶

対値 r と位相 θ という２つの実数の組で表される
数理科学

NO. 523, XXX 2007

方程式は次の形となります。

1

i!

∂
|Ψ(t)" = H(t)|Ψ(t)"∂t(1)

ここで H は古典的にはエネルギーに対応するハ
ミルトニアン演算子です。もし H(t) が時間に依存しないときには波動関数に |Ψ(t)" = e−iωt |ψ"

時間依存性を変数分離した形を仮定するとシュレ
ディンガー方程式は H|ψ" = E|ψ", (E = !ω)図 1 断熱定理：断熱ポテンシャル中を急いで動かせばどこかに行ってしまうかごの鳥もゆっ
くり動かせばかごの中から逃げられません。

という固有値方程式に帰着します。自明な時間依存性をもつこの波動関数 |Ψ(t)" により記述される状態は定常状態とよばれます。以下、束縛状態の束縛状態にあるとすれば、H(Rt )|Ψ(t)" = E(Rt )|Ψ(t)".

(3)とよばれる空間的に局在した粒子的な状態に関して考察を進めてみます。この時、エネルギーはE = E1 , E2 , · · · と離散的な値をとります。ここ( R0 = R(0), Rt = R(t)) と時刻 t でもスナップ
ショットハミルトニアンの固有状態となります。

で少し条件をゆるめて物理系は時間に依存するパラメターの組 R(t) = (x1 (t), x2 (t), · · · ) により

3. ベリー位相とゲージ変換

しないがパラメター R(t) を経由して H(R(t)) と式 (3) は、断熱過程においては波動関数 |Ψ(t)"指定されハミルトニアンは明示的には時間に依存
時間変化するとしましょう。ここでシュレディン

が一般には無関係なはずのスナップショットハミ
ルトニアン H(Rt ) の固有値 E(Rt ) の固有状態

ガー方程式による時間発展とは独立にであることを意味します。固有方程式は斉次の線

H(R)|R" = E(R)|R"(2)という各時間（スナップショット）ごとの固有値問題を考えることができます。ここでパラメター

R 依存性を明示的に示しました。また各固有状態
は !R|R" = 1 と規格化されているとします。

一方、波動関数の時間発展はこの固有値問題とは一般には無関係ですが、断熱定理 とよばれる重要な定理はつぎの事実を主張します。3) 「スナップショットハミルトニアンの特定の束縛状態から時間発展を開始し、各固有値が孤立したままであり、
なおかつパラメターの変化が十分にゆっくりしているとき、シュレディンガー方程式の解である波動関数は各時刻ごとのハミルトニアンの固有状態であり続ける。∗1）」
つまり、この断熱過程において |Ψ(0)" = |R0 "と初期状態（時刻 0 とします）で固有値 E(R0 )*1） 特異性のないハミルトニアンに対して、注目する束縛状態と近接する状態のエネルギー間隔 ∆、パラメター変化の典
型的時間 T として T ∆ >> ! が必要である。正確には文
献3) 。

2

形方程式なので解の定数倍もまた解ですから、たとえ規格化の条件 !Ψ(t)|Ψ(t)" = 1 を要求してもその位相は断熱定理だけでは決定できません。しかしシュレディンガー方程式の解を次の形で探してもよいこととなります。

|Ψ(t)" = |Rt "U (t)(4)

あらかじめ準備しておいた |R" からの位相のずれ
が U というわけです。ここで規格化の条件から|U |2 = 1 が必要です。これをシュレディンガー方
程式に代入することでこの位相因子 U = eiΓ は次の方程式を満たすことがわかります。∗2）

U̇U −1 = −!Rt |Ṙt " − iE(Rt )/!

(5)

これを Γ について解けば、Γ(t) = γ(C) + Γ0 (t),!t
Γ0 (t) = − 0 dt E(Rt )"γR (C) = −iAµ dxµ(6)CAµ = !R|∂µ |R"`´*2） i! |Ṙt #U + |Rt #U̇ = E|Rt #U(7)

∂となります。ここで微分は ∂µ = ∂Rというパµ

ラメター空間内での微分を表し、系は時刻０からt までにパラメター空間上の経路 C を動くとします。Γ0 (t) は定常状態の位相変化に対応するよく知られた項ですが、これ以外の経路 C に依存する位相である γ(C) をその発見者にちなんでベリー位相と呼びます4) 。
ここでいくつかの注意が必要です。まず、固有値問題 (2) は斉次の線形方程式ですから、上述のよ
うに |R" 自体の位相も完全に任意であり不定です。
そこで |R" の代わりに、位相が ω = eiθ , (|ω|= 1)

異なる、次の |R$ " を使うこともできます。

|R" = |R" ω = |R" e
$

(8)

$ iθ

図 2 風向きの観測データをあつめた気象官はこの経路の中に低気圧が一つあると判断できるのでしょうか？

ここでの議論において極めて本質的であることに注意しましょう。例えば (7) などで微分の概念が出てきますが、定義に立ち返れば∂µ |R" = lim

h→0|R(x + hx̂µ )" − |R(x)"(11)
hと考えるわけですが、勝手な２点での位相はまったく不定なのでこの極限は単純には存在せず微分すら定義されていないわけです。微分の定義のためには、経路上でなめらかな固有関数 |R" を構成しておかなければなりません。上で述べたように固有関数の位相を変えることがゲージ変換をも

これを直接代入するとたらしますから、これはなめらかなゲージを一つAµ = A$µ + ω −1 ∂µ ω = A$µ + i∂µ θ(9)と書けますが、これは Maxwell 方程式系でのゲー
ジ変換とまったく同型の変換です。つまり異なる

固定することに対応します。つまり、なめらかなゲージ固定ごとにベリー位相 γ は確定し、ゲージに依存してベリー位相は 2π の整数倍の不定性を位相の固有関数を用いることは Aµ に関してゲーもつこととなります。この不定性は離散的ですの

ジ変換を引き起こすのです。この Aµ は電磁場のベで、連続変形できるゲージ変換で移りあうことは

クトルポテンシャルとは直接関係ないものですが、なく、この不定性をもたらすゲージ変換は「大き逆にベクトルポテンシャルの拡張としてのゲージな (Large) ゲージ変換」と呼ばれます。具体例にポテンシャルとして新しい意義をもちます。(R は対して数値的にこの微分を計算することを考えれ電磁場のように４次元である必要はなく、いくらば、ここでのゲージ固定の重要性は明らかでしょでも高次元のものが考えられます) このゲージ変う5) 。ゲージ固定により Aµ を確定したことはパ

換 (8)(9) の下でベリー位相 γR は γR = γR! + ∆ラメター空間上異なる２点間の相互関係である一と ∆ だけ変更を受けます。

∆ = −i"ω −1 ∂µ ωdxµ =

C

"

dθ

(10)

θ(C)

このように γR は固有状態の位相の選び方に依存するので一般には意味を持ちませんが、
経路 C がパラメター空間で閉曲線をつくる、
つまりパラメターの値がもとに戻る時には、変換（行列）ω = eiθ の一価性 の要求の下で ∆ は2π の整数倍となることに注意しましょう。4, 5) つまり閉曲線 C に対しては ベリー位相そのものγR(C)には 2π の整数倍の不定性がありますが、位相因子として eiγR は確定し、実験による量子的な観測量となり得ます。更に、この固有関数の不定性は
数理科学

NO. 523, XXX 2007

つの「接続」が与えられていると考えられますのでこの Aµ はベリー接続とよばれることもあります。6)基準 (ゲージ）の取り方は全く任意ですが、この基準に特異性があればその分だけ観測に不定性が生じるわけです。この状況を異なる例ですが、類似のわかりやすい例をとりあげて説明してみましょう。地球上図 (2) のように閉曲線上で風向きを測定することから閉曲線内に低気圧があるかどうかを決定することを考えてみましょう。この例
では中に１つ風の吸い込みとしての低気圧がある
ことはほぼ自明のように思われますが、この思考過程ではだれでも隣り合った地点での風向きを比

3

N
N

S

S

N
N

N

S

S

N

S
N
N

S

N

N

S

N

S

N

N

S

S

N

S

S

N

N

S

S

図 4 北極近くで同じ観測をやったらどうでしょう？この時、どこで観測したのかを知らな
い（新米の）予報官は低気圧がこの経路内にあることに気づかないかもしれません。

図 3 方位磁石を頼りに風向きを測定したデータを聞いた日本の（新米の）予報官は自信を持って低気圧がこの経路内に１つあると判断するでしょう。

整数倍だけ異なることになるわけです。ベリー位
相の 2π の不定性はこのアナロジーにより理解することができます。

較していることに注意しましょう。比較するため

固有ベクトル |R" は位相の自由度、つまりゲーには隣接２点間の風向きを示す矢印を重ねるわけジの自由度の分不定ですが、固有空間への射影演ですが、よく考えると重ねるための平行移動の過算子 P = |R"!R| はゲージ不変量となりますので∗3）

程に任意性があることに気づきます。

このようなアプリオリな平行移動の概念なしに低気圧の存在を確定するには例えば方位磁石を基準にすればよいことに気がつきます。つまり大局的に与えられた基準である方位磁石（ゲージ固定）と比較す

これとパラメター空間内の経路上で一意的な状態|φ"（例えば固定された状態）を用いて−1/2|R"φ = P |φ"Nφ, Nφ = |!R|φ"|2 (12)を用いれば、 Nφ $= 0 の条件下でゲージ固定されることで船の上で風向きのデータを集めるわけでたベリー接続が定義できることとなります。5) こす。例えば、北極から十分離れた日本近海で観測の Nφ = 0 を与える点が低気圧の例での「北極」した場合には、図 (3) のようになります。この時、に対応し、特定の |φ" によるゲージが特異的にな

方位磁石を基準とした閉曲線上での風向きデータ

る点を与えます。実は、この特定のゲージが特異だけから（大局的考察なしに）この閉曲線のなかになる点はチャーン数とよばれる有名な位相不変に低気圧が１つは最低限存在すると判断すること量を特徴づけることともなります4, 5, 8, 9) 。

ができるでしょう。つぎに図 (4) のように観測地
点が北極近辺で閉曲線が北極を囲んでいる場合を

4. 幾何学的位相としてのベリー位相

考えてみましょう。この時の方位磁石による風向きの観測データは、常にほぼ北風が吹いているこ

ベリーは断熱過程を考察することによりパラメとを示しますから、そのデータのみを使えば、低ター空間内の経路のみに依存するいわば幾何学的気圧の存在には気がつかないかもしれません。

な起源を持つ位相を発見しました4)。しかし歴史を

このように、大域的な基準であるゲージ固定の振り返れば、同様な起源を持つ位相は他にもあり、方法によって、風向きが定義する回転数は 2π の量子世界においては極めて基本的な意義があることが近年認識されるようになりました。∗4） これら*3） 実際の観測は決して連続した点すべてで行うことはでき

ず、離散的な点でしか行えませんが、その離散的な点を用いても同様な考察ができることが格子ゲージ理論との関連から議論されています。5, 7)

4*4） より広く光学および古典論においても類似の概念が重要な役割を果たすの現象もあることがベリー位相の発見以降、再認識されるようになりました。10, 11)

は総称して 幾何学的位相 と呼ばれ、広くベリー位相と呼ばれることもあります。例えば、「アハロノ
フ・ボーム効果」はベリー位相として幾何学的位相の典型例と見なすことができます。より古くはディラックの磁気単極子の発見12) , ならびに Wu-Yang
による磁気単極子のゲージ理論としての幾何学的意義ならびに数学的構造（ファイバーバンドル）

図 5 断熱ポテンシャル中をゆっくり動かせば３羽の鳥は逃げられないわけですが、かごの中を動いて入れ替わることは許されるでしょう。

の記述13) がベリー位相の発見以前から知られていました。さらには一般のゲージ場ならびに接続の的には内部自由度と呼ばれます。ディラックによ理論は幾何学的観点から数学のみならず物理としるスピン自由度の導入、通常のの電磁場の拡張とてもベリーの発見時に既に一般的に整理された状しての非可換ゲージ場の導入がその典型的な例と況にありました。14) ただしベリー位相の発見によ見なせるでしょう。

り、このゲージ構造は量子論至る所に極めて普遍

今まで U (1) とよばれる位相の自由度に関する的に存在することが広く認識されるようになった

ゲージ構造がベリー接続として束縛状態の断熱的のだとは言えるでしょう。また、ベリー位相の発時間発展に現れることをみてきましたが、断熱過見と同時期に量子ホール効果においてホール伝導程を少し拡張した状況で考えるとより一般に非可度のきわめて高い量子化はホール伝導度の幾何学

換ゲージ場としての構造がとてもわかりやすい形

ベリー

で現れることがベリーの議論の直後に示されまし位相の定義の際、仮想的なベクトルポテンシャル

た19) 。今まではスナップショットハミルトニアン

である Aµ が定義されたのですから、現実の電磁の束縛状態には縮退がないことを仮定してきまし

場の例を考えれば対応する磁場 B = rot A を考えたが、この状態が M 重に縮退することをゆるしることは自然でしょう。この磁場から定義されるてみましょう。この時、断熱定理はパラメターのパラメーター空間全体での全磁束が前述のチャー

変化が十分ゆっくりであればこの縮退した固有空ン数と呼ばれる整数の位相不変量となり、この整間内から出発した状態はシュレディンガー方程式

数値が量子化されたホール伝導度を与えるのです。

による時間発展によってもこの固有空間から外にこの発見はベリー位相の発見とともにその後の幾でることはないことを主張します。これによって

何学的位相及びトポロジカルな物理量が物性論にパラメター空間内での閉曲線上を系が動くとき、おいて重要な役割をはたす種々の現象の鋳型とし系の状態は位相変化を受けるだけでなく M 個のて位置づけることができます。例えば異方的超伝

縮退した状態の適当な線形結合の状態に変化する

導体、スピン軌道相互作用が重要な半導体におけことが許されるわけです。つまり縮退のないとき

るスピンホール伝導度の議論はその典型的な応用

は位相変化を受けるだけでしたが、多成分の場合

例と理解できます16〜18) 。

は固有状態の組はユニタリ変換の作用を受けるこ

8, 15)的意義にあることが発見されました。

ととなります。例えれば図 (5) のように断熱的にゆっくり鳥かごを動かせば鳥は逃げられないがお5. 非可換ベリー位相

互いの場所を入れ替えることはできるというとこ物理学においてはある理論が特定の現象を記述

ろでしょう。この時、M 個の状態をまとめてするとしても普通は必ずその適用限界があり、その限界を乗り越えるために理論の拡張が要求されます。そこでの理論的整合性を最も簡単に保つ方
法は一成分の物理量を多成分とすることで、物理
数理科学

NO. 523, XXX 2007

|R" = (|R"1 , · · · , |R"M )と書いて (4) の |Ψ" としても縮退した固有状態を一緒に考えるのが便利です。つまり波動関数とし5

て多成分のものを考えることになります。この縮過程における時間発展現象のみならず、様々な物

退した空間を内部自由度とみなすわけです。対応

理現象、物理系に対して幾何学的位相の観点からして (5) での時間発展を記述する U (t) がユニタリ

普遍的な新しい概念と解釈を与えることがわかっ行列となります。行列は一般に非可換ですから逐てきました。その応用の範囲は極めて広いのです

次近似で積分する必要がありますが、経路順序積

が、いくつか例を挙げれば、平行スピン対を許す∗5）

P を使って U は次のように書き下せます。
"#$
U = eiΓ0 P exp − i
Aµ dxµC

Aµ = !R|∂µ |R"
!R1 |∂µ |R1 "

..
=
.


非ユニタリ状態等の一般的な異方的超伝導相の分類20) 、グラフェンと呼ばれる炭素２次元単層系でのディラック電子の量子ホール効果21) 、
（分数）量

(13)

!R1 |∂µ |RM "

..

.

!R1 |∂µ |RM "子ホール状態、スピン液体状態等、量子液体相におここで縮退の無いときと同様にスナップショット

し詳しく説明しましょう。ハニカム格子をつくるベの固有状態の選び方は任意であることに注意しま

ンゼン環を２次元平面に無限に敷きつめた系はグしょう。この任意性は縮退した状態間のユニタリ

ラフェンと呼ばれ、近年実験的に構成されたため変換に対応しますから (8) の ω がユニタリ行列

多くの興味を集めています。この系のフェルミエネ

(ω = ω

ルギー近傍の電子構造はエネルギー分散が運動量

···
..
.

!RM |∂µ |R1 " · · ·†−1) に拡張されることとなります。よって Aµ の変換は∗6）

ける新しい秩序であるトポロジカル秩序、量子秩序等の議論22) 等において、この（非可換）ベリー
接続は興味深い役割を果たすこととなります5)。
最後に非可換ベリー接続の適用例を一つだけ少変化に対して線形のものとなり、通常の固体中でAµ = ω −1 A$µ ω + ω−1 ∂µ ω

(14)

極めて有効な有効質量近似がそのままでは適用できずそのためゼロギャップ半導体とよばれる特異なとゲージ理論でよく知られたものとなりこの A µ が系となります。これはその低エネルギーの有効理

非可換のゲージポテンシャルを定義することになり論が、自由電子もしくはホールにより記述されず、ます。特にその U (1) 部分と呼ばれる Tr Aµ につある種の相対論的なディラック電子となることをいては iθ = det ω として Tr Aµ = Tr A$µ + i∂µ θ

意味します。この系に磁場をかけるとランダウ縮と (9) と同じ変換性を持ち、縮退のない場合とほ退は通常と同様に生じますが、そのエネルギー間

ぼ同様の意義をもちます。 この U (1) 部分だけを

隔は等間隔とはなりません。さらにその量子ホー議論する立場からは、完全な縮退は必ずしも必要

ル効果を考えるとホール伝導度に 特異な量子化則

でなく M 個の状態をまとめて考えることのみがが生じ (図 6)、また特異なエッジ状態が重要な役割

本質となります。鳥かごは適当に広くてもよいとを果たします。この系における量子ホール効果のいうわけです。断熱性の仮定はこの場合、M 個の研究の際に非可換ベリー位相、ならびにベリー接状態が他の状態からエネルギー的に分離されてい続が有効に使われました21)。グラフェンではそのることだけを要求し M 状態間の相対的なエネル

ハニカム格子に起因する格子の効果が本質的であ

ギーの縮退、分離は無関係となります。

り、実験で観測されるような弱磁場極限であって5)

近年の研究でこの非可換ベリー接続は単に断熱*5） 0 から t までの時間を順に t1 , · · · , tN と N 等分してR

P e−i C Aµ dxµ

Rt= T e−i 0 Aµ ẋµ dt= lim e−i∆tẋµ Aµ (x(tN )) · · · e−i∆tẋµ Aµ (x(t1 ))
N →∞

!

!

*6） ∂µ |R# = (∂µ |R# )ω + |R# ∂µ ω

6も、起因するエネルギー準位に極めて小さなエネ
ルギーギャップが多数存在します。そのため、個々の一電子状態ごとのホール伝導度を求めることには困難となりますが、フェルミ球全体としての多電子系のベリー接続を扱うことにより物理的なエ
ネルギーギャップが仮定できれば理論的取り扱い

30
20
10
-3

-2

-1

1

2

3

-10
-20
-30

図 6 グラフェンにおける量子化されたホール伝
導度を e2 /h を単位としてフェルミエネル
ギーの関数としてしめしたもの。青線部分は
電子またはホールによる通常の量子化ホール伝導度を示し、赤線部分のゼロエネルギー近傍ではホール伝導度の量子化則は特異なディラック電子型となります21) 。

がきわめて容易になります。技術的にはフェルミ
エネルギー以下の全ての一電子軌道からなる非可換ベリー接続を構成し、その U (1) 部分を計算することで多電子系のホール伝導度を求めることができることとなります。
非可換ベリー位相の解説を目標にその基礎から

8) D. J. Thouless, M. Kohmoto, P. Nightingale, andM. den Nijs. Phys. Rev. Lett. 49, p. 405, 1982.
9) M. Kohmoto. Ann. Phys. (N.Y.)160, p. 343, 1985.
10) S. Pancharatnam. Proc. Ind. Acad. Sci. A44, p.
247, 1956.
11) J. H. Hannay. J.Phys.A:Math.Gen.18, p. 221,
1984.
12) P. A. M. Dirac. Proc. R. Soc., A 133, p. 60, 1931.
13) T. T. Wu and C. N. Yang. Phys. Rev. D 12, p.
3845, 1975.
14) T. Eguchi, P.B.Gilkey, and A.J.Hanson. Phys.
Rep.66, p. 213, 1980.
15) トポロジカルな観点からのレビューとして例えば
Y. Hatsugai. J. Phys. C, Condens. Matter9,
p. 2507, 1997.
16) Y. Morita and Y. Hatsugai. Phys. Rev. B62, p. 99,
2000.
17) S. Murakami, N. Nagaosa, and S.-C. Zhang. Science5, p. 1348, 2003.
18) J. Sinova and et al. Phys. Rev. Lett.92, p. 126603,
2004.
19) F. Wilczek and A. Zee. Phys. Rev. Lett.52, p.
2111, 1984.
20) Y. Hatsugai, S. Ryu, and M. Kohmoto. Phys. Rev.
B70, p. 054502, 2004.
21) Y. Hatsugai, T. Fukui, and H. Aoki. Phys. Rev.
B74, p. 205414, 2006.
22) X. G. Wen. Phys. Rev. B40, p. 7387, 1989.
（東京大学大学院工学系研究科）

最近の研究例までご紹介してきましたが、パラメ
ターを含む量子系における固有値問題という極めて一般的であり量子物理学の至る所にあらわれる
基本的な設定のなかにゲージ構造がとてもわかりやすい形で存在することが理解していただければ、幸いです。今後ベリー位相等、幾何学的位相の現代の物理における重要な意義はより広範囲の現象の中で明らかとなっていくことは間違いないと思
われます。
参考文献
1) Y. Aharonov and D. Bohm. Phys. Rev.115, p.
485, 1959.
2) A.Tonomura and et al. Phys. Rev. Lett.48, p.
1443, 1982.
3) A. Messiah. Quantum Mechanics. Dover, 1999.
4) M. V. Berry. Proc. R. Soc., A392, p. 45, 1984.
5) Y. Hatsugai. J. Phys. Soc. Jpn., 73, p. 2604,
2004, 74, p. 1374, 2005, 75, p. 123601, 2006.
6) 例えば Y.-S. Wu and H.-Z. Li. Phys. Rev. B38, p.
11907, 1988.
7) T. Fukui, Y. Hatsugai, and H. Suzuki. J. Phys.
Soc. Jpn.74, p. 1674, 2005.

数理科学

NO. 523, XXX 2007

7