ベレジンスキー＝コステリッツ＝サウレス転移（ベレジンスキー＝コステリッツ＝サウレスてんい、BKT転移、英: Berezinskii-Kosterlitz-Thouless transition）、または、コステリッツ＝サウレス転移（KT転移）とは、統計力学の2次元XY模型において起こる相転移である。1971年にベレジンスキー[1][2]、1973年にジョン・M・コステリッツとデイヴィッド・J・サウレス[3]によって理論的に提案され、1978年にヘリウム4の超流動薄膜において実験的に観測された[4]。
概要

通常の相転移は、低温相では長距離秩序が存在し、温度上昇に伴い秩序が壊れて無秩序相へと移行する。対称性の言葉で言えば、低温相では系の持つ対称性が自発的に破れており、高温相では対称性が回復している。統計力学において相転移現象を記述するために用いられる模型の一つであるXY模型の場合、3次元空間においては通常の二次相転移が起こる。しかし、2次元空間においては、マーミン＝ワグナーの定理から、対称性が自発的に破れず、長距離秩序を持つ相は存在し得ないことが示されている。空間の次元が低くなるとゆらぎが大きくなり、秩序相が不安定となって相転移が起こらなくなるというのは一般的に知られている性質であり、例えば、2次元イジング模型や3次元ハイゼンベルク模型は二次相転移を起こすが、1次元イジング模型や2次元ハイゼンベルク模型は相転移を起こさない。しかし、これらとは異なる特殊な例として、2次元XY模型は低温相において通常の長距離秩序を持たない代わりに、特殊な秩序を持つことで相転移を起こす。これが、ベレジンスキー＝コステリッツ＝サウレス転移と呼ばれる相転移である。

一般的な相転移では、低温相における相関関数は一定の値を持つが、BKT転移の場合には距離に対してべき的に減衰する。このような減衰は本来は臨界点において起こるはずの挙動である。つまりこれは、BKT転移の転移温度以下の低温相においては、有限の温度領域に渡って臨界点としての状態が続いていることを意味している。この状態の秩序は準長距離秩序(quasi long range order)と呼ばれる。

BKT転移は、実験的には、ヘリウム4の超流動薄膜[4]、2次元超伝導体[5]、ジョセフソン接合列[6]、冷却原子気体のボース＝アインシュタイン凝縮[7]などにおいて観測されている。
準長距離秩序

2次元XY模型のハミルトニアンを H = − J ∑ ⟨ i , j ⟩ S i ⋅ S j = − J ∑ ⟨ i , j ⟩ cos ⁡ ( θ i − θ j ) {\displaystyle H=-J\sum _{\langle i,j\rangle }\mathbf {S} _{i}\cdot \mathbf {S} _{j}=-J\sum _{\langle i,j\rangle }\cos(\theta _{i}-\theta _{j})}

とする。ここで、Jは相互作用の強さを表す定数、 S i = ( cos ⁡ θ i , sin ⁡ θ i ) {\displaystyle {\boldsymbol {S}}_{i}=(\cos {\theta _{i}},\sin {\theta _{i}})} は2成分のスピンベクトル、θiは2次元空間でスピンの向きを指定する角度である。

この模型において、ゆらぎを持つスピン S i {\displaystyle {\boldsymbol {S}}_{i}} と S j {\displaystyle {\boldsymbol {S}}_{j}} についての相関関数は以下のように定義される。 G ( r i − r j ) = ⟨ S i ⋅ S j ⟩ = ⟨ cos ⁡ ( θ i − θ j ) ⟩ = ⟨ e i ( θ i − θ j ) ⟩ {\displaystyle G({\boldsymbol {r}}_{i}-{\boldsymbol {r}}_{j})=\langle {\boldsymbol {S}}_{i}\cdot {\boldsymbol {S}}_{j}\rangle =\langle \cos(\theta _{i}-\theta _{j})\rangle =\langle e^{i(\theta _{i}-\theta _{j})}\rangle }

ここで、 r i {\displaystyle {\boldsymbol {r}}_{i}} はi番目のスピンの位置である。最後の等式は、虚数部分の熱平均がゼロとなることを用いている。

高温領域の振る舞いは、高温展開を用いて調べることができる。このとき、相関関数はスピン間の距離rの関数として指数関数的に減少することが示される。 G ( r i − r j ) ∼ exp ⁡ ( − r log ⁡ 2 β J ) {\displaystyle G({\boldsymbol {r}}_{i}-{\boldsymbol {r}}_{j})\sim \exp \left(-r\log {\frac {2}{\beta J}}\right)}

ここで、βは逆温度である。

一方、低温領域については、スピン波近似を用いて計算される。系が十分低温である場合、各スピンのゆらぎは小さくなり、隣接するスピン同士の向きはほとんど同じとみなすことができる。このような近似のもとで、相関関数は G ( r i − r j ) ∼ ( 1 r ) 1 / 2 π β J {\displaystyle G({\boldsymbol {r}}_{i}-{\boldsymbol {r}}_{j})\sim \left({\frac {1}{r}}\right)^{1/2\pi \beta J}}

となる。この式は距離rに反比例する減少関数となっている。

通常の相転移では、秩序相であれば長距離秩序が存在し相関関数は減衰せず、無秩序相であれば相関関数は指数関数的に減衰する。つまり2次元XY模型の低温領域における相関は、距離が離れるにつれてゆっくりと減衰するような特殊な秩序を持っている。このような秩序は準長距離秩序と呼ばれる。
量子渦のダイナミクス
渦の役割

2次元XY模型の相転移においては、量子渦のダイナミクスが重要となる。

BKT転移を示す多くの系では、低温においては渦とその渦度が逆符号の渦とが束縛されたペア（渦対）となって存在し、温度上昇に伴い無秩序相へ移ると、渦対は束縛されていない単独の渦2個へと解離する。これは、低温においては単独の渦の存在は不安定となっているが、渦対として存在することは可能であり、逆に高温では、単独の渦の存在が安定となるためである。これより、転移温度より低温のBKT相では、束縛された渦対のみが存在する。ただし、低温相における渦対はスピン波による準長距離秩序にほとんど寄与しない。一方で、転移温度より高温では、単独の渦が大量に発生することで、スピン間の相関は指数関数的に抑制され、秩序は失われる。
単独の渦のエネルギー

BKT転移のメカニズムを解明するために、渦生成に必要なヘルムホルツの自由エネルギー（つまり、エネルギーとエントロピーとの競合）について考える。

単独の渦1個（簡単のため、渦度を1とする）が持つエネルギーは、スピン波近似のハミルトニアンから、 E ∼ J 2 ∫ d 2 r ( ∇ θ ) 2 = J 2 ∫ a L d 2 r 2 π d r 1 r 2 = π J log ⁡ ( L a ) {\displaystyle E\sim {\frac {J}{2}}\int {\mathrm {d} }^{2}r(\nabla \theta )^{2}={\frac {J}{2}}\int _{a}^{L}d^{2}r2\pi {\mathrm {d} }r{\frac {1}{r^{2}}}=\pi J\log \left({\frac {L}{a}}\right)}

と表せる。第2式の動径方向の積分は、格子間隔（渦芯の半径）aを下限として系全体の半径Lまでの範囲で行っている。