相関 を考慮すること

相関を考慮すること (1) 相関がない時 ( バンド理論たとえば密度汎関数理論など ) の 一電子エネルギーおよび波動関数 相関がない時のスペクトル ピークの位置 非摂動グリーン関数

相関を考慮すること (2) 相関を考慮すること ＝ グリーン関数論の文脈では「自己エネルギー補正」 を考慮すること 自己エネルギー

相関を考慮すること (3) 実部 FF FF | 虚部 | 自己エネルギーは基本的にこんな感じ 複素数！ 注意！マイナス∵傾き負

相関を考慮すること (4) ピークの位置

相関を考慮すること (5) 占有バンドの場合 FF 非占有バンドの場合 FF FF 自己エネルギー補正を 考慮するとピークが E F 近傍に集まってくる バンドが縮む ピークの位置 : 注意！分母 >1

相関を考慮すること (6) 実部 | 虚部 | もともとのバンド幅 こういう場合どうなるか ? この場合は直接計算するしかない ( 強相関 )

相関を考慮すること (7) 自己エネルギーの虚部 = “ 電子の散乱 ” “ 素励起 ” のエネルギースケール “ 素励起 ” = フォノン、ポーラロン、ポラリトン、プラズモン、 マグノン、ソリトン、スピノン、ホロン、ダブロン, フェイゾン, etc | 虚部 |

相関を考慮すること (8) ハバードモデル 1963 年にジョン・ハバードによって提出された、電子相関の効果の強い固体 中の電子の振る舞いを量子論的に記述するモデルである。 元々は、遷移金属 の様に最外殻電子が d 軌道や f 軌道にあり、電子の波動関数の広がりが大きく、 電子同士の波動関数の重なりのために生じる電子相関が大きな固体中の電子 を良く記述するモデルとして提出されたものである。 ハバードモデルは非常 に単純なハミルトニアンを持つモデルであるにも関わらず、非常に多様な電 子の振る舞いを表現できる。 この様な電子の振る舞いの多様さは電子同士の 相互作用 ( 電子相関 ) によってもたらされていると考えられている。 局所 ( オンサイト ) クーロン相互作用

“ 素励起 ” = フォノン、ポーラロン、ポラリトン、プラズモン、 マグノン、ソリトン、スピノン、ホロン、ダブロン, フェイゾン, etc “ 素励起 ” のエネルギースケール | 虚部 | 相関を考慮すること (9) 自己エネルギーの虚部 = “ 電子の散乱 ”

プラズマ振動 ある場所の電子集団が局所的に動くとそこで電気的中性が破れて電荷密度を 生じ、電子を引き戻す方向に電場を生ずる. イオンは電子より質量がはるかに 大きい ので、電場によって加速されるのは電子だけである. こうしてその電場により 電子群が動いて、電気的中性を取り戻す. しかし、電子には慣性があるので、 中性 を取り戻した時点では止まらず、逆の方向に行き過ぎる. そこでまた中性が 破れて電場が生じまた電子群が引き戻される。かくして電子群の往復運動、 すなわ ち振動が起こる. これは巨視的には電荷密度の波動となる。これが プラズマ振動である 長距離クーロン相互作用 相関を考慮すること (10)

低エネルギープラズマロン状態の GW 解析 中村和磨 ( 九工大 ) 共同研究者 : 酒井志朗 ( 東大 ), 吉本芳英 ( 鳥取大 ), 野原善郎 (MPI), 野村悠祐 ( 東大 ), 有田亮太郎 ( 東大 ), 黒木和彦 ( 阪大 ) PHYSICAL REVIEW B 88, (2013) 有機導体 TMTSF 遷移金属酸化物 (SrVO 3 )

研究目的 : 強相関物質の特徴 ・ 狭バンド幅 ・ バンド幅に匹敵する局所（オンサイト） クーロン相互作用 ・ 相互作用による自己エネルギー補正大 有機導体 : EtMe 3 Sb[Pd(dmit) 2 ] 2 W = 0.44 eV U = 0.61 eV V = 0.20 eV PRB 86, (2012) ・ 電子状態変化 ( 金属絶縁体転移、 高温超伝導, etc)

SrVO 3 K 3 C 60 Na 封入ゼオライト W = 0.85eV U = 2.71eV V = 0.61eV W = 0.62 eV U = 0.82 eV U’= 0.76 eV J = eV V = 0.25 eV 研究目的 : W = 2.55 eV U = 3.48eV U’= 2.37eV J = 0.51eV V = 0.79eV PRB 80, (2009) PRB 85, (2012)JPSJ 77, (2008) 孤立バンド性

物質 ω ｐ l (bare) W U cRPA V cRPA U-V 実験 SrVO 金属 Na-SOD AF 絶縁体 K 3 C 金属・超伝導 TMTSF 金属・ SDW Al 金属 1.25 (q x →0) 0.41 (q y →0) unit: eV 研究目的： Plasmon excitation

SrVO 3 反射率 EELS ○ 理論 Expt. 1.4eV Theory 1.9eV 反射率・ EELS はプラズモン励起を観測する実験手段 Energy (eV)

K 3 C 60 ○ 理論 q=(2π/5a)(1,0,0) Energy (eV) 反射率 反射率・ EELS はプラズモン励起を観測する実験手段

y x (TMTSF) 2 PF 6 E//x E//y ○ 理論 × 実験 1000 cm -1 = 0.12 eV 反射率 反射率・ EELS はプラズモン励起を観測する実験手段

孤立狭バンド系では、 バンド幅 ～ プラズモン励起エネルギー ～ 相互作用 であり、エネルギー的に拮抗状態にある プラズモン励起は実験でも実際に観測されている この問題を第一原理計算の観点から検討したい 研究目的 ゆえに、この励起は低エネルギー物性に本質的 (“relevant” という ) な影響を与えうる

プラズモン励起による自己エネルギー G Σ W GW 自己エネルギー

GW calculation detail Fock exchange RPA correlation

do q = 1, N k ( MPI) do n = 1, N i (MPI) do  = 1, N o do k = 1, N k (omp) call FFT enddo enddo n enndo q do k = 1, N k do  = 1, N o do  = 1, N o do ω = 1, N ω enddo w do ωj = 1, Nω enddo w enddo a enddo b enddo k do G=1, N PW do G’=1, N PW (omp) Enddo enddo GW do q = 1, N k ( MPI ) do  = 1, N pair ( MPI+om p) do k = 1, N k call FFT enddo call TETRAHEDRON do G=1, N PW do G’=1, N PW do  = 1,N  do k=1, N k Enddo enddo enndo enddo RPA 並列化

京速コンピュータ

y x (TMTSF) 2 PF 6 PRB 88, (2013)

y x 反射率： TMTSF E//x E//y ○ 理論 × 実験 1000 cm -1 = 0.12 eV 1eV

Spectral function DFT PRB 88, (2013) Red: Kohn-Sham band GW

Origin of poles B. I. Lundqvist, Phys. Kondens. Mater. 6, 193 (1967). D. C. Langreth, Phys. Rev. B 1, 471 (1970).

因果関係 ImΣ(k,w) A(k,w) Red: Kohn-Sham band

Plasma frequency 1000 cm -1 = 0.12 eV ω p ～ (( )/2)*0.12 ～ 0.8eV この「 0.8eV 」というエネルギースケールの「外側」で スペクトル関数が明るくなっている A(k,w)

グリーン関数の精度 δ=0.01eVδ=0.10eV Y  X M  Red: Kohn-Sham band

まとめ : TMTSF 有機化合物 (TMTSF) 2 PF 6 の反射率を計算したところ、 電場偏光を a 軸に取った場合のプラズマ周波数は 1 eV 、 ab 面内でかつ a 軸に垂直な場合は 0.2 eV であり、実験 を再現した. 「低エネルギープラズモン励起」の電子構造への効果を GW 計算コードを用いて調べたところ、以下が確認され た : (i) 占有 / 非占有バンドの約 0.5eV 下 / 上に、プラズモン励 起に由来する新たな状態が出現する ( プラズマロン状 態 ). (ii)X-M 線に沿った電子占有領域において特に大きな電 子散乱が生じている. (iii) 温度効果（グリーン関数に導入されるボケ ) が顕著で ある. 他の物性 ( 磁性・超伝導など ) への影響・競合効果の検討 が 将来課題.

SrVO 3

反射率 EELS ○ 理論 Expt. 1.4eV Theory 1.9eV

Energy (eV) FF 電子構造 : ARPES: SrVO3 Sr O V 1.4eV 付近に状態がある

ARPES and LDA+DMFT: SrVO3 LDA+DMFT Expt LDAonly Lower Hubbard band Nekrasov, PRB 72, (2005) (U=5.5eV)

ω (eV) GW+DMFT: SrVO3 Sakuma-Werner-Aryasetiawan, PRB 88, (2013) (U=5.5eV) ReΣ(k=Γ,ω) ImΣ(k=Γ,ω)

 =0.01eV  =0.05eV  =0.10eV  =0.20eV ReΣ(k=Γ,ω)ImΣ(k=Γ,ω) 自分の計算 : SrVO3

 dependence of the ImΣ DOS  =0.01eV  =0.05eV  =0.10eV  =0.20eV ImΣ(ω) EFEF

A(kw): δ=0.20eV, sht=2.20eV (EF=8.14eV) KS (EF=8.13eV)

A(kw): δ=0.10eV, sht=2.45eV (EF=8.13eV) KS (EF=8.13eV)

A(kw): δ=0.05eV, sht=2.55eV (EF=8.14eV) KS (EF=8.13eV)

A(kw): δ=0.01eV, sht=2.64eV (EF=8.13eV) KS (EF=8.13eV)

0.2eV0.1eV 0.05eV 0.01eV  dependence of A(kw)

 dependence of the GW band KS  =0.01eV  =0.05eV  =0.10eV  =0.20eV

 dependence of the GW DOS KS  =0.01eV  =0.05eV  =0.10eV  =0.20eV

Photoemission & inverse photoemission: SrVO 3 PhotoemissionInverse photoemission  =0.01eV  =0.05eV  =0.10eV  =0.20eV

Photoemission: SrVO 3

まとめ : SrVO3 「低エネルギープラズモン励起」の電子構造 への効果を GW 計算コードを用いて調べた (i) プラズマロン状態 (ii)δ 依存性に鋭敏 ( 準粒子バンド幅 ) (iii) 強相関効果との競合

ImΣ(kw): dlt=0.20eV, sht=2.20eV (EF=8.14eV) KS (EF=8.13eV)

ImΣ(kw): dlt=0.10eV, sht=2.45eV (EF=8.13eV) KS (EF=8.13eV)

ImΣ(kw): dlt=0.05eV, sht=2.55eV (EF=8.14eV) KS (EF=8.13eV)

ImΣ(kw): dlt=0.01eV, sht=2.64eV (EF=8.13eV) KS (EF=8.13eV)

Progress report : SrVO3, 11x11x11, n=1/24 Kazuma Nakamura 計算条件： - SrVO3, 11x11x11, GGA, - 1/24 filling - Ecut = 10 Ry for χ, Ecut = 49 Ry for ψ, - Nb=50 (docc:12, pocc=t 2g =3, uocc=35), - χ grid ( 点 ), - Σ grid -22 ～ 38 eV (6000 点 ), - δ=0.01eV

n=1/6 (t2g=1 electron)n=1/24 (t2g=0.25 electron) Band structure

n=1/6 (t2g=1 electron)n=1/24 (t2g=0.25 electron) Inverse of dielectric function q=(2π/11a)(1,0,0)

EF n=1/6 n=1/24 ImΣ DOS

GW-DOS sht=0.0eV (EF=7.86eV) KS (EF=6.846eV)

GW-DOS sht=1.0eV (EF=6.80eV) KS (EF=6.846eV)

GW-DOS sht=2.0eV (EF=3.04eV) KS (EF=6.846eV)

GW-DOS sht=3.0eV (EF=1.45eV) KS (EF=6.846eV)

GW-DOS sht=4.0eV (EF=0.33eV) KS (EF=6.846eV)

フェルミ準位、どうもうまく決まらない。 しょうがないから、シフト =+3 eV について A(kω) と ImΣ(kω) を描いてみた。

A(kw): sht=3eV (EF=1.45eV) KS (EF=6.846eV)

ImΣ(kw): sht=3eV (EF=1.45eV) KS (EF=6.846eV)

Progress report : SrVO3, 11x11x11, n=1/3 Kazuma Nakamura 計算条件： - SrVO3, 11x11x11, GGA, - 1/3 filling - Ecut = 10 Ry for χ, Ecut = 49 Ry for ψ, - Nb=50 (docc:12, pocc=t 2g =3, uocc=35), - χ grid ( 点 ), - Σ grid -22 ～ 38 eV (6000 点 ), - δ=0.01eV

n=1/6 (t2g=1 electron)n=1/3 (t2g=2 electron) Band structure

n=1/6 (t2g=1 electron)n=1/3 (t2g=2 electron) Inverse of dielectric function q=(2π/11a)(1,0,0)

EF n=1/6 n=1/24 ImΣ DOS EF n=1/3

GW-DOS sht=0.0eV (EF=11.34eV) KS (EF=9.55eV)

GW-DOS sht=1.0eV (EF=10.76eV) KS (EF=9.55eV)

GW-DOS sht=2.0eV (EF=10.17eV) KS (EF=9.55eV)

GW-DOS sht=2.98eV (EF=9.55eV) KS (EF=9.55eV)

GW-DOS sht=3.0eV (EF=9.54eV) KS (EF=9.55eV)

GW-DOS sht=4.0eV (EF=8.85eV) KS (EF=9.55eV)

A(kw): sht=2.98eV (EF=9.55eV) KS (EF=9.55eV)

ImΣ(kw): sht=2.98eV (EF=9.55eV) KS (EF=9.55eV)