Web会議ツールZoomを利用したオンライン物性理論セミナーです。

登録する際のメールアドレスは、できるだけ大学もしくは研究機関のものをご使用ください。

• ZoomのミーティングURLおよびパスワードは、登録されたメールアドレスに送信されます。転載・転送は控えてください。
• 参加時の表示名は「登録時の名前＠登録した機関名」に設定してください。
• 質問する際は、ハウリング防止のため、ヘッドフォン着用をお願いします。

世話人:　高三和晃（東大）、下川統久朗（OIST）、永井佑紀（原子力機構）、品岡寛（埼玉大学)

• 第108回: 副島 智大 (ハーバード大)「対称性の破れから生まれるトポロジー:異常量子ホール結晶の最近の展開」
• 第107回: 岩崎 舜平 (慶應義塾大学)「バルクな奇周波数超伝導のペア特性と核スピン格子緩和率1/T_1　」
• 第106回: Mateo Fontaine (Keio Univ.)「Phase diagram of a spin-1/2 XXZ ladder with chirality-chirality interaction」
• 第105回: 髙橋 秀顕 (トリノ大学) 「テンソルトレイン形式の階層型運動方程式：散逸ホルスタイン模型による有機半導体中の電荷移動の研究」
• 第104回: 𠮷井 真央 (東大物工) 「トポロジカル輸送と準周期構造の関係」
• 第103回: 玉谷 知裕 (東大物性研) 「高次高調波発生を例にした非線形光学の摂動論から非摂動論への拡張」
• 第102回: 宮﨑 優希 (青山学院大学) 「量子相転移における密度行列の置換不変部分に関するレニーエントロピーの振る舞い」
• 第101回: 杉本 高大 (慶應義塾大学) 「DC電場による一次元モット絶縁体の励起スペクトルの離散化」
• 第100回: 福井 毅勇 (東大) 「模型の拡張によるキタエフスピン液体の実現可能性の探索: 汎関数繰り込み群による研究」
• 第99回: 金子 隆威 (上智大学) 「量子多体系の実時間発展のテンソルネットワーク法による計算と機械学習手法による予測」
• 第98回: 鈴木 遼太郎 (ベルリン自由大学) 「観測ランダム回路における量子複雑性相転移」
• 第97回: 山岸 愛 (東大理) 「量子ウォークの高次元化と量子アクティブ粒子の提案」
• 第96回: 深井 康平 (東大物性研) 「一次元Hubbard模型における局所保存量の表式」
• 第95回: Rico Pohle (Keio Univ.) 「Ground state of the S=1/2 pyrochlore Heisenberg antiferromagnet:A quantum spin liquid emergent from dimensional reduction」
• 第94回: 金賀 穂 (千葉大) 「電子系・スピン系における非線形光学応答の制御法についての理論解析」
• 第93回: 大熊 信之 (九工大) 「2粒子トポロジーを用いた分数チャーン絶縁体の特徴付け
• 第92回: 鏡原 大地 (近大) 「理想Bose粒子系および強く相互作用するBose粒子系におけるRényiエンタングルメントエントロピーのクエンチダイナミクス」
• 第91回: 新城 一矢 (理研) 「光励起したモット絶縁体の吸収スペクトル」
• 第90回: 大門 俊介 (QST) 「深層学習AIによる量子伝導現象の解読」
• 第89回: 水田 郁 (東大) 「時間周期系/準周期系の最適な量⼦計算」
• 第88回: 越智 正之 (阪大) 「波動関数理論に基づく固体の第一原理計算手法の開発」
• 第87回: 幸城 秀彦 (理研) 「多段磁化プラトーを実現するフラストレートスピンラダーの理論的研究」
• 第86回: 小沢 耀弘 (東大) 「強磁性積層カゴメ系ワイル半金属Co3Sn2S2の有効模型を用いた磁気秩序とスピン電荷輸送の理論解析」
• 第85回: 佐藤 駿丞 (筑波大) 「第一原理電子ダイナミクス計算とアト秒光科学への応用」
• 第84回: 仲澤 一輝（東大) 「量子輸送現象によって検出する物質中のトポロジカル構造」
• 第83回: 中西 優馬（東工大) 「PT対称性由来の散逸時間結晶 」
• 第82回: 大岩 陸人（明治大) 「対称性適合多極子基底を用いた有効模型生成」
• 第81回: 野垣 康介（京大) 「局所空間反転対称性の破れが誘起する新奇超伝導相におけるトポロジーと強相関効果」
• 第80回: Kimberly Remund（OIST) 「Spin-1 Magnet - a u(3) Formalism」
• 第79回: 川野 雅敬（ミュンヘン工科大学) 「サイン二乗変形を用いた長周期量子相の効率的記述」
• 第78回: 佐藤 年裕（ヴュルツブルク大学) 「フラストレート量子スピン系への補助場量子モンテカルロアプローチ」
• 第77回: 足立 景亮（理研) 「斥力相互作用するアクティブマターモデルの非平衡協同現象」
• 第76回: 桂 ゆかり（東大・NIMS) 「論文からの大規模実験データの収集と理論物性科学への活用」
• 第75回: 三澤 貴宏（BAQIS) 「トポロジカル絶縁体におけるグリーン関数のゼロ -トポロジカル状態の簡便な可視化方法」
• 第74回: 乾 幸地（東大) 「自動微分を用いた目的の性質を満たすハミルトニアンの逆設計」
• 第73回: 竹森 那由多（阪大) 「動的平均場近似を用いた超伝導ハイパーマテリアルの解析」
• 第72回: 荻野 卓啓（東大物性研) 「スピン1/2XXZ梯子系における量子相転移と相構造の解明」
• 第71回: 上田 篤（東大) 「テンソルネットワークによる二次元古典模型の繰り込みフローの視覚化」
• 第70回: 森田 克洋（東京理科大) 「フラストレート量子スピン系の有限温度特性と磁気熱量効果」
• 第69回: 大同 暁人（京大) 「内因的超伝導ダイオード効果」
• 第68回: 中 惇（東京電機大学)「有機反強磁性体のスピン・電荷輸送」
• 第67回: 曽根 和樹（東大）「非線形振動子のトポロジカルな同期現象」
• 第66回: 増木 亮太（東大）「非調和フォノン理論に基づく有限温度における構造最適化」
• 第65回: 安田 憲司（MIT）「積層制御による二次元強誘電体の設計」
• 第64回: 山田 昌彦（阪大）「ハニカム格子上のトポロジカルZ_4スピン液体」
• 第63回: 南谷 英美（分子研）「第一原理計算によるエネルギー損失過程の解析」
• 第62回: 幸城 秀彦（東大）「多重磁化プラトーを実現するフラストレートスピンラダーの理論的解明」(延期)
• 第61回: 吉永 敦紀（東大）「弱横磁場下のイジング模型におけるHilbert Space Fragmentation」
• 第60回: 播木 敦（大阪府立大）「内殻X線分光を用いた超伝導層状ニッケル酸化物の電子状態の理論研究」
• 第59回: 水野 竜太（京大）「DMFTの不純物問題におけるフルバーテックスの簡略化と非局所相関の取り込みのためのその応用」
• 第58回: 田財 里奈（名古屋大）「汎関数くりこみ群/密度波方程式による非局所相転移の理論」
• 第57回: 近藤 寛記（東大）「対称性に保護されたマグノンのディラック表面状態と反強磁性体におけるマグノンの非線形スピンネルンスト効果」
• 第56回: 高橋 雅大（阪大）「キタエフ磁性体におけるトポロジカルネマティック相転移と点欠陥に束縛されたマヨラナゼロモードの確立」
• 第55回: 伏屋 雄紀（電通大）「ナノスケールのチューリング・パターン：単層ビスマスの形態形成」
• 第54回: 秦 徳郎（東工大）「量子液体における三体相関の実験的検出」
• 第53回: 吉川 尚孝（東大）「電荷密度波の振幅モード励起による絶縁体的状態」
• 第52回: 竹内 宏光（大阪市立大）「量子楕円渦：U(1)対称性が破れた超流動相における巻き数１の量子渦の自発的“分裂”」
• 第51回: 中村 真（中央大）「ゲージ・重力対応を用いた非平衡相転移の解析と電流駆動型非平衡三重臨界点の発見」
• 第50回: 吉岡 信行（東大）「ニューラルネットワークによる固体量子物性の第一原理計算」
• 第49回: 金杉 翔太（京大）「多バンド系におけるアナポール超伝導の理論」(第42回の延期分)
• 第48回: 水上 雄太（東大）「キタエフ量子スピン液体におけるマヨラナ励起の磁場角度依存性」
• 第47回: 那須 譲治（横浜国立大）「キタエフ量子スピン系の熱輸送特性とその磁場角度依存性」
• 第46回: 関 和弘（理研）「強相関格子模型に対する変分量子計算法の古典計算機シミュレーション」
• 第45回: 北折 暁（東大）「創発インダクター：材料開拓とこれからの課題」
• 第44回: 星野 晋太郎（埼玉大）「ボゴリウボフ準粒子に対する相互作用の効果」
• 第43回: 中川 大也（東大）「散逸Hubbard模型における磁性と超流動ペアリング」
• 第42回: 金杉 翔太（京大）「多バンド系におけるアナポール超伝導の理論」(延期)
• 第41回: 國見 昌哉（分子研）「冷却原子系における準古典近似を用いた非平衡ダイナミクスの研究」
• 第40回: 篠原 康「固体からの高次高調波発生の時間依存スペクトルの半古典ダイナミクスによる記述と、時間依存Hartree-Fockへの拡張」
• 第39回: 内藤 智也「物性理論は原子核理論の夢を見るか? ―原子核理論と物性理論のクロスオーバー」
• 第38回: 池田 暁彦「超強磁場に誘起されるLaCoO3の励起子凝縮 〜ひずみ計測による研究〜」
• 第37回: 後藤 慎平「冷却原子気体の観測誘起転移」
• 第36回: 石塚 淳「ウラン系超伝導体UTe2の最近の進展と第一原理計算によるトポロジカル超伝導予測」
• 第35回: 金子 竜也「η-pairing状態の光誘起理論」
• 第34回: 石塚 大晃「テラヘルツ光によるスピン流・軌道流の整流効果」
• 第33回: 古谷 峻介「量子スピン系の磁化プラトーとアノマリー」
• 第32回: 関孝一「ランダム量子スピン系での実空間数値くりこみ群法」
• 第31回: 工藤耕司「量子ホール効果の断熱接続における一般化ストレーダ公式」
• 第30回: 望月健「フロケ系や開放系特有のトポロジカルな現象：量子ウォークを舞台として」
• 第29回: 永井瞭「機械学習手法を用いた交換相関汎関数の構築とその展望」
• 第28回: 高橋惇「脱閉じ込め量子臨界現象に対する“de-signer”ハミルトニアンによるアプローチ」
• 第27回: 磯部大樹「ファンホーブ特異性のもたらす電子相関物性：超伝導・絶縁体状態および非フェルミ流体」
• 第26回: 岩﨑祐昂「バンドエンジニアリングによる近似結晶半導体の探索と熱電材料への応用」
• 第25回: 山本大輔「磁場中SU(3) Heisenberg模型の量子相と隠れたネマティック性を持つスピンダイポール秩序」
• 第24回: 兎子尾理貴「空間反転対称性が破れた高伝導金属における流体力学的異常輸送現象」
• 第23回: 村上雄太「励起子絶縁相の非線形応答と非平衡誘起」
• 第22回: 井戸康太「励起状態のための変分モンテカルロ法」
• 第21回: 小野淳「反強磁性ディラック半金属におけるネールベクトルの高速光制御」
• 第20回: 永井恒平「高強度光駆動下における固体電⼦系の動的対称性とその実験的検証」
• 第19回: 田村駿「カイラル対称性を持つ系における奇周波数クーパー対の理論」
• 第18回: 藤原理賀「カムチャツカ半島産低次元磁性体における量子スピン状態」
• 第17回: 小野清志郎「対称性指標に基づくトポロジカル超伝導相の分類と物質探索への応用」
• 第16回: 野村悠祐「機械学習手法を用いた量子多体系の研究 ー手法拡張を含む最近の進展」
• 第15回: 渡邉光「磁気的なパリティの破れに誘起される円偏光ガルバニック効果」
• 第14回: 原口祐哉「新規ハニカム格子化合物における量子物性の探索」
• 第13回: 青山和司「磁性絶縁体のスピン伝導とトポロジカル転移」
• 第12回: 岩木惇司「熱的量子純粋行列積状態」※曜日がいつもと違います
• 第11回: 野本拓也「第一原理計算に基づくGd系スカーミオン物質の研究」
• 第10回: 松下太樹「線ノード半金属におけるトポロジカルな圧電効果」
• 第9回: 藤代有絵子「高密度なトポロジカル磁気構造における巨大創発磁場と量子輸送現象」
• 第8回: 明石遼介「第一原理計算に基づく高圧下における水素化物超伝導体の研究」
• 第7回: 道下佳寛「強相関電子系における量子開放系としての性質及び非エルミート性」
• 第6回: 藤本純治「動的ねじれによるマグノン流生成」
• 第5回: 山田昌彦「アンダーソン・キタエフスピン液体」
• 第4回: 北谷基治「動的バーテックス近似(DΓA)を用いた非従来型超伝導の研究」
• 第3回: 奥村駿「空間反転対称性の破れた金属における磁気ヘッジホッグ格子の安定性とトポロジカル転移」
• 第2回: 植木輝「第二種超伝導体の渦状態におけるホール効果」
• 第1回: 光元亨汰「凍結した乱れのない系におけるスピングラス転移」

第108回

日時: 2024年11月26日(火)11:00-
発表者: 副島 智大 (ハーバード大)
発表題目: 対称性の破れから生まれるトポロジー:異常量子ホール結晶の最近の展開
Registration form

近年、グラフェンなどの二次元物質を組み合わせることによってモアレ模様を作り出し、それによって物性を制御するという試みが大きな成功を収めてきた。これらの物質群のことをモアレ物質と呼び、超伝導に代表されるような新奇物性が多く実現されている。特に、チャーン絶縁体に代表されるようなトポロジカルに非自明な物質の実現は多くの注目を集めている。 今回のトークでは、その中でも最近特に注目を集めているRhombohedral pentalayer grapheneの物理を紹介する。特に、この系で近年実験的に観測された整数チャーン絶縁体が、対称性の破れによってトポロジーが生まれた「異常量子ホール結晶」として理解できることを議論する。

第107回

日時: 2024年11月12日(火)11:00-
発表者: 岩崎 舜平 (慶應義塾大学)
発表題目: バルクな奇周波数超伝導のペア特性と核スピン格子緩和率1/T_1
Registration form

超伝導は、2つの電子がクーパーペアと呼ばれる分子を形成し、この”分子ボソン”がある種のボース凝縮を起こした状態と解釈出来る。通常の超伝導は同時刻の電子同士のペア形成が可能な偶周波数超伝導であるのに対し、本研究が対象とする奇周波数超伝導は、異なる時刻にある電子の間でしかペア形成が生じないという著しい特徴を有している。このような異なる時刻の2電子が組むペアは、素朴な分子の描像とは大きく異なるため、“分子ボソン”として振る舞うのか、或いは、低温においてボース凝縮を起こすのかどうか明らかではない。

また、近年、重い電子系反強磁性超伝導体CeRh_{0.5}Ir_{0.5}In_{5}の圧力誘起磁気量子臨界点近傍における核スピン-格子緩和率1/T_1の異常な温度依存性が観測されており[1]、これに基づき、奇周波数p波超伝導状態の実現が指摘されている[1]。しかし、従来のバルクな奇周波数超伝導の理論は、1粒子励起スペクトルが負になるという問題を抱えており、実験と定量的に比較することが困難であった[2]。そのため、奇周波数p波状態が、1/T_1の振る舞いを本当に説明できるかどうかは明らかではなかった。

そこで、本講演ではバルクな奇周波数超伝導におけるクーパーペアの性質[2]、及び、その核スピン格子緩和率1/T_1に対する効果[3]について理論的に議論する。

まず、我々が考案した1粒子励起スペクトルの正定値性を満たす奇周波数超伝導のモデル[2]に対して、近年提案された経路積分形式の理論[4-6]とBCS-Eagless-Leggett理論[7,8]を組み合わせて適用し、絶対零度T=0における奇周波数ペア波動関数の時空間構造を計算する。加えて、NozieresとSchmitt-Rink により開発された強結合理論[9]で超伝導相転移温度T_cを評価する。これにより、強結合領域（従来の偶周波数s波の場合では理想Bose分子気体によって良く記述される領域）において、奇周波数ペアが分子ボソン的に振る舞うかどうかを明らかにする[2]。

さらに、奇周波数p波超伝導状態の1/T_1の温度依存性を、平均場近似[4-6]の範囲で明示的に計算し、CeRh_{0.5}Ir_{0.5}In_{5}において観測された1/T_1の温度依存性をどの程度説明できるかを検討する。これに加えて、奇周波数s波の場合についても、従来の偶周波数s波BCS状態でT_c直下に現れるコヒーレンスピークが、奇周波数Cooper対形成によってどのように影響を受けるかを解明する[3]。

[1] S. Kawasaki, T. Oka, A. Sorime, Y. Kogame, K. Uemoto, K. Matano, J. Guo, S. Cai, L. Sun, J. L. Sarrao, J. D. Thompson and G.-q. Zheng, Commun. Phys. 3, 148 (2020).
[2] S. Iwasaki, T. Kawamura, K. Manabe, and Y.Ohashi, Phys. Rev. A. 109, 063309, (2024).
[3] S. Iwasaki, and Y. Ohashi, arXiv:2409.10894 (2024).
[4] D. Solenov, I. Martin, and D. Mozyrsky, Phys. Rev. B 79, 132502 (2009).
[5] H. Kusunose, Y. Fuseya, and K. Miyake, J. Phys. Soc. Jpn. 80 (2011) 044711.
[6] Y. V. Fominov, Y. Tanaka, Y. Asano, and M. Eschrig, Phys. Rev. B 91, 144514 (2015).
[7] D. M. Eagles, Phys. Rev. 186, 456, (1969).
[8] A. J. Leggett, in Modern Trends in the Theory of Condensed Matter, edited by A. Pekalski and J. Przystawa (Springer-Verlag, Berlin, 1980), pp. 13-27.
[9] P. Nozieres and S. Schmitt-Rink, J. Low Temp. Phys. 59, 195 (1985).

第106回

日時: 2024年10月29日(火)11:00- This talk will be given in English
発表者: Mateo Fontaine (Keio Univ.)
発表題目: Phase diagram of a spin-1/2 XXZ ladder with chirality-chirality interaction
Registration form

We investigate the ground-state phase diagram of a spin-1/2 XXZ model with chirality-chirality interaction (CCI) on a two-leg ladder, considering both isotropic [1] and weakly anisotropic [2] cases. This model provides a minimal setup to explore the interplay between spin and chirality degrees of freedom and is potentially relevant for understanding the effects of four-spin ring exchange [3,4] and phonon-mediated interactions on a ladder.

Spin-chirality duality transformations [4] allow us to relate the regimes of weak and strong CCIs. By applying Abelian bosonization as in Refs. [5,6] in combination with two types of spin-chirality dualities, we derive a comprehensive phase diagram. Numerical simulations are also performed to confirm the predicted phase structure and critical properties.

In the isotropic (Heisenberg) model, the phase diagram consists of the rung singlet (RS) phase, the Haldane phase, and dimer phases, and a scalar chiral phase [1]. We find that the RS and Haldane phases consist of multiple regions with different dominant correlations in terms of spin, dimer, or chirality. In particular, the dominance of staggered and uniform vector chiral correlations in the RS phase for strong positive and negative CCIs can be understood as dual counterparts to the conventional Néel correlation.

The introduction of weak XXZ anisotropy further enriches the phase structure within the region of positive CCI and antiferromagnetic rung interaction [2]. Specifically, Néel and vector chiral orders emerge for easy-axis anisotropy, while two distinct symmetry-protected topological (SPT) phases emerge for easy-plane anisotropy. These two SPT phases can be viewed as twisted variants of the Haldane phase [6,7]. We demonstrate that the two SPT phases, along with a trivial phase, can be distinguished by topological indices in the presence of certain symmetries [8].

[1] M. Fontaine and S. Furukawa, arXiv:2408.11731 (2024).
[2] M. Fontaine, K. Sugimoto, and S. Furukawa, Phys. Rev. B 109, 134413 (2024).
[3] A. Läuchli, G.Schmid, and M. Troyer, Phys. Rev. B 67, 100409 (2003).
[4] T. Hikihara, T. Momoi, and X. Hu, Phys. Rev. Lett. 90, 087204 (2003).
[5] S. Takayoshi and M. Sato, Phys. Rev. B 82, 214420 (2010).
[6] T. Ogino, S. Furukawa, R. Kaneko, S. Morita, and N. Kawashima, Phys. Rev. B 104, 075135 (2021).
[7] Z.-X. Liu, Z.-B. Yang, Y.-J. Han, W. Yi, and X.-G.Wen, Phys. Rev. B 86, 195122 (2012).
[8] F. Pollmann and A. M. Turner, Phys. Rev. B 86, 125441 (2012).

第105回

日時: 2024年10月8日(火)16:30- 通常と開催時刻が異なります
発表者: 髙橋 秀顕 (トリノ大学)
発表題目: テンソルトレイン形式の階層型運動方程式：散逸ホルスタイン模型による有機半導体中の電荷移動の研究
Registration form

高移動度をもつ有機半導体（OSC）の開発は、有機太陽電池や有機電界効果トランジスタの性能向上において重要な課題である。OSCにおける電荷移動度の予測には近似理論が広く用いられている[1]。しかし、バンド伝導からホッピング伝導への移行や、輸送中の電荷の過渡的局在化などの電荷輸送機構の解明において、分子振動との相互作用の効果を正確に扱える理論が必要不可欠である。また、分子間のファンデルワールス相互作用に起因する低周波数の分子間振動が、過渡的局在化において重要な役割を果たすという機構が提案されているが[2]、局所的および非局所的な相互作用の相対的な寄与については、依然として議論が残されている。

このような課題に対処するためにいくつかの手法が開発されているが[3,4]、本研究[9]ではテンソルトレイン（TT）形式の階層型運動方程式 (HEOM-TT)[5,6]を採用した。散逸ポーラロン系のような従来のHEOMでは取り扱いが困難な大規模系に対してHEOM-TTが有効であることを示し、さらにTT形式に対する異なる時間発展アルゴリズム[7,8]の比較を行った。また、モデルとしては分子内振動に起因する局所的なホルスタイン型の相互作用に焦点を当て、低周波の局所的な相互作用が電荷移動に与える影響を調べた[9]。

[1] S. Fratini, D. Mayou, and S. Ciuchi, Adv. Funct. Mater. 26, 2292–2315 (2016).
[2] S. Ciuchi and S. Fratini, Phys. Rev. B 86, 245201 (2012).
[3] M. Lian, Y. Wang, Y. Ke and Y. Zhao, J. Chem. Phys. 151, 044115 (2019).
[4] W. Li, J. Ren, and Z. Shuai, Nat Commun 12, 4260 (2021).
[5] Q. Shi, Y. Xu, Y. Yan, and M. Xu, J. Chem. Phys. 148, 174102 (2018).
[6] R. Borrelli, J. Chem. Phys. 150, 234102 (2019).
[7] C. Lubich, I. Oseledets, and B. Vandereycken, SIAM J. Numer. Anal. 53, 917–941 (2015).
[8] S. V. Dolgov, Comput. Methods Appl. Math. 19, 23–38 (2019).
[9] H. Takahashi and R. Borrelli, J. Chem. Theory Comput. 20, 7052-7064 (2024).

第104回

日時: 2024年7月16日(火)11:00-
発表者: 𠮷井　真央 (東大物工)
発表題目: トポロジカル輸送と準周期構造の関係
Registration form

周期系と乱雑系の中間状態の1つである準周期系はTwisted bilayer grapheneなどを代表とする薄膜積層系[1]や非従来型の超伝導状態[2]や異常臨界現象[3]などが見つかった準結晶系を含むことから近年注目を浴びている。これらは実験や応用などの側面でも重要なだけでなく、理論研究の側面でも様々な面白い現象が報告されている。本講演ではその一側面としてトポロジカル現象の研究を紹介する。

結晶の電子構造やBlochの定理を導入する際によく使うNearly-Free electron近似では自由電子系に周期ポテンシャルを印加すると逆格子ベクトル分だけ波数がシフトした波動関数たちが混成し、Bragg面でエネルギーギャップを開ける。この描像により周期系ではBrillouin Zoneなどの概念が導入されるが、準周期系でもこの描像を使う事が可能でありGap labelling定理という定理につながる[4,5]。これは周期系ではほぼ自明な定理だが、準周期系では準周期性を特徴づける無理数によってエネルギーギャップがトポロジカルに保護されるという非自明な定理になる。本研究ではこのGap labelling定理の導入を行い、これを利用したトポロジカル電荷輸送の研究を薄膜積層系[6]と準結晶系[7]で紹介する。特に準結晶系ではフラクタル性によって多バンドでのトポロジカル電荷輸送が生じ、さらにはこれらが繰り込み群でつながっているフラクタルな電荷輸送である事を説明する。

[1] A. Geim and I. Grigorieva, Nature 499, 419–425 (2013).
[2] K. Kamiya, T. Takeuchi, N. Kabeya, N. Wada, T. Ishimasa, A. Ochiai, K. Deguchi, K. Imura, and N. K. Sato, Nature Communications 9 (2018).
[3] K. Deguchi, S. Matsukawa, N. Sato, et al.. Nature Materials 11, 1013–1016 (2012).
[4] J. Bellissard, From Number Theory to Physics, Springer Proc. Physics 47 (1990)
[5] M. Yoshii, S. Kitamura, T. Morimoto, Physical Review B 107, 064201 (2023).
[6] M. Koshino and H. Oka, Physical Review Research 4, 013028 (2022).
[7] M. Yoshii, S. Kitamura and T. Morimoto, Physical Review B 104, 155126 (2021).

第103回

日時: 2024年7月9日(火)11:00-
発表者: 玉谷 知裕 (東大物性研)
発表題目: 高次高調波発生を例にした非線形光学の摂動論から非摂動論への拡張
Registration form

近年、高強度テラヘルツ光源の開発に伴い、半導体を用いた高強度光物理学の分野が大きな注目を集めている。実際、これらの光源を用いることで、強い光と物質の相互作用に起因する非摂動論的な非線形光学現象が観測されるようになった。現在、これらの現象がブロッホ電子のバンド内遷移（電場による加速運動）によって引き起こされているのか、それともバンド間遷移（光による励起過程）に起因するものなのかという問題を巡って様々な議論がなされているが、未だ信頼性のある解答は得られていない。

そこで本講演では、これまでに講演者が高強度光物理現象を取り扱うために構築してきた理論的枠組みについて紹介する[1-3]。本物理現象の最大の特徴は、高強度光で駆動されたブロッホ電子が常に価電子帯と伝導帯の重ね合わせ状態（量子純粋状態）にあるという点である。このブロッホ電子の量子純粋状態は、非線形光学にとって中心的な概念である電気分極を生み出すと共に、光の周波数で超高速に変化する。そのため、非摂動論的な非線形光学現象はブロッホ電子の単純な電場による加速運動だけで記述することはできず、バンド内遷移とバンド間遷移の双方の同時プロセスを考慮することで初めて理解可能になる。講演者はこの理論を非線形光学現象の一種である高次高調波発生に適用し、実験結果と比較することでその妥当性を明らかにした[4,5]。さらに本理論における考察を深化させることで、高次高調波発生と他の高強度光物理現象との関係性を明らかにし[6]、高強度領域における非線形光学現象を統一的に理解するための物理的描像を示した。本講演ではこれらの詳細な説明を行うと共に、従来の摂動論的な領域における非線形光学を非摂動論的な領域にまで拡張するために必要な考え方について述べる。

[1] T. Tamaya, A. Ishikawa, T. Ogawa, and K. Tanaka, Phys. Rev. Lett, 116, 016601 (2016).
[2] T. Tamaya, A. Ishikawa, T. Ogawa, and K. Tanaka, Phys. Rev. B 94, 241107(R) (2016).
[3] T. Tamaya, A. Akiyama and T. Kato, Phys. Rev. B 107, L081405 (2023).
[4] N. Yoshikawa, T. Tamaya, and K. Tanaka, Science 356, 736 (2017).
[5] P. Xia, T. Tamaya, et al., Phys. Rev. B 104, L121202 (2021). (＊equal contribution)
[6] T. Tamaya and T. Kato, Phys. Rev. B 100, 081203(R) (2019).

第102回

日時: 2024年6月25日(火)11:00-
発表者: 宮﨑 優希 (青山学院大学)
発表題目: 量子相転移における密度行列の置換不変部分に関するレニーエントロピーの振る舞い
Registration form

量子情報理論に基づくアプローチは、量子多体系の示す物理現象の解析において有効に用いられている。例えば、エンタングルメントエントロピーのスケーリング則は、従来の「対称性の破れ」では説明できないトポロジカル秩序相を特徴づけることが知られている[1, 2]。他にもエンタングルメントスペクトラムの構造は、磁気秩序相と量子スピン液体相の区別を可能とすることが知られている[3, 4]。このように今日の量子物性物理学では、量子情報に関係した量の理論的解析、及び実験的な測定の需要が高まっている。一方、光格子中の冷却原子気体系[5]、イオントラップ系[6]などをはじめとする人工量子系では、その高いクリーンネスと制御性を活かした量子多体系のシミュレーションが盛んに行われている。近年ではこのような人工量子系を用いて、エンタングルメントエントロピーの測定[7]だけでなく、様々な物理量の測定を繰り返し収集したデータにより密度行列そのものを再構築する、量子状態トモグラフィー（以下QST）[8]が行われている。

本研究では系の密度行列の置換不変な部分（以下PIDM）に注目し、量子多体系におけるPIDMの持つ役割を探索した。PIDMは、密度行列に対し全通りの粒子の置換を施したものを平均することで得られる。QSTの文脈では、完全な密度行列の再構築に必要な測定回数が粒子数に関して指数関数的に増加するのに対し、PIDMの構築は多項式回数で可能であることが知られている。我々は、このスケーラブルに取得可能な量が量子相転移における基底状態の変化をどのように反映するかを探索した[9]。特に横磁場イジング鎖を例に、基底状態に対してPIDMの2次のレニーエントロピーを密度行列繰り込み群法を用いて計算した。結果として、PIDMの2次のレニーエントロピーは各相によって異なるシステムサイズ依存性を示した。本講演ではこれらの結果を紹介すると共に、その振る舞いの要因について議論する。

[1] A. Kitaev and J. Preskill, Physical Review Letters 96, 110404 (2006).
[2] M. Levin and X.-G. Wen, Physical Review Letters 96, 110405 (2006).
[3] F. Pollmann et al., Physical Review B 81, 064439 (2010).
[4] A. M. Turner et al., Physical Review B 83, 075102 (2011).
[5] I. Bloch, Nature Physics 1, 23 (2005).
[6] R. Blatt and C. F. Roos, Nature Physics 8, 277 (2012).
[7] R. Islam et al., Nature 528, 77 (2015).
[8] H. Häffner et al., Nature 438, 643 (2005).
[9] Y. Miyazaki et al., arXiv:2404.08389.

第101回

日時: 2024年6月11日(火)11:00-
発表者: 杉本 高大 (慶應義塾大学)
発表題目: DC電場による一次元モット絶縁体の励起スペクトルの離散化
Registration form

電磁場による系のエネルギー準位の分裂に関する研究は、古くは磁場下でのランダウ準位や電場下でのシュタルク効果から始まり、今なお数多の議論がなされている。1960年にはワニエが傾斜ポテンシャルによって周期ポテンシャル中の波動関数が局在化すると同時にエネルギー準位が離散化すること（ワニエ・シュタルク効果）を予言し、人工的に作られた半導体超格子を使った観測も行われた[1]。近年では放射光による高強度テラヘルツ光を用いたポンプ・プローブ分光の技術が大きく発展したことで、周期ポテンシャルとして固体結晶を使ったワニエ・シュタルク効果の観測が可能となりつつあり、実際に半導体を使った研究も報告されている[2]。

自由電子系でのワニエ・シュタルク効果はよく理解されているが、それでは量子多体効果が本質的となる強相関電子系においてはワニエ・シュタルク効果がどのように現れるだろうか。本研究の目的は、モット絶縁体におけるワニエ・シュタルク効果を理解し、どのように観測されるのかを理論的に明らかにすることである。理論計算にあたっては半充填の一次元ハバード模型を取り扱い、時間発展シミュレーションにiTEBDを用いることで、モット絶縁体に静電場を印加したときの光励起スペクトルを調べた。本講演では、特に光学伝導度[3]および一粒子励起スペクトル[4]に現れるワニエ・シュタルク効果について紹介する。

[1] M. Glück, A. R. Kolovsky, and H. J. Korsch, Phys. Rep. 366, 103 (2002).
[2] C. Schmidt et al., Nat. Commun. 9, 2890 (2018).
[3] M. Udono, T. Kaneko, and K. Sugimoto, Phys. Rev. B 108, L081304 (2023).
[4] K. Sugimoto, arXiv:2401.17466.

第100回

日時: 2024年5月28日(火)11:00-
発表者: 福井 毅勇 (東大)
発表題目: 模型の拡張によるキタエフスピン液体の実現可能性の探索: 汎関数繰り込み群による研究
Registration form

キタエフ模型 [1]は、量子スピン液体状態が基底状態として厳密に求まることや、その分数励起のトポロジカル量子計算への応用可能性などから、理論と実験の双方から精力的な研究がなされてきた。その候補物質として、α-RuCl3をはじめとした多彩な物質が研究されているが、不可避的に存在するハイゼンベルグ相互作用などの非キタエフ型相互作用によって、そのほぼ全てにおいて、低温で磁気秩序が生じてしまう。したがって、非キタエフ型相互作用の効果の解明と、キタエフスピン液体を基底状態に持つ新たなプラットフォームの開拓がキタエフスピン液体実現のために重要である。しかしながら、これらの問題に理論的に取り組むためにはキタエフ模型に含まれない効果を議論する必要がある。これらは模型の可解性を損なうため、フラストレートした量子スピン系の問題となるが、数値的にも解析的にも万能の手法は知られておらず、キタエフ模型の文脈に関わらず長年の課題である。

我々は、上述の研究のための有望な手法の一つとしてPFFRG (pseudofermion functional renormalization group) 法 [2, 3]を採用し、冷却極性分子系、高スピン系、そして、３次元格子系の３通りのキタエフ量子スピン液体の実現可能性の研究を行ってきた。極性分子系については、先行研究で提案された長距離相互作用を持つ模型について調べ、相互作用の異方性に依らず基底状態が磁気秩序を持つこと、等方的な場合にスピン液体実現に最も近いことを解明した [4]。また、高スピン物質のミニマルな模型の一つとしてスピンSのキタエフ-ハイゼンベルグ模型について、スピンの長さと相互作用の比の双方を変化させながら基底状態相図を解明し、S<2の場合には有限の幅のスピン液体相が得られた [5]。最後に、3次元格子系として、候補物質も豊富であるハイパーハニカム格子上のキタエフ-ハイゼンベルグ模型を考え、別の手法による先行研究とは異なり、模型の満たすべき対称性を満たした信頼性の高い基底状態相図を解明した [6]。 セミナー当日は、手法について簡単に紹介し、次に上述の研究をそれぞれ紹介する。時間に余裕があれば、発表者の行った最近の研究 [7]について触れる。

[1] A. Kitaev, Ann. Phys. 321, 2 (2006).
[2] J. Reuther and P. Wölfle, Phys. Rev. B 81, 144410 (2010).
[3] T. Müller et al., Rep. Prog. Phys. 87, 036501 (2024).
[4] K. Fukui, K. Kato, J. Nasu, and Y. Motome, Phys. Rev. B 106, 014419 (2022).
[5] K. Fukui, K. Kato, J. Nasu, and Y. Motome, Phys. Rev. B 106, 174416 (2022).
[6] K. Fukui, K. Kato, and Y. Motome, J. Phys. Soc. Jpn. 92, 064708 (2023).
[7] K. Fukui, K. Kato, and Y. Motome, arXiv:2402.05516.

第99回

日時: 2024年5月7日(火)11:00-
発表者: 金子 隆威 (上智大学)
発表題目: 量子多体系の実時間発展のテンソルネットワーク法による計算と機械学習手法による予測
Registration form

一般に、量子多体系の実時間発展を長時間にわたって計算することは非常に困難な課題である。時間に依存したSchrödinger方程式を解いて各時刻での波動関数を愚直に求める場合は系の全固有状態の情報が必要で、Hilbert空間の次元がシステムサイズの指数関数で増大するため小さな系の時間発展しか計算できない。量子モンテカルロ法は実時間経路積分の計算で負符号問題が生じるため、一般には適用できない。エンタングルメントの面積則を満たす状態を用いたテンソルネットワーク法は、現状では空間1次元系や短時間のダイナミクスへの適用に限られる。時間依存変分モンテカルロ法は大きな系の長時間のダイナミクスへも適用できる可能性があるが、変分試行関数の最適化のコストが高い。

本講演では、量子多体系の実時間発展を2通りの方法で計算する。1つ目は、テンソルネットワーク法で各時刻の量子状態を求める方法である。空間2次元の孤立量子系のクエンチダイナミクスに着目し、ブルートフォースな方法でどこまで長時間のダイナミクスが追えるのかを検証する。ここで得られた結果が、近年注目されているアナログ量子シミュレータのクエンチダイナミクスの実験にどう役立つかも議論する[1,2]。2つ目は、物理量の正確な時系列データが既に得られている場合に、それをインプットにして機械学習手法で長時間予測する方法である。とくに、流体解析の分野で知られていた動的モード分解と呼ばれるデータ駆動型の手法を取り上げる。長時間予測が困難であると予想される量子多体系で臨界的なふるまいを示す時系列データに適用した場合でも、インプットデータの時間よりも1桁長い時間まで物理量の時間発展を正確に予測できることを事例研究を通して紹介する[3]。この手法は波動関数の情報を必要とせず、インプットデータの詳細によらず適用できるため、今後、様々な量子多体系の長時間ダイナミクスの予測に適用できると期待される。

[1] R. Kaneko and I. Danshita, Commun. Phys. 5, 65 (2022).
[2] R. Kaneko and I. Danshita, Phys. Rev. A 108, 023301 (2023).
[3] R. Kaneko, M. Imada, Y. Kabashima, T. Ohtsuki, arXiv:2403.19947.

第98回

日時: 2024年4月9日(火)16:30- 通常と開催時刻が異なります
発表者: 鈴木　遼太郎 (ベルリン自由大学)
発表題目: 観測ランダム回路における量子複雑性相転移
Registration form

ユニタリ発展と量子測定の両方から成る、量子系の観測されたダイナミクスが、主に測定誘起相転移と呼ばれるある種の相転移現象を示すことにより近年は注目を集めている[1,2]。一方で、量子状態の複雑性を定量的に求めることは、量子計算、量子多体ダイナミクス、高エネルギー物理学における現象を理解する上で重要な役割を果たすことが期待されている[3]。ここである量子状態の量子複雑性とは、その状態を直積状態から量子計算機を用いて生成するために必要な局所ユニタリ演算子の最小数である。

本発表では、観測ランダム回路において量子複雑性がどう振る舞うかを議論する。具体的に我々は、多数の量子ビットがランダムユニタリ回路に従って時間発展し、かつ各時間ステップにある頻度で個々の量子ビットが測定されるような場合における、量子複雑性のダイナミクスを考察した。そして、測定頻度を変えると量子複雑性の時間発展の様子がある臨界測定頻度を境に大きく変わることを厳密に示した。本結果は、観測ランダム回路が量子複雑性の相転移を示す興味深いトイモデルであることを示唆している。本発表は論文[4]に基づく。

[1] Brian Skinner, Jonathan Ruhman, and Adam Nahum. “Measurement-induced phase transitions in the dynamics of entanglement.” Physical Review X 9.3 (2019). [2] Yaodong Li , Xiao Chen, and Matthew PA Fisher. “Measurement-driven entanglement transition in hybrid quantum circuits.” Physical Review B 100.13 (2019). [3] Scott Aaronson. “The complexity of quantum states and transformations: from quantum money to black holes.” arXiv preprint arXiv:1607.05256 (2016). [4] Ryotaro Suzuki, Jonas Haferkamp, Jens Eisert, Philippe Faist. “Quantum complexity phase transitions in monitored random circuits.” arXiv preprint arXiv:2305.15475 (2023).

第97回

日時: 2024年2月27日(火)11:00-
発表者: 山岸 愛 (東大理)
発表題目: 量子ウォークの高次元化と量子アクティブ粒子の提案
Registration form

本講演では、(i)量子アクティブマターの最小理論モデルの提案と、(ii)離散時間量子ウォークの高次元系への拡張の二点について講演する。

前半のトピックであるアクティブマターとは、外界から取り込んだエネルギーを自身の内部に蓄え、その一部を運動エネルギーに変換して動き回る個体や、そのような個体の集団のことを指す。古典力学系では多くの研究がなされており、従来の物理学では扱いづらかったエネルギー・運動量が非保存の系の統一的な理解が進んでいる。量子力学系への拡張は発展途上であるが、足立らの論文[1]を皮切りに、研究が盛り上がりつつある[2-5]。

本講演では、非ユニタリ量子ウォークを用いた、量子アクティブマターの最小理論モデル[2]を提案する。この量子アクティブ粒子のモデルで、古典アクティブBrown粒子[6]との類似点と、量子性を同時に観測した。類似点としては、古典アクティブBrown粒子のエネルギー流入のような、非エルミート性のパラメータを大きくすると非自明に運動がよりアクティブになる様子を観測した。量子性としては、ピークの弾道的な伝搬（1次元）、等エネルギー面上での運動（2次元）、基底状態と励起状態の間での共鳴起源の振動を観測した。

後半では、離散時間量子ウォークを高次元の格子系に拡張する手法[7]について紹介する。二次元系における量子ウォークのモデルはいくつか知られているが、1次元系のように、連続極限の有効ハミルトニアンがDiracハミルトニアンと合致するような離散時間量子ウォーク[8]のモデルはこれまでなかった。そこで、2次元系でも同様の極限が取れるモデル[7]を提案し、その振る舞いとトポロジカルな性質（エッジ状態、コーナー状態）について議論する。

[1] K. Adachi, K. Takasan and K. Kawaguchi, Phys. Rev. Research 4, 013194 (2022).
[2] M. Yamagishi, N. Hatano and H. Obuse, arXiv:2305.15319 (2023).
[3] Y. Zheng and H. Löwen, arXiv:2305.16131 (2023).
[4] R. Khasseh et al., arXiv:2308.01603 (2023).
[5] K. Takasan, K. Adachi and K. Kawaguchi, arXiv:2308.04382 (2023).
[6] F. Schweitzer, W. Ebeling and B. Tilch, Phys. Rev. Lett. 80, 5044 (1998).
[7] M. Yamagishi, N. Hatano, K.-I. Imura and H. Obuse, Phys. Rev. A 107, 042206 (2023).
[8] F. W. Strauch, Phys. Rev. A 73, 054302 (2006).

第96回

日時: 2024年2月6日(火)11:00-
発表者: 深井 康平 (東大物性研)
発表題目: 一次元Hubbard模型における局所保存量の表式
Registration form

量子可積分系はBethe仮設を用いた厳密解を有するが、その背景には無限個の局所保存量の存在がある。量子逆散乱法における転送行列の議論より、局所保存量の存在自体は保証されるが、それらの具体的な表式を導くことは困難である。代表的な量子可積分系であるspin-1/2 XYZ鎖と一次元Hubbard模型について、これまで特にそれらの表式が調べられてきた。XYZ鎖については、Boost operatorを用いた局所保存量の漸次的構成法が知られ[1]、最近全ての局所保存量の表式が陽に得られた[2]。XXZの場合については、Temperley-Lieb代数を用いた構成法も最近発見された[3]。

一方、一次元Hubbard模型の局所保存量の一般的な構造は未解明であった。Hubbard模型の場合の難しさとして、Boost operatorを用いた漸次的構成法が適用できないという点がある。これはHubbard模型のR行列がspectral parameterの差分の依存性を持たないことに起因する。これまで、低次の3つの局所保存量の表式が得られていたが、より高次の具体的な表式については知られていなかった[1]。

本研究では、一次元Hubbard模型の全ての局所保存量の表式を求めた[4,5]。それらはある特殊なXX鎖の局所保存密度の積の線型結合から構成され、高次の保存量では非自明な係数が線型結合に登場し、その係数を計算する漸化式を導出した。一部の係数はXXX鎖の局所保存量にも登場するCatalan数[1]であることが分かった。今回の結果はBoost operatorが適用できない量子可積分系において、局所保存量の一般的な表式が求まった最初の例である[4]。今回求まった局所保存量と独立な局所保存量が存在しないことも示した[5]。

[1] M. P. Grabowski and P. Mathieu, Ann. Phys. 243, 299 (1995).
[2] Y. Nozawa and K. Fukai, Phys. Rev. Lett. 125, 90602 (2020).
[3] B. Nienhuis and O. E. Huijgen, J. Phys. A: Math. Theor. 54, 304001 (2021).
[4] K. Fukai, Phys. Rev. Lett. 131, 256704 (2023).
[5] K. Fukai, arXiv:2309.09354 (2023).

第95回

日時: 2024年1月30日(火)11:00-This talk will be given in English
発表者: Rico Pohle (Keio Univ.)
発表題目: Ground state of the S=1/2 pyrochlore Heisenberg antiferromagnet: A quantum spin liquid emergent from dimensional reduction

Registration form

The quest for quantum spin liquids (QSLs) has posed challenges to both experimental and theoretical disciplines for approximately 50 years since their conceptual proposal on the triangular lattice [1]. Two-dimensional frustrated magnets, such as Kitaev, triangular, or kagome materials, have proven to be excellent platforms for the exploration of QSLs [2]. In three dimensions, the pyrochlore lattice stands out as a prime example, offering a diverse spectrum of spin liquids, including spin ice, quantum spin ice, and higher-rank spin liquids [3,4]. These have been comprehensively understood through classical approaches and their perturbative expansions.

The ground state of the fully quantum S=1/2 Heisenberg antiferromagnet (HAF) on the pyrochlore lattice has been a subject of intense debate for over three decades, with conflicting arguments suggesting both the presence and the absence of a QSL ground state. Recent state-of-the-art numerical investigations propose a symmetry-broken and essentially classical ground state [5]. However, to achieve a convincing and conclusive understanding, theoretical work must incorporate a comprehensive finite-size analysis to make reliable estimates in the thermodynamic limit. Without such an analysis, the nature of the ground state remains an open question.

In this talk, we demonstrate that the ground state of the S=1/2 pyrochlore HAF is indeed a QSL [6]. By employing the state-of-the-art “many-variable variational Monte Carlo” method, we solve the full quantum many-body Hamiltonian and unveil an emergent dimensional reduction in the ground state with dominant algebraic correlations in a confined 2D subspace. By examining an effective 2D lattice model, we validate the algebraic decay of spin and singlet correlations in large system sizes and observe a vanishing excitation gap in the thermodynamic limit. The wave-function structure supports the fractionalization of spins into spinons. Importantly, this QSL proves to be robust against spin-orbit interactions, expanding the possibilities of realizing QSLs in real materials such as iridium pyrochlore oxides.

[1] P. Anderson, Mater. Res. Bull 8, 153 (1973).
[2] L. Balents, Nature 464, 199 (2010).
[3] J. G. Rau and M. J. Gingras, Annu. Rev. Condens. Matter Phys. 10, 357 (2019).
[4] H. Yan, et al., PRL 124, 127203 (2020).
[5] I. Hagymási, et al., PRL. 126, 117204 (2021), N. Astrakhantsev, et al., PRX 11, 041021 (2021), R. Schäfer and B. Placke et al., Phys. Rev. Lett. 131, 096702 (2023).
[6] R. Pohle, Y. Yamaji, M. Imada, arXiv:2311.11561

第94回

日時: 2023年12月19日(火)11:00-
発表者: 金賀 穂 (千葉大)
発表題目: 電子系・スピン系における非線形光学応答の制御法についての理論解析

Registration form

固体電子系（グラフェンや半導体・モット絶縁体など）やスピン系（反強磁性体・スピン鎖・キタエフ磁性体など）における高次高調波発生（HHG）が精力的に研究されている[1]。HHGを始めとする非線形光学応答の性質は結晶の対称性に大きく左右され、特に空間反転対称性を持つ系では偶数次高調波と光整流の発生が禁止される。空間反転対称な系において偶数次高調波や光整流を生成するには何らかの外因的な要素によって動的対称性[2]を破る必要がある。この外因要素の導入は、最初から空間反転対称性の破れた系を考えるより複雑なセットアップを要するが、一方、物質を変えることなく偶数次高調波の発生をコントロールできるという利点をもたらす。

反転対称性で禁止された応答を引き出すための有望な「外因的な要素」として、(i) DC電流(電場)の印加[3-5]、あるいは、(ii) 2色レーザーの印加[6,7]が考えられる。本講演では(i), (ii)に関連し、量子マスター方程式[8]を用いることで散逸の効果を考慮した数値解析の結果について紹介する[4,5,7]。まず固体電子系において、グラフェンにおけるDC電流による偶数次高調波の誘起と2色レーザーによる光電流生成について解析した[5,7]。広いパラメータ空間で摂動的領域と非摂動的領域の両方に同時にアクセスできる数値解析のメリットを活かして、特徴的なスペクトル構造の変化を捉えた。次に、スピン系、特に電気磁気結合を持つキタエフスピン液体[9]におけるHHGを解析した[4]。磁歪型電気磁気結合の導入による対称性の低下から第二高調波が発生することや、スペクトルの特徴的な磁場依存性を見出した。講演ではこれらの解析結果の詳細について報告する。

[1] O. Schubert, et al., Nat. Photon. 8, 119 (2014); N. Yoshikawa, et al., Science 356, 736 (2017); Z.-Y. Zhang, et al., Nat. Commun. 14, 1795 (2023), etc.
[2] O. Neufeld, D. Podolsky, and O. Cohen, Nat. Commun. 10, 1 (2019).
[3] S. Wu, et al., Nano Lett. 12, 2032 (2012); K. Takasan, et al., Phys. Rev. B 104, L161202 (2021); J. L. Cheng, et al., Opt. Express 22, 15868 (2014).
[4] M. Kanega, T. N. Ikeda, and M. Sato, Phys. Rev. Res. 3, L032024 (2021).
[5] M. Kanega and M. Sato, in preparation.
[6] Y. Ikeda, S. Kitamura, and T. Morimoto, Phys. Rev. Lett. 131, 096301 (2023).
[7] M. Kanega and M. Sato, in preparation.
[8] H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems (Oxford University Press, 2007); T. N. Ikeda and M. Sato, Phys. Rev. B 100, 214424 (2019).
[9] A. Kitaev, Ann. Phys. 321, 2 (2006); Y. Motome and J. Nasu, J. Phys. Soc. Jpn. 89, 012002 (2020); Y. Tokura, S. Seki, and N. Nagaosa, Rep. Prog. Phys. 77, 076501 (2014).

第93回

日時: 2023年11月28日(火)11:00-
発表者: 大熊 信之 (九工大)
発表題目: 2粒子トポロジーを用いた分数チャーン絶縁体の特徴付け

Registration form

強相関系におけるトポロジー効果に関する研究は、自由粒子系の場合に比べ多くの課題を残している。 そのような状況にあって「分数量子ホール効果」の物理は、長年の研究蓄積のある数少ない研究領域の一つである。 分数量子ホール系は、強磁場下で形成されるランダウ準位と強相関効果の協奏で生じる事が知られ、エニオン励起状態や基底縮退の観点から、量子コンピュータへの応用が期待される物性系の一つである。 また理論面では、ランダウ準位の持つ解析的性質から、2次元強相関系において厳密な議論が出来る数少ない例として知られている。

強磁場を必要とする分数量子ホール系に対し、注目するバンドがランダウ準位と同じトポロジカル構造を持つ「分数チャーン絶縁体 [1]」は、強磁場が不要であるという実験的な利点を持つ。 一方理論面からみると、ランダウ準位の持つ解析的な性質を失う事から、その取り扱いの難しさは典型的な強相関系のそれとなる。 厳密な議論の出来ない分数チャーン絶縁体の理論的探索は従来、ランダウ準位と似た性質を持つチャーンバンドを探す方向性で進展してきた。 特に、バンド構造の平坦条件や量子幾何テンソルの理想的条件 [2]といった、1粒子の性質に着目した研究が盛んに行われている。 本発表ではそれらの試みをレビューしつつ、我々の行った「2粒子束縛状態のトポロジー」に注目した研究 [3]について紹介する。

[1] T. Neupert, L. Santos, C. Chamon, and C. Mudry, PRL 106, 236804 (2011).
[2] R. Roy, Phys. Rev. B 90, 165139 (2014).
[3] NO and T. Mizoguchi, Phys. Rev. Res. 5, 013112 (2023).

第92回

日時: 2023年11月7日(火)11:00-
発表者: 鏡原 大地 (近大)
発表題目: 理想Bose粒子系および強く相互作用するBose粒子系におけるRényiエンタングルメントエントロピーのクエンチダイナミクス

Registration form

エンタングルメントは量子力学特有の非局所的な相関のことであり、現在では量子多体系の理解に欠かせない概念となっている。エンタングルメントは典型的にはエンタングルメントエントロピーを用いて定量化され、その着目系のサイズ依存性の振る舞いは様々な有用な情報を提供する。近年、光格子中の1次元Bose粒子系において、Rényiエンタングルメントエントロピーの実験測定が実現し注目を集めている[1]。我々は、実験状況に即し光格子中の1次元Bose粒子系を調べ、相互作用のない極限[2]および強く相互作用する領域[3]の2つの領域において、解析的にRényiエンタングルメントエントロピーを導出した。

論文[2]では、Mott絶縁体状態や電荷密度波状態などの積状態の初期状態から出発し相互作用のないハミルトニアンへクエンチした場合のダイナミクスを調べた。この場合、Rényiエンタングルメントエントロピーは1粒子相関関数からなる行列の永久式(permanent)で特徴付けられることが分かった。得られた関係式を利用し、エンタングルメントエントロピーが着目系のサイズに比例する体積則の条件を得るとともに、先行研究より大きなシステムサイズの系のシミュレーションを実行した。

論文[3]では、初期状態をMott絶縁状態とし、ホッピング振幅のわずかな増加によって誘導されるダイナミクスを調べた。この場合の低エネルギー励起状態は、クエンチによってエンタングルされた対として励起される、doublonとholonというフェルミオン準粒子によってよく記述される。Rényiエンタングルメントエントロピーをこれらの準粒子の相関関数を用いて解析的に表せることが分かった。また、得られた表式に基づいて、Rényiエンタングルメントエントロピーの振る舞いの物理的な解釈を与えた。

[1] R. Islam, R. Ma, P. M. Preiss, M. Eric Tai, A. Lukin, M. Rispoli, and M. Greiner, Nature (London) 528, 77 (2015); A. M. Kaufman, M. E. Tai, A. Lukin, M. Rispoli, R. Schittko, P. M. Preiss, and M. Greiner, Science 353, 794 (2016).
[2] D. Kagamihara, R. Kaneko, S. Yamashika, K. Sugiyama, R. Yoshii, S. Tsuchiya, and I. Danshita, Phys. Rev. A 107, 033305 (2023).
[3] S. Yamashika, D. Kagamihara, R. Yoshii, and S. Tsuchiya, arXiv:2209.13340 (accepted in Phys. Rev. Res.).

第91回

日時: 2023年10月24日(火)11:00-
発表者: 新城 一矢 (理研)
発表題目: 光励起したモット絶縁体の吸収スペクトル

Registration form

モット絶縁体は銅酸化物高温超伝導体の母物質であることもあり，これまでに膨大な研究が重ねられてきた．そのため，できることはやり尽くされた分野のイメージもあるかもしれない．しかし，最近のパルスレーザーを駆使した実験に目を向けると，サブフェムト秒という電子の運動を刻々と追えるような超高速分光が行われ始め，高強度テラヘルツ電場を印加した実験も活発化し始めた，新しい研究領域が現在進行形で広がっているホットな分野である．高精度で制御されたパルス電場はモット絶縁体の複雑な電子状態の情報を引き出すのに効果的なだけでなく，新奇量子多体状態を生み出す手段としても優れた役割を果たす．

このような光励起状態を理論的に理解するには，量子多体系の非摂動的数値計算による取り扱いが，しばしば必須になる．そのための主要なアプローチとして，これまで時間依存数値対角化法や非平衡動的平均場理論が広く使われてきた．本セミナーでは，これらの相補的な手法として期待されている時間依存密度行列繰り込み群法を用いて，光励起状態を解析した結果を紹介する [1-6]．線形吸収スペクトルは，モット絶縁体の光励起ダイナミクスを理解するための最も基礎的な情報を含む．我々は，時間依存密度行列繰り込み群法を用いることで，一次元だけでなく，二次元ハバード模型でも実験と比較可能な線形吸収スペクトルを計算することができた [1-3]．その結果，次元性に依存したスピンと電荷の結合・分離の様子が吸収スペクトルに色濃く現れることがわかった．さらに，ポンプパルスでハバード模型を励起した状態の過渡吸収スペクトルを計算すると，元々のモット絶縁体の性質が大きく変化することがわかった [4-6]．特に，系の次元性だけでなく，ポンプパルスの周波数によっても，吸収の構造は大きく変化することがわかったので，その詳細について明らかになったことを紹介する．

[1] S. Ohmura, A. Takahashi, K. Iwano, T. Yamaguchi, K. Shinjo, T. Tohyama, S. Sota, and H. Okamoto, Phys. Rev. B 100, 235134 (2019).
[2] K. Shinjo, S. Sota, and T. Tohyama, Phys. Rev. B 103, 035141 (2021).
[3] K. Shinjo, Y. Tamaki, S. Sota, and T. Tohyama, Phys. Rev. B 104, 205123 (2021).
[4] K. Shinjo, S. Sota, and T. Tohyama, Phys. Rev. Res. 4, L032019 (2022).
[5] K. Shinjo, S. Sota, S. Yunoki, and T. Tohyama, Phys. Rev. B 107, 195103 (2023).
[6] T. Tohyama, K. Shinjo, S. Sota, and S. Yunoki, Phys. Rev. B 108, 035113 (2023).

第90回

日時: 2023年7月25日(火)11:00-
発表者: 大門 俊介 (QST)
発表題目: 深層学習AIによる量子伝導現象の解読

Registration form

ナノ量子系の電気伝導は、一般に周期性をもたない不規則な磁場依存性を示すことが知られている。この磁気伝導度ゆらぎは、試料内で電子が複雑に散乱・干渉することで引き起こされ、試料固有の振動が観測されることから「量子指紋」と呼ばれる。一方で、不規則に振動する量子指紋を解読することは従来の解析手法では困難であり、ミクロな内部情報を抽出することはできなかった。そこで本研究では、量子指紋を解読可能な深層生成ネットワークを開発し、複雑な磁気伝導度ゆらぎからミクロな内部情報の抽出に挑戦した。講演では、開発ネットワークによる試料のミクロ構造および量子干渉情報の生成と、実験データへの応用について報告する。

第89回

日時: 2023年6月27日(火)11:00-
発表者: 水田 郁 (東大)
発表題目: 時間周期系/準周期系の最適な量⼦計算

Registration form

量⼦多体系の時間発展に対する量⼦計算は、物性物理・量⼦化学への幅広い応⽤が期待される量⼦計算の最も基本的かつ重要な対象である。それ故に、時間発展をどこまで⾼速かつ正確に計算できるか、ということが量⼦計算初期以来の重⼤な問いである。それが近年、時間⾮依存の量⼦多体系の時間発展に対しては、Hamiltonianの⼀般的な⾏列多項式を印加する量⼦特異値変換アルゴリズムという⼿法により時間発展の時刻・許容誤差について理論上最も効率よく計算できることが⽰された[1]。量⼦特異値変換は素因数分解,探索問題など様々な量⼦アルゴリズムを包含し”量⼦アルゴリズムの⼤統⼀理論”と呼ばれるほどに汎⽤性が⾼い⼿法であるが、時間依存する系は扱えない。時間依存系はDyson級数を含むなどの時間発展演算⼦の複雑さもあり、時間⾮依存系ほど効率的な⼿法は確⽴されていなかった。

そこで我々は時間依存系の中でもH(t)=H(t+T)という時間周期性を持つ時間周期系に着⽬し、その時間発展を計算する量⼦アルゴリズムを構築した[2]。具体的には、Floquet理論・Lieb-Robinson限界といった⾮平衡系の⼿法と、量⼦特異値変換・振幅増幅などの量⼦計算の⼿法を組み合わせ、時間周期系に対しても時間⾮依存系とほぼ同じ計算コストで時間発展を計算できることを⽰した。時間周期系は、光照射下の物質など物性物理・量⼦化学で最も関⼼のある⾮平衡量⼦系のクラスであるため、我々の結果は量⼦計算のそれらへの応⽤性を広めるものである。同時に、本⼿法は時間準周期系などへも拡張可能であり[3], より広範な時間依存系に対する量⼦アルゴリズムの探索にも繋がると期待される。本講演では、量⼦特異値変換などの量⼦計算のレビューから始め、我々の結果である時間周期系/準周期系の時間発展の最適な量⼦計算⼿法について議論する。

[1] J. M. Martyn, Z. M. Rossi, A. K. Tan, and I. L. Chuang, PRX Quantum 2, 040203 (2021).
[2] K. Mizuta and K. Fujii, Quantum 7, 962 (2023).
[3] K. Mizuta, arXiv: 2301.06232.

第88回

日時: 2023年6月13日(火)11:00-
発表者: 越智 正之 (阪大)
発表題目: 波動関数理論に基づく固体の第一原理計算手法の開発

Registration form

固体の第一原理電子状態計算は多くの場合、密度汎関数理論に基づいて行われるが、近年は多体波動関数をあらわに用いる、波動関数理論に基づいた計算も活発に行われている。波動関数理論には、Hartree-Fock（HF）波動関数（Slater行列式）を出発点として、多数のSlater行列式の線型結合を考える理論、電子間相互作用の効果を摂動的に取り扱う理論、Jastrow-Slater型の波動関数を用いる理論など様々なものが存在するが、一般には複雑な波動関数を用いるほどその計算コストも劇的に増加してしまう。

本講演では、Jastrow相関因子によってハミルトニアンを相似変換することで電子相関効果を取り込むトランスコリレイティッド（TC）法[1]の紹介を行う。その特徴のひとつは、Jastrow-Slater型の波動関数を扱いつつも、その計算時間が（粒子数について）HF法と同等のオーダーで済むことである[2]。また、相似変換ハミルトニアンについてpost HF理論を適用できるなど、他の波動関数理論との組み合わせもできる。最近、固体に対する計算コードTC++を公開した[3]のでその紹介も行う。

[1] S. F. Boys and N. C. Handy, Proc. R. Soc. London Ser. A 309, 209 (1969).
[2] M. Ochi et al., J. Chem. Phys. 136, 094108 (2012).
[3] https://github.com/masaochi/TC

第87回

日時: 2023年5月30日(火)11:00-
発表者: 幸城 秀彦 (理研)
発表題目: 多段磁化プラトーを実現するフラストレートスピンラダーの理論的研究

Registration form

S=1スピンラダー系が良い有効模型だと考えられている有機磁性体BIP-TENOの磁化測定において複数の磁化プラトーが観測された[1]. 磁化プラトーは磁性体に磁場を印加した際に現れる, 磁場の変化に対して磁化が変化しない現象であり, エネルギーギャップの空いた安定状態の存在を示唆する. BIP-TENOの磁化過程で特徴的なのは, 磁場の発生方法によって観測される磁化プラトーの数が異なることである. ミリ秒程度の時間スケールで印加磁場が変化する場合では全磁化の1/4の大きさを持つプラトーのみが観測されている. 一方, マイクロ秒程度の時間スケールで変化する場合では, 全磁化の1/4の大きさを持つプラトーに加えて, 全磁化の1/3, 1/2の大きさを持つプラトーが観測されている. 磁化プラトーの数が磁場の発生方法によって増減する機構や, これらの磁化プラトーを実現する磁気構造に関しては未解明である.

この現象を説明する単純な仮説をとして, 磁化過程の間にスピンは十分緩和しており, 格子は高速磁場に追従できず緩和できていないという仮説が考えられる. この仮説に基づき異なる2つのスピン模型の基底状態を考えることで異なる2つの磁化過程を説明できる可能性を検討した. そのために, まず観測されている2つの磁化過程を再現する有効模型を考える必要があるが, これまで1/4, 1/3, 1/2プラトーを基底状態に持つS=1スピンラダー模型が知られていなかった. 我々はこれら複数の磁化プラトーを実現する可解なS=1スピンラダー模型を発見し, 解析的に磁化プラトー状態が基底状態となることを示した[2]. 本講演ではこの可解なS=1スピンラダー模型の基底状態の構成方法について詳しく述べる. 時間が許せば解析的に解くための条件を緩めた模型での数値計算結果や実験結果との関連性についても述べる.

[1] Kazuya Nomura, Yasuhiro H. Matsuda, Akihiko Ikeda, Yoshimitsu Kohama, Hiroshi Tsuda, Naoki Amaya, Toshio Ono, and Yuko Hosokoshi, Phys. Rev. B 105, 214430 (2022).
[2] Hidehiko Kohshiro, Ryui Kaneko, Satoshi Morita, Hosho Katsura, and Naoki Kawashima, Phys. Rev. B 104, 214409 (2021).

第86回

日時: 2023年5月9日(火)11:00-
発表者: 小沢 耀弘 (東大)
発表題目: 強磁性積層カゴメ系ワイル半金属Co3Sn2S2の有効模型を用いた磁気秩序とスピン電荷輸送の理論解析

Registration form

ワイル半金属とは、ワイル点と呼ばれる縮退を持たないギャップレス点を有するトポロジカル半金属であり、2011年にその報告がなされた[1][2]。 特に、磁気秩序によって時間反転対称性が破れた磁性ワイル半金属は、異常ホール効果に代表される特異なスピン依存伝導を示す。 積層カゴメ格子強磁性体であるCo3Sn2S2は、2018年に報告された磁性ワイル半金属の有力な候補物質の1つである[3]。 この系は、ワイル点のベリー曲率に起因した巨大な異常ホール効果を示し、フェルミ準位に掛かる金属バンドも比較的少ない。 つまり、磁性ワイル半金属特有の電磁応答を調べる条件が揃った系といえる。 磁性ワイル半金属Co3Sn2S2における電磁応答を理論的に理解・予測するためには、少ない自由度を用いた単純な模型による解析が必要である。

この講演では、磁性ワイル半金属Co3Sn2S2の有効低エネルギー模型を提案し、それを用いた磁気秩序とスピン/電荷輸送の解析を紹介する。 この模型は、カゴメ格子を構成するCoのd軌道と、それらに挟まれた三角格子を構成するSnのp軌道を用いた2軌道dpモデルである[4]。 スピン自由度を含めると8x8行列の比較的小さなハミルトニアンで表されるが、第一原理計算で得られたフェルミ準位付近のワイル点構造や異常ホール伝導率の極大を再現する。 また、電子間相互作用を平均場近似により導入すると、非ドープでの面直強磁性秩序に加え、電子数の変化による常磁性および反強磁性が生じることを示した[5]。 特に反強磁性状態では、正味の磁化はゼロである一方、異常ホール伝導率と軌道磁化は有限となることを見出した。 最後に、非磁性状態 (ディラック半金属相)および強磁性状態 (ワイル半金属相)のスピンホール効果について議論し、時間が許せば実験グループとの共同研究についても紹介する[6][7]。

[1] X. Wan, A. M. Turner, A. Vishwanath, and S. Y. Savrasov, Phys. Rev. B 83, 205101 (2011) [2] A. A. Burkov and L. Balents, Phys. Rev. Lett. 107, 127205 (2011). [3] E. Liu et al., Nat. Phys. 14, 1125 (2018) [4] A. Ozawa and K. Nomura, J. Phys. Soc. Jpn. 88, 123703 (2019). [5] A. Ozawa and K. Nomura, Phys. Rev. Materials 6, 024202 (2022). [6] Y. C. Lau, J. Ikeda, K. Fujiwara, A. Ozawa, T. Seki, K. Nomura, A. Tsukazaki and K. Takanashi, arXiv:2203.02356 (2022). [7] T. Seki, Y.-C. Lau, J. Ikeda, K. Fujiwara, A. Ozawa, S. Iihama, K. Nomura and A. Tsukazaki, Phys. Rev. Research 5, 013222 (2023).

第85回

日時: 2023年4月25日(火)11:00-
発表者: 佐藤 駿丞 (筑波大)
発表題目: 第一原理電子ダイナミクス計算とアト秒光科学への応用

Registration form

近年の光科学技術の発展により、光が駆動する固体中の電子ダイナミクスをアト秒の時間分解能で観測することが出来るようになってきた。このようなアト秒分光技術の固体物理への応用は、量子多体系の非線形・非平衡なダイナミクスに関する重要な知見を与える。その一方で、固体の複雑な電子構造を反映した実験結果を直接解釈することは難しいのが現状である。時間依存密度班関数理論(TDDFT)[1]に基づく第一原理電子ダイナミクス計算は、このような複雑な電子ダイナミクスを微視的に記述・理解するための有用なツールである。

本セミナーではまず、固体中の電子ダイナミクスを解析する手法として実時間TDDFTシミュレーションを概観する。その後TDDFTシミュレーションの応用として、光によって駆動された系の非線形・非平衡な応答を調べるための第一原理ポンプ・プローブ計算について解説する[2]。また、このポンプ・プローブ計算の応用例として、アト秒過渡吸収分光の第一原理シミュレーションを金属Ti[3]や磁性体Co[4]に適用し、これらの金属内部に光駆動される電子やスピンのダイナミクスについて解析した結果についても紹介する。

[1] Erich Runge and E. K. U. Gross, Phys. Rev. Lett. 52, 997 (1984)
[2] Shunsuke A. Sato, Comput. Mater. Sci. 194, 110274 (2021)
[3] M. Volkov, S. A. Sato, F. Schlaepfer, L. Kasmi, N. Hartmann, M. Lucchini, L. Gallmann, A. Ru-bio, U. Keller, Nature Physics 15, 1145 (2019)
[4] Shunsuke A. Sato, Electron. Struct. 4, 014007 (2022)

第84回

日時: 2023年4月11日(火)11:00-
発表者: 仲澤 一輝（東大)
発表題目: 量子輸送現象によって検出する物質中のトポロジカル構造

Registration form

　物質中で創発するトポロジカルな構造は，量子Hall効果やBKT転移などを端緒として，精力的に調べられてきた．前者は電子バンドが，後者はスピンがつくるトポロジカル構造であり，その後トポロジカル絶縁体，Weyl半金属，そして磁気スカーミオンといったコンセプトへと発展した．本講演ではこのような物質中のトポロジカルな構造を量子輸送現象によって検出することを念頭に，異なる2つのテーマについて議論する．

　一つ目に，カゴメ強磁性体Co3Sn2S2の原子層薄膜のトポロジカル物性について報告する．この物質では，自発的な時間反転対称性の破れによるWeyl点に起因する巨大な異常Hall効果や負の磁気抵抗効果などが報告されている[1]．一方で，Coカゴメ単層物質においてChern絶縁体状態が報告されており[2]，この実験的観測が待たれる．ところが，3次元空間のトポロジカルな特異点であるWeyl点が2次元物質におけるChern絶縁体状態にどのように反映されているのかは非自明であり，これを理解することは極薄膜作製の上でも重要である．我々は第一原理計算に基づき，Coカゴメ単層・2層・3層構造について，表面終端原子がSnの場合とSの場合を考え，トポロジカル物性の原子層数依存性を調べた．その結果，表面終端原子がSnの場合には強磁性，バンドトポロジー，異常Hall伝導度などの輸送特性が薄膜からbulkにかけて系統的に接続する一方，Sの場合にはこれらの性質が層数に依存し，原子層薄膜では表面終端原子によって物理的特性が大きく異なることを見出した[3]．

　二つ目に，反強磁性金属におけるスピンHall効果について議論する．磁化のつくるスカラーカイラリティーに誘起されるトポロジカルHall効果はよく知られているが，近年ではNéelベクトルのそれが存在するときのトポロジカルスピンHall効果(TSHE)が議論されており，反強磁性スカーミオン系などで調べられている[4]．ところが，これまでのスカラーカイラリティーによる解釈に疑問が残るほか，強結合から弱結合にわたる包括的な理解が得られていない．我々はTSHEがNéelベクトルのスカラーカイラリティーではなく，むしろベクトルカイラリティーに誘起されること，反強磁性スカーミオンでは現れず，反強磁性メロン（半スカーミオン）で現れることを見出した．また，反強磁性金属特有の横スピン緩和(dephasing)を考慮することで，スピンHall伝導度の強結合から弱結合にかけての振る舞いを議論し，特に弱結合の場合にスピンHall伝導度が増大することを示した[5]．

[1] E. Liu et al., Nat. Phys. 14, 1125 (2018); Q. Wang et al., Nat. Commun. 9, 3681 (2018).
[2] L. Muechler et al., Phys. Rev. B 101, 115106 (2020).
[3] K. Nakazawa, Y. Kato, and Y. Motome, arXiv:2212.09026.
[4] B. Göbel, A. Mook, J. Henk, and I. Mertig, Phys. Rev. B 96, 060406(R) (2017); C. A. Akosa, O. A. Tretiakov, G. Tatara, and A. Manchon, Phys. Rev. Lett. 121, 097204 (2018).
[5] K. Nakazawa, K. Hoshi, J. J. Nakane, J. Ohe, and H. Kohno, in preparation.

第83回

日時: 2023年2月28日(火)11:00-
発表者: 中西優馬（東工大)
発表題目: PT対称性由来の散逸時間結晶

Registration form

時間結晶は、「連続時間並進対称性を自発的に破り、離散時間並進対称性を保った非自明な状態」として特徴づけられ、ノーベル物理学者であるWilczekによって導入された[1]。しかし、孤立系の基底状態や熱平衡状態では、広いクラスの物理系において時間結晶は存在しないことが示された[2]。その後時間結晶は、散逸系やフロケ系などの非平衡系に拡張され、理論的および実験的に活発に研究されている。特に、環境と相互作用したGKSL方程式で記述される集団散逸スピン系において、熱力学極限にのみ現れる時間結晶は境界時間結晶(BTC)と呼ばれている[3]。

本セミナーでは、BTCがパリティ・時間反転（PT）対称性由来の散逸時間結晶と解釈できる根拠をいくつか提示する[4,5]。まず、BTCを引き起こす最も基本的なモデルである、１スピンBTCモデル[3]がGKSL方程式のPT対称性を満たし、BTC相では、定常状態のPT対称性が保持されることを解析的に示す。さらに、他のBTCモデルについても同様のことが成立していることを数値的に示す。最後に、弱い散逸下での摂動を用いて、ゲインとロスのある集団１スピン系のクラスにおいて、ゲインとロスが釣り合った場合のみ、BTCが現れることを示す。これらの結果は、BTCがPT対称性由来の時間結晶であることを強く示唆している。

[1] F. Wilczek, Phys. Rev. Lett. 109, 160401 (2012).
[2] H. Watanabe and M. Oshikawa, Phys. Rev. Lett. 114, 251603 (2015). [3] F. Iemini et al., Phys. Rev. Lett. 121, 35301 (2018).
[4] YN and T. Sasamoto, Phys. Rev. A 105, 022219 (2022). [5] YN and T. Sasamoto, Phys. Rev. A 107, L010201 (2023).

第82回

日時: 2023年2月7日(火)11:00-
発表者: 大岩 陸人（明治大)
発表題目: 対称性適合多極子基底を用いた有効模型生成

Registration form

近年、ワニエ関数を用いた密度汎関数理論（DFT）に基づく有効模型計算法が発展しており[1]、各種物性の定量的解析に加え、物質が示す隠れた秩序を特徴付ける秩序変数の特定などに応用されている。しかし、ワニエ関数を用いた有効模型は、系の対称性を完全には満たさず、遠方の飛び移り積分を多数含むため、模型が複雑になり、電子格子相互作用等の計算コストが増大するという課題がある。また、隠れた秩序に対応する演算子表現と関連物性の微視的発現機構を系統的に明らかにするためには、新たな理論的アプローチを開発する必要がある。

そこで我々は、物質の電子状態を詳細に考慮できるDFT計算法と、電子自由度を系統的に記述できる電子多極子表現論[2,3]を組み合わせることで、汎用的な有効模型計算法を開発した。本手法で求まる有効ハミルトニアンは、対称性適合多極子基底の線形結合で表されるため、結晶の対称性を完全に満たす。このとき、各基底の重みはモデルパラメータに対応する。DFT計算で求まるバンド分散を再現するように重みを最適化することで、低エネルギーの有効模型を構築できる。また、本手法を用いることで物質中の隠れた多極子自由度とその演算子表現を明らかにでき、それらの寄与を定量的に評価することも可能である。本講演では、グラフェンやSrVO3、単層MoS2などの典型物質に加え、カイラル結晶の単体テルルに対して本手法を適用した結果を示す[4]。

[1] N. Marzari, A. A. Mostofi, J. R. Yates, I. Souza, and D. Vanderbilt, Rev. Mod. Phys. 84, 1419 (2012).
[2] S. Hayami, M. Yatsushiro, Y. Yanagi, and H. Kusunose, Phys. Rev. B 98, 165110 (2018).
[3] H. Kusunose, RO, and S. Hayami, J. Phys. Soc. Jpn. 89, 104704 (2020).
[4] RO and H. Kusunose, Phys. Rev. Lett. 129, 116401 (2022).

第81回

日時: 2023年1月24日(火)11:00-
発表者: 野垣 康介（京大)
発表題目: 局所空間反転対称性の破れが誘起する新奇超伝導相におけるトポロジーと強相関効果

Registration form

奇パリティ超伝導状態は、トポロジカル超伝導実現の良い舞台であり、長年にわたり精力的な研究が続けられてきた。電子スピンのみを考慮する超伝導理論では、スピン三重項超伝導が奇パリティ超伝導の候補となる。一方で、スピン以外の内部自由度を考慮することで新しい奇パリティ超伝導を実現することができる。特に、局所空間反転対称性の破れにより活性化する副格子自由度に由来した奇パリティ超伝導の可能性が理論的に指摘された[1]。

2021年に、局所空間反転対称性が破れた新奇重い電子系超伝導体CeRh2As2が発見され、大きな注目を集めている[2]。特にその温度-磁場相図は理論的に指摘されていた相図と良い一致を見せており、新たな奇パリティ超伝導状態の候補物質となり得る。

本セミナーでは、CeRh2As2におけるトポロジカル超伝導の可能性とその相図における強相関効果に関する最近の進展について紹介する。講演の前半では、群論に基づいてフェルミ面公式を導出し、第一原理計算の結果と組み合わせることでトポロジカル超伝導であることを予言する[3]。後半では、局所空間反転対称性が破れた強相関超伝導体のミニマルモデルを解析し、CeRh2As2の相図を強相関効果の視点から理解する[4]。

[1] T. Yoshida, M. Sigrist, and Y. Yanase, Phys. Rev. B 86, 134514 (2012).
[2] S. Khim et al., Science, 373, 1012 (2021).
[3] K. Nogaki, A. Daido, J. Ishizuka, and Y. Yanase, Phys. Rev. Research 3, L032071 (2021).
[4[ K. Nogaki and Y. Yanase, Phys. Rev. B 106, L100504 (2022).

第80回

日時: 2022年12月13日(火)11:00- This talk will be given in English
発表者: Kimberly Remund（OIST)
発表題目: Spin-1 Magnet - a u(3) Formalism

Registration form

Systems based on spin-1 moments exhibit many fascinating properties, as demonstrated by spin nematics [1,2], Fe-based superconductors [3], and cold atom systems [4]. However, because spin-1 moments beside displaying dipoles can also support quadrupoles at the level of a single site, they are not properly described by O(3) vectors, even in the classical limit.

In this talk, we explore a new approach which correctly describes the (semi-)classical physics of spin-1 moments [5]. This approach, based on the algebra u(3) [6], provides a convenient starting point for Monte Carlo and molecular dynamics simulations, as well as for analytic theory. We benchmark our method by applying it to the ferroquadrupolar phase of the spin-1 bilinear-biquadratic model on the triangular lattice [5], reproducing and extending known results [1,2].

These results open the prospect of exploring dynamical properties of spin-1 magnets which lie outside the spin-wave paradigm, including spin-liquid phases, and the dynamics of topological excitations.

[1] H. Tsunetsugu and M. Arikawa, J. Phys. Soc. Jpn 75, 083701 (2006)
[2] A. Laeuchli et al, Phys. Rev. Lett. 97, 087205 (2006)
[3] R. M. Fernandes et al., Nature Physics 10, 97 EP (2014)
[4] E. Demler and F. Zhou, Phys. Rev. Lett. 88, 163001 (2002)
[5] K. Remund et al., Phys. Rev. Research 4 033106 (2022)
[6] N. Papanicolaou, Nuclear Physics B 305, 367 (1988)

第79回

日時: 2022年12月6日(火)17:00- 通常と開催時刻が異なります
発表者: 川野 雅敬（ミュンヘン工科大学)
発表題目: サイン二乗変形を用いた長周期量子相の効率的記述

登録フォーム

我々が近年考案した、長周期秩序を示す量子相を効率的に記述する数値計算手法を紹介する。 こういった相は磁気スカーミオン相、高温超伝導体候補物質における電荷・スピン密度波相、ねじれ二重層グラフェンなど枚挙にいとまがなく、 その性質の解明が急務の状況である。 しかしながら、その長周期性のため特に有限系の数値計算による記述が非常に困難であった。 今回我々はサイン二乗変形および変形Fourier変換によりこの困難を解消した [1]。 サイン二乗変形は系のエネルギースケールを空間的に変調し、系の中央付近で熱力学極限の周期をもった秩序状態が実現する。 本研究で新たに導入した変形Fourier変換は、サイン二乗変形下において波数空間の物理量計算を可能にし、 それらは小さいシステムサイズで既に熱力学極限へと収束する。 本講演ではこの手法のベンチマーク計算およびスピン軌道相互作用を有する強相関電子系への応用 [2]を紹介したい。

[1] M. Kawano and C. Hotta, Phys. Rev. Research 4, L012033 (2022).
[2] M. Kawano and C. Hotta, arXiv:2208.09902 (2022).

第78回

日時: 2022年11月15日(火)17:00- 通常と開催時刻が異なります
発表者: 佐藤 年裕（ヴュルツブルク大学)
発表題目: フラストレート量子スピン系への補助場量子モンテカルロアプローチ

登録フォーム

量子モンテカルロ計算は、与えられた格子サイズと温度に対して、すべての熱力学量を近似なしで得ることができ、動的特性や輸送特性を調べることに優れているため、対象モデルの固有の物理的性質の探究やモデル計算と量子物質の多くの実験データとの比較研究を通して、科学的に新たな知見を提供する。しかしながら、フラストレート量子スピン系などの一般的な問題には、量子モンテカルロ計算に特有の負符号問題に悩まされるが、この問題の深刻さは定式化に依存するため、最適化戦略を打ち出すことができる。本講演では、多体フェルミオン系をシミュレートする方法の一つである補助場量子モンテカルロ法を用いて、フラストレート量子スピン系の負符号問題の深刻さを緩和する位相ピンニングアプローチを紹介する[1]。フラストレート量子スピン系の一例として、実験的にも理論的にも継続的に注目されているハニカム格子上の一般化されたキタエフモデルを挙げる。このモデルは、多くの磁気秩序状態と同様に$Z_2$スピン液体状態をサポートすると考えられているだけでなく、電子相関とスピン軌道相互作用が競合する層状ハニカム型イリジウム酸化物や$\alpha$-RuCl$_3$との関係から特に人気がある。このアプローチを用いた我々の最近の取り組みとして、$\alpha$-RuCl$_3$における磁場中の熱力学的性質の研究を紹介する。また、高い対称性を持つ一般化されたキタエフモデルに対して、負符号問題が生じない補助場量子モンテカルロ法の定式化についても紹介する[2]。

[1] T. Sato and F. F. Assaad, Phys. Rev. B. 104, L081106 (2021).
[2] T. Sato and F. F. Assaad, Phys. Rev. B. 106, 155110 (2022).

第77回

日時: 2022年10月25日(火)11:00-
発表者: 足立 景亮（理研)
発表題目: 斥力相互作用するアクティブマターモデルの非平衡協同現象

登録フォーム

超伝導や液液相分離など、多数の要素が集まることで起こる協同現象は固体や液体といった凝縮系で観測され、背後にある相転移のメカニズムや系の詳細によらない普遍的性質の研究はシンプルな理論モデルを使って進められている。一方、生物個体や細胞のようにエネルギーを使って自ら動く要素の集まりはアクティブマターと呼ばれ、魚の群れのように一方向に進む集団運動など、特有の協同現象を示す。近年、アクティブマターを一種の凝縮系とみなし、協同現象の特徴を調べる実験・理論研究が盛んに行われている[1,2]。

本セミナーでは、バクテリアの凝集のようなアクティブマターの凝集現象に注目し、シンプルな粒子モデルを使ってそのメカニズムや普遍的性質を議論する。まず、平衡系の気液相分離や液液相分離とは対照的に、粒子同士を引きつける引力相互作用がなくても、粒子が自ら動く性質と斥力相互作用に起因して凝集が起こりうることを説明する。次に、基板の配向性などによって系が空間的異方性をもつと、パラメータの詳細によらずに粒子数密度の長距離相関が生じ、凝集状態への相転移には外場駆動系と呼ばれる非平衡系と共通する特徴が現れることを解説する[3]。

[1] D. Needleman and Z. Dogic, Nat. Rev. Mater. 2, 1 (2017). [2] G. Gompper et al., J. Phys. Condens. Matter 32, 193001 (2020). [3] K. Adachi, K. Takasan, and K. Kawaguchi, Phys. Rev. Research 4, 013194 (2022).

第76回

日時: 2022年9月6日(火)11:00-
発表者: 桂 ゆかり（東大・NIMS)
発表題目: 論文からの大規模実験データの収集と理論物性科学への活用

登録フォーム

第一原理計算の結果や 理論物性科学で生まれたモデルの妥当性を示すには、実験データと合致することを示すことが重要である。だが論文中で報告されている実験データの多くはグラフ画像に埋もれたままでデジタル化されておらず、簡単には使えない状態であった。

本研究者らが取り組んできたStarrydataプロジェクトでは、論文中で頻出するグラフ画像から実験データを収集することオープンデータベースを作成してきた。そのプロトタイプである熱電特性データベースでは、2016年から現在まで、7200本以上の論文から、33000枚以上のグラフ画像に掲載された、38000試料以上の熱電特性の温度依存性のデータを集めてきた。専属のデータ収集者らが論文を読解して、自動データ点検出機能の補助を借りつつ、グラフを丸ごとトレースして数値化してきた。

これらの大規模実験データを活用することで、第一原理計算結果の検証や、理論モデルの妥当性の検証、機械学習による「研究者の勘」の再現を行った例を共有する。

第75回

日時: 2022年8月30日(火)11:00-
発表者: 三澤 貴宏（BAQIS)
発表題目: トポロジカル絶縁体におけるグリーン関数のゼロ -トポロジカル状態の簡便な可視化方法

登録フォーム

トポロジカル絶縁体は非自明なトポロジカル不変量をもつ絶縁体として特徴づけられる。C. L. KaneとE. J. Meleによる時間反転対称性に保護された2次元のトポロジカル絶縁体の発見[1]を契機として、爆発的な勢いで研究が進展してきた。特に、トポロジカル絶縁体の系統的な分類と、それを指導原理とした新しいトポロジカル絶縁体の発見は、近年の凝縮系物理学のハイライトの一つである[2]。しかしながら、与えられた系がトポロジカル絶縁体であるかを判定するのは一般には簡単ではなく、簡便な判定法の開発が続けられている[3,4]。

最近、我々はグリーン関数の対角成分のゼロが、無相関のトポロジカル絶縁体を検出できる有用な量であることを見出した[5]。特に、トポロジカル相では、バンド反転によって、グリーン関数のゼロがバンドギャップを横断することを示した。ゼロがバンドギャップを横断するかどうかによって、トポロジカル状態を「目で見て」、簡単に判別できるのかこの方法の大きな特徴である。また、具体的な模型でグリーン関数のゼロを計算することで、６個の基本的なクラスのトポロジカル絶縁体のみならず、最近注目を集めている高次トポロジカル絶縁体も判定できることを示した。 面白いことに、解析の鍵となったのは、最近再発見された、線形代数における固有ベクトルと固有値の間に成り立つ関係式「eigenvector-eigenvalue identity」である[6]。長らく見過ごされてきた、線形代数における基本的ともいえるこの関係式がグリーン関数のゼロを通じて、トポロジカル絶縁体の判定に役に立つことを紹介する。時間が許せば、グリーン関数のゼロとエッジ状態の関係、グリーン関数のゼロを用いた現実的な系のトポロジカル相の判定についても解説する。

[1] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 (2005).
[2] 初期のレビューとして、M. Z Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010), X.-L.Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 (2011).
[3] J. Kruthoff et al., Phys. Rev. X 7, 041069 (2017).
[4] H. C. Po, A. Vishwanath, and H. Watanabe, Nat. Commun. 8, 50 (2017).
[5] T. Misawa and Y. Yamaji, Phys. Rev. Research 4, 023177 (2022).
[6] P. Denton, S. Parke, T. Tao, and X. Zhang, Bull. Am. Math. Soc. 59, 31 (2022).

第74回

日時: 2022年8月9日(火)11:00-
発表者: 乾 幸地（東大)
発表題目: 自動微分を用いた目的の性質を満たすハミルトニアンの逆設計

登録フォーム

物質科学の究極の目的は、望みの性質を持つ物質や材料を自在に作り出すことである。理論研究においては、第一原理計算、経験や直感などに基づいて候補となるモデルを構築し、そこに含まれるパラメタを変化させながら目的の物理量を計算することによって、最適なハミルトニアンを探索することがしばしば行われる。しかしながら、こうしたアプローチでは、質的に新しい物理に到達することは難しい。このような状況を打破するために、目的となる物理的な性質からハミルトニアンを推定するという逆問題を解く形のアプローチが注目され[1]、ベイズ推定や生成モデル、摂動論などを用いた研究が行われている[2-5]が、計算コストや汎化性能、計算可能な物理量に制約があるといった課題が残されている。

そこで我々は機械学習の分野で活用されている自動微分の手法を用いて、目的とする物理的な性質を示すハミルトニアンを自動で構築する汎用的なフレームワークを開発した。発表では、量子異常ホール効果に適用することで、このフレームワークがHaldaneモデルを再発見するだけでなく、その6倍大きな異常ホール効果を示す新しいハミルトニアンを自動で生成できることを紹介する。また、光起電力効果[6]への応用では、太陽光の照射によって電子が非共面的なスピン配置の上を動くことによって、約900A/m2の光電流が発生することを示す[7]。コードはgithubに公開してある[8]。

[1] E. Lander et al., Materials Genome Initiative Strategic Plan (National Science and Technology  Council, 2021).
[2] G. L. W. Hart et al., Nature Materials 4, 391 (2005).
[3] R. Tamura and K. Hukushima, Phys. Rev. B 95, 064407 (2017).
[4] B. Sanchez-Lengeling and A. Aspuru-Guzik, Science, 361, 6400 (2018).
[5] H. Fujita et al., Phys. Rev. B 97, 075114 (2018).
[6] R. C. Miller, Phys. Rev. 134, A1313 (1964).
[7] K. Inui and Y. Motome, arXiv:2203.07157 (2022).
[8] https://github.com/koji-inui/automatic-hamiltonian-design

第73回

日時: 2022年7月26日(火)11:00-
発表者: 竹森 那由多（阪大)
発表題目: 動的平均場近似を用いた超伝導ハイパーマテリアルの解析

登録フォーム

準結晶における相関効果は、強相関準結晶物質AuAlYbの合成および準結晶と局所的な構造が類似している近似結晶においては現れない量子臨界性の発見[1]により脚光を浴びている。また、2018年に、Al-Mg-Zn準結晶でバルクの超伝導が発現することが発見[2]され、電子比熱の温度依存性から、電子間に働く有効相互作用は弱結合的であると結論付けられた。しかしながら、並進対称性のない系における弱結合超伝導は非自明な問題であり、超伝導状態の発現機構および特性の解明が求められている。

本発表では、準結晶の構造モデルであるペンローズタイリングを取り上げ、電子間引力相互作用を仮定した引力ハバード模型を用いて、実空間動的平均場近似やBdG方程式を用いたs波超伝導状態の解析の結果を紹介する。超伝導転移点はすべての格子点で唯一に決まる一方で、秩序変数が空間変化する超伝導状態が発現する。特に、弱結合領域においては、各々のクーパー対の重心運動量が0とならない非BCS的な弱結合超伝導状態であること[3]を明らかにした。さらに、実験観測量にもBCS超伝導体との違いが現れる[4]。比熱の跳びはBCS理論値より10-20%小さくなるという結果はAl-Mg-Zn準結晶において観測された実験値と矛盾せず、非一様な秩序変数の構造を反映していると考えられる。この結果は、周期系の超伝導と準結晶の超伝導を区別する指標となりうるため、今後の実験測定が期待される。以上の結果から、準結晶と周期近似を施した近似結晶では、比熱の跳びに定量的な違いが現れることが推測される。同様の解析により、ペンローズ近似結晶の比熱を計算し、ペンローズ準結晶や正方格子の場合の結果と比較する。

[1] K. Deguchi, et al., Nat. Mater. 11 1013-1016 (2012).
[2] K. Kamiya, et al., Nature Communications, 9 154 (2018).
[3] S. Sakai, N. Takemori, A. Koga, and R. Arita, Phys. Rev. B 95 024509 (2017).
[4] N. Takemori, R. Arita and S. Sakai, Phys. Rev. B 102 115108 (2020).

第72回

日時: 2022年7月19日(火)11:00-
発表者: 荻野 卓啓（東大物性研)
発表題目: スピン1/2XXZ梯子系における量子相転移と相構造の解明

登録フォーム

正方格子上のJ1-J2模型は、フラストレート量子磁性体の原型として精力的に研究が行われてきた。 交替Neel相とストライプNeel相の間にスピン液体相やvalence bond solid (VBS)相 などが現れることが提案された。 正方格子上のJ1-J2模型から一次元的に切り出したフラストレートXXZ梯子模型は、 類似性を保ちながら、場の理論や密度行列繰り込み群を活用し詳細に研究することのできる模型として興味深い。 特に、等方的な模型(Δ=1)では、 繰り込み群の流れの中でダイマー間の実効的引力相互作用が発生し、 Haldane相とrung singlet (RS)相の間にcolumnar dimer (CD)相が現れることが、 ボソン化解析と数値計算により明らかになった[1,2]。

我々は、四体相互作用のあるXXZ梯子模型[3,4]、および、フラストレートXXZ梯子模型[5]の基底状態相図を解析した。 四体相互作用を直接入れることで、ダイマー間相互作用の効果をより明確にすることができる。 小さなXXZ異方性(Δ～1)と四体相互作用が組み合わされることで、種々の秩序相と featulessギャップ相 (RS、Haldane、RS、Haldane)からなる豊かな相図が得られることを明らかにした。 ここで、RS* (Haldane*)相は、通常のRS (Haldane)状態と一方の鎖のスピンのz軸周りπ回転によって結びつくものである。 featuless相間、秩序相間ではc=1ガウシアン転移が、 featuless相と秩序相の間ではc=1/2イジング転移が起こり、それらの転移線がΔ=1で交差することを示した。 フラストレートXXZ梯子模型でも、同様の相構造を得た。 以上の結果は、ボソン化による有効場の理論と行列積状態を用いた数値計算によって得られたものである。 特に、4つのfeatureless相がトポロジカル指数によって明確に区別されることを数値的に実証した。

時間が許せば、フラストレートXXZ梯子模型について、 上述の等方的な場合の近傍(Δ～1)だけでなく、 容易面領域(Δ～0)におけるギャップ相とギャップレス相のなす相構造についても解説する。

[References]
[1] O. A. Starykh and L. Balents, Phys. Rev. Lett. 93, 127202 (2004).
[2] T. Hikihara and O. A. Starykh, Phys. Rev. B 81, 064432 (2010).
[3] T. Ogino, S. Furukawa, R. Kaneko, S. Morita, and N. Kawashima, Phys. Rev. B 104, 075135 (2021).
[4] T. Ogino, R. Kaneko, S. Morita, S. Furukawa, and N. Kawashima, Phys. Rev. B 103, 085117 (2021).
[5] T. Ogino, R. Kaneko, S. Morita, and S. Furukawa, arXiv:2206.13048.

第71回

日時: 2022年7月12日(火)11:00-
発表者: 上田 篤（東大)
発表題目: テンソルネットワークによる二次元古典模型の繰り込みフローの視覚化

登録フォーム

近年、モンテカルロ法と肩を並べる量子系/古典系を調べる数値的手段としてテンソルネットワークが人気です。 特に量子一次元系においては、共形場理論との相性の良さから臨界現象やトポロジカル相の分類に広く貢献してきました。[1-2] 同様に、二次元古典臨界模型においても、テンソル繰り込みによって得られたfixed point tensorを用いてスケーリング次元やセントラルチャージが計算できることが示され、”古い問題を解く新しい手段”として注目されてきました。[3]

本公演では、臨界点周りの相図である繰り込み群のフロー(RG flow)もテンソル繰り込みを用いて定量的に計算できることを示します。[4] 実例として、古典XY模型から実際に各スケールでの場の理論における結合定数を計算しそれがKosterlitzのRG flowに一致することをみせます。 時間があれば、より複雑な古典二次元系であるRP^2模型への応用などを紹介します。[5]

[1] P. Calabrese and J. Cardy, J. Phys. A: Math. Theor. 42, 504005 (2009).
[2] F. Pollmann et al, Phys. Rev. B 81, 064439 (2010).
[3] ZC. Gu and XG. Wen, Phys. Rev. B 80, 155131 (2009).
[4] A. Ueda and M. Oshikawa, Phys. Rev. B 104, 165132 (2021).
[5] A. Ueda and M. Oshikawa, arXiv:2202.07042 (2022).

第70回

日時: 2022年6月28日(火)11:00-
発表者: 森田　克洋（東京理科大)
発表題目: フラストレート量子スピン系の有限温度特性と磁気熱量効果

登録フォーム

これまで多くの磁性体の強磁場磁化測定が精力的に行われてきた。理論ではT＝０の磁化曲線の研究が主流であるが、実験では必ずしも十分低温で実施されているわけではない。そのため、有限温度下での磁化過程の理論的研究は重要である。また、近年、数百テスラを超えるパルス磁場発生装置による磁性体の磁化測定も実施されている[1]。そのパルス幅は数μ秒であり、その磁化過程は等温過程ではなく、断熱過程だと思われている。しかしながら、理論的な断熱磁化過程の研究はまだ十分にされているとは言えない状況である。さらに近年、実験的に磁性体の高磁場下における磁気熱量効果を求めることにも成功している。この磁気熱量効果を詳細に調べることにより磁気相転移を検知できることが実験的に明らかになっている[2]。

以上から我々はモデル物質が多数存在する二等辺三角格子S=1/2 Heisenberg模型の等温及び断熱磁化過程の調査を行った[3]。今回我々が用いた手法は有限温度ランチョス法である。その結果、等温磁化過程では1/3プラトーが確認されたが、断熱過程では1/3プラトーは傾斜を持つことが明らかになった。また磁気熱量効果も同時に調べた。これらの結果は今後の実験と比較することが可能であると考えている。

講演では有限温度ランチョス法の原理やその改良方法、本模型の等温及び断熱磁化曲線、磁気熱量効果について詳しく述べる。また時間があれば同じく有限温度ランチョス法を用いた研究であるカゴメ格子（及び三角格子）Kitaev-Heisenberg模型の磁気特性についても紹介したい[4]。

[1] N. Miura, T. Osada, and S. Takeyama, J. Low Temp. Phys. 133, 139 (2003).
[2] Y. Kohama, et al., Proc. Natl. Acad. Sci. USA 116, 10686 (2019).
[3] K. Morita, Phys. Rev. B 105, 064428 (2022).
[4] K. Morita and T. Tohyama, Phys. Rev. Research 2, 013205 (2020).

第69回

日時: 2022年6月7日(火)11:00-
発表者: 大同 暁人（京大)
発表題目: 内因的超伝導ダイオード効果

登録フォーム

近年，物質の新しい機能性として非相反な輸送現象が注目を集めている[1,2]．非相反性とは右に動くものと左に動くものが異なる振る舞いを示すことであり，例えば電気抵抗の非相反性が半導体から超伝導体にいたるまでまで様々な物質で観測されている．一方，超伝導体では超伝導性によって一方向の電気抵抗がゼロになり，もう一方は有限に留まるといった状況が生じうる．このような超伝導ダイオード効果(SDE)が，反転中心のないNb/V/Ta超格子において最近観測された[3]．SDEは省エネルギー電気回路にも応用が期待できるほか，基礎物理学的にも興味深い現象である．しかしSDEが発生するメカニズムは解明されておらず，理論提案が待たれていた．

本講演ではバルク物質のSDEに関する理論研究の第一歩として，対破壊機構による臨界電流の非相反性を議論し，SDEを引き起こす内因的なメカニズムを提案する[4]．平均場近似によるミクロな計算により，超伝導転移温度近傍でのダイオード効果の温度スケーリングや，低温での著しい増大が明らかになった．特に低温高磁場での振る舞いは反転対称性の破れた超伝導体で提案されていたエキゾチック超伝導相（ヘリカル超伝導）と深く関連しており，超伝導ダイオード効果がヘリカル超伝導のバルクプローブとして重要な役割を果たしうることを議論する．

[1] Y. Tokura and N. Nagaosa, Nat. Commun. 9, 3740 (2018). [2] T. Ideue and Y. Iwasa, Annu. Rev. Condens. Matter Phys. 12, 201 (2021). [3] F. Ando et al., Nature 584, 373 (2020). [4] A. Daido, Y. Ikeda, and Y. Yanase, Phys. Rev. Lett. 128, 037001 (2022).

第68回

日時: 2022年5月24日(火)11:00-
発表者: 中 惇（東京電機大学)
発表題目: 有機反強磁性体のスピン・電荷輸送

登録フォーム

近年、スピン流やベリー位相に立脚した新しい量子輸送現象の開拓が大きく進展している。スピンホール効果[1]と異常ホール効果[2]はこの代表例であり、これらは強い相対論的効果や強磁性磁化、スピンカイラリティなどをその起源とすることから、研究の舞台はスピン軌道結合が強い重金属や半導体、多彩な磁気構造を示す無機磁性体に限られてきた。これに対して本研究では、スピン軌道結合が小さい典型的な有機反強磁性体であるκ-(ET)2Xにおいて、これらと本質的に異なる新しいスピン・電荷輸送現象が発現することを見出した。

κ-(ET)2Xは、ETと呼ばれる分子が配向の異なる2種類のダイマーを形成し、三角格子上に配列した構造を持つ。ダイマー上のホール間に働く強いクーロン相互作用によりモット転移や反強磁性秩序を示す典型的な強相関系である。これらの性質は古くから分子配向を無視した単純なハバードモデルでよく説明できるとされてきた[3]。これに対して我々は、この系の電子状態を分子配向に着目して再考察することで、反強磁性状態のバンドに特異なスピン分裂が存在すること、さらに電場や熱勾配によりスピン流生成が可能となることを明らかにした[4]。また、スピン軌道結合を取り入れた解析により異常ホール効果が発現することを見出した[5]。これは磁化やスピンカイラリティを必要としない新しいタイプの異常ホール効果である。

[1] S. Murakami, N. Nagaosa, and S. C. Zhang, Science 301, 1348 (2003).
[2] K. Ohgushi, S. Murakami, and N. Nagaosa, Phys. Rev. B 62, R6065 (2000).
[3] H. Kino and H. Fukuyama, J. Phys. Soc. Jpn. 65, 2158 (1996).
[4] M. Naka et al., Nat. Commun. 10, 4305 (2019).
[5] M Naka et al., Phys. Rev. B 102, 075112 (2020).

第67回

日時: 2022年5月10日(火)11:00-
発表者: 曽根 和樹（東大)
発表題目: 非線形振動子のトポロジカルな同期現象

登録フォーム

トポロジカル絶縁体などで議論されるバンド構造のトポロジーは、物質のロバストな性質を実現するための指針として多く研究されている。従来このようなバンドトポロジーは線形系のハミルトニアンに基づいて議論されてきたが、近年の研究で非線形系においてもその対応物が存在することが明らかになってきた[1]。非線形方程式は生物から量子に至る様々な系を（実効的に）記述するため、非線形系へのトポロジーの拡張は基礎理論だけでなく応用の観点からも重要と考えられる。しかしながら、どのようなクラスの非線形系まで拡張が可能か、拡張した場合のバルクエッジ対応はどのように定義されるかといった問題が未解明である。

本講演では我々の論文[2]に基づき、散逸のある非線形系にトポロジーの概念を拡張することで得られたトポロジカル同期について紹介する。トポロジカル同期状態では、バルクの非線形振動子はカオスな振る舞いを見せる一方で、エッジに存在する振動子は同期し周波数を揃えて振動する。これは非線形振動子で従来研究されてきたものとは異なり、振動子間の結合のトポロジーに対応して出現する局在同期状態であり、非線形系へのエッジモードの拡張となっている。これらの結果は非線形系をロバストに制御する手法を与えるとともに、非線形トポロジーの基礎理論を解明するための足掛かりになると期待される。

[1] D. Smirnova, D. Leykam, Y. Chong, and Y. Kivshar, Appl. Phys. Rev. 7, 021306 (2020).
[2] K. Sone, Y. Ashida, and T. Sagawa, arXiv: 2012.09479 (2020).

第66回

日時: 2022年4月19日(火)11:00-
発表者: 増木 亮太（東大)
発表題目: 非調和フォノン理論に基づく有限温度における構造最適化

登録フォーム

結晶構造は物質の性質や対称性を決める基本的な要素であり、非経験的にこれを決定することは計算物性物理学の重要課題の一つである。結晶構造の第一原理計算においては、密度汎関数法(DFT)による構造最適化が最も頻繁に用いられる。これは原子核を古典粒子とみなした基底状態エネルギーの最適化であり、原子核の量子性(フォノンの零点振動)と格子振動による有限温度効果が無視されている。一方で、熱膨張や構造最適化に見てとれるように、物質の結晶構造においてこれらの効果は無視できない影響を持つ。

近年、計算機性能の向上により、従来は困難だったフォノン間相互作用の第一原理計算の研究が盛んに行われており、フォノン振動数の温度依存性や熱伝導率などが標準的に計算されるようになりつつある[1,2]。しかし、その結晶構造への影響を調べることは容易ではなく、多くの物質に適用できるような効率的な手法は存在しなかった。

そこで、本研究では、逆空間自己無撞着フォノン理論(momentum-space self-consistent phonon theory)[3]を拡張し、変分自由エネルギーの最適化によって有限温度の結晶構造を効率的に計算する手法開発を行った。開発した手法は、非調和フォノン計算の公開パッケージであるALAMODE[3,4]に実装した。我々はこれを代表的な強誘電体であるBaTiO3に適用し、その3段階の構造相転移を再現することに成功した。

[1] T. Tadano and S. Tsuneyuki, J. Phys. Soc. Jpn. 87, 041015 (2018).
[2] L. Monacelli et al., J. Phys.: Condens. Matter 33, 363001 (2021).
[3] T. Tadano and S. Tsuneyuki, Phys. Rev. B 92, 054301 (2015).
[4] Y. Oba et al., Phys. Rev. Mater. 3, 033601 (2019).

第65回

日時: 2022年4月12日(火)11:00-
発表者: 安田 憲司（MIT）
発表題目: 積層制御による二次元強誘電体の設計

登録フォーム

ファンデルワールス層状物質からなる二次元強誘電体は、原子層レベルまで強誘電性が保持されるため、強誘電メモリへの応用の観点から注目を集めている。従来、これらの二次元強誘電体はバルクの層状化合物を剥離することによるトップダウンの手法によって得られてきた。一方、これらの本手法では極性の層状結晶が必要となるため、得られる強誘電体には大きく制限が課されていた。本講演では、近年発展してきたファンデルワールス積層技術を応用することで、結晶の対称性を制御し、強誘電性を持たない母物質から強誘電体を設計できることを示す。 バルクの六方晶窒化ホウ素(BN)は隣接する層が反平行した積層構造のため、バルクBNは空間反転対称性を有する。一方、２枚のBNを回転せずに積層すると、結晶は強誘電体となる[1]。平行積層されたBNは、ABまたはBA stackingとなり、面間のスライディングに伴う面直分極反転が生じる。本研究では、単層BNを切断、積層することで平行積層二層BNを作製し、圧電応答力顕微鏡並びにグラフェンをプローブとして分極を測定した[2]。本設計指針は幅広い二次元物質に適用可能であり、二次元強誘電体の種類を大幅に拡大する可能性を有している。実際我々は、同様のアイディアを遷移金属ダイカルコゲナイド(TMD)にも適用することで、MoS2, MoSe2, WS2, WSe2という4つの半導体TMDにおいても強誘電性を実現した[3]。２枚のBNに小角捻りを加えることで、分極が周期的に配置したモアレ強誘電体が生じること、またその基板としての利用についても紹介したい。

[1] L. Li, M. Wu ACS Nano. 11, 6382–6388 (2017).　[2] K. Yasuda et al. Science 372, 1458–1462 (2021).　[3] X. Wang, K. Yasuda et al. Nature Nanotech. (2022).

第64回

日時: 2022年3月1日(火)11:00-
発表者: 山田 昌彦（阪大)
発表題目: ハニカム格子上のトポロジカルZ_4スピン液体

登録フォーム

新奇な量子物質を発見し特定することは昨今大変注目を集めている。 我々は今回創発SU(4)対称性[1,2]という、SU(2)よりも遥に高い対称性を 実現することが期待されている、スピン軌道相互作用の強いd^1物質に注目した。 実際、α-ZrCl_3は強いフラストレーションを持つ、ハニカム格子上の SU(4)ハイゼンベルク模型[3]を実現することが予想され、これらのSU(4)物質は ギャップトZ_4スピン液体[4]という新しいタイプのスピン液体を実現する可能性が ある。我々は新たにこのハニカム格子上SU(4)ハイゼンベルク模型に対して 密度行列くりこみ群(DMRG)の研究を行い、基底状態が実際にギャップトで全ての 対称性を保持することを示し、この系がギャップトスピン液体状態にあることを 示唆する結果を得た[5]。トポロジカルエンタングルメントエントロピーは Z_4スピン液体に期待される値に近い値をとり、この状態は熱力学極限で16の 縮退基底状態を持つことが期待される。

[1] M.G.Y., M. Oshikawa, and G. Jackeli, Phys. Rev. Lett. 121, 097201 (2018). [2] M.G.Y., M. Oshikawa, and G. Jackeli, Phys. Rev. B 104, 224436 (2021). [3] P. Corboz et al., Phys. Rev. X 2, 041013 (2012). [4] M.G.Y., and S. Fujimoto, arXiv:2111.14470 (2021). [5] M.G.Y., K. Penc, and F. Pollmann, to appear.

第63回

日時: 2022年2月8日(火)11:00-
発表者: 南谷 英美（分子研)
発表題目: 第一原理計算によるエネルギー損失過程の解析

登録フォーム

電子とフォノンの相互作用は電気抵抗、ポーラロン形成、電荷密度波、熱電効果から超伝導に至るまで固体物理の随所にあらわれる最も基本的な相互作用の一つである。電子とフォノンが相互作用する仕組み自体は、格子が揺れると、電子が感じるポテンシャルが変わり、それによって電子状態が変調することで直感的に理解できる。しかし、相互作用強度の定量的な見積もりは容易ではない。それを可能にする一つの方法が、密度汎関数摂動理論に基づく第一原理計算を応用した電子フォノン相互作用の評価である。近年の計算機能力や、各種計算パッケージの発達により、今では各フォノンモードが各電子状態と相互作用する強度を直接的に見積もることができるようになっている[1,2]。 電子フォノン相互作用を精密に計算することは、主に超伝導転移温度の推定に使われているが、それ以外にも固体中での電子のエネルギー損失の定量化に直結する。本講演ではその例として、非弾性電子トンネル過程や、半導体におけるジュール熱の素過程についての第一原理解析を紹介する[3,4]。

[1] S. Baroni et al., Rev. Mod. Phys. 73, 515 (2001). [2] F. Giustino, Rev. Mod. Phys. 89, 015003 (2017). [3] E. Minamitani et al., Phys. Rev. B 94, 085411 (2016). [4] E. Minamitani, Phys. Rev. B 104, 085101 (2021).

第62回

日時: 2022年1月25日(火)11:00- ※講演者の都合により延期になりました
発表者: 幸城 秀彦（東大)
発表題目: 多重磁化プラトーを実現するフラストレートスピンラダーの理論的解明

登録フォーム

第61回

日時: 2022年1月18日(火)11:00-
発表者: 吉永 敦紀（東大)
発表題目: 弱横磁場下のイジング模型におけるHilbert Space Fragmentation

登録フォーム

量子系が熱浴から孤立していても、多体系で物理量が熱平衡値へ緩和することは古くから議論され、冷却原子系などを用いた多くの研究もこれを示唆している。この孤立系の熱平衡化を説明する試みとして固有状態熱化仮説が提案され、可積分性や強い乱雑性を持たない模型でその有効性が広く検証されている。一方で近年、初期状態に依存して長時間の振動を示す量子系[1]の発見に伴い、固有状態熱化仮説を破る非可積分系が注目され、特にHilbert Space Fragmentation (HSF)と呼ばれる機構[2]が注目を集めている。HSFは、ダイナミクスに対する制約条件が原因で、保存量を指定した部分空間内であってもHilbert空間がさらに指数関数的な個数の部分空間に分離する現象である。本研究では、2次元以上の正方格子上の横磁場イジング模型において、その弱磁場下の有効模型が HSF を示すことを発見した[3]。講演では、このHSFがドメイン壁の保存に起因することを説明し、またクエンチ後の時間発展で予想される初期状態依存の多様なダイナミクスについても紹介する。 [1] H. Bernien, et al., Nature 551, 579 (2017). [2] P. Sala, et al., Phys. Rev. X 10, 011047 (2020). [3] A. Yoshinaga, H. Hakoshima, T. Imoto, Y. Matsuzaki, R. Hamazaki, arXiv:2111.05586

第60回

日時: 2021年12月14日(火)11:00-
発表者: 播木 敦（大阪府立大)
発表題目: 内殻X線分光を用いた超伝導層状ニッケル酸化物の電子状態の理論研究

登録フォーム

2019年に層状Ni酸化物(Nd0.8Sr0.2NiO2)で超伝導が発見され、現在盛んに議論されている。高温超伝導銅酸化物との違いや類似性を理解するために、Niイオンの価数や正孔が注入される軌道のキャラクタの決定が望まれる。本研究では、最近報告されたNi 2p内殻X線光電子分光(XPS)、内殻X線吸収分光(XAS)、共鳴非弾性X線散乱(RIXS)の実験スペクトル[1][2]を密度汎関数理論(DFT)＋動的平均理論(DMFT)[3]を用いて解析した結果を報告する [4]。この手法は内殻スペクトルの解析に標準的に用いられるクラスタ模型では考慮できない価電子帯との混成効果や正孔ドープに伴う価電子帯の再構成を考慮した解析が可能である。Ni 2p XPSの解析から、電荷移動エネルギーΔなどの固体パラメータを決定し、Niの価数がほぼ1価であることがわかった。この結果は高温超伝導銅酸化物と同様の3d9基底状態が実現していることを示唆するが、Δの値は銅酸化物より2eV程度大きく、よりMott-Hubbard型の電子構造を有するという違いがある。今回のセミナーでは、相補的な情報をもたらすXASやRIXS過程の解析結果も示し、層状Ni酸化物Nd1-xSrxNiO2の電子状態を多角的に議論する。

[1] Y. Fu et al., arXiv:1911.03177 (2019)
[2] M. Hepting et al., Nat. Mat. 19, 381 (2019)
[3] A. Hariki et al., Phys. Rev. Lett. 121,126403 (2018)
[4] K. Higashi et al., Phys. Rev. X 11, 041009 (2021)

第59回

日時: 2021年11月30日(火)11:00-
発表者: 水野　竜太（京大）
発表題目: DMFTの不純物問題におけるフルバーテックスの簡略化と非局所相関の取り込みのためのその応用

登録フォーム

動的平均場理論(DMFT)[1]を超えて非局所相関を取り込む拡張のうち、dual fermion method [2]や Dynamical vertex approximation (DΓA) [3]のようなダイアグラム展開の方法を 用いるタイプの拡張において、2 体の局所バーテックスは重要な役割を持つ。しかし、連続時間量子モンテカルロ法(CT-QMC)[4]のような厳密な手法でこれを求めることは大きな計算リソースを必要する困難な作業である。これまで、このような困難を回避・軽減するための様々な努力がなされてきたが [5-7]、2 体の局所バーテックスを求めることは未だ容易ではない。 そこで、本研究では、厳密なフルバーテックスが持つ周波数構造に着目し、2 体の局所フルバーテックスの簡略形を提案した。その簡略化されたフルバーテックスを用いて、 DMFTで得られた自己エネルギーから簡略化形のフルバーテックスを見積もる方法と、 dual fermion method の計算を低い数値コストで行うための方法を開発した。本公演では、フルバーテックスの簡略化、および、開発した2 つの手法の詳細を説明し、いくつかのモデルへの適用結果を紹介する。

Reference
[1] A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg: Rev. Mod. Phys. 68, 13 (1996).
[2] J. Otsuki, H. Hafermann, and A. I. Lichtenstein, Phys. Rev. B, 90, 235132 (2014).
[3] Motoharu Kitatani, Thomas Schäfer, Hideo Aoki, and Karsten Held, Phys. Rev. B, 99, 041115 (2019).
[4] A. N. Rubtsov, V. V. Savkin, and A. I. Lichtenstein: Phys. Rev. B 72, 035122 (2005).
[5] H. Shinaoka, D. Geffroy, M. Wallerberger, J. Otsuki, K. Yoshimi, E. Gull, and J. Kunes, SciPost Phys., 8, 12 (2020).
[6] J. Kaufmann, P. Gunacker, and K. Held, Phys. Rev. B, 96, 035114 (2017).
[7] N. Wentzell, G. Li, A. Tagliavini, C. Taranto, G. Rohringer, K. Held, A. Toschi, and S. Andergassen, Phys. Rev. B, 102, 085106 (2020).
[8] R, Mizuno, M. Ochi, K. Kuroki: arXiv:2101.04466.

第58回

日時: 2021年11月16日(火)11:00-
発表者: 田財 里奈（名古屋大）
発表題目: 汎関数くりこみ群/密度波方程式による非局所相転移の理論

登録フォーム

電子相関に由来する自発的対称性の破れは、長年の中心課題である。最近そのバラエティーが拡大し、多体効果によって電子の飛び移り積分が変調する「ボンド秩序」や、古典には実現しない永久自発カレントが生まれる「電流秩序」が、銅酸化物、Ir系、フラストレートV化合物[1]において相次いで報告された。これらの新規相転移の特徴は、その秩序変数が、複数の格子点にまたがって定義されるという“非局所性“である。しかし、「これらの非局所転移が、なぜ通常のスピン秩序や電荷秩序よりも安定化しやすいのか？」や、「非局所秩序が発現する条件は？」といった根本的な疑問は未解明のままである。

本セミナーでは初めに、非局所秩序を統一的に記述するため、構造因子という関数を導入する。これによって対称性に注目した見通しの良い議論が可能となる。次に、多体効果を考慮して構造因子を決定する理論手法として、(i)変分法と併せた汎関数くりこみ群法、 (ii) Luttinger-Ward理論に基づいて熱力学ポテンシャルを最小化する密度波方程式、という２つの手法を紹介する。これらの手法を用いて得られた、幾何学フラストレーションに誘起される電荷ループ秩序の理論、さらに最近の話題として、新規カゴメ金属AV3Sb5 (A=Cs, Rb, K)におけるDavid-Star型ボンド秩序の機構について紹介する[2]。特に、AV3Sb5の超伝導相の転移温度は、電荷秩序の量子臨界点近傍で最大となり、ボンド秩序の量子揺らぎが媒介する新規超伝導相が発現している可能性がある[3]。本研究では、ボンド揺らぎによって、s波超伝導や、トリプレットp 波超伝導が実現する可能性を指摘したい。

[1] B. R. Ortiz et al., Phys. Rev. Materials 3, 094407 (2019).
[2] R. Tazai, et al.,arXiv:2107.05372 (2021).
[3] K. Y. Chen, et al.,Phys. Rev. Lett. 126, 247001 (2021).

第57回

日時: 2021年11月2日(火)11:00-
発表者: 近藤 寛記（東大）
発表題目: 対称性に保護されたマグノンのディラック表面状態と反強磁性体におけるマグノンの非線形スピンネルンスト効果

登録フォーム

トポロジカル物性の研究はマグノン系でも進んでおり、量子ホール絶縁体[1]や量子スピンホール絶縁体[2]の対応物などが調べられている。3次元トポロジカル絶縁体などがもつディラック表面状態に対応するマグノン状態では多くの新奇物性が期待されるが、電子系の場合でクラマース対としてトポロジカルに保護された表面状態を生じていた対称性を、マグノン系でどう考えればよいかは明らかではなかった。 そこで我々は、時間反転（スピンの向きの反転）と半並進を組み合わせた操作に関する対称性がマグノンのクラマース対を生じることに着目し、この対称性に保護されたマグノンのディラック表面状態を有する磁性体のモデルを構築した[3]。この系において、電場をかけたときのジャロシンスキー・守谷相互作用の変化により、向かい合った2つの面のディラック分散はそれぞれ上昇/下降するが、それに伴うマグノン流も議論する。上記の表面状態は、ファンデルワールス磁性体CrI3において現れることが期待される。

また、マグノンのトポロジカルな非線形応答に関する研究も紹介する。近年電子系において、電場に関する2次のホールカレントがベリー曲率の双極子モーメント（ベリー曲率双極子）で記述されることが明らかになっている[4]。我々はマグノンの温度勾配に関する2次の応答を考え、ベリー曲率双極子のような物理量を含んだ温度勾配とは垂直なカレントの表式を導いた[5]。特にネール反強磁性体では、PT対称性によって温度勾配とは垂直なカレントは上/下スピン由来のマグノンで逆方向に流れるため、期待できるのはスピン流である。このような非線形スピンネルンスト効果は、中心対称性を破った多くの単純なネール反強磁性体で起きることを議論する。

[1] H. Katsura, N. Nagaosa, and P. A. Lee, Phys. Rev. Lett. 104, 066403 (2010). [2] V. A. Zyuzin and A. A. Kovalev, Phys. Rev. Lett. 117, 217203 (2016). [3] H. Kondo and Y. Akagi, arXiv:2012.02034. (Accepted in Physical Review Letters) [4] I. Sodemann and L. Fu, Phys. Rev. Lett. 115, 216806 (2015). [5] H. Kondo and Y. Akagi, arXiv:2109.09464.

第56回

日時: 2021年10月19日(火)11:00-
発表者: 高橋　雅大（阪大）
発表題目: キタエフ磁性体におけるトポロジカルネマティック相転移と点欠陥に束縛されたマヨラナゼロモードの確立

登録フォーム

強い量子ゆらぎにより絶対零度極限においても長距離磁気秩序を有しない物質は量子スピン液体と呼ばれ、様々な模型やその候補磁性体において精力的に研究されている。二次元蜂の巣格子の各頂点にスピン1/2 が配置され、ボンド方向によって三種類のイジング相互作用が働くキタエフ模型[1]はその代表例であり、本模型は驚くべきことに、各スピン自由度を遍歴マヨラナ準粒子と局所Z2ゲージ場に分数化することで厳密に解くことができる。また、スピン軌道相互作用の強い遷移金属を有するキタエフ候補物質群の発見(例えば[2])により、さまざまなスピン相互作用を含んだキタエフ磁性体に関する研究は、新奇なマヨラナ物理現象を発現させる舞台として、さらには誤り耐性のあるトポロジカル量子計算への応用も視野に、非常に注目を集めている。

本講演では、キタエフ模型とキタエフ磁性体について述べた後、遍歴マヨラナ準粒子の多体相関によって引き起こされる一次相転移[3]について説明する。この相転移により、系は回転対称性の破れを伴いながらトポロジカル相を変化させることがわかっており、このトポロジカルネマティック相転移と近年観測されたα-RuCl3でのネマティック相転移[4]との対応についても触れていきたい。またもう一つ別の話として、最近の研究からキタエフ模型の点欠陥に束縛されたマヨラナゼロモードと、その存在を反映すると期待できる非局所なスピン-スピン相関についても併せて発表したい。

[1] A. Kitaev, Ann. Phys. 321, 2 (2006).
[2] G. Jackeli and G. Khaliullin, Phys. Rev. Lett. 102, 017205 (2009).
[3] M.O.T. et al, Phys. Rev. Research 3, 023189 (2021).
[4] O. Tanaka et al., arXiv:2007.06757 (2020).

第55回

日時: 2021年10月5日(火)11:00-
発表者: 伏屋 雄紀（電通大基盤理工）
発表題目: ナノスケールのチューリング・パターン：単層ビスマスの形態形成

登録フォーム

生物における形態形成の過程では，無秩序から秩序が生まれる．チューリングは，異なる拡散速度を持つ活性因子と抑制因子間の相互作用により，形態形成の過程が簡単な連立微分方程式（反応拡散方程式）で記述できることを発見した[1]．チューリングの理論は，生物学や化学など，幅広い分野で極めて大きな影響を及ぼした．生物パターンの典型的な長さスケールはcmからmmで，純粋な化学系ではsub-mm程度である．しかし，μm以下のチューリング・パターンは自然界においては極めて稀であり，ほとんど調べられてこなかった．また，チューリング・パターンが生成される舞台はソフトマターに限られていた．

本研究では，NbSe2上に成長したビスマス単原子膜の形態形成を理論的に調べた．我々は，三種の原子間相互作用からなる模型を構築し，その時間発展方程式を解くことで，ビスマス原子位置の安定解を求めた．その結果，原子５個分（約2nm）の周期を持つチューリング・パターンが生成されることが分かった．計算によって得られたストライプ構造，ドメイン壁およびそれらのY字型接合は，何れも実験で観測されたパターンと非常によい一致を示す[2,3]．我々の結果は，“ハード”マターにおいて原子スケールのチューリング・パターンが得られた最初の実例となる．

[1] A. Turing, Phil. Trans. R. Soc. Lond. B, 237, 37 (1952).
[2] A. Fang, et al., Sci. Adv., 4, eaaq0330 (2018).
[3] Y. Fuseya, H. Katsuno, K. Behnia, A. Kapitulnik, Nature Physics, (2021), https://doi.org/10.1038/s41567-021-01288-y

第54回

日時: 2021年9月28日(火)11:00-
発表者: 秦 徳郎（東工大）
発表題目: 量子液体における三体相関の実験的検出

登録フォーム

量子ドットにおける近藤効果は、量子液体における平衡状態から非平衡状態までの連続的な振る舞いを扱う最適な舞台である[1]。その線形応答領域は、平衡状態における電子間の二体相関に基づく揺動散逸関係で記述される。この「二体相関」は帯磁率に相当する量であり、これまでの研究によって実験的にも理論的にも確立されている。しかし、非線形非平衡領域を扱うには、「三体相関」といった高次の相関の役割を明らかにする必要がある。

本研究では、ユニタリ極限に近い理想的な近藤効果を実現し、その非線形非平衡領域を解析することで三体相関を評価した[2]。磁場が印加されていない状態では三体相関はゼロで、磁場を印加すると三体相関が増大することがわかった。この結果は、近年発表された非平衡領域にある量子液体の理論（局所フェルミ液体論）[3]に定量的に合致した。本成果は、非平衡状態における量子液体において、実験的に初めて三体相関を検出することに成功したものであり、他の非平衡現象を定量的に理解するためのよい試金石となる。

[1] D. Goldhaber-Gordon et al., Nature 391, 156–159 (1998). [2] T. Hata et al., Nat. Commun. 12, 3233 (2021). [3] A. Oguri and A. C. Hewson, Phys. Rev. Lett. 120, 126802 (2018).

第53回

日時: 2021年9月7日(火)11:00-
発表者: 吉川 尚孝（東大）
発表題目: 電荷密度波の振幅モード励起による絶縁体的状態

登録フォーム

光によって物質の電気伝導性や磁性などの性質を自在に制御する試みは兼ねてから盛んに研究されてきた。最近、テラヘルツ(THz)波や中赤外光などの光を用いることで、熱化を抑えながら特定のモードを励起することが可能となり、これによる低温相や隠れた秩序相への相転移現象が報告されている[1]。我々は秩序相を示す物質の秩序変数を光で直接操作することによる物性制御の実現を目的とし、テラヘルツ周波数帯に電荷密度波の秩序変数の振幅モードを示す２次元物質3R-Ta1+xSe2に着目して高強度テラヘルツ波を用いたポンププローブ分光を行った。高強度THz波を用いてこの物質の電荷密度波の振幅モードを二光子過程で駆動することに成功し、さらに振幅モード励起に伴って金属的な応答を示す平衡状態のCDW相から平衡状態にない絶縁体的な状態へ変化することを示した[2]。

[1] X. Li et al., Science 364, 1079 (2019). [2] N. Yoshikawa et al., Nat. Phys. (2021). (https://doi.org/10.1038/s41567-021-01267-3)

第52回

日時: 2021年8月24日(火)11:00-
発表者: 竹内 宏光 (大阪市立大理学部)
発表題目: 量子楕円渦：U(1)対称性が破れた超流動相における巻き数１の量子渦の自発的“分裂“

登録フォーム

自発的対称性の破れが起こる系でどのような形状の位相欠陥が生じるかという問題は その対称性の破れ方を位相幾何学で分類することで理解されてきた． このような位相欠陥に対する従来の認識を否定する物理現象が最近観測された[1]． U(1)対称性が自発的に破れた従来の超流体では量子渦と呼ばれる“糸状“の位相欠陥が実現し， その断面は台風の目の様に回転対称である． ところが，この実験系の超流動相では同様の対称性の破れが起こるにも関わらず， そこで観測された位相欠陥の形状はそれとは似ても似つかないものであった． 本講演ではこの奇妙な位相欠陥が量子楕円渦と呼ばれる回転対称性を破った渦状態として定量的に理解され， “糸状“ではなく“帯状“の渦として存在することを示す[2]．

[1] Seji Kang, Sang Won Seo, Hiromitsu Takeuchi, and Yong-il Shin Observation of Wall-Vortex Composite Defects in a Spinor Bose-Einstein Condensate Phys. Rev. Lett. 122, 095301 (2019) [2] Hiromitsu Takeuchi Quantum Elliptic Vortex in a Nematic-Spin Bose-Einstein Condensate Phys. Rev. Lett. 126, 195302 (2021)

第51回

日時: 2021年8月3日(火)11:00-
発表者: 中村 真（中央大）
発表題目: ゲージ・重力対応を用いた非平衡相転移の解析と電流駆動型非平衡三重臨界点の発見

登録フォーム

ゲージ・重力対応とは、強く相互作用するゲージ粒子の理論を、高次元の曲がっ た時空における一般相対性理論の古典論に置き変える対応関係である。この対応 関係は微視的理論のレベルで成立しており、系のマクロな状態が平衡状態か非平 衡状態であるかは問わない。同時に、対応するゲージ理論を粗視化してマクロな 物理を表現する機構が重力理論側には備わっている。これらを総合すると、非平 衡状態に対応するバックグラウンドを重力理論側で適切に設定することで、非平 衡状態にあるゲージ粒子多体系のマクロな性質を、重力の古典論から読み取るこ とができるものと考えられる。講演者らは、ゲージ・重力対応のこのような性質 を利用し、定常電流の流れる非平衡定常状態において、ゲージ粒子多体系が示す 相転移を調べてきた。特に本研究では、カイラル対称性の自発的破れに対する定 常電流の影響を解析し、定常電流の存在により新奇な三重臨界点が現れる場合が あることを発見した[1]。本講演では、ゲージ・重力対応の概要を説明するとと もに、今回発見した三重臨界点の臨界指数など、詳細を報告する予定である。

[1] T. Imaizumi, M. Matsumoto and Shin Nakamura, ``Current Driven Tricritical Point in Large- Nc Gauge Theory’’ Phys. Rev. Lett. 124 (2020) 19, 191603, [arXiv:1911.06262 [hep-th]]

第50回

日時: 2021年7月20日(火)11:00-
発表者: 吉岡 信行（東大）
発表題目: ニューラルネットワークによる固体量子物性の第一原理計算

登録フォーム

画像認識や自動翻訳などの情報処理タスクにおいて、大きな威力を発揮しているニューラルネットワークは、非常に柔軟な表現能力を持つことが知られている。近年、その表現能力をうまく活用することで、量子多体系をはじめとした、複雑な相関を持つ状態のシミュレーションを効率的に行えることがわかってきた。これは、波動関数法による多体電子状態の第一原理計算のように、変分関数の表現能力に難点を抱えてきた数値手法を大幅に改善させる可能性を示唆している。

本講演では、ニューラルネットワークを用いて、固体系の電子状態を計算する手法を導入する。特に、変分モンテカルロ法を用いた基底状態計算および部分空間展開によるバンド計算を行う手法を議論する。密な原子構造・高い対称性に起因して電子相関が強まり、結合クラスター法などの従来手法が破綻するような領域においても、ニューラルネットワークにより電子状態を精密に表現できることを示す。本講演は論文[1]に基づく。

[1] N. Yoshioka, W. Mizukami, and F. Nori, Communications Physics 4, 106 (2021).

第49回

日時: 2021年7月13日(火)11:00-
発表者: 金杉翔太（京大）
発表題目: 多バンド系におけるアナポール超伝導の理論

登録フォーム

Cooper対が内部自由度を持つ非従来型超伝導は、U(1)ゲージ対称性だけでなく様々な対称性の破れを引き起こす。特に、量子系において重要な時間反転対称性・パリティ対称性を破る場合は活発に研究されている。例えば時間反転対称性を破るカイラル超伝導では、そのトポロジカルな性質が注目を集めている。空間反転対称性の破れた系では、パリティ混成した超伝導が精力的に研究されてきた。さらに最近では、空間反転対称性を保つ系でもパリティ混成超伝導の秩序化が議論されている。奇パリティ多極子秩序の揺らぎが超伝導を媒介する系では、偶パリティと奇パリティの対形成引力の拮抗が指摘されている[1]。近年発見されたスピン三重項超伝導体の候補物質UTe2では、結晶構造は空間反転対称性を保つにも関わらず、圧力下でパリティ混成した超伝導状態の安定化が提案されている[2]。上記のようなパリティの異なる超伝導不安定性が拮抗する系では、時間反転対称性を破るパリティ混成超伝導の安定化が示されているが[3]、その詳しい性質は明らかにされていない。

そこで本講演では、時間反転対称性を破るパリティ混成超伝導の性質を、新たに多バンド性を取り入れて解析した結果[4]を紹介する。まず、超伝導が時間反転対称性とパリティ対称性の両方を破ることに起因して、Bogoliubov準粒子のスペクトルに非対称性が現れることを明らかにする。この非対称性はバンド間ペアリングが担っており、多バンド性を考慮することで初めて見える現象である。さらに超伝導の対称性が極性を持つ場合には、Bogoliubovスペクトルの非対称性が有効的にアナポールモーメントを生じ、その方向に空間変調した超伝導（アナポール超伝導）が実現することを示す。最後に以上の結果をUTe2に適用し、結晶構造の局所的な空間反転対称性の破れに起因してアナポール超伝導が実現することを議論する。

[1] V. Kozii and L. Fu, Phys. Rev. Lett. 115, 207002 (2015). [2] J. Ishizuka and Y. Yanase, Phys. Rev. B 103, 094504 (2021). [3] Y. Wang and L. Fu, Phys. Rev. Lett. 119, 187003 (2017). [4] S. Kanasugi and Y. Yanase, submitted.

第48回

日時: 2021年6月29日(火)11:00-
発表者: 水上 雄太（東大）
発表題目: キタエフ量子スピン液体におけるマヨラナ励起の磁場角度依存性

キタエフ模型は、蜂の巣格子上のS=1/2スピンがボンド方向に依存するイジング相互作用を持つ量子スピン模型である[1]。その相互作用に由来したフラストレーションにより、キタエフ模型の基底状態としては量子スピン液体が厳密解として得られる。この量子スピン液体状態では、スピンが分数化して二種類のマヨラナフェルミオンによる低エネルギー励起が生ずることが理論的に知られている。近年、キタエフ模型が現実の物質で実現すると理論的に予測され[2]、強いスピン軌道相互作用を持つモット絶縁体において盛んに研究が行われている[3,4]。キタエフ量子スピン液体状態においては、系のトポロジカルな性質によりマヨラナフェルミオンによるカイラルエッジ流が生じ、半整数の熱量子ホール効果が期待されているが[4]、近年それがキタエフ物質α-RuCl3において観測されている[5]。一方、そのバルク状態においては、磁場中ではマヨラナフェルミオンの励起スペクトルはギャップΔを持ち、磁場𝐻=(ℎ𝑥,ℎ𝑦,ℎ𝑧)に対してΔ∝|ℎ𝑥ℎ𝑦ℎ𝑧|の形の強い磁場方向依存性を持つことが理論的に予言される。

本講演では、我々が行っているα-RuCl3単結晶試料に対する磁場角度回転比熱測定を紹介する。蜂の巣格子面内において磁場角度を回転させた際の比熱は顕著な六回振動を示す。これは、系の励起ギャップの磁場角度依存性に六回振動が存在していることに起因する。比熱から得られる励起ギャップの磁場角度依存性と磁場強度依存性は、キタエフ模型で予言されるものと非常に良く一致しており、マヨラナ励起のエッジ流に対するバルク状態を観測していると考えられる[6]。

[1] A. Kitaev, Ann. Phys. 321, 2 (2006). [2] G. Jackeli and G. Khaliullin, Phys. Rev. Lett. 102, 017205 (2009). [3] H. Takagi et al., Nat. Rev. Phys. 1, 264 (2019). [4] Y. Motome, and J. Nasu, J. Phys. Soc. Jpn. 89, 012002 (2020), 第47回講演 [5] Y. Kasahara et al., Nature 559, 227 (2018). [6] O. Tanaka, Y. Mizukami et al., arXiv:2007.06757 (2020).

登録フォーム

第47回

日時: 2021年6月22日(火)11:00-
発表者: 那須 譲治（横浜国立大）
発表題目: キタエフ量子スピン系の熱輸送特性とその磁場角度依存性

登録フォーム

キタエフ量子スピン模型は、極低温まで磁気秩序を示さない量子スピン液体状態を基底状態に持つ可解な模型として注目を集めている。さらに、スピン軌道相互作用の強い強相関電子系において、その模型の実現可能性が指摘されたことで、候補物質の様々な物理量に対する実験が行われている。キタエフ量子スピン液体においては、量子スピンの分数化により、創発マヨラナ準粒子が生じ、それが磁場によってトポロジカル状態になることで、マヨラナ端状態と非可換エニオンが実現すると考えられている。最近のキタエフ候補物質α-RuCl3における実験で、マヨラナ粒子の特徴である電子系の場合の半分の値に熱ホール係数が量子化する振る舞いが観測されており[1]、精力的な研究が行われている。

本講演では、まず、キタエフ量子スピン液体に対してこれまで我々が行ってきたマヨラナ描像でのモンテカルロ法を用いて計算した熱ホール係数の数値計算結果を紹介する[2]。さらに、ハイゼンベルク相互作用やΓ相互作用など付加的な相互作用の効果を議論するために、マグノン描像において熱ホール係数を評価し、その磁場角度依存性を計算した結果も紹介する[3]。当日は、それらの計算結果と実験結果[4]との関係も議論する予定である。

[1] Y. Kasahara et al., Nature 559, 227 (2018), M. Yamashita et. al., Phys. Rev. B 102, 220404(R) (2020), J.A.N. Bruin et. al., arXiv:2104.12184. [2] J. Nasu, J. Yoshitake, and Y. Motome, Phys. Rev. Lett. 119, 127204 (2017). [3] S. Koyama, J. Nasu, arXiv:2105.03113. [4] T. Yokoi et. al., arXiv:2001.01899, O. Tanaka et. al., arXiv:2007.06757, 第48回講演.

第46回

日時: 2021年6月8日(火)11:00-
発表者: 関 和弘（理研）
発表題目: 強相関格子模型に対する変分量子計算法の古典計算機シミュレーション

近年利用可能な量子計算機は、誤り訂正を行わないものであり、NISQ (noisy intermediate-scale quantum computers)と呼ばれている。NISQを用いた計算では、 量子ビット数や量子ゲート数のなるべく小さい回路を設計し、古典計算機で実行できる部分は古典計算機で行うことで、エラーの影響を抑えて意味ある結果を得ようという量子-古典ハイブリッ ド法が用いられる。NISQを量子多体系のハミルトニアンの固有値問題へ応用する枠組みのうち代表的なものとして、 variational quantum eigensolver (VQE)法が提案された[1]。VQE法は量子多体系の適当な変分波動関数を量子計算機で用意してその変分波動関数に対するハミルトニアン期待値を計算し、ハミルトニアン期待値を最小化するように古典計算機で変分パラメータを最適化する変分法である。また、元の固有値問題を、複数の波動関数で張られる適当な部分空間で一般化固有値問題として定式化して解くquantum subspace expansion (QSE)法も提案された[2]。QSE法は部分空間におけるハミルトニアン行列と重なり行列の要素を量子計算機で計算し、一般化固有値問題を古典計算機で解くという量子と古典の役割分担がある。VQE法やQSE法では、標的とする波動関数に応じて適切な変分波動関数や適切な部分空間を選ぶことが課題となっている。

本講演では、我々の最近の取り組みとして、VQE法の派生として射影演算子を作用させた変分波動関数を用いる方法[3]や、適当な波動関数とハミルトニアンから生成されるクリロフ部分空間で計算を行う方法[4]について紹介する。

[1] A. Peruzzo, J. McClean, P. Shadbolt, M.-H. Yung, X.-Q. Zhou, P. J. Love, A. Aspuru-Guzik, and J. L. O’Brien, Nat. Commun. 5, 4213 (2014). [2] J. McClean, M. E. Kimchi-Schwartz, J. Carter, and W. A. de Jong, Phys. Rev. A 95, 042308 (2017). [3] K. Seki, T. Shirakawa, and S. Yunoki, Phys. Rev. A 101, 052340 (2020). [4] K. Seki and S. Yunoki, PRX Quantum 2, 010333 (2021).

登録フォーム

第45回

日時: 2021年5月25日(火)11:00- ※日程が変更になりました
発表者: 北折 暁（東大）
発表題目: 創発インダクター：材料開拓とこれからの課題

登録フォーム

磁気スキルミオンを始めとした非共面磁気構造においては、仮想的な磁場である創発磁場が伝導電子に作用する。創発磁場はこれに由来するトポロジカルホール効果を通じて、これまで多数の物質においてその静的な側面が実験的に観察されてきた。近年では創発電磁場の動的な側面に注目した研究も開拓されつつあり、その一環として提唱されたのが創発インダクター[1]である。創発インダクターはらせん型の磁気構造を電流駆動した際に生じる創発電場を検出・利用した素子で、素子の微細化が容易であることから応用面にも期待できる。既に本現象はGd3Ru4Al12において実証実験が行われている[2]が、この物質は転移温度が20K以下と低いことが問題であった。

本講演ではYMn6Sn6を用いた室温における創発インダクタンスの観測[3]を解説する。この実験結果は単に転移温度が上がっただけでなく、予期しなかった複数の特徴や示唆を含んでおり、その起源や特性に関しては理論的解明が待たれる状況である。さらに講演の最後には今後残された課題を列挙していき、本研究領域の発展を願って布教を行う。

[1] N. Nagaosa, Emergent inductor by spiral magnets. Jpn. J. Appl. Phys. 58, 12090 (2019). [2] T. Yokouchi et al. Emergent electromagnetic induction in a helical-spin magnet. Nature 586, 232-236 (2020). [3] A. Kitaori et al. Emergent electromagnetic induction beyond room temperature. arXiv:2103.02814

第44回

日時: 2021年4月20日(火)11:00- 注:曜日が4月から変更になります
発表者: 星野 晋太郎（埼玉大）
発表題目: ボゴリウボフ準粒子に対する相互作用の効果

登録フォーム

フェルミ面を持つ金属は、電子間相互作用によって多くの興味深い現象を引き起こす。その代表例である超伝導状態では、フェルミ面近傍の電子状態が再構成され、常伝導状態よりも低いエネルギーをもつ状態が実現する。通常、状態密度はフェルミ準位上でゼロになるためフェルミ面は消失している。ところが、ある種の超伝導体では、転移温度以下においてもフェルミ面が残るという可能性が指摘されている。このとき、フェルミ面を構成する粒子は電子そのものではなく、電子と正孔が量子力学的に重ねあわされたボゴリウボフ準粒子(bogolon) である点が特徴的であり、ボゴリウボフ・フェルミ面と呼ばれる。実験的にも、Fe(Se,S) の超伝導状態における比熱や熱伝導度の振る舞いから、ボゴリウボフ・フェルミ面が実現している可能性が報告されている。このような系は、通常の電子から成るフェルミ面とは定性的に異なる振る舞いを示すことが予想される。

本講演ではbogolon間の相互作用の効果が、通常の電子系とは異なる振る舞いを引き起こすことを紹介する。まず、bogolonに対する低エネルギー有効モデルを解析すると、奇周波数クーパー対という異時間でのみペアを組む特異な量子状態が形成されることがわかる[1]。また、相互作用はフェルミ面を不安定化させる効果もある。その一つの帰結が自発的対称性の破れを伴う偶周波数（同時刻）boglonクーパー対の形成である[2]。講演では、これらの性質を元の電子描像から調べた結果についても議論する。

[1] T. Miki, S.-T. Tamura, S. Iimura, and S. Hoshino, arXiv:2103.02251 (2021).
[2] S.-T. Tamura, S. Iimura, and S. Hoshino, Phys. Rev. B 102, 024505 (2020).

第43回

日時: 2021年4月6日(火)11:00- 注:曜日が4月から変更になります
発表者: 中川大也（東大）
発表題目: 散逸Hubbard模型における磁性と超流動ペアリング

登録フォーム

量子系が環境と相互作用することによる散逸の効果は、系の量子コヒーレンスを破壊する主要因であるため基礎・応用の両面において重要な問題である。このような量子開放系の問題は少数自由度の量子系においても重要であるが、近年冷却原子気体を用いることにより強相関・大自由度の量子多体系に対して実験的に制御可能な形で散逸の効果を調べることが可能になってきた。デコヒーレンスやダンピングといったネガティブな効果を想起させる散逸だが、このような散逸が多体物理と絡み合うことによって平衡状態では現れない秩序や転移が起こることが明らかとなってきた。本講演では、強相関多体系の代表例であるFermi-Hubbard模型において散逸が引き起こす秩序状態について議論する。まず、冷却原子気体における典型的な散逸である非弾性散乱の存在するHubbard模型において、スピン交換相互作用の効果が通常の平衡系とは逆転し、非平衡定常状態で強磁性が発現することを示す[1,2]。次に、散逸として光の自然放出による原子状態の遷移を考えると、非BCS型の非平衡ペアリング機構によりηペアリングという超流動状態が安定化することを示す[3]。これら2つの例は、斯波変換によって数学的に裏表の関係にある。最後に、[3]の研究過程で判明した、平衡・非平衡を問わない任意のフェルミオン系に対して成り立つCooperペアの運動量分布の普遍的な上限についても議論したい。

[1] M. Nakagawa, N. Tsuji, N. Kawakami, and M. Ueda, Phys. Rev. Lett. 124, 147203 (2020). [2] M. Nakagawa, N. Kawakami, and M. Ueda, Phys. Rev. Lett. 126, 110404 (2021). [3] M. Nakagawa, N. Tsuji, N. Kawakami, and M. Ueda, arXiv:2103.13624.

第42回

日時: 2021年3月29日(月)11:00- ※講演者の都合により延期になりました
発表者: 金杉翔太（京大）
発表題目: 多バンド系におけるアナポール超伝導の理論

第41回

日時: 2021年3月8日(月)11:00-
発表者: 國見昌哉（分子研）
発表題目: 冷却原子系における準古典近似を用いた非平衡ダイナミクスの研究

登録フォーム

冷却原子気体では、レーザー等の実験技術を駆使することにより、制御可能かつクリーンな量子多体系を実験的に実現することが可能である。このことを利用し、数値計算が困難である大自由度の量子多体ハミルトニアンを実験的にできるだけ忠実に実装し、”実験的に”多体シュレーディンガー方程式を解いてもらう、というのがファインマンが提唱した量子シミュレーションの概念である。実際に、冷却原子系、イオントラップ系、Rydberg原子系等を舞台とした様々な量子シミュレーションが行われている。 実験的にハミルトニアンがきれいに実装できるとはいうものの、より信頼性のある結果を得るためには可能な範囲内で理論計算との定量的な比較が重要である。本講演では実際の量子シミュレーションとの比較に用いられ、実験の定量的な再現に成功しているTruncated-Wigner近似(TWA)[1,2]を用いた研究結果を紹介する。TWAは一種の準古典近似であり、実時間ユニタリー時間発展を記述する手法である。特長としては、平均場近似を超えた近似法であること、大きなサイズ、高次元系への適用が可能などが挙げられる。 本講演の前半では、3次元Bose-Hubbard模型におけるクエンチダイナミクスをTWAを用いて計算した結果を紹介する。また、その結果と京都大学高橋グループが行った3次元光格子中のボース気体の実験結果を比較し、両者がフィッティングパラメータ無しでよく一致することを示す[3]。講演の後半では、量子スピン系に特化したTWAであるdiscrete TWA(DTWA)の計算結果を紹介する[4]。本研究ではDTWAを用いてRényiエンタングルメントエントロピーを計算する新しい手法の開発をした。また、これを利用したDTWAによる近似が妥当な時間スケールを評価する手法の提案をしたのであわせて紹介する。 [1] P. B. Blakie, et al., Adv. Phys. 57, 363 (2008). [2] A. Polkovnikov, Ann. Phys. 325, 1790 (2010). [3] K. Nagao, MK, Y. Takasu, Y. Takahashi, and I. Danshita, Phys. Rev. A 99, 023622 (2019). [4] MK, K. Nagao, S. Goto, and I. Danshita, Phys. Rev. Research, 3, 013060 (2021).

第40回

日時: 2021年3月1日(月)11:00-
発表者: 篠原 康（東大）
発表題目: 固体からの高次高調波発生の時間依存スペクトルの半古典ダイナミクスによる記述と、時間依存Hartree-Fockへの拡張

登録フォーム

10 V/nm級の光電場を物質に照射した際に、入射光の数十倍の光子エネルギーを持った電磁波が輻射される高次高調波発生と呼ばれる現象は気相(主に希ガス)による実験で確認された。この現象は、3ステップモデルと呼ばれる、1) 電子のイオン化, 2)電場による自由粒子様イオン化電子の加速, 3)親イオンへの再衝突という簡便な古典模型/描像でよく理解できる事が知られている[1]。 2011年、Ghimire等がZnO結晶から高次高調波発生の観測を報告した[2]。そのスペクトルのいくつかの特徴が気相におけるものと違い、議論を呼んだ。結晶は各々独自の電子状態をもっていて、強い電場で電子が駆動された際に、電子がどのように動き、そしてその運動の情報を高調波のスペクトルにどう焼き付けるかも結晶に応じて変わりうる。これは希ガスのような球対称系とは対照的である。結晶においては、電子間の相互作用を無視してもなお、特定の電子状態がどのような高調波を生成するかを説明する簡便な描像の決定版は未だ得られていない。 我々は、強い電場で駆動された結晶中電子のダイナミクスと、それがどのように輻射スペクトルに反映されるかを説明する簡便な描像を得るために、固体版の3ステップモデルともいうべき解析法を考案して、独立電子の一次元結晶における量子シミュレーションから得られる結果と比較した[3]。我々の解析法では3ステップモデルを固体版に翻訳したもので、i) 結晶ではエネルギー準位に禁制帯があるのでそこでの遷移を考慮すること ii) 価電子帯はエネルギーに幅があるので、全ての価電子帯から出発した電子の軌跡か考慮する の二点が原版との大きな違いである。我々の古典的な解析法は、量子シミュレーションで評価した時間分解スペクトルの形状をよく記述していることを示した。我々はこの理解を基に、一次元結晶の時間依存Hartree-Fock方程式の結果の解釈にとりくんだ[4]。その結果、平均場のレベルで考慮された電子・正孔対の間の引力が、電子・正孔対生成を促進する兆候を見出した。 セミナーでは、これらの詳細と、これらの模型をより詳細に理解するための最近の取り組みについて紹介する。

[1] P. B. Corkum, Phys. Rev. Lett. 71 (1993) 1994. [2] Shambhu Ghimire, Anthony D. DiChiara, Emily Sistrunk, Pierre Agostini, Louis F. DiMauro & David A. Reis, Nature Physics 7 (2011) 138.[3] Takuya Ikemachi, Yasushi Shinohara, Takeshi Sato, Junji Yumoto, Makoto Kuwata-Gonokami, Kenichi L. Ishikawa, Phys. Rev. A 95 (2017) 043416. [4] Takuya Ikemachi, Yasushi Shinohara, Takeshi Sato, Junji Yumoto, Makoto Kuwata-Gonokami, and Kenichi L. Ishikawa, Phys. Rev. A 98 (2018) 023415.

第39回

日時: 2021年2月22日(月)11:00-
発表者: 内藤 智也（東大）
発表題目: 物性理論は原子核理論の夢を見るか? ―原子核理論と物性理論のクロスオーバー

登録フォーム

物性物理 (原子, 分子, 固体) は電子の多体問題, 原子核物理は陽子と中性子 (核子) の多体問題と捉えられる. 電子, 陽子, 中性子は全てフェルミオンであることから, 原子核物理と物性物理は全てフェルミオンの多体問題であると捉えられ, 魔法数 (閉殻構造) などの共通した性質を持つ [1, 2] 一方, 相互作用の違いに由来する異なる性質も多く知られている. また, 原子核もフェルミオン多体系であることから, 密度汎関数理論や波動関数理論など, 電子系と共通した微視的な計算手法が広く用いられている [3]. 本講演では, まず, 原子核構造について簡単のレビューののち, 電子系と原子核の共通点・相違点に着目した原子核変形の起源に関する研究 [4] と, クーロン力の原子核構造への影響に関する近年の研究 [5, 6, 7] について簡単に紹介する. [1] K. Hagino and Y. Maeno. Found. Chem. 22, 267 (2020). [2] Y. Maeno, K. Hagion, and T. Ishiguro. Found. Chem. (In press). [3] J. A. Maruhn, P.-G. Reinhard, and E. Suraud. “Simple Models of Many-Fermion Systems” (Springer, 2010). [4] T. Naito, S. Endo, K. Hagino, and Y. Tanimura. arXiv:2009.05955 [physics.atom-ph] (2020). [5] T. Naito, R. Akashi, and H. Liang. Phys. Rev. C 97, 044319 (2018). [6] T. Naito, X. Roca-Maza, G. Colò, and H. Liang. Phys. Rev. C 99, 024309 (2019). [7] T. Naito, X. Roca-Maza, G. Colò, and H. Liang. Phys. Rev. C 9101 0263011 (2020).

第38回

日時: 2021年2月15日(月)11:00-
発表者: 池田 暁彦（ISSP）
発表題目: 超強磁場に誘起されるLaCoO3の励起子凝縮 〜ひずみ計測による研究〜

登録フォーム

相関電子系ではせめぎ合う2相の境界領域に新奇相が発現する。近年、固体中の励起子自由度から多軌道相関電子系を理解する試みがなされた結果、モット絶縁体とバンド絶縁体の境界領域に励起子絶縁体相が現れることが予言され着目されている 。実験ではスピンシングレットの励起子凝縮が遷移金属カルコゲナイド系で、スピントリプレット励起子の凝縮がコバルト酸化物系で候補とされている。なかでも磁場中のLaCoO3がトリプレット励起子凝縮の有力候補として報告されたが [1, 2]、実験的実証は進んでいなかった。これは実験的困難として、超強磁場を要することと、その中での観察手段が限られることに起因していた。近年、物性研では1000 T発生装置の稼働に成功し [3]、その中で利用できる磁歪測定装置の開発が進んだ [4]。そこで我々は、超強磁場中でLaCoO_3_のスピン状態変化（励起子数の変化）をひずみを通じて測定することを試みた。

これまでに実験により600 Tまでの磁歪計測を、初期温度5 , 78 , 108 Kの条件で得ることに成功した。この結果、超強磁場中に複数の相が存在することがわかった。まず低温領域にある磁場に対してギャップを持つ相（α相）[5]は600 Tまで安定であることがわかった。30 K以上100 T以上で発現するβ相はこれまでの実験でギャップを持つと考えていたが [5]、今回磁場範囲が大きく広がったことでギャップレスであることがわかった。さらに > 380 Tの領域に新たなギャップレス相（γ相）が発現することがわかった[6]。我々はα相の起源として励起子の超格子相（Solid）、βとγ相の起源として励起子凝縮相（BEC）の可能性を検討している。セミナーではこれらの現状を紹介する。

[1] T. Tatsuno et al., J. Phys. Soc. Jpn. 85, 083706 (2016). [2] A. Sotnikov and J. Kunes, Sci. Rep. 6, 30510 (2016). [3] A. Ikeda et al., Rev. Sci. Instrum. 88, 083906 (2017). [4] D. Nakamura, A. Ikeda et al., Rev. Sci. Instrum. 89, 095106 (2018). [5] A. Ikeda et al., Phys. Rev. Lett. 125, 177202 (2020). [6] A. Ikeda et al., unpublished.

第37回

日時: 2021年2月8日(月)11:00-
発表者: 後藤 慎平（近畿大）
発表題目: 冷却原子気体の観測誘起転移

登録フォーム

量子系のユニタリ時間発展は系の持つ量子性の指標であるエンタングルメントエントロピーを一般に増大させる。一方で環境との相互作用の結果である観測は、非ユニタリな影響を系に及ぼしてエンタングルメントエントロピーを減少させることがある。そしてこれらユニタリ時間発展と観測の効果の競合によって、古典シミュレーション可能な程度のエンタングルメントしか持たない面積則状態と一般に古典シミュレーション不可能な体積則状態の間に転移が生じることが最近の量子回路模型の研究から明らかになった [1-3]。我々はこの観測によって誘起される転移が冷却原子系で実現しうるか、またどのように検出できるかについて研究を行ったので [4] その結果について解説する。冷却原子系はユニタリ時間発展が観測可能なほど良い孤立系であり、また外部レーザーによって観測を導入することでその影響さえ制御可能である [5]。我々は散逸下の空間一次元Bose-Hubbard模型を数値的に解析することにより、冷却原子系においても観測誘起転移が生じることを示した。またエンタングルメントエントロピーは散逸下では実験的に測定することが非常に困難であるため、状態を特徴づける代わりの測定量が必要となる。そこで実験で測定可能な運動量分布とダイナミクスを計算してみたところ、系の状態を反映した振る舞いを見せることが分かった。

[1] Y. Li, X. Chen, and M. P. A. Fisher, Phys. Rev. B 98, 205136 (2018). [2] A. Chan et al., Phys. Rev. B 99, 224307 (2019). [3] B. Skinner, J. Ruhman, and N. Adam, Phys. Rev. X 9, 031009 (2019). [4] S. Goto and I. Danshita, Phys. Rev. A 102, 033316 (2020). [5] T. Tomita et al., Sci. Adv. 3, e1701513 (2017).

第36回

日時: 2021年2月1日(月)11:00-
発表者: 石塚 淳（京大）
発表題目: ウラン系超伝導体UTe2の最近の進展と第一原理計算によるトポロジカル超伝導予測

登録フォーム

2018年末にウラン系化合物UTe2において超伝導が報告され、強磁性臨界点近傍のスピン三重項超伝導体の可能性が議論されており大きな注目を集めている[1, 2]。他のウラン系強磁性超伝導体URhGeなどとは異なり、常圧で強磁性が起こっていないことは実験研究にとって極めて有利であり、スピン三重項超伝導の特性を調べるのに理想的な系であるといえる。また、時間反転対称性が保たれていれば、これまでとは異なるクラスのトポロジカル超伝導となっている可能性がある。本講演では、これまでに明らかとなったUTe2の電子状態や超伝導特性について解説する。UTe2の第一原理バンド計算は、絶縁体を予言するため実験の金属的な振舞いと大きな矛盾がある。そこでまず我々は、電子相関を取り入れるためGGA+U法を用いてUTe2の電子状態を調べた[3]。結果、クーロン相互作用U = 1 eV付近で絶縁体金属転移を観測し、金属的な電子状態を得た。このフェルミ面に基づいて超伝導状態のトポロジカルな性質を解明した。その結果、空間群Immmの全ての可能な奇パリティ超伝導は、Uが中間的な領域でトポロジカルに非自明であることを示した。さらに、有効的な周期アンダーソン模型を構築し、磁気揺らぎによって誘起される超伝導について調べたところ、スピン三重項状態(B3u状態とAu状態)とスピン一重項状態(Ag状態)が安定となることを示した[4]。

[1] S. Ran, et al., Science 365, 684 (2019). [2] D. Aoki, et al., J. Phys. Soc. Jpn. 88, 043702 (2019). [3] J. Ishizuka, S. Sumita, A. Daido, and Y. Yanase, Phys. Rev. Lett. 123, 217001 (2019). [4] J. Ishizuka and Y. Yanase, arXiv:2008.01945 (2020).

第35回

日時: 2021年1月25日(月)11:00-
発表者: 金子 竜也（コロンビア大）
発表題目: η-pairing状態の光誘起理論

登録フォーム

光を用いたポンプ-プローブ法によって相関電子系の非平衡物性を捉えることが可能となったことを背景に、非平衡過程を経て誘起される量子状態の研究が進んでいる[1]。本講演では、Mott絶縁体の励起状態であるη-pairing状態に着目し、その光誘起理論について解説する。η-pairing状態は約30年前に数理物理学者の楊振寧 (C. N. Yang) によって指摘されたHubbard模型の厳密な固有励起状態で、長距離に発達したペア密度波的な超伝導相関を有する状態である[2]。η-pairing状態は基底状態で実現しなかったためあまり注目されてこなかったが、我々は厳密対角化法などを駆使し、η-pairing状態が光照射によって誘起できることを理論的に示した[3]。特に、外場を印加したHubbard模型におけるη-pairing相関の上昇が、Hubbard模型の対称性とそれに基づく光学的選択則から説明できることを明らかにした[3]。また、我々は光誘起状態の伝道特性がη-pairing相関と深く結びついていることも示した[4]。今回はη-pairingの数学的な背景もふまえ、対応する数値計算の結果も示しながら、光誘起η-pairingの理論について紹介する予定である。

[1] C. Giannetti, M. Capone, D. Fausti, M. Fabrizio, F. Parmigiani, and D. Mihailovic, Adv. Phys. 65, 58 (2016). [2] C. N. Yang, Phys. Rev. Lett. 63, 2144 (1989). [3] T. Kaneko, T. Shirakawa, S. Sorella, and S. Yunoki, Phys. Rev. Lett. 122, 077002 (2019). [4] T. Kaneko, S. Yunoki, and A. J. Millis, Phys. Rev.Research 2, 032027 (2020). [5] 金子竜也, 白川知功, 柚木清司, 固体物理 55(1), 21-29 (2020).

第34回

日時: 2021年1月18日(月)11:00-
発表者: 石塚 大晃（東工大）
発表題目: テラヘルツ光によるスピン流・軌道流の整流効果

登録フォーム

光を用いた磁性体の制御（光スピントロニクス）は，高輝度のレーザーの発達により現実的な研究課題となっている．これらの研究では磁荷の制御に主眼が置かれており，スピン流に関してはスピンポンプなど拡散流に限られており，太陽電池のような整流効果による流れの直接制御については最近まで研究されてこなかった． 本講演では，テラヘルツ光を用いた磁性体や半導体における整流効果について解説する．講演の前半では磁性体の磁気励起を用いたスピン流の整流について，一次元スピン鎖におけるスピノン流の整流[1]と，反強磁性秩序におけるマグノンによる整流効果[2]に関する非線形応答理論を用いた我々の最近の研究について紹介する．具体的なモデルを用いた計算からは，現在得られるレーザー強度でも観測可能なスピン流が得られることが解った．また，様々なセットアップを用いた解析から，反転対称性の破れた磁性体であれば，一般にスピン流の生成ができることを確認した．これらの結果は反転対称性の破れた系を対象としているが，反転対称なトポロジカル絶縁体においても光による軌道流の整流ができることを最近見出した[3]．時間があれば，本結果についても紹介する．

[1] H. Ishizuka & M. Sato, Phys. Rev. Lett. 122, 197702 (2019). [2] H. Ishizuka & M. Sato, Phys. Rev. B 100, 224411 (2019). [3] M. Davydova et al., in preparation.

第33回

日時: 2021年1月12日(火)11:00- 注:曜日がいつもと違います
発表者: 古谷 峻介（茨城大）
発表題目: 量子スピン系の磁化プラトーとアノマリー

登録フォーム

Lieb-Schultz-Mattis (LSM) 定理は最近、様々な系・文脈において研究されているが、元の定理はS=1/2 XXZ量子スピン鎖がuniqueかつgappedな基底状態を持たないことを主張するものであった。LSM定理は量子多体系のU(1) fluxへの応答という観点で定式化することができる。LSM定理は時間反転対称性存在下で議論されることが多いが、量子スピン系においてこのU(1) fluxによる定式化は、磁場存在下の状況に容易に拡張でき、量子スピン系の磁化プラトー上の基底状態縮退度（の下限）についての制限を課す。この条件はOshikawa-Yamanaka-Affleck条件[1]として知られ、2次元や3次元の量子スピン系でも成り立つことが経験的に知られている。最近ではLSM定理とその関連する現象が、量子相に対応する場の理論の持つ’t Hooftアノマリーの観点から調べられている[2]。他方で、磁化プラトーの性質のアノマリーによる研究はあまり行われていないのが現状である。 本講演ではLSM定理とOshikawa-Yamanaka-Affleck条件、そしてアノマリーという観点から、フラストレート磁性体の磁化プラトー上の基底状態の性質を議論する。まず、チェッカーボード上の反強磁性Heisenberg模型の磁化プラトー上の基底状態のアノマリーを、空間的なtwistを伴う境界条件下で調べる[3]。次にカゴメ格子や三角スピンチューブ上の量子スピン系における1/3プラトー上の基底状態の共通点と違いについて、SU(3) Wess-Zumino-Witten理論のアノマリー[4]を使って議論する[5]。

[1] M. Oshikawa, M. Yamanaka, and I. Affleck, Phys. Rev. Lett. 78, 1984 (1997). [2] S. C. Furuya and M. Oshikawa, Phys. Rev. Lett. 118, 021601 (2017). [3] S. C. Furuya and Y. Horinouchi, Phys. Rev. B 100, 174435 (2019). [4] Y. Yao, C.-T. Hsieh, and M. Oshikawa, Phys. Rev. Lett. 123, 180201 (2019). [5] S. C. Furuya, Y. Horinouchi, and T. Momoi, arXiv:2011.13095.

第32回

日時: 2020年12月21日(月)11:00-
発表者: 関孝一（新潟大）
発表題目: ランダム量子スピン系での実空間数値くりこみ群法

登録フォーム

相互作用にランダムネスを含んだ二次元量子フラストレートスピン系では、スピン液体のような物理量の振る舞いを示す可能性が指摘されており[1,2]、現在も活発な研究が行われている。これまでの研究では、厳密対角化法や量子モンテカルロ法を用いた数値計算が行われてきた。しかし、これらの手法にはシステムサイズや扱える模型に制限があるため、二次元ランダム量子スピン系を効率よく計算する新しい手法の開発は重要なテーマである。一方で、一次元ランダムスピン系においては、Maらによる実空間くりこみ群法[3]および、スピンブロック間のエネルギーに着目して状態を選択する改良[4]によって、定量的に良い基底状態が得られることが示されている。

そこで我々は、スピンブロック間のエネルギースケールに着目することで、実空間数値くりこみ群法に対するさらなる改良を行い、二次元ランダム量子スピン系の数値計算への適用と性能評価を行った。特に、反強磁性相互作用にボックスタイプのランダムネスが存在するS=1/2ハイゼンベルグ模型に対する数値計算を行い、ランダムネスの強さを変えて基底状態を求めた。セミナーでは、実空間数値くりこみ群法に対する改良点および、正方格子・三角格子に対する計算結果や、モデルと改良型数値繰り込み群法の相性の詳細等について紹介する。

[1] T. Shimokawa, K. Watanabe, and H. Kawamura, Phys. Rev. B 92, 134407 (2015). [2] K. Uematsu and H. Kawamura, J. Phys. Soc. Jpn. 86, 044704 (2017). [3] S.-k. Ma, C. Dasgupta, and C.-k. Hu, Phys. Rev. Lett. 43, 1434 (1979). [4] T. Hikihara, A. Furusaki, and M. Sigrist, Phys. Rev. B 60, 12116 (1999).

第31回

日時: 2020年12月14日(月)11:00-
発表者: 工藤耕司（筑波大）
発表題目: 量子ホール効果の断熱接続における一般化ストレーダ公式

登録フォーム

断熱過程におけるトポロジカル数の不変性は、近年のトポロジカル相の理論の基礎である。その源流の一つとして、adiabatic heuristic argument[1]が挙げられる。これによると分数量子ホール(FQH)状態は、粒子の統計性を変化させることで、整数量子ホール(IQH)状態に断熱的に変形できる。この議論はFQH状態のギャップの起源を説明する上で大変有用だが、統計性を規定する組紐群との整合性に注意する必要がある。例えば、バルクギャップを調べる上で重要な周期系では、組紐群の代数的制約から統計性の連続変形は許されない。さらに周期系で生じるトポロジカル縮退は、この断熱変形に一層複雑な構造をもたらすため、トーラス上におけるその正当性はより非自明な問題となる。

本研究で我々は、トーラス上のadiabatic heuristic argumentの正当性を数値的に示すとともに、「トポロジカル不変性の本質は状態の連続性ではなくギャップの連続性にある」ことを示唆する一般化ストレーダ公式を見出した[2]。まず我々は、多体状態の連続変形はできないもののIQH・FQH間のギャップは連続かつ断熱的につながることを数値的に示した。その過程においてトポロジカル縮退度は不規則に変化するが、多体チャーン数は断熱不変量となった。さらに不変量としての多体チャーン数が、縮退度の変化を記述することを解析的に示した。この関係式は、ホフスタッターの蝶におけるストレーダ公式と類似した構造を示唆している。

[1] M. Greiter and F. Wilczek, Modern Physics Letters B 04, 1063 (1990); M. Greiter and F. Wilczek, Nuclear Physics B 370, 577 (1992). [2] K. Kudo and Y. Hatsugai, Phys. Rev. B 102, 125108 (2020)

第30回

日時: 2020年12月07日(月)11:00-
発表者: 望月健（北大）
発表題目: フロケ系や開放系特有のトポロジカルな現象：量子ウォークを舞台として

登録フォーム

過去数十年に渡って、対称性に守られたトポロジカル相に関する膨大な研究が行われ、電子系における新奇現象の発見や既知現象の再解釈が成されてきた。近年では、電子系のみならず光学系においても非自明なトポロジカル相に由来する現象が議論されている[1]。

本講演では、そのような光学系の一種である量子ウォークに着目し、トポロジカルな現象を調べる。量子ウォークでは系のパラメータを通じてトポロジカル数を制御できるため、光子の存在確率から実空間で直接的にエッジ状態を観測する事が可能である[2]。そのため、量子ウォークの様々なトポロジカル相及びそれに由来するエッジ状態が調べられてきた。量子ウォークのダイナミクスとシュレーディンガー方程式もしくはディラック方程式のダイナミクスは類似性が高いため、量子ウォークは時間非依存ハミルトニアンを生成子とするユニタリダイナミクスの“シミュレータ”と称される事が多い。実際、量子ウォークのトポロジカル相に関する多くの性質は、時間非依存ハミルトニアンで記述される電子系との類推から理解できる。しかし、実験的制約から現時点ではそのような“シミュレーション”は決して効率的ではない。また、電子系の“模倣”は新しい光学デバイスなどの応用可能性を秘めている一方で、本質的に新しい現象の発見及び理解には繋がり辛いと考えられる。

本研究では、時間非依存ハミルトニアンで記述される孤立系との類推からは理解できない現象を探索するため、量子ウォークの持つ二つの特徴に着目した。一つ目は、光の増幅減衰効果の制御性である。増幅減衰効果が存在する量子ウォークは、非ユニタリな時間発展演算子で記述される開放系となる。二つ目は、量子ウォークはフロケ系の一種であり離散的な時間発展をする事である。その結果我々は、開放フロケ系特有の長寿命エッジ状態が出現する事を明らかにした[3]。この時、エッジ状態の存在はカイラル対称性に守られ、寿命の延びは空間時間反転対称性の破れに起因する。我々の理論予想は単一光子を用いた系で実証され、それは開放量子系における時間空間反転対称性の破れに対する初の観測実験となった[4]。さらに時間があれば、時間グライド対称性というフロケ系特有の対称性に起因するトポロジカル相についても紹介したい[5]。

[1] T. Ozawa, H. M. Price, A. Amo, N. Goldman, M. Hafezi, L. Lu, M. C. Rechtsman, D. Schuster, J. Simon, O. Zilberberg, and I. Carusotto. Review of Modern Physics, 91 015006 (2019). [2] T. Kitagawa, M. A. Broome, A. Fedrizzi, M. S. Rudner, E. Berg, I. Kassal, A. Aspuru-Guzik, E. Demler, and A. G. White. Nature communications, 3 882, (2012) [3] Ken Mochizuki, D. Kim, N. Kawakami, and H. Obuse, arXiv:1609.09650. [4] L. Xiao, X. Zhan, Z. H. Bian, K. K. Wang, X. Zhang, X. P. Wang, J. Li, Ken Mochizuki, D. Kim, N. Kawakami, W. Yi, H. Obuse, B. C. Sanders, and P. Xue, Nature Physics 13, 1117 (2017). [5] Ken Mochizuki, T. Bessho, M. Sato, and H. Obuse, Physical Review B 102, 035418 (2020).

第29回

日時: 2020年11月30日(月)11:00-
発表者: 永井瞭（ISSP）
発表題目: 機械学習手法を用いた交換相関汎関数の構築とその展望

登録フォーム

密度汎関数理論（DFT）は、穏やかな計算量で比較的高精度な計算を可能とするため、電子状態計算における標準的な手法として用いられている。しかし、DFTの電子間相互作用を書き表す部分の中に、厳密に書き下せない交換相関汎関数という部分が存在する。

我々は機械学習の手法を用いて交換相関汎関数を系統的に構築する手法を研究している。まず、単純な１次元系において手法を適用し、DFTの数式操作の中に機械学習で構築された汎関数が有効であることを示す。その後、実際の物質への手法の拡張について紹介する。小分子の厳密な電子状態データをもとに機械学習を行うことで、幅広い分子系に用いることができる汎関数を構築した事例を示す。 また、交換相関汎関数のような立式が困難な理論に対して機械学習を用いることによって期待される改善や新たな知見について議論する予定である。

[1] R. Nagai, R. Akashi, S. Sasaki, and S. Tsuneyuki: Neural-network kohn-sham exchange-correlation potential and its out-of-training transferability. J. Chem. Phys. 148, 241737 (2018). [2] R. Nagai, R. Akashi and O. Sugino: Completing density functional theory by machine learning hidden messages from molecules. npj Comput. Mater. 6, 43(2020).

第28回

日時: 2020年11月24日(火)11:00- 注:曜日がいつもと違います
発表者: 高橋惇（中国科学院 物理研究所）
発表題目: 脱閉じ込め量子臨界現象に対する“de-signer”ハミルトニアンによるアプローチ

登録フォーム

反強磁性的な基底状態として古典的にはstaggered磁化を持つNéel状態が代表的だが、量子系特有の状態としてスピンのSU(2)などの対称性を保ったまま空間対称性のみを自発的に破るvalence-bond固体(VBS)状態も複数種類存在する。これらの間のNéel-VBS転移は最も「基本的」な量子反強磁性体における量子相転移と言えるが、その詳しい性質は未解明である。素朴なGinzburg-Landau描像に立てば二つの相は全く関係のない対称性を自発的に破るため、直接転移は一般的には不連続転移になると予想される。しかし、二つの相を秩序の欠陥が別の様態で閉じ込められている状態とみなし、いくつかの理論的仮定を置くことで転移点では欠陥が脱閉じ込めを起こし一般に臨界的(連続転移)になるという、場の理論によるDeconfined Quantum Criticality (DQC)のシナリオも提案されている[1]。DQCの成立を巡っては理論的[2]、数値的[3]、そして実験的[4]な研究も盛んに行われているが、未だに決着はついていない。

そこで本研究では負符号問題を避けつつ目標の相図を持つハミルトニアンをデザインするという“de-signer”ハミルトニアンの思想[5]に基づいて、実際のNéel-VBS転移の様子を量子モンテカルロ法によって調べるアプローチを取った。プラケット型VBS(PVBS)という従来実現されていなかったVBS状態からのPVBS-Néel転移を起こす模型を構成し、その結果今まで調べられていた柱状VBS(CVBS)と異なり、弱い不連続転移を示すことが分かった[6]。興味深いことに、転移点直上ではO(3)のNéel秩序とO(2)のVBS秩序が非自明に混合し創発的なSO(5)対称性を生じるという理論的予言[2]を部分的に実現していた。またPVBS秩序に残るZ2対称性を更に自発的に破る相転移も発見し、PVBS相が「破れうる対称性を途中までしか自発的に破らない」中間相であると判明した。このような多段階の離散的対称性の破れを手軽に整理し、可能な相図のトポロジーを議論するためのグラフ理論的な枠組みについても紹介する。

時間が許せば、ある意味でCVBSとPVBSの中間とみなせるような、今まで議論されていなかったヘリカルvalence-bond相へのNéel相からの直接転移を起こす模型をde-signerハミルトニアンの枠組みで発見／構成した研究[7]についても紹介したい。

[1] T. Senthil, et al., Science 303 1490 (2004). [2] C. Wang, et al., Phys. Rev. X 7 031051 (2017). [3] H. Shao, et al., Science 352 213 (2016). [4] J. Guo, et al., Phys. Rev. Lett. 124 206602 (2020). [5] R. K. Kaul, et al., Annu. Rev. Condens. Matter Phys. 4 179 (2013). [6] J. Takahashi and A. W. Sandvik, Phys. Rev. Research 2 033459 (2020). [7] B. Zhao, J. Takahashi, and A. W. Sandvik, arXiv:2005.10184 (2020).

第27回

日時: 2020年11月16日(月)11:00-
発表者: 磯部大樹（MIT）
発表題目: ファンホーブ特異性のもたらす電子相関物性：超伝導・絶縁体状態および非フェルミ流体

登録フォーム

グラフェン2層の超格子構造における絶縁体および超伝導状態の実験的観測に端を発し，ファンデルワールス物質からなる超格子構造における電子相関効果の研究が実験と理論の両面から盛んに行われている．このような物質系では，積層の自由度に制御された幅の狭いエネルギーバンドにおいて電子相関効果が強く現れる．本セミナーでは，まず2層グラフェンにおける絶縁体および超伝導状態をファンホーブ特異点とフェルミ面のネスティングに着目して理論的に説明することを試みる[1]．解析は繰り込み群の手法を用いる．続いて，2層グラフェンにおける魔法角をファンホーブ特異性の観点から考察する[2]．すなわちグラフェン2層の相対角度の制御により，ファンホーブ特異点における状態密度の発散がエネルギーの関数としてべき的になることを示し，それに対応する角度を魔法角と見なせることを説明する．このような状態を，一般に高次ファンホーブ特異性と呼ぶこととする．最後にこの高次ファンホーブ特異点における電子相関効果について議論する[3]．繰り込み群による解析により，非フェルミ流体が安定な固定点として生じることを示す．

[1] H. Isobe, N. F. Q. Yuan, and L. Fu, Phys. Rev. X 8, 041041 (2018). [2] N. F. Q. Yuan, H. Isobe, and L. Fu, Nat. Commun. 10, 5769 (2019). [3] H. Isobe and L. Fu, Phys. Rev. Research 1, 033206 (2019).

第26回

日時: 2020年11月9日(月)11:00-
発表者: 岩﨑祐昂（東大）
発表題目: バンドエンジニアリングによる近似結晶半導体の探索と熱電材料への応用

登録フォーム

準結晶は並進対称性と両立しない回転対称性を持った固体であり、これまでに70種類以上の熱力学的に安定な準結晶が確認されてきた。しかし、それらの電子状態はいずれも金属であり、原子スケールで半導体や絶縁体的な性質を持ったものは見つかっていない。準結晶半導体が存在するかどうかは固体物理学の基本的な問題の一つであると同時に高性能熱電材料への応用も期待されている[1]。準結晶と同じクラスターを持ち，それが周期的に配列した構造を有する近似結晶は，準結晶の物性を理解する際に重要な概念であり，両者を比較するアプローチによって準結晶の理解が深められてきた。周期性を持たない準結晶のバンド構造を計算し、半導体を見つけるのは困難であるが、近似結晶であればそれは可能であるため、近似結晶半導体の探索は準結晶半導体を発見するための良いアプローチとなり得る。本発表では、第一原理計算を援用したアルミ系近似結晶半導体の探索についての一連の研究について紹介する。

近年、Mihalkovičらによる計算の結果、Al-Ir系1/0近似結晶が特定の組成と結晶構造を満たしたときに、40 meV程度のバンドギャップを有する半導体となる可能性が指摘された[2]。しかし、実験的には、Alの原子欠損により半導体化は実現せず、高圧合成[3]や元素置換[4]による欠損の抑制を試みたがいずれも半導体化には至らなかった。そこでバンドギャップが確認されているAl-Ir系1/0近似結晶を基に、第一原理計算により、ギャップの大きさを制御する指針を確立し、その指針に基づいて近似結晶半導体を探索した。第一原理計算によって価電子帯上端と伝導体下端を構成する軌道を可視化することで、Alの一部をSiで置換し、IrをRuで置換した、Al-Si-Ru系が0.26 eV程度のギャップをもつ半導体になる可能性を見出した。この予想に基づいて、Al-Si-Ru系で近似結晶の合成を試み、Al67.6Si8.9Ru23.5近傍の組成で、1/0近似結晶の単相試料が作製できることが分かった。単相試料の熱電物性測定から、Al-Si-Ru系近似結晶が半導体的な特性を示すことが分かった。第一原理計算と実験値の比較から、バンドギャップは0.15 eV程度であると考えられ、初めて近似結晶半導体を実験的に作製することに成功した[5]。

時間に余裕があれば、Al-Si-Ru系近似結晶の熱電特性の最適化や高次近似結晶の探索に関する最近の研究についても紹介したい。

[1] Y. Takagiwa and K. Kimura, Science and Technology of Advanced Materials 15, 044802 (2014). [2] Yutaka Iwasaki, Koichi Kitahara and Kaoru Kimura, Journal of Alloys and Compounds 763, 78 (2018). [3] 岩﨑祐昂，北原功一，木村 薫, 日本熱電学会誌, 16, 139 (2020). [4] Y. Iwasaki, Koichi Kitahara, and Kaoru Kimura, Journal of Alloys and Compounds 851 156904 (2021). [5] Yutaka Iwasaki, Koichi Kitahara, and Kaoru Kimura, Physical Review Materials 3, 061601 (2019).

第25回

日時: 2020年11月2日(月)11:00-
発表者: 山本大輔（青学）
発表題目: 磁場中SU(3) Heisenberg模型の量子相と隠れたネマティック性を持つスピンダイポール秩序

登録フォーム

固体物質系では通常スピン自由度σ=↑,↓を持つ電子が主役となるため、Hubbard模型やHeisenberg模型といった代表的な統計模型はN=2次の特殊ユニタリ群SU(2)の対称性（もしくはそれ以下）を持つ。より高次の対称性（N>2）はスピン液晶物質や遷移金属酸化物などの模型において、複数の異なる相互作用の大きさがたまたま一致するときにしか得られない。したがってN>2の高次対称性を固体物質で実現するためには系のパラメータの偶発的なfine-tuningか、何らかのエキゾチックな発現機構による創発的なSU(N>2)相としての出現が必要となる。

一方で近年、レーザー光の対向波のよるポテンシャル（光格子）中に閉殻構造を持つ173Yb原子や87Sr原子などのアルカリ土類（様）原子気体を充填することで、直接的にSU(N)の拡張Hubbard模型やHeisenberg模型を作成する試みが盛んに為されている。例えば173Yb原子は（閉殻構造のために）電子スピン由来のスピン自由度は持たないが、核スピンI=5/2に由来して６成分の自由度を持つ。この成分間の散乱長はすべて等しく等価であり、さらに光ポンピングの操作によって６成分のうちのいくつかのみに偏らせることもできる。実際、既にSU(3), SU(4)やSU(6)などのHubbard系におけるMott絶縁体の作成[1,2]や超格子構造を用いた短距離磁気相関の観測[3]などが実験的に達成されている。

本研究ではSU(3)系に対する強磁場効果を議論する。強磁場の印加は磁性体研究の基本的な物性計測実験であり、格子幾何やトポロジー、量子・熱ゆらぎとの協奏効果によって磁化プラトーやネマティック状態、磁場誘起スピン液体などの非自明な磁気状態を生み出すことが知られている。ここでは三角格子上のSU(3)反強磁性Heisenberg模型に対する磁場効果を調べた[4]。密度行列繰り込み群をソルバーとした大規模クラスター平均場＋スケーリング法[5]を用いて絶対零度の量子磁性を、直積状態近似波動関数を用いた準古典Monte Carlo法を用いて有限温度の相転移現象を議論した。その結果、磁化およびスカラーネマティック秩序変数の磁場依存性における非自明なプラトーの形成、副格子スピン（ダイポール）モーメントが有限な強磁場相における「隠されたネマティック秩序」を持つ“ネマティック-スピンダイポール相”、渦度1/2を持つ分数渦対励起が誘起するトポロジカル相転移などの様々な興味深い物性現象を見出した[4]。SU(3)対称点から若干ずれた場合（スピン-1 bilinear-biquadratic模型に対応）にこれらの性質が保たれるか否かに関しても議論したい。

[1] S. Taie, R. Yamazaki, S. Sugawa, and Y. Takahashi, Nat. Phys. 8, 825 (2012). [2] C. Hofrichter et al., Phys. Rev. X 6, 021030 (2016). [3] H. Ozawa, S. Taie, Y. Takasu, and Y. Takahashi, Phys. Rev. Lett. 121, 225303 (2018). [4] D. Yamamoto et al., Phys. Rev. Lett. 125, 057204 (2020). [5] D. Yamamoto, G. Marmorini, and I. Danshita, Phys. Rev. Lett. 112, 127203 (2014).

第24回

日時: 2020年10月26日(月)11:00-
発表者: 兎子尾理貴（京大）
発表題目: 空間反転対称性が破れた高伝導金属における流体力学的異常輸送現象

登録フォーム

近年、graphene, PdCoO2, WP2等を代表とする、様々な高伝導金属・半導体において、「流体力学領域（hydrodynamic regime）」と呼ばれる、新しいクラスの非平衡領域が実現していることが明らかになってきた[1-3]。この領域では、不純物・フォノン散乱などの運動量緩和散乱が強く抑制された結果、電子-電子散乱が最も支配的な散乱プロセスとなり、電子系のダイナミクスは流体力学理論によって有効的に記述されることが理論的に予言される。そして実際にこのような兆候が近年、メゾ系での「負の非局所抵抗・磁気抵抗」[3,4]や電荷中性点近傍のgrapheneでの「Wiedemann-Franz 則の破れ」[5]などを通して実験的に確認され、実験と理論の両面から多くの関心を集め始めている。特に、近年実現したNoncentrosymmetricな電子流体物質（bilayer-graphene, WP2等）は空間反転対称の破れに伴う有限のBerry 曲率によって、これまでにない異常な流体力学的応答を実現すると期待されるため、これら結晶対称性や幾何学的性質の効果を流体力学と統合することは今日、極めて重要な課題になっていると言える。

そこで、本研究[6] では空間反転対称が破れた金属における流体力学理論を定式化し、空間反転対称性の破れ・Berry曲率がもたらす非従来的な流体力学的輸送応答を明らかにした。この拡張された流体力学理論では、ミクロなBloch波動関数が示す幾何学的効果が流体スケールにおいて付加的な駆動力として記述され、電子流体に特殊なflow（異常エッジ電流、非対称性Poiseuille flowなど）を引き起こす。特に興味深いのは、これらの系においてカイラル渦効果（渦誘起の異常電流）を結晶系に一般化した新奇な異常輸送現象（Generalized vortical effect: GVE）が発現する点であり、これは「空間反転対称が破れた系における電子流体」と「真空中のカイラル流体」との間にある種のアナロジーが成り立っていることを示している。本講演では、微視的モデルに基づいた定量的評価を与えるとともに、上記の輸送応答を観測する実験的セットアップも提案したい。

[1] P. J. W. Moll, et al., Science 351 1061 (2016). [2] R. K. Kumar, et al., Nat. Phys. 13 1182 (2017). [3] J. Gooth, et al., Nat. Comm. 9 4093 (2018). [4] P. S. Alekseev, Phys. Rev. Lett. 117 166601 (2016). [5] J. Crossno, et al., Science. 351 6277 (2016). [6] R. Toshio, K. Takasan and N. Kawakami, Phys. Rev. Research 2, 032021(R).

第23回

日時: 2020年10月19日(月)11:00-
発表者: 村上雄太（東工大）
発表題目: 励起子絶縁相の非線形応答と非平衡誘起

登録フォーム

励起子絶縁相とは、電子とホールのペアである励起子が自発的に凝縮することで生じる秩序相である。 予言から半世紀以上経った現在でも励起子絶縁相を示す物質は見つかっていないが、 平衡状態と非平衡状態の両方のセットアップで出現する可能性が指摘されており、その実現/発見に向けて精力的に研究が進められている。 本発表では、我々が最近行なった、i) 平衡状態における励起子絶縁相の集団励起と非線形応答に関する研究[1][2]と、 ii) 強相関電子系における非平衡状態として現れる励起子絶縁相に関する研究[3]、を紹介する。

平衡系においては、近年Ta2NiSe5などの遷移金属カルコゲナイト系の物質が励起子絶縁相を示す候補物質としてあげられ、注目を浴びている[4]。 しかし、現実物質では電子格子相互作用によっても似たような物性を示す秩序相が出現しうるため、励起子絶縁相の同定が困難である。 異なる秩序相の起源は、物質の集団励起の性質に反映されることが期待される。 そこで、我々は、格子自由度と結合した励起子絶縁相を考え、その集団励起モードの性質を明らかにした[1]。 有限の電子格子相互作用により位相モードは有限のギャップを持ち、線形応答の範囲ではフォノンモードとの区別が難しくなるが、 これらのモードを共鳴的に励起した時の応答が、位相モードとフォノンモードで大きく異なることがわかった。 講演では、このような非線形応答の違いとそれを用いたモードの同定に関する議論を行うとともに、 実空間実時間分解ポンププローブ実験を用いた集団励起モード観測の可能性[2]に関しても議論する。

一方、励起子絶縁相は、光励起された半導体の非平衡定常状態においても実現すると考えられている。 これまで、非平衡定常状態による励起子絶縁相の出現は報告されていないものの、光励起中の半導体において 励起子絶縁相に類似したコヒーレントな過渡状態が報告された[5]。このような状態は系が熱化する前の prethermalな状態と捉えることができるが、このような状態は電子間相互作用や電子格子交互作用による散乱 が弱い場合に長寿命になる。 我々は、励起された強相関電子系においても長寿命のprethermalな状態が実現しうることに注目し、 スピン遷移近傍の多軌道ハバード模型の非平衡ダイナミクスを動的平均場理論を用いて調べた[3]。 これにより、多軌道強相関系においてもprethermal状態において励起子絶縁相が実現することを 明らかにした。講演では、ダイナミクスの詳細に加え、有効模型を用いた解析を紹介し、 過渡的な励起子絶縁相の出現が局所的なエントロピー保存とエネルギー保存に由来することを議論する。

[1] Y. Murakami, D. Golež, T. Kaneko, A. Koga, A. J Millis, P. Werner, PRB(Editors’ suggestion) 101, 195118(2020) ; D. Golez, Z. Sun, Y. Murakami, A. Georges, A. J Millis, arXiv:2007.09749 [2] P. Andrich, H. Bretscher, Y. Murakami,et al., arXiv:2003.10799. [3] P. Werner and Y. Murakami, arXiv:2006.02096. [4] Y. Wakisaka et al., PRL 103, 026402 (2009). [5] Y. Murotani, et al., PRL 123, 197401 (2019).

第22回

日時: 2020年10月12日(月)11:00-
発表者: 井戸康太（東大物性研）
発表題目: 励起状態のための変分モンテカルロ法

登録フォーム

変分モンテカルロ法(VMC)は様々な量子多体系を負符号問題なく解析できる強力な数値計算手法である。VMCでは多数の変分パラメータを試行波動関数に導入することで様々な量子状態を柔軟に表現できるが、その適用の多くは、エネルギーや秩序変数などの静的物理量の測定に基づいた基底・低励起状態の解析に留まっていた。近年、VMCの枠組みで励起状態を求める研究が活発に行われている。本講演では、我々が開発してきた励起状態のためのVMCに関連した研究を2つ紹介する。

まず、有限温度の物理量を計算するためのVMCを紹介する。近年、熱的純粋状態(TPQ)などの発展により、各温度の混合状態を再現できる純粋状態が生成できるようになってきた[1]。この考えを元に、我々は純粋状態に基づいた有限温度VMCを開発した[2]。Hubbard模型におけるベンチマークの結果、エネルギーなどの物理量の温度依存性がTPQによる結果をよく再現することを示した。

二つ目に、動的構造因子やスペクトル関数などの動的物理量の計算方法を紹介する。VMCを用いた動的物理量の計算手法として、スピンの動的なゆらぎの情報を含んだスピン動的構造因子を計算する手法が提案されている[3]。この手法は量子スピン鎖のスピン励起を高精度に計算できることが示されており[4]、中性子散乱による実験結果との比較にも使われてきている[5]。このようにVMCを用いたスピンダイナミクスに関する研究は発展しているものの、共鳴非弾性X線散乱などで測定される電荷ダイナミクスは未だに計算できるよう拡張されていない。我々は、VMCに基づいた遍歴電子系における電荷動的構造因子を計算する手法を開発した[6]。Hubbard模型におけるベンチマークの結果、Mottギャップを生み出すと考えられている複合粒子が電荷ゆらぎを高精度に捉えるうえで重要な役割を担っていることが明らかになった。

[1] S. Sugiura and A. Shimizu, Phys. Rev. Lett. 108, 240401 (2012). M. Imada and M. Takahashi, J. Phys. Soc. Jpn. 55, 3354 (1986). J. Jaklič and P. Prelovšek, Phys. Rev. B 49, 5065 (1994). A. Hams and H. De Raedt, Phys. Rev. E 62, 4365 (2000). [2] K. Takai, K. Ido, T. Misawa, Y. Yamaji and M. Imada, J. Phys. Soc. Jpn. 85, 034601 (2016). [3] T. Li and F. Yang, Phys. Rev. B 81, 214509 (2010). [4] F. Ferrari, A. Parola, S. Sorella, and F. Becca, Phys. Rev. B 97, 235103 (2018). [5] B. Dalla Piazza, et al., Nat. Phys. 11, 62 (2015). [6] K. Ido, M. Imada, T. Misawa, Phys. Rev. B 101, 075124 (2020).

第21回

日時: 2020年10月5日(月)11:00-
発表者: 小野淳（東北大）
発表題目: 反強磁性ディラック半金属におけるネールベクトルの高速光制御

登録フォーム

反強磁性体は外部磁場の擾乱に対して頑強であることや強磁性体に比して高いエネルギースケールを持つことから，スピントロニクス分野においても精力的に研究が進められている．空間・時間反転対称性を持つ反強磁性状態における交替磁化（ネールベクトル）を制御する原理の一つとして，スピン軌道相互作用を介して定常電流から誘起される交替的な非平衡磁化（ネールスピン軌道トルク）を用いるものが理論的に提案され，実際にネールベクトル方向の電気的なスイッチングが可能であることが報告されている[1]．また，ある種の非共型な空間群対称性によりディラック点や線ノードが保護される場合があることが知られているが，これによりネールベクトルを介してトポロジカルな電子構造を操作するというアイディアが提案され[2]，トポロジカル反強磁性スピントロニクスという新たな分野が生まれている[3]．一方で，反強磁性体の典型的なエネルギースケールはテラヘルツ程度であり，これに対応するピコ秒程度の高速なダイナミクスは実験と理論共にほとんど明らかになっていない．

本研究では反強磁性ディラック半金属に注目し，直流電場ならびにテラヘルツパルスによって誘起される実時間ダイナミクスを解析した．局在的な古典スピン系のランダウ・リフシッツ・ギルバート方程式に加えて遍歴的な電子系のシュレーディンガー方程式を連立して解くことで，磁気構造と電子構造の両者を取り扱う解析を行った．閾値を超える外部電場の印加により，ピコ秒程度に相当する時間スケールでネールベクトルが回転し，これに伴ってディラック点のギャップが変化することが見出された．ここではネールスピン軌道トルクのみならず磁気異方性も重要な役割を果たしている．本セミナーではこれらの詳細について説明するとともに，磁気光学効果を用いたネールベクトル方向の高速時間分解測定の可能性についても議論する．

[1] J. Železný et al., Phys. Rev. Lett. 113, 157201 (2014); P. Wadley et al., Science 351, 587 (2016). [2] L. Šmejkal, J. Železný, J. Sinova, and T. Jungwirth, Phys. Rev. Lett. 118, 106402 (2017). [3] L. Šmejkal, Y. Mokrousov, B. Yan, and A. H. MacDonald, Nat. Phys. 14, 242 (2018).

第20回

日時: 2020年9月28日(月)11:00-
発表者: 永井恒平（京大）
発表題目: 高強度光駆動下における固体電⼦系の動的対称性とその実験的検証

登録フォーム

高強度レーザー光電場に周期的に駆動される固体中の電子状態はフロケ状態の観点から議論され、実験ではフロケ・ブロッホバンドの観測[1]や光誘起異常ホール効果の観測[2]などによってフロケ状態の存在が確認されている。フロケ状態の生成および制御は光による超高速な物性制御の実現に向けて期待されており、対称性はそのために最も重要な概念である。中でも動的対称性は時間並進、反転操作と空間対象操作を組み合わせたもので、非摂動論的な光と物質の相互作用の下で起こる高次高調波発生（HHG）の理解に有用であるとして1993年に元々のアイデアが発案され[3]、2019年にフロケ群論として定式化が行われた[4]。しかし、これまでほぼすべての動的対称性の研究がHHGの偏光選択則の理解を目的としていたものであり、また固体における動的対称性の実験的な検証も乏しかった。

我々は固体における非摂動論的な光学過程である高次サイドバンド発生(HSG)の偏光選択則の理解に対して動的対称性を応用し、その実験的検証を行った[5]。HSGは半導体に対してバンド間遷移に近共鳴な弱強度のパルス光と高強度赤外パルス光を同時照射することで生じる高次の和周波発生過程であり、2色の光を制御することでHHGより詳細に固体中の電子の励起状態を調べられるとして注目されている。これまでHSGは高強度電場下での電子の実時間発展に注目して研究されてきた[6]。我々はHSGをフロケ状態におけるラマン散乱過程として捉え直すことで、動的対称性によるHSGの偏光選択則の理解を可能にした。また、我々は実験的にHSGの偏光選択則を調べ、理論と比較した。対象の半導体として遷移金属ダイカルコゲナイドの単一原子層薄膜を用いることでHSGにおいて本質でない伝搬効果を避けて実験を行った。偏光選択則を系統的に調べた結果が理論的に予想された選択則と完全に一致し、固体における動的対称性を実証できた。

[1] Wang, Y. H., Steinberg, H., Jarillo-Herrero, P. & Gedik, N. Observation of Floquet–Bloch states on the surface of a topological insulator. Science 342, 453–457 (2013). [2] McIver, J. W. et al. Light-induced anomalous Hall effect in graphene. Nat. Phys. 16,38–41 (2020). [3] Ben-Tal, N., Moiseyev, N. & Beswick, A. The effect of Hamiltonian symmetry on generation of odd and even harmonics. J. Phys. B 26, 3017–3024 (1993). [4] Neufeld, O., Podolsky, D. & Cohen, O. Floquet group theory and its application to selection rules in harmonic generation. Nat Commun 10, 405 (2019). [5] Nagai, K., Uchida, K., Yoshikawa, N. et al. Dynamical symmetry of strongly light-driven electronic system in crystalline solids. Commun Phys 3, 137 (2020). [6] Zaks, B., Liu, R. B. & Sherwin, M. S. Experimental observation of electron–hole recollisions. Nature 483, 580–583 (2012).

第19回

日時: 2020年9月23日(水)11:00- 注:曜日がいつもと違います
発表者: 田村駿（名大工）
発表題目: カイラル対称性を持つ系における奇周波数クーパー対の理論

登録フォーム

超伝導体のクーパー対はクーパー対を成す２電子の入れ替えに関する対称性によって分類される．通常の超伝導体ではスピン一重項s波やスピン三重項p波などが考えられてきたが，2電子の時間に関する入れ替えも考慮すると，これらのクーパー対は時間の入れ替えに関して偶関数になっている．時間の入れ替えに関して奇関数になる場合についても考慮するとスピン一重項p波やスピン三重項s波などがあり得る．時間の入れ替えに関して奇関数になるようなクーパー対は奇周波数クーパー対と呼ばれ，偶周波数クーパー対の示す磁場を排斥するマイスナー効果とは逆の常磁性マイスナー効果を示すなどの特徴を持つ．奇周波数クーパー対は一般的にバルクでは不安定であるが，超伝導体の表面などには広く存在することが理論研究により示されており[1]，μSRによる実験でその存在が確認されている[2,3]．近年，トポロジカル超伝導体においてバルクの偶周波数クーパー対を超える量の奇周波数クーパー対が表面に誘起されることが理論的に示唆されており，トポロジカル数と奇周波数クーパー対の間に何らかの関係があることが予想されていた．

本セミナーでは奇周波数クーパー対に関して簡単な紹介をした後，多体相互作用のない一次元のカイラル対称性がある系に関して奇周波数クーパー対とトポロジカル数の関係について紹介する[4,5]．

[1]Y. Tanaka and N. Nagaosa, J. Phys. Soc. Jpn. 81, 011013 (2012). [2]A. Di Bernardo, et al., Phys. Rev. X 5, 041021 (2015). [3]J. A. Krieger, et al., Phys. Rev. Lett. 125, 026802 (2020). [4]S. Tamura, S. Hoshino and Y. Tanaka, Phys. Rev. B 99, 184512 (2019). [5]A. Daido and Y. Yanase, Phys. Rev. B 100, 174512 (2019).

第18回

日時: 2020年9月14日(月)11:00-
発表者: 藤原理賀（東京理科大）
発表題目: カムチャツカ半島産低次元磁性体における量子スピン状態

登録フォーム

量子スピン液体状態は，スピンフラストレーションを内包する低次元磁性体での観測が期待されている．しかしスピンフラストレーションは，格子の僅かな歪みや弱い鎖間・層間相互作用であっても解消されてしまうため，モデル物質の探索が難航する場合が多い．

無機化合物の場合，構成元素と結晶構造を自在に操作することは難しい．よって結晶構造が既知である鉱物は，低次元磁性体の良い探索領域であると言える．カムチャツカ半島には多くの火山が存在する．加えてこの地帯には銅鉱床が存在するため，新種の銅鉱物が多く発見され，その数は現在も増え続けている [1]．カムチャツカ半島産銅鉱物の結晶構造を調査したところ，スピンフラストレーションを内包し得る低次元量子磁性体が多く存在することが判明し，これらの鉱物が生成された自然環境を模倣することで，数種類の新低次元量子磁性体の開発に成功した．

本セミナーでは，正方カゴメ格子反強磁性体の初のモデル物質であるアトラス鉱 KCu6AlBiO4(SO4)5Cl [2]、辺共有四面体量子スピンクラスタ鎖と名付け研究を続けてきたフェドトフ鉱 K2Cu3O(SO4)3 [3] を中心に、それらのスピン状態を紹介したい．これらの鉱物では，極低温においても磁気秩序が形成されないため，理想的な低次元量子磁性体であるといえる．J-PARCおよびANSTOに設置されている複数の分光器を利用して行われた中性子非弾性散乱実験の結果，前者では量子スピン液体状態、後者ではハルデン状態の形成を示す実験結果が得られた。詳細は当日紹介するが、発表者が開発したその他の低次元磁性体の研究結果についても，時間の許す限り紹介したい [4, 5]。

[1] Scientific Investigations Report 2010-5090-W. [2] M. Fujihala et al., Nat. Commun. 11, 3429 (2020). [3] M. Fujihala et al., Phys. Rev. Lett. 120, 077201 (2018). [4] M. Fujihala et al., Sci. Rep. 7, 16785 (2017). [5] M. Fujihala et al., Phys. Rev. B 101, 024410 (2020).

第17回

日時: 2020年9月7日(月)11:00-
発表者: 小野清志郎（東大）
発表題目: 対称性指標に基づくトポロジカル超伝導相の分類と物質探索への応用

登録フォーム

マヨラナ粒子のプラットフォームであるトポロジカル超伝導体は、今なお精力的に研究されている。超伝導体のトポロジーを研究するもう1つの意義は、トポロジーと超伝導相の性質に密接な関係があることだ。例えば、安定な超伝導ノード構造がトポロジーによって系統的に理解されている[1]。ここ数年、結晶対称性とトポロジーの研究には急速な進展があった。これまでは、トポロジカル不変量をブロッホ関数の空間反転対称性の固有値で計算できるFu-Kane公式[2]がその簡便さから重宝されていた。近年その一般化として、結晶対称性の既約表現からトポロジカル相を判定する対称性指標[3]が提案された。さらにこの理論を密度汎関数計算と組み合わせることで、データベースのトポロジカル(結晶)絶縁体探索が行われトポロジカル絶縁体・半金属データベースが完成した[4]。最近の研究[5]により、超伝導体の対称性指標は既約表現だけでなく、Pfaffian不変量が重要な役割を果たすことが明らかになった。しかしPfaffian不変量は0、または1に値を取るため、従来の方法で系統的に計算することが困難であった。実際、これまではごくわずかの対称性に対してのみ個別に計算されていた。

そこで我々は数学的に等価な補助問題を導入することで、Pfaffian不変量を含めて対称性指標を系統的に計算できるように定式化し、磁気空間群の1次元表現に属する秩序変数を持つ超伝導体全てに対して対称性指標を計算した[6]。本セミナーでは、絶縁体の対称性指標をレビューしたのちに、トポロジカルに自明な超伝導体を再考する[7]。最後に、Pfaffian不変量を含めた対称性指標の定式化するための補助問題を説明し、物質への応用を紹介する。

[1] 例えば、S. Kobayashi, et al., Phys. Rev. B 97, 180504 (2018)., S. Sumita, et al., Phys. Rev. B 99, 134513 (2019). [2] L. Fu and C. L. Kane, Phys. Rev. B. 76, 045302 (2007). [3] H. C. Po, A. Vishwanath, and H. Watanabe, Nat. Commun. 8, 50 (2017). (See also B. Bradlyn et al., Nature, 566, 298-305 (2017)) [4] T. Zhang, et al., Nature, 566, 475-479 (2019)., M. G. Vergoniry, et al., Nature, 566, 480-485 (2019)., F. Tang, et al., Nature, 566, 486-489 (2019). [5] K. Shiozaki, arXiv:1907. 13632., M. Geier, et al., Phys. Rev. B 101, 245128 (2020). [6] S. Ono, H. C. Po, and K. Shiozaki, arXiv:2008.05499. [7] A. Skurativska, T. Neupert, and M. H. Fischer, Phys. Rev. Research 1, 013012 (2019)., S. Ono, H. C. Po, and H. Watanabe, Sci. Adv. 6, eaaz8367 (2020).

第16回

日時: 2020年8月24日(月)11:00-
発表者: 野村悠祐（理研）
発表題目: 機械学習手法を用いた量子多体系の研究 ー手法拡張を含む最近の進展

登録フォーム

指数関数的に大きな次元を持つ量子多体系の波動関数を有限個のパラメータで精度よく表すことは、物性物理のみならず素粒子、原子核、量子化学などに共通するグランドチャレンジである。本講演では機械学習で用いられるボルツマンマシンが量子多体波動関数表現に有用であることを紹介する。

前半では、隠れ層が二層ある深層ボルツマンマシン(DBM)を用いると、基底状態の波動関数を任意の精度で解析的に表現することが可能であることを示す[1]。これは経路積分を包含するより一般的な量子古典対応のフレームワークを提供するが、物理量計算に隠れ層自由度をモンテカルロサンプリングする必要性があるため負符号問題が生じる場合がある[1]。

一方で、隠れ層が一層に制限された制限ボルツマンマシン(RBM)では、DBMのように解析的な表現はできないものの、数値的にパラメータを最適化することで、量子多体波動関数の表現が可能になる。RBMの良いところは解析的に隠れ層自由度をトレースアウトできるところであり、量子系の数値計算を行う上では負符号問題を回避できる点でDBMより利点がある場合がある。RBMを用いた手法は、フェルミオン系、フェルミオン-ボソン結合系、フラストレーションのあるスピン系な度への適用が進んでおり、講演の後半部分でそれら一連の研究について紹介する[2-4]。

[1] G. Carleo, Y. Nomura, and M. Imada, Nat. Commun. 9, 5322 (2018) [2] Y. Nomura, A. S. Darmawan, Y. Yamaji, and M. Imada, Phys. Rev. B 96, 205152 (2017) [3] Y. Nomura, J. Phys. Soc. Jpn. 89, 054706 (2020) [Editor’s choice] [4] Y. Nomura and M. Imada, arXiv:2005.14142

第15回

日時: 2020年8月17日(月)11:00-
発表者: 渡邉光（京大）
発表題目: 磁気的なパリティの破れに誘起される円偏光ガルバニック効果

登録フォーム

光によって直流電流を生み出す光電流応答は太陽電池や光検出器などに応用されており、 我々の実生活にとって重要な物性応答といえる。交流外場によって直流電流を取り出すため、 応答には反転対称性の破れが必要とされる。典型例としての半導体pn接合では作り付けの内部電界が その役割を担っているが、一方で結晶のミクロな構造も同様の役割を果たす。この、ミクロな 反転対称性の破れに誘起される光電流、いわゆるバルク光電流の理解は急速に進められており、 強誘電体などを題材に多くの研究がなされている。また、TaAsなどのトポロジカル物質と組み合わせる ことで巨大な光電流応答が得られることが分かっており、近年のトポロジカル物性の発展と 相まって、光電流応答の研究は更なる展開を見せている [1] 。

これら目覚ましい発展がある一方、光電流を調べるセットアップは非磁性体がほとんどであった。 マルチフェロイックスや反強磁性スピントロニクス分野で議論されるように [2] 、 反転対称性は 磁気秩序によっても破られるが、そのような磁性誘起の光電流応答に着目した研究は少ない。 非磁性的・磁気的な反転対称性の破れが系の電子状態に与える影響は定性的に全く異なるから、 光電流応答にもその違いが色濃く反映されるものと期待できる。

そこで我々は、非磁性的・磁気的な反転対称性の破れを持つ系の基本的な対称性を利用し、 光電流応答を系統的に分類した [3] 。その結果として、これら二つのタイプの破れが対照的な役割を果たす ことを明らかにすることができた。分類結果は先行研究を包括するだけでなく、磁性体に特有な新しい光電流の クラスを明らかにしている。本セミナーでは新たな応答として、円偏光によって誘起される光電流 について議論する。この応答は、非磁性体に特徴的ないわゆるシフトカレントと双対的な関係にある。 時間があれば、この新奇な光電流応答がトポロジカル物質において大きく増強される可能性についても議論したい。

[1] 例えば. J. W. McIver, D. Hsieh, H. Steinberg, P. Jarillo-Herrero, & N. Gedik, Nat. Nanotech. 7, 96 (2011). [2] M. Fiebig, J. Phys. D. Appl. Phys. 38, R123 (2005); V. Baltz, A. Manchon, M. Tsoi, T. Moriyama, T. Ono, and Y. Tserkovnyak, Rev. Mod. Phys. 90, 015005 (2018). [3] HW & Y. Yanase, arXiv:2006.06908; 同様の結果として J. Ahn & N. Nagaosa, arXiv:2006.06709

第14回

日時: 2020年8月3日(月)11:00-
発表者: 原口祐哉（東京農工大）
発表題目: 新規ハニカム格子化合物における量子物性の探索

登録フォーム

ハニカム格子は、最近接・次近接相互作用の拮抗(J_1-J_2ハニカム)やボンド方向に依存したイジング的な相互作用(キタエフ相互作用)に由来するフラストレート物理の舞台として近年精力的に研究が行われている。これらのモデル物質となる化合物は少ないため、完全な物性解明には至っていないのが現状である。無機化合物の探索には固相反応法は用いることがほとんどであるが、単純な固相反応法ではすでに莫大な物質が調べられている上に熱力学的に安定な典型的な構造を取る。一方、低温メタセシス反応は非平衡プロセスが多用されているため、相図の枠組みを超えた準安定相を合成することができるという利点をもつ。

本公演では、登壇者が開発した層状岩塩型化合物からイルメナイト型への低温メタセシス反応を用いて合成した5つの新規ハニカム格子磁性体MgMnO3, ZnMnO3, MgIrO3, ZnIrO3, CdIrO3における合成および実験的物性研究について紹介する。これらの化合物はいずれも500°C程度で分解する準安定相であることが確かめられた。マンガンハニカムについて、MgMnO3ではJ_1およびJ_2がいずれも反強磁性である一方ZnMnO3ではJ_1が強磁性・J_2が反強磁性となり、いずれもJ_1-J_2フラストレートハニカム格子になっていることがわかった。またイリジウムハニカムでは、三種類の化合物の結晶構造と磁性の変化の系統性から、Irイオン周りの三方晶歪の増大にともないキタエフ磁性から遠ざかっていることが明らかになった。これらの結果から、イルメナイト型酸化物における理想的なJ_1-J_2フラストレート磁性およびキタエフ磁性の実現可能性とそれに向けた物質設計の指針について提案したい。

第13回

日時: 2020年7月27日(月)11:00-
発表者: 青山和司（阪大）
発表題目: 磁性絶縁体のスピン伝導とトポロジカル転移

登録フォーム

2次元の磁性体において容易面(XY)的磁気異方性がある場合には、トポロジカル欠陥である渦励起の存在に伴いKosterlitz-Thouless(KT) 転移が生じる。一方、2次元のフラストレート磁性体の 典型例である三角格子上の古典Heisenberg反強磁性体においては、スピンの120°構造に起因したZ2渦と呼ばれる特殊渦の出現とそのKTタイプのトポロジカル転移（渦対-単独渦転移）の存在が川村・ 宮下によって理論的に提案されている[1]。今回我々は、これらのトポロジカル転移が系の輸送現象にどのような影響を与えるのかを調べた。

線形応答理論に基づく数値解析により、スピン流の伝導率がトポロジカル転移温度において発散的な増大を示すことが明らかとなった[2]。Z2渦転移は比熱や帯磁率などに明確な異常が現れないため検出が困難であったが、我々の結果は、スピン流測定がトポロジカル転移を検出する有力なプローブになることを示唆している。

[1] H. Kawamura and S. Miyashita, J. Phys. Soc. Jpn. 53, 4138 (1984). [2] K. Aoyama and H. Kawamura, Phys. Rev. Lett. 124, 047202 (2020); Phys. Rev. B 100, 144416 (2019).

第12回

日時: 2020年7月21日(火)11:00- 注:曜日がいつもと違います
発表者: 岩木惇司（東大）
発表題目: 熱的量子純粋行列積状態

登録フォーム

比熱や磁化率の温度依存性は、系の低エネルギー状態の情報を含む基本的な物理量であり、実験との比較を通して理論モデルを検証する上でも重要な役割を果たす。しかし、相互作用の強い量子多体系において、これらを計算する汎用的で簡便な方法がないことが長年の課題である。例えば、有限温度対角化の手法はシステムサイズが高々30程度までに限られ、量子モンテカルロ法には負符号問題がある。困難の主な原因は、有限温度の状態がギブズ状態と呼ばれる量子混合状態であることにある。

混合状態のアンサンブル平均を取ることを回避する方向性の一つとして、典型性という考え方を元にギブズ状態を熱的純粋量子（TPQ）状態と呼ばれる純粋状態で代表させる方法があり、TPQ状態を計算機上で比較的容易に生成する有効な方法論がいくつか提案されている[1]。ところがこの方法では純粋状態を表現するために扱う変数の数がシステムサイズに対して指数関数的に大きくなるため計算コストの節約にはならず、やはりシステムサイズが高々30程度しか扱えない。システムサイズを大きくするためには何らかの近似が必要であるが、有限温度の純粋状態のエンタングルメントエントロピーが体積則に従うことがこれを難しくしている。システムサイズに対して多項式的な数の変数で純粋状態を記述する行列積状態（MPS）は、歴史的には密度行列繰り込み群の変分関数として導入され[2]、エンタングルメントの小さい基底状態の物理量を簡便な計算で求めることに成功している。

本研究[3]では、あえてMPSを用いてTPQ状態を構築することを試みた。MPS単独では有限温度量子純粋状態を表現できないだろうというナイーブな予想に反し、今回、MPSに補助系（auxiliary）を付けることにより、系が保持するべき基底の数を確保し、ほぼ基底状態と同じ計算コストで諸物理量を得ることに成功したので、その詳細について述べたい。

[1]A. Hams and H. De Raedt, Phys. Rev. E 62, 4365 (2000); S. Sugiura and A. Shimizu, Phys. Rev. Lett. 108, 240401 (2012) [2]J. Dukelsky, M. A. Martín-Delgado, T. Nishino and G. Sierra, Europhys. Lett. 43, 457 (1998) [3]A. Iwaki, A. Shimizu and C. Hotta, arXiv:2005.06829

第11回

日時: 2020年7月13日(月)11:00-
発表者: 野本拓也（東大）
発表題目: 第一原理計算に基づくGd系スカーミオン物質の研究

登録フォーム

近年Gd2PdSi3[1]、Gd3Ru4Al12[2]およびGdRu2Si2[3]などのGd系化合物でスカーミオン相が観測され、注目を集めている。これらの系は常磁性相で空間反転対称性を保ち、実現するスカーミオン格子が数ナノメートルと小さな格子間隔をもつ点で、他のスカーミオン物質と一線を画している。現象論的な模型に基づく解析により、磁気フラストレーション[4]や高次スピン間相互作用[5]がこのようなスカーミオン格子を安定化させることが知られているが、格子間隔の小さな高密度スカーミオンは応用上も重要であり、第一原理計算に基づく形成機構の微視的な解明が望まれる。

本講演では、我々が最近行ったDFT+Lichtenstein法による、GdRu2Si2およびGd2PdSi3の磁気構造解析の結果について紹介する[6]。我々の計算は、磁気スパイラル相に関して、実験とコンシステントな転移温度および磁気格子間隔を与えるが、その発現機構は従来考えられてきたRKKY機構とは大きく異なっている。我々の解析によれば、有限波数の磁気構造は強磁性的なGd-5d間相互作用と反強磁性的なGd-4f間相互作用の競合で理解され、通常のRKKY相互作用は無視できるほど小さい。講演では第一原理計算に基づく磁気構造の計算手法と我々の研究の詳細について紹介する。

[1] T. Kurumaji et al., Science 365, 914 (2019). [2] M. Hirschberger et al., Nat. Commun. 10, 5831 (2019). [3] N. D. Khanh et al., Nat. Nanotechnol. 15, 444 (2020). [4] A. O. Leonov and M. Mostovoy, Nat. Commun. 6, 8275 (2015). [5] S. Hayami, R. Ozawa, and Y. Motome, Phys. Rev. B 95, 224424 (2017). [6] T. Nomoto, T. Koretsune, and R. Arita, arXiv:2003.13167.

第10回

日時: 2020年7月6日(月)11:00-
発表者: 松下太樹（阪大）
発表題目: 線ノード半金属におけるトポロジカルな圧電効果

登録フォーム

線ノード半金属は１次元のFermi線を有するトポロジカル半金属である。この物質固有の物性として、Drumhead表面状態や半量子化した電気分極、軌道磁化が知られており、精力的に研究されてきた[1,2]。しかし、この物質を特徴付ける応答や輸送現象は知られておらず、この物質のトポロジカル性がどのように物理現象に顔を出すかは未解決な問題として残されていた。 　最近、電場応答で線ノード半金属のトポロジカル性が現れるかが検討された[3]。解析の結果、トポロジカルカレントは、Fermi線を構成する各二次元Dirac点の寄与を足し上げるとキャンセルされてしまうことが示されている。我々は、電気応答ではなくAxial電気応答を考えることで、トポロジカルカレントの相殺が避けられることに着目し研究を行った[4]。このセミナーではその成果について紹介する。Axial電場を得るための具体的な処方箋を示し、動的な格子歪みによってトポロジカルカレントが誘起できることを示す。また、動的な格子歪みによって誘起されるトポロジカルカレントは、線ノード半金属に存在する半量子化値をとる電気分極の運ぶ分極電流であることについても議論する。 [1] A. A. Burkov, M. D. Hook, and L. Balents, Phys. Rev. B 84, 235126 (2011). [2] S. T. Ramamurthy and T. L. Hughes, Phys. Rev. B 95, 075138 (2017). [3] W. B. Rui, Y. X. Zhao, and A. P. Schnyder, Phys. Rev. B 97, 161113 (2018). [4] TM, S. Fujimoto and A. P. Schnyder : arXiv:2004.11014(2020).

第9回

日時: 2020年6月29日(月)11:00-
発表者: 藤代有絵子 （東大）
発表題目: 高密度なトポロジカル磁気構造における巨大創発磁場と量子輸送現象

登録フォーム

トポロジカルに非自明な磁気構造は、伝導電子との結合によりゲージ場である創発磁場を生み出し、様々な電磁気応答をもたらす。特に最近では、磁気周期の短い（< 3 nm）高密度なスキルミオン格子・ヘッジホッグ格子が多くの物質で発見されており、巨大な創発磁場の舞台として注目されている。一方、2次元的なスキルミオン格子と３次元的なヘッジホッグ格子では、創発磁場の空間分布や外部磁場への応答が異なるため、劇的に異なる物性がこれまでに観測されてきた[1]。

本セミナーでは、カイラル磁性体MnGeにおける高密度なヘッジホッグ格子に着目し、巨大な創発磁場[2]とそのダイナミクス[3]がもたらす特異な量子輸送現象について、これまでの実験結果を紹介する。特に、ヘッジホッグ格子が磁場によって消失する強磁性転移付近では、既存の概念では説明のできないような、非自明な輸送現象が数多く観測されている。例えば、巨大な磁気ゼーベック効果[4]やトポロジカルホール効果の符号反転に加え、最近では既存の強磁性体の中でも最大級の異常ホール効果（ホール伝導度40,000 Ω-1cm-1、ホール角20 %）が出現することを発見した。強磁場測定を中心とした様々な検証実験と、関連する理論研究[5]の内容を紹介しつつ、ヘッジホッグ格子の消失に伴う巨大な磁気揺らぎという視点から、これらの実験結果を議論したい。

[1] Y. Fujishiro, N. Kanazawa, and Y. Tokura, Applied Physics Letters 116, 090501 (2020) [2] N. Kanazawa et al., Physical Review Letters 106, 156603 (2011) [3] N. Kanazawa et al., Nature Communications 7, 11622 (2016) [4] Y. Fujishiro et al., Nature Communications 9, 408 (2018) [5] H. Ishizuka and N. Nagaosa, arXiv:1906.06501

第8回

日時: 2020年6月22日(月)11:00-
発表者: 明石遼介 （東大）
発表題目: 第一原理計算に基づく高圧下における水素化物超伝導体の研究

登録フォーム

2014年末に報告された200ケルビンにおける硫化水素の超伝導転移を契機として、超高圧下における水素化合物の超伝導相の探索が近年盛んに進められている[1,2]。BCS理論によれば、フォノン機構による典型的な超伝導体の転移温度は電子間引力を媒介するフォノンの振動数に比例する。ゆえに最も軽い元素である水素を含む化合物を金属化できれば高温超伝導が実現する、という仮説は古くから提唱されており[3]、近年の発見はこの仮説を肯定的に証明したものとみられる。しかし、実際に圧力セルの中でどのような結晶相が形成されているかを確かめることは難しい。光や中性子などを用いた回折法により圧力セル内の水素の位置を決定することが困難なためである。一方、第一原理計算ならば任意の圧力下で熱力学的に安定な構造を探索し、電子・格子物性をシミュレートすることが可能である。どのような構造ならば測定された物性値を再現できるか、というアプローチから、第一原理シミュレーションはセル内の超伝導相について重要な示唆を与えてきた。

本講演では、我々が第一原理計算に基づき解明した高圧下水素化物の特殊な性質について紹介する。主要なトピックとしては、硫化水素金属相の(i)安定な結晶構造に類似した無数の準安定構造系列の存在 [4]、および(ii)電子状態がフェルミレベル近傍に特異的に集中し、ペアリングを増強する仕組み[5]について述べる予定である。

[1] A. P. Drozdov et al., Nature 525, 73 (2015). [2] J. A. Flores-Livas et al., Phys. Rep. 856, 1 (2020), for a review. [3] N. W. Ashcroft, Phys. Rev. Lett. 21, 1748; 92, 187002 (2004). [4] R. Akashi et al., Phys. Rev. Lett. 117, 075503 (2016). [5] R. Akashi, Phys. Rev. B 101, 075126 (2020).

第7回

日時: 2020年6月15日(月)11:00-
発表者: 道下佳寛 (京大）
発表題目: 強相関電子系における量子開放系としての性質及び非エルミート性

登録フォーム

近年、非エルミート有効ハミルトニアンで記述される新たな物性が注目を集めている。 例えば特異な臨界現象[1]、センサー感度の向上[2]、光の透過の非相反性の創出[3]、新 奇トポロジカル相[4,5], skin effect[6]などがある。非エルミート物性は主に二つの文脈、量子開放系と強相関電子系、において研究がなさ れている。量子開放系においては、ポストセレクションという観測結果に対する制限下で の密度演算子のダイナミクスを非エルミートハミルトニアンで記述する事ができ、また強 相関電子系においては、一体の Green 関数を記述するものとして特に条件を貸す事なく自 然に導入される。この二つの文脈で独立に研究が進んでいく一方で、両者での非エルミー トハミルトニアンの導入のされ方、導入される際の条件は異なっており、そもそも両者で 本当に同じものを考えているのか明らかではなかった。

そこで、本研究[7]では、強相関電子系を量子開放系として記述する事で、同じ模型に おける両者の文脈から得られる非エルミートハミルトニアンを比較することを目指す。そ の際、量子マスター方程式のダイナミクスは、Green 関数法におけるダイアグラムの手法 を用いて計算できることを示し、また結果として、同じ模型から出発した際に、量子開放 系の文脈で非マルコフ過程のもとで得られた非エルミートハミルトニアンと、強相関電子 系の文脈(Green 関数法の下)で得られたものは一致する事を示す。また、固体物質中での 電子のマルコフ的なダイナミクスと非マルコフ的なダイナミクスの、エネルギー分散(ス ペクトル関数)への寄与も明らかにする。

[1] Y. Ashida et al. , Nature communications 8, 15791 (2017). [2] W. Cheng et al. , Nature 548, 192 (2017). [3] A. Regensburger, et al, Nature 488, 167-171 (2012). [4] H. Shen et al. , PRL 120, 146402 (2018). [5] K. Kawabata et al. ,PRX 9, 041015 (2019). [6] S. Yao and Z. Wang PRL 121,086803 (2018). [7] Y. Michishita and R. Peters PRL 124,196401 (2020)

第6回

日時: 2020年6月8日(月)11:00-
発表者: 藤本純治 (中国科学院大学Kavli理論科学研究所）
発表題目: 動的ねじれによるマグノン流生成

登録フォーム

磁気と力学運動の結合に対する研究の歴史は長く，Einstein-de Haas効果などの磁気機械効果は現在も盛んに研究されている．特に，科学技術の進歩によってナノスケールの力学装置とスピンの結合に関心が向けられるようになった[1]．これまでは半導体や金属などのナノメカニカルな装置が主流であったが，最近，磁性絶縁体のカンチレバーがサブミクロンの大きさで作製されるようになり[2]，その磁気機械効果[3]に注目が集まっている． そこで，我々はそのようなナノメカニカルな強磁性絶縁体カンチレバーを考え，外的にねじれ振動が誘起されたときマグノン流が生成されるという新たな磁気機械結合を示した[4]．具体的には，交換相互作用と容易軸異方性のみで構成されるシンプルな局在スピンのラグランジアンを仮定し，力学的なねじれの効果を容易軸の変調とみなすことで取り入れた．そしてスピンに対する局所的な回転変換を施すことで，力学的なねじれがスピンゲージ場として記述され，Dzyaloshinskii-Moriya(DM)相互作用を生むことを示した．このDM相互作用を外場としてマグノン流の線形応答を計算することで，ねじれ振動によるマグノン流生成を記述した．また，生成されたマグノン流は，逆スピンHall効果で検出可能な大きさであることを確認した．

[1] P. Mohanty et al., Phys. Rev. B 70, 195301 (2004); G. Zolfagharkhani et al., Nat. Nanotechnol. 3, 720 (2008). [2] Y.-J. Seo et al., Appl. Phys. Lett. 110, 132409 (2017). [3] K. Harii et al., Nat. Commun. 10, 2616 (2019). [4] J. Fujimoto, and M. Matsuo, arXiv:2004.14707.

第5回

日時: 2020年6月1日(月)11:00-
発表者: 山田昌彦 (阪大)
発表題目: アンダーソン・キタエフスピン液体

登録フォーム

半世紀以上にわたってアンダーソン局在の研究は行われてきたが、それは主に低次元の電子系やフェルミオン系においてである。磁気秩序を持った通常の磁性体においてはフェルミオンの励起が存在せず、アンダーソン局在とは違った枠組みで研究されていた。一方で、基礎励起がマヨラナフェルミオンで記述されるキタエフスピン液体においては状況が違っている。私はボンド乱れのあるキタエフ模型を研究し、キタエフスピン液体が乱れの増大に伴いアンダーソン局在を起こすことを発見した。カーネルポリノミアル法に基づく基底状態計算によると、少ない量の乱れでもマヨラナフェルミオンを局在化させるのに十分であり、量子化された熱ホール効果が次第に消え去ることがわかった。局在のメカニズムそのものが今までのアンダーソン局在とは違ったタイプのものであり、このことについても詳しく議論していきたい。

[1] M.G.Y., https://arxiv.org/abs/2004.06257

第4回

日時: 2020年5月25日(月)11:00-
発表者: 北谷基治 (理研)
発表題目: 動的バーテックス近似(DΓA)を用いた非従来型超伝導の研究

登録フォーム

動的平均場理論(DMFT)は相関電子系を扱う強力な多体手法であるが、自己エネルギーの波数依存性を無視してしまうために擬ギャップや異方的超伝導などの物理を記述できない。この問題点を克服するために、非局所的な揺らぎの効果を取り込む拡張が近年盛んに議論されている。

今回はDMFTのダイアグラム展開に基づく拡張[1]の一種である動的バーテックス近似(DΓA)を用いた非従来型超伝導に関する最近の我々の研究を紹介する。２次元Hubbard模型での電子スペクトルや超伝導転移温度の計算を通して、DMFTで得られる局所的なバーテックス関数の動的な構造(振動数依存性)が非局所的な揺らぎを通してペアリング相互作用にどのような構造をもたらすかを議論する[2]。また、最近発見されたニッケル酸化物超伝導に関する計算結果[3]も紹介する。

[1] For a review, see G. Rohringer et al., RMP 90 025003 (2018). [2] M. Kitatani et al., PRB 99 041115(R) (2019). [3] M. Kitatani et al., arXiv:2002.12230 (2020).

第3回

日時: 2020年5月18日(月)11:00-
発表者: 奥村駿 (東大)
発表題目: 空間反転対称性の破れた金属における磁気ヘッジホッグ格子の安定性とトポロジカル転移

登録フォーム

近年、空間反転対称性が破れた系に現れるカイラルな磁気テクスチャが注目を集めている。とりわけ、B20型の結晶構造をもつMnSi1-xGexでは、2次元的な磁気スキルミオン結晶や3次元的な磁気ヘッジホッグ格子といったトポロジカルに特徴的な磁気テクスチャが実現することが知られている[1]。特に、磁気ヘッジホッグ格子では、非共面的な磁気テクスチャの生み出す創発磁場が磁気モノポールと反モノポールを形成し、それに伴う巨大異常ホール効果や熱電効果などの特異な量子輸送現象を引き起こすことが観測されている。 こうした磁気ヘッジホッグ格子に関する理論研究としては、局在スピン模型を用いた計算が行われているが[2,3]、実験で見られるように低温かつゼロ磁場で磁気ヘッジホッグ格子が安定に存在することを説明できていない。さらに、連続体近似において磁気構造を仮定することで、磁気モノポールと反モノポールの磁場中での振る舞いも調べられているが[4]、離散的な格子の影響や、磁場中でのエネルギー最適化を考慮した解析が必要である。

本研究では、磁気ヘッジホッグ格子の安定化機構や外場に対する応答を調べることを目的とした研究を行ってきた。対象とする物質が金属的な性質を持つことに着目し、遍歴電子の効果とスピン軌道相互作用を取り入れた有効スピン模型の解析を行った。変分計算と焼きなまし法による数値シミュレーションの結果、ゼロ磁場での基底状態において磁気ヘッジホッグ格子が発現しうることを見出した。基底状態相図を解明することにより、短周期の磁気ヘッジホッグ格子を安定化させるためには、スピン電荷結合による双二次相互作用とスピン軌道結合によるDzyaloshinskii-Moriya相互作用の両方が重要であることを明らかにした。さらに、これらの磁気ヘッジホッグ格子に磁場を印加すると、磁気モノポールと反モノポールの対消滅による多段的なトポロジカル転移が起こりうることも明らかにした[5]。

[1] Y. Fujishiro et al., Appl. Phys. Lett. 116, 090501 (2020). [2] J. Park et al., Phys. Rev. B 83, 184406 (2011). [3] S. Yang et al., Phys. Rev. B 94, 054420 (2016). [4] X. Zhang et al., Phys. Rev. B 94, 174428 (2016). [5] S. Okumura et al., Phys. Rev. B 101, 144416 (2020).

第2回

日時: 2020年5月11日(月)11:00-
発表者: 植木輝 (北大)
発表題目: 第二種超伝導体の渦状態におけるホール効果

登録フォーム

銅酸化物超伝導体における熱ホール係数の符号反転が観測されて以来、第二種超伝導体の渦状態におけるホール効果の研究は実験的にも理論的にも盛んに行われてきた。しかし現在においても、私たちの渦状態におけるホール効果に関する理解が十分であるとは言い難い。その原因のひとつとして、通常、磁場中の超伝導体を研究するために用いられる、アイレンバーガー方程式（準古典方程式）にホール効果を記述するための項が欠落していることがあげられる。最近、私たちはゴルコフ方程式から出発し、ホール効果を記述する力を含むようにアイレンバーガー方程式を拡張した[1]。この拡張準古典方程式は銅酸化物超伝導体やCeCoIn5などでみられる、熱ホール角の増大についても研究できる。今回のセミナーでは、拡張準古典方程式を解くことによって得られた次の結果、(1)s波超伝導体における渦の帯電効果[1]、(2)渦の運動が伴う抵抗状態におけるホール効果[2]、(3)ピン止めされた渦をもつd波超伝導体におけるホール角の増大[3]、について報告する。

[1] H. Ueki, M. Ohuchi, and T. Kita J. Phys. Soc. Jpn. 87, 044704 (2018). [2] H. Ueki, Doctral thesis, Hokkaido University, 2017. [3] H. Ueki, H. Morita, M. Ohuchi, and T. Kita, arXiv:1903.10733.

第1回

日時: 2020年4月27日(月)11:00-
発表者: 光元亨汰 (豊田理研)
発表題目: 凍結した乱れのない系におけるスピングラス転移

スピングラス転移は、”quenched disorder”と、そのdisorderによる”フラストレーション”の効果によって引き起こされると考えられている。しかし一方で、Y2Mo2O7に代表される、化学的な乱れのない幾何学的フラストレート磁性体においても、スピングラス転移の存在が実験的に報告されており[1]、このquenched disorderなしでのスピングラス転移の起源は数十年に渡り未解決問題として残っていた。近年、Y2Mo2O7において、MoイオンのJahn-Teller歪みの存在が実験的に示唆されている[2]。我々は、この実験結果に基づき、格子歪みを動的自由度として捉えたスピン模型を考案した。MCシミュレーションの結果、この模型において、スピンと格子が同時に凍結することを見出した[3]。さらにはこの模型を無限大次元に拡張した平均場模型においても同時ガラス転移が起こる得ることがわかったので、これらについて議論する。

[1] M. J. P. Gingras et al. PRL (1997) [2] P. M. M. Thygesen et al. PRL (2017) [3] K. Mitsumoto, C. Hotta, H. Yoshino PRL (2020)