藤堂研究室へようこそ

日本国内においても、高性能な物質科学シミュレーションソフトウェアが数多く開発・公開されているが、その知名度は必ずしも高くない。また、ドキュメントの作成やユーザサポートにも問題が多く、普及の妨げとなっている。物質科学アプリケーションのさらなる公開・普及を目指し、物質科学シミュレーションのポータルサイト「MateriApps」の整備を行っている。また、気軽にシミュレーションを始めることのできる環境構築を目指し、仮想Linuxシステム「MateriApps LIVE!」、MateriAppsアプリケーションのインストールスクリプト集「MateriApps Installer」の開発・公開も進めている。

• MateriApps - 物質科学シミュレーションのポータルサイト, MateriApps LIVE!, MateriApps Installer
• 井戸康太, 加藤岳生, 三澤貴宏, 藤堂眞治, ソフトウェア紹介「物質科学シミュレーションのポータルサイトMateriApps ― 分子動力学シミュレーションのはじめの一歩, 分子シミュレーション学会誌 “アンサンブル” 22, 74-80 (2020).
• 本山裕一, 三澤貴宏, 加藤岳生, 藤堂眞治, 物質科学シミュレーションのポータルサイト MateriApps, 固体物理 52, 743-755 (2017).
• Yusuke Konishi, Ryo Igarashi, Shunsuke Kasamatsu, Takeo Kato, Naoki Kawashima, Tsutomu Kawatsu, Hikaru Kouta, Masashi Noda, Shoichi Sasaki, Yayoi Terada, Syng e Todo, Shigehiro Tsuchida, Kazuyoshi Yoshimi, Kanako Yoshizawa, MateriApps --- a Portal Site of Materials Science Simulation, JPS Conf. Proc. 5, 011007 (2015).
• 平成31年度文部科学大臣表彰 http://www.mext.go.jp/b_menu/houdou/31/04/1415044.htm
• 理学系研究科のページ https://www.s.u-tokyo.ac.jp/ja/info/6354/

モンテカルロ法では和をとるべき「状態」に制約があり、制約を満たしながらサンプリングすることが難しい場合があります。（例として、充足可能性問題。）そのような状況でうまく状態を変えてサンプリングする方法として代表的なのが、ワームアルゴリズム（worm algorithm）と呼ばれる手法です。ワームとは制約を破るキンク（点）のことで、ワームアルゴリズムのアイデアは「制約を常に満たすのは難しいから、いったん制約を破ってしまって後でつじつまを合わせよう」というものです。このモンテカルロ法では制約を破るワーム（キンク）を導入して、ワームを確率的に動かすことでサンプリングを行います。従来の方法では、ワームをほぼ完全にランダムに動かすようにしていました。そこで我々は効率の良い計算をするためのワームの動かし方に関する指針——できるだけ前へ進めという指針——を提案しました。これを可能にするためには、ワームが動く確率を最適化する必要があるのですが、以前我々が開発した確率最適化アルゴリズム Phys. Rev. Lett. 105, 120603 (2010) を使うと簡単に実装できます。そうして最適化したワームアルゴリズムを制約つき問題に書き直したイジングモデルに応用し、従来のワームアルゴリズムと比較して計算効率を２５倍改善しました。これはイジングモデルに対して最も効率的だと信じられているクラスターアップデートよりもさらに効率的です。

• Hidemaro Suwa, Geometric allocation approach to accelerating directed worm algorithm, Phys. Rev. E 103, 013308 (2021). (preprint: arXiv:1703.03136)

近年注目を浴びている 5d 軌道の電子系では電子間相互作用とスピン軌道相互作用の両方が重要になります。我々はスピンと軌道が強く相互作用する場合、電子間のクーロン相互作用が励起子（電子と正孔の結合状態）を安定化させボーズ・アインシュタイン凝縮を引き起こすことを示しました。このように励起子のボーズ・アインシュタイン凝縮により生じる絶縁体相は励起子絶縁体と呼ばれ、現実の物質ではこれまでほとんど見つかっていませんでした。今回提案した理論を 5d 軌道の電子系であるイリジウム酸化物に適用し、２層系が長い間探されてきた励起子絶縁体であることを予言しました。またこの２層系の物質は、励起子がボーズ・アインシュタイン凝縮を起こす量子相転移のすぐ近くにあることを理論的に見出しました。相対論の効果でスピンと軌道が強く相互作用する 5d 電子系では、電子間のクーロン相互作用と電子の運動エネルギーが同程度であることが本質的に重要となることを示しました。

• Hidemaro Suwa, Shang-Shun Zhang, Cristian D. Batista, Exciton condensation in bilayer spin-orbit insulator, Phys. Rev. Research 3, 013224 (2021).

テンソルネットワークを用いた実空間くりこみ群法は、多体スピン系に対する数値計算手法として近年広く用いられてきている。テンソルくりこみ群法では大規模な古典系/量子系の物理量を効率的に計算することができる。しかしながら、Tensor Renormalization Group (TRG)やHigher-Order Tensor Renormalization Group (HOTRG)といった既存の手法では、高次元になるほど計算量が膨大になってしまうという問題点があった。この問題を解決するために、既存手法に比べて計算量が低い手法や、同じ計算時間でより精度の高い複数の手法を開発してきた。我々の提案したAnisotropic Tensor Renormalization Group (ATRG)では、3次元量子系など高次元系における計算量を劇的に減らすことができる。さらに、テンソルネットワークをテンソルを頂点上だけでなく頂点を結ぶ線上にも置いた形に拡張することで、実空間くりこみ法の精度を、既存の同程度の計算時間を要する数値計算手法に比べて100倍程度高めることに成功した。この手法は一般のテンソルネットワークに対して適用可能である。また、ランダム系やフラストレート古典スピン系へのテンソルネットワークの応用についても研究を進めている

• Daiki Adachi, Tsuyoshi Okubo, Synge Todo, Anisotropic tensor renormalization group, Phys. Rev. B 102, 054432 (7pp) (2020). (preprint: arXiv:1906.02007)
• Daiki Adachi, Tsuyoshi Okubo, Synge Todo, Bond-weighted Tensor Renormalization Group, preprint: arXiv:2011.01679

可積分系は物理と数学において豊富な知見を提供してきた。特に古典可積分系は古くから研究されながらも、今なお多くの新たな発見が報告されている。これまで、古典可積分系は、偶然あるいはひらめきによって、あるいはある種の仮設の範囲でのみ発見・構成されてきた。我々はこれに対し、深層学習を用いることで、これまでにない新たな構成法を提案した。具体的には、「作用角変数」と呼ばれる特殊な正準座標で与えられたハミルトニアンを与え、それに対する自然なハミルトニアンを深層学習によって発見する。この手法では特に、(i)~ニューラルネットで全単射な座標変換を表現できる。(ii)~時間発展の自動微分を随伴法により効率化できる点に着目した。実際に古典可積分系の代表である戸田格子に対して我々の手法を適用し、戸田格子のポテンシャルを見つけ出すことができることを示した

• Fumihiro Ishikawa, Hidemaro Suwa, Synge Todo, Neural Network Approach to Construction of Classical Integrable Systems, preprint: arXiv:2103.00372.