藤堂研究室へようこそ

シミュレーションで探る量子多体現象

物質の状態を知るには、多体のシュレディンガー方程式を解き、統計力学の分配関数を求めればよい。しかしながら、現代のスーパーコンピュータの計算能力をもってしても、完全な解を求めることはできない。そこで、もとの方程式の中に含まれる、物理的に重要な性質を失うことなく、シミュレーションを実行しやすい形へ表現しなおすことが、計算物理における重要な鍵となる。藤堂研究室では、モンテカルロ法に代表される確率的なシミュレーション、経路積分に基づく量子ゆらぎの表現、特異値分解やテンソルネットワークによる情報圧縮、統計的機械学習の手法などを駆使し、古典/量子スピン系やボーズハバード系などに代表される強相関多体系における新奇な状態や相転移現象の探索・解明を目指している。また、最先端のスーパーコンピュータの能力を活かすための並列化手法の研究、次世代シミュレーションのためのオープンソースソフトウェアの開発・公開も進めている。

セミナー


最近の研究より

古典調和振動子模型の非エルゴード性

ergodicity.png記憶を含んだランジュバン方程式である一般化ランジュバン方程式は、通常のランジュバン方程式と異なり、その記憶効果によって特異な拡散現象を示すことが知られている。近年、異常拡散の解析で活発に利用され、またこの異常拡散の示す非エルゴード性が注目を集めている。しかしながら、従来の研究ではエルゴード性の定義について混乱があり、また具体的な物理的モデルを用いた解析は少なく、その物理的意味は明らかではなかった。我々は、非エルゴード性を示す古典調和振動子系を提案し、系における注目粒子の振る舞いを、分子動力学計算、厳密対角化、解析計算を用いて調べた。また同時に、エルゴード性の定義について再考した。その結果、本モデルにおける非エルゴード性の起源は注目粒子の周りに励起される孤立局在モードであることが明らかとなった。

共振器系での動的な協力現象

bistabilityレーザー照射下で現れる光双安定性は、レーザー強度に対して共振器中のフォトン数が双安定な状態を示し、またその間を不連続に跳ぶといった一次相転移現象に似た振る舞いの現れる相転移現象である。このような非平衡系で現れる動的な協力現象を、大規模な数値計算によって量子力学的な微視的模型から解析を行った。具体的には多数のフォトンと多数の二準位原子からなる量子マスター方程式を近似なく解く並列計算を実行した。量子マスター方程式の時間発展演算子の固有値・固有状態から、定常状態でのフォトン数分布関数のサイズ依存性、および緩和時間のサイズ依存性を調べ、平衡系での一次相転移に対応する結果を得た。従来の研究と比べフォトン数密度の小さな領域では、準安定状態のレーザー周波数依存性が定性的に異なることを明らかにした。また、レーザー強度を時間周期的に変調させることで、その周期に対し動的な相転移現象が現れることを明らかにした。

実験データと第一原理計算を組み合わせた結晶構造決定

xray結晶構造推定は非常に難しい問題として古くから知られており、様々な推定方法が開発されてきた。特に最近では、実験データがある場合にはエネルギーの最適化と実験データの再現を同時に行うことによって、結晶構造推定の成功率を上げられることが知られてきている。その際に用いられる方法は、実験データの再現度とエネルギーを足し合わせた新しい評価関数を用いて最適化を行うと言った方法である。しかしながら、2つの評価関数を 足し合わせてしまっているため、それぞれの評価関数の情報が失われてしまうといった欠点がある。我々は2つの評価関数の同時最適点を探る方法として重ね合わせ最適化法を開 発し、その性能について調査を行なっている。例えばSiO2系は、エネルギーの局所最 適点が多く存在するため結晶構造を決定するのが難しい系であるが、我々の方法を用いることで結晶構造の推定精度が大幅に上昇することことを確認した。

連続空間における経路積分モンテカルロ法

pimc.pngHe4の2次元系は低温で並進対称性とゲージ対称性が同時に破れ、超固体と呼ばれる量子相へ転移することが予想されている。有限温度での2次元系He4の平衡分布を調べ超固体の予想を裏付けるため、連続空間における経路積分モンテカルロ法の改良を進めている。Event-chainモンテカルロの手法と、worm algorithm の手法を援用し、詳細釣り合い条件を破るアルゴリズムを開発した。このアルゴリズムを相互作用のない理想Bose粒子に適用し、従来の方法より分布の収束が速く、トロッター数に対する時間計算量が改善した結果を得た。

非一様な系に対するテンソルネットワーク繰り込み

テンソルネットワークは、厳密対角化といった従来の指数的にコストがかかる計算手法に比べ、ベキ的な計算量のコストですむ計算手法として近年注目されている。代表的な手法にTRGや、HOTRGといったものがあるが、それ以外にProjectorによる手法も注目されている。これまでテンソルネットワークは一様な系に対して用いられてきたが、我々は、HOTRGとProjectorによる計算手法を非一様系に拡張し、ボンド希釈を含むイジングモデルに適用し、非一様系においてProjectorによる手法が収束が速いことを確認した