0
大学数学基礎解説
文献あり

超伝導体におけるvortex解

172
0

はじめに

最近LK-99という常温常圧超伝導体の候補が話題になりました。ネット上で大きな盛り上がりを見せ、関連株の上昇も起こったほどなのですが、他グループによる検証結果は尽く否定的です[1]

LK-99が常温・常圧下で磁石に浮く動画がいくつか違うグループから公開されています。超伝導体が磁石の上に安定して浮く現象は、マイスナー効果と「ピン留め効果」の結果であることはよく知られています。一方LK-99の動画を見ると、完全には浮いておらず、一部は必ず接地しているように見受けられます。よって反発力は働くが安定はしておらず、少なくともピン留めは起こっていないようです。(※ピン留めにも見える"LK-99"の浮上動画も1つだけありますが、果たして...)

磁束の侵入・常伝導スポットとピン留め効果

このピン留め効果とはなんでしょうか[2]

大前提として、超伝導体は、その状況下で常伝導よりも超伝導であるほうがエネルギーが低いために超伝導状態になっています。

超伝導体を磁場の下に置きます。マイスナー効果により超伝導体は磁場を排斥しますが、排斥するにもエネルギーコストがかかるので、ある程度大きな磁場にさらされると常伝導になります。第1種超伝導体(Type I)は、臨界磁場に達すると全体が完全に常伝導になります(図1上右)。一方第2種超伝導体(Type II)には下部臨界磁場というものが存在します。これを超えると磁場の一部が超伝導体内に紐状・渦状で侵入し、侵入した部分だけ常伝導になります(図1下右)。そしてそれ以外の部分は超伝導を保ちます。これによりType IIは比較的強い磁場でも安定してバルクとしての超伝導状態を保ちます。

Type IとIIの違い。Type Iは臨界磁場を超えると全体が常伝導体になる。Type IIは下部臨界磁場を超えると、磁場が侵入した部分のみが常伝導体となり(青い部分)、それ以外の部分は超伝導体(オレンジ色の部分)を保つ。 Type IとIIの違い。Type Iは臨界磁場を超えると全体が常伝導体になる。Type IIは下部臨界磁場を超えると、磁場が侵入した部分のみが常伝導体となり(青い部分)、それ以外の部分は超伝導体(オレンジ色の部分)を保つ。

ところがType IIに電流を流すと、磁束と電流の間にローレンツ力が働き動いてしまいます。これは電流に対する抵抗となって現れ、せっかくのマイスナー効果よる抵抗ゼロの性質が破れてしまいます。

ここでType IIに常伝導のスポットが点々と存在しているとします。磁場が侵入した部分は常伝導になり、超伝導よりもエネルギーの高い状態になります。それなら磁場はなるべく常伝導のスポットを通ったほうが全体のエネルギーが低くなります。これにより、磁束線は常伝導スポットに「ピン留め」され、電流を流しても磁束は動かず抵抗ゼロの状態が保たれます。

ピン留めの概念図。青い部分は常伝導状態。エネルギーを下げるため、磁束(黒い線)は常伝導スポットを通る。 ピン留めの概念図。青い部分は常伝導状態。エネルギーを下げるため、磁束(黒い線)は常伝導スポットを通る。

さらにはこの効果により、超伝導体が磁石に反発しても安定的に浮くことができます(注1)。

本記事の内容:Nielsen-Olesen vortexの導出

本記事では、超伝導体内に侵入しする磁束のモデルであるNielesen-Olesen vortexに関して説明します。(※ピン留めには言及しません)

Type IIに磁場が侵入する機構は、Ginzburg-Landau理論に電磁場を結合させた理論により記述できます。これは"Higgs場"にU(1)ゲージ場が結合する、Abelian-Higgs模型と呼ばれる系に等価です。この系には渦状の解が存在し、Nielsen-Olesen vortexと呼ばれます。この解に関する位相の条件から、磁束の量子化が示せます。

これは第2種超伝導体内に磁場が渦として侵入し、またその磁束が量子化されることのモデルとなっています。

本記事はRefs.[3][4][5]を参考にしています。

Abelian-Higgs模型

Ginzburg-Landau理論(GL理論)というものがあります。この理論では、系の対称性が自発的に破れているとき、系の自由エネルギー等を対称性の破れのオーダーパラメータのべきで記述します。そして破れが小さいとしてその低次の項のみを取り出し解析します。非常に広範に適用できる理論・考え方です。物性物理・素粒子物理・宇宙論等様々な分野で有用です。

これから考える系はAbelian Higgs模型(AH模型)と言われる以下のような系です:

Abelian Higgs模型
L=14FμνFμν+12|Dμϕ|2V(ϕ),   V(ϕ)=λ2(|ϕ|2ϕ02)2,Dμ:=μiqAμ
ここで
Fμν:=μAννAμ,  ϕ:スカラー場

この模型はU(1)ゲージ場Aμにスカラー場ϕが結合している系です。これはGL理論とみなすことができ、その場合ϕがオーダーパラメータに対応します。超伝導の場合ϕはU(1)対称性の自発的破れのオーダーパラメータである、電子のCooper対の期待値に対応します。Higgsの名前が冠されているのは、ϕがHiggs場と同様の役割を果たすからです。

でもとりあえず物理的なことは忘れ、この系の解を探します。

Abelian-Higgs模型におけるvortex

AH模型における渦状の解、vortex解に関して考察します。

運動方程式と軸対称なAnsatz

Euler-Lagrange方程式
{μL(μϕ)Lϕ=0,μL(μAν)LAν=0
に上記Langangianを代入して計算すると以下のようになります:
{DμDμϕ=2λϕ(|ϕ|2ϕ02)     (ϕ02:=μ2/λ)νFμν=Jμ,   Jμ:=12iq(ϕμϕϕμϕ)+q2Aμ|ϕ|2
Jμはカレントです。以下静的な渦状の解を求めることにして、t依存性は考えません。

渦は軸対称なのでz軸対称かつx-y平面の円対称解を考えます。そのため円柱座標(ρ,φ,z)に移ります(ρ=x2+y2,φ=tan1(y/x))。次のようなansatzをおきます:
{A=eφA(ρ),ϕ=f(ρ)einφ
すなわちゲージ場はeφ方向の成分のみを持ちかつρにのみ依存。ϕρ依存性の他にφに関する位相をもつとします。またnは、φφ+2πに対するϕの一価性より整数です。

EoMを円柱座標で書き直し、上記のansatzを代入すると以下を得ます:
EoM for ϕ:(2ρ2+1ρρ)f+(nρqA)2f+2λ(f2ϕ02)f=0,(1)EoM for Aμ:ρ[1ρρ(ρA)]+(qnρq2A)f2=0

境界条件

ρが大きい領域での解の振る舞いに関して考えます。この領域で解はエネルギー最小の状態=真空に近づかなければいけません。fに関してはϕVを最小にする条件から
f(φ)ϕ0
になります。一方Aμに関してはゼロなら電場も磁場もゼロでありエネルギー最小です。しかし実際にはAμ=0は強すぎる条件であり、Aμ=0とゲージ同値であるpure gauge form
(2)1qμχ
に近づけば電磁場はゼロになります。ここでχ=nφとすると、:=(/x,/y)eφφ方向の基底として
1qχ=nqρeφ
であるから、このとき
(3)A=nqρ
となります。よってこのAはpure gauge formでありエネルギーを最小にします。そしてEq.(1)をみると
(4)A=nqρ,   f=ϕ0
はEoMの解となっています。よってEq.(4)はエネルギー最小のEoMの解、すなわち真空解です。

さらにeinφは原点でsingularなので、f(r=0)=0を課すことで解を原点でもなめらかにします。

磁束の量子化

ある局面Sを貫く磁束Φを計算します。dSμνμ,ν方向の微小面積要素とします。するとStokesの定理より
Φ=SFμνdSμν=SAμdxμ
のようにΦをゲージ場の線積分で書くことができます。Aρの小さい領域での振る舞いはわかりませんが、Sρが十分大きいところにとれば、Φを真空解で評価できます。Sの経路をeρ,ezに直交するようにz=const.の円にします。すると
Φ=SAμdxμ=S1qμχdxμ=1q[χ(φ)]02π=1q[nφ]02π=2πnq
を得ます。ϕ2π回転に対する一価性からnは整数なので、これは磁束が量子化されていることを示します。これはトポロジカルな不変量です。

BPS limit

Eq.(1)の解を一般的に解析的に求めることは無理なのですが、BPS limitと呼ばれる極限では運動方程式が簡単になります。またこの極限は、AH系に限らず物理的に重要なことが多いです。そこでここではBPS limitにおいて運動方程式を導き、数値的に解を計算します。

そのために系のエネルギーを計算しておきます。z方向には並進不変性があるので、ここでは改めて空間2次元の場合を考えます。静的な配位を考え場は時間に依存しないとして時間微分は落とし、また後に上記の場のansatzを入れることを考慮しA0はゼロとします。すると
E=d2x(12F122+12|Dkϕ|2+V(ϕ))
ここでF12は磁場に相当します:B=F12。考察対象の2次元空間は、3次元中の2次元だと考え、磁場はz方向を向いているとします。

これを変形し「平方完成」すると以下のようになります(Bogomol'nyi completionと呼ばれます)
E=d2x{|(D1+iD2)ϕ|2+[F12+q(|ϕ|2ϕ02)]2}+2πnϕ02+λq22d2x(|ϕ|2ϕ02)2
2πnϕ02は量子化された磁束に関する項であり、nは整数です。この式からλq2ならば、エネルギーに関して下限が存在します:
E2πnϕ02
このようなエネルギーの下限をBogomol'nyi boundと呼びます。

ここで
λ=q2
の場合を考えます。このとき最後の項は消え、エネルギーは上記の下限を実現し得ます。エネルギーの下限を実現することを「Bogomol'nyi boundをsaturateする」と言い、このときEoMを満たす解はBPS解と呼ばれます。

BPS解を求めます。このとき第1項・第2項がゼロになるので
{|(D1+iD2)ϕ|2=0,F12+e(|ϕ|2ϕ02)=0
が成立します。極座標に移り((x,y)(ρ,φ))、前章の解のansatzを代入すると
A=1q(nρff)f=f2fffρ+q2(f2ϕ02)
を得ます。fρ微分です。上の式においてAf,fで書かれており、下の方程式はfのみで構成されています。下の方程式を解いてfを求め、その解を上の方程式に入ればAが求まります。前述したように、f(r0)=0f(r)=ϕ0が境界条件です。

数値計算と解に関する考察

数値的にこれを解き、プロットしたのが図3です。
BPS limitでの!FORMULA[86][1119053195][0]の振幅!FORMULA[87][37794][0]と磁場!FORMULA[88][36678][0]。パラメータは!FORMULA[89][-65762220][0]にしている。 BPS limitでのϕの振幅fと磁場B。パラメータはλ=q2=ϕ0=1にしている。

磁場Bρが小さい範囲のみで存在しており、紐状に絞られていることがわかります。逆に"Higgs場"ϕの値はこの領域で小さく、ρが大きい領域でϕ0=1に近づきます。

ここでAH系の対称性について述べておきます。AH系のLagrangianはU(1)のゲージ変換
AμAμ+1qμχ(x),   ϕeiχ(x)ϕ(x)
に関して不変です。しかしρが大きいところでHiggs場が真空の配位に近づくと、量子論的にはϕからその真空での値ϕ0を差し引いて、ポテンシャルVの最小の場所の周りで場を展開しなおさないといけません。よって新たにϕ~:=ϕϕ0という場を導入します。このときϕ~Aμの結合項により、ϕ~で書かれたLagrangianにおいて、ゲージ場に質量項が生じ、上記のU(1)対称性が破れます。そしてゲージ場が質量を持つことで磁場が排斥される真空が実現します。これがマイスナー効果に対応します(注2)。対称性が破れている相をHiggs phase、回復している相をCoulomb phaseと呼びます。このようなメカニズムをHiggs機構と呼びます。

詳しいことは省略しますが、磁場Bが減少する場所は、Higgs機構で生じたゲージ場の質量の逆mV1=1/(2qϕ0)付近です[5]。一方Higgs場ϕが増大するのは、Higgs場の質量の逆mH1=1/(2λϕ0)付近です。これらの値はBPS limit:λ=q2で等しくなります。図3で用いたパラメータでは、これらはどちらも1/2になります。そして実際数値計算を見るとその付近でB,fは共に減少・増大しています。

複数のvortexの相互作用はmV,mHの大小で定まります:

  • mV>mH: vortex間の相互作用は引力。超伝導体で言えばtype Iに相当する。複数のvortexは安定には存在できない。
  • mV<mH: vortex間の相互作用は斥力。type IIに相当。複数のvortexが安定して存在。
  • mV=mH: BPS limit。vortex間には相互作用が働かない。これはBogomol'nyi boundをsaturateしていることによる。

type IIの場合、vortexはその斥力により格子を形成します。これはAbrikosov latticeと呼ばれます。

ただしAH系は相対論的不変性を持つので、その点は実際の超伝導とは異なります。実際の超伝導体のような非相対論的な場合、type Iの状況ではそもそもvortex解が存在しません。

そのほかのこと

2つほどコメントです。

  1. 一般にBPS解の性質は有用です。例えばDブレーン(超弦理論の低エネルギー有効理論である超重力理論におけるソリトン解)のRamond-Ramond chargeという量を計算する際には、BPS状態のDブレーンの性質が利用されます[6]。またSeiberg-Witten理論において、 N=2超対称Yang-Mills理論における低エネルギー有効理論はBPS状態のゼロ質量粒子(モノポールとダイオン)で書けるという予想がされる際にも、BPS状態の性質が使われます。
  2. vortexの端点にはmagnetic monopoleが付随することができます。このmonopole-vortex系のエネルギーは、vortexの長さに比例します。このため有限のエネルギーではmonopoleを無限に引き離すことができません。これは有限のエネルギーで電荷・物体を無限遠まで引き離せる電磁気力・重力とは著しく異なります。このような現象を「閉じ込め」と呼びます。クォークはハドロン内に閉じ込められていますが、それを強い相互作用の真空の超伝導性で説明する理論が存在します。これを「双対超伝導描像」と呼びます[3]

まとめ

第2種超伝導体内に侵入した磁束のモデルであるNielsen-Olesen vortexに関して述べました。Abelian Higgs模型はGinzburg-Landau理論から導かれる有効作用です。解に関して軸対称のansatzをおき、遠方で真空の配位に近づく紐状・渦状のvortex解を探します。このvortexは磁束を持ち、それは量子化されることがわかります。

BPS limit:λ=q2では運動方程式が簡単になります。これを数値的に解き、実際にvortex解を求めました。ρが小さいところではBは有限、ϕはゼロ、ρが大きいところではその逆になります。Higgs機構によりAH系が持つU(1)対称性は自発的に破れ、ゲージ場は質量mVを獲得しますが、Bが減少する位置はこの質量の逆ρmV1=1/(2qϕ0)付近です。一方Higgs場ϕはその質量mHの逆ρmH1=1/(2λϕ0)のあたりで増大します。BPS limitではmV=mHです。

vortex間の相互作用はmV,mHの大小で変化します。mV<mHなら斥力であり、これがtype IIに相当します。一方mV>mHの場合引力であり、このとき複数のvortexは安定には存在できず、type Iに相当します。ただし実際の超伝導体のように非相対論的な場合、type Iに相当するケースではそもそもvortex解自体存在しません。

おしまい。




  • 注1: ピン留めという言葉、またそのイメージから推測すると、磁石の上に浮上している超伝導体は回転もできないように思えます。しかしこれは間違っていて、超伝導体は磁石上で非常になめらかに回転します。また磁石のレールの上を超伝導体がなめらかに走ることも可能です。
  • 注2: マイスナー効果のもうひとつの重要な事実に「電気抵抗がゼロになる」があります。これには、超伝導状態における準粒子の励起にCooper対の凝縮により定まるエネルギーギャップが必要であることが関わります。このエネルギーギャップにより、電子の流れが変化するような散乱のエネルギーが高くなり、散乱が起きなくなることで抵抗がゼロになります。

参考文献

[2]
村上雅人, 超伝導のピン留めって何? (「間違いだらけの物理概念」P39-54), 丸善株式会社, 1993
[3]
近藤慶一, ゲージ場の量子論入門 −質量ギャップとクォーク閉じ込めの解決に向けて−, SGCライブラリー45, サイエンス社, 2005, 188-190
[4]
Năstase, Horaţiu, Classical Field Theory, Cambridge University Press, 2019, 220-222
[5]
Shifman, Mikhail, Advanced Topics in Quantum Field Theory − Lecture Note −, Cambridge University Press, 2022, 93-99
[6]
太田信義, 超弦理論・ブレイン・M理論, Springer-Verlag Tokyo Modern Theoretical Physics Series, 丸善出版, 2002, 144-158
投稿日:2023810
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

bisaitama
bisaitama
143
68396

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. はじめに
  2. 磁束の侵入・常伝導スポットとピン留め効果
  3. 本記事の内容:Nielsen-Olesen vortexの導出
  4. Abelian-Higgs模型
  5. Abelian-Higgs模型におけるvortex
  6. 運動方程式と軸対称なAnsatz
  7. 境界条件
  8. 磁束の量子化
  9. BPS limit
  10. 数値計算と解に関する考察
  11. そのほかのこと
  12. まとめ
  13. 参考文献