Finite-state automata for exact matrix product operators of tight-binding Hamiltonians: fractals, quasicrystals, trees and hyperbolic lattices

物質中を移動する電子の挙動を理解しようとしている場面を想像してみてください。最も単純なモデルでは、科学者たちは電子が、地図上の都市間を移動する旅行者のように、格子状に隣の原子へと飛び移っていく様子を描きます。この「タイトバインディング（強結合）」モデルは、電気がどのように流れるか、物質がどのように熱を伝達するか、そしてなぜ一部の物質が磁性を持つのかを理解するための基本的なツールです。数十年にわたり、このモデルは、タイル張りの床のような規則的で繰り返されるパターンに対して最もよく機能してきました。しかし、自然界はこれほど整然としていることは稀です。電子は、奇妙で、繰り返されているものの決して全く同じではないパターンを持つ物質や、海岸線や雪の結晶のように、あらゆるスケールで自己相似的に繰り返される形状の中も移動します。これらの複雑な幾一、すなわちフラクタル、準結晶、樹木状の構造は、通常の格子では説明できないエキゾチックな電子挙動を生み出します。しかし、これらを研究することはコンピュータにとって悪夢でした。シミュレーションにおける原子の数が増えるにつれ、電子が通り得るあらゆる経路をマッピングするために必要なメモリが爆発的に増加し、最も強力なスーパーコンピュータでさえもすぐに限界に達してしまうのです。

研究チームは今、この複雑さを制御する方法を見出しました。彼らは、コンピュータが言語を処理する方法から借りてきた巧妙なトリックを用いて、これらの複雑な電子マップを記述する新しい手法を開発しました。原子間の膨大な詳細リストをすべて保存しようとする代わりに、彼らは材料の構造全体を極めて小さなデジタル空間に捉える、ルールに基づいたコンパクトなシステムを構築しました。このシステムにより、彼らは数十億個の原子を持つ材料をシミュレートすることが可能になりました。これは、以前は厳密な計算では到達不可能だったスケールです。この新しいツールと、エネルギー準位を分析するための標準的な手法を組み合わせることで、彼らはフラクタルの自己相似的なパターンから、双曲空間の曲がった無限の広がりまで、さまざまな環境における電子の挙動を高い精度で算出することができました。

この画期的な成果の核心は、研究者が材料の形状に関する情報をどのように整理したかにあります。標準的なコンピュータシミュレーションでは、システムを記述するために必要なメモリは通常、原子の数に正比例して増加します。材料のサイズを2倍にすれば、必要なメモリもおよそ2倍になります。しかし、ここで研究対象となっている複雑な形状の場合、必要なメモリは通常あまりにも急速に増大するため、大規模なシステムのシミュレーションには地球上に存在するすべてのストレージ量が必要になってしまいます。研究者たちは、これらの特定の形状——フラクタルの、準結晶、および樹木状の構造——が、再帰的なルールを用いて構築されていることに気づきました。つまり、構造全体が、予測可能な方法で積み重ねられた、より小さな自身のコピーによって構成されているのです。彼らはこれらのルールを、通過する車に応じて状態を変化させる信号機のような、シンプルな機械へとエンコードしました。ただし、この場合の機械は、材料の構築レイヤーを移動する電子の位置を追跡します。

著者たちが「有限状態オートマトン」と呼ぶこの機械は、設計図（ブループント）として機能します。それはすべての原子をリストアップするのではなく、電子が構造の次のレベルへと移動する際に、コンピュータに原子間の接続を生成させるための一連の指示を保持しています。これらの形状を構築するためのルールは有限かつ反復的であるため、この機械は、システムがいかに大きくなっても、構造全体を記述するために固定された少数の内部状態しか必要としません。研究者たちは、この機械を「行列積演算子（matrix product operator）」として知られる数学的構造へと翻訳しました。平たく言えば、これは電子の移動のルールを、小さく連結されたブロックの連鎖として書き出す方法です。これらのブロックのサイズは、形状のルールの複雑さのみに依存し、総原子数には依存しません。これは、たとえシミュレーションに数十億のサイトが含まれていても、コンピュータが電子の移動や相互作用を規定する数学的対象であるハミルトニアンを記述するために、管理可能なデータ量のみを保存すればよいことを意味します。

彼らはこの手法をテストするために、幅広い困難な幾何学形状に適用しました。まず、手法が正しく機能することを確認するために、馴染みのある正方形や三角形の格子から始め、次に、正方形を9つの小さな正方形に繰り返し分割し、中央の1つを取り除いた有名なフラクタル形状であるシェルピンスキー・カーペットへと進みました。また、決して繰り返さないものの厳格なルールに従う一次元のフィボナッチ鎖、およびその二次元的な拡張についてもシミュレーションを行いました。最後に、中心から外側へと枝分かれしていく構造であるケイリー木、および中心から離れるにつれて接続の数が指数関数的に増加する、鞍型のように内側に湾曲した双曲格子に取り組みました。あらゆるケースにおいて、この手法は、すべての接続の全リストを書き出すことなく、システムのルールに関する正確な記述を生み出しました。

ルールがエンコードされると、研究者たちは「カーネル多項式法（kernel polynomial method）」と呼ばれる手法を用いて、材料のスペクトル特性を計算しました。このアプローチにより、システム全体を一度に解くことなく、電子が利用可能なエネルギー準位や、材料全体における電子の分布を決定することができます。シェルピンスキー・カーペットについては、20万サイトを超えるシステムをシミュレートし、電子波が構造全体に広がっており、一箇所に留まる兆候がないことを発見しました。また、20万近いサイトを持つフィボナッチ鎖についても観察を行い、準結晶に特有の、電子密度が複雑で非反復的な方法で変化する、独特で複雑なエネルギーパターンを確認しました。双曲格子については、n = 10の世代の3-Cayley木としてモデル化し、材料の無秩序さが電子の流れにどのように影響するかをマッピングし、構造の対称性が混沌とした環境下でも維持されていることを示しました。

この研究の意義は、これまでシミュレーションが極めて困難であった幅広いクラスの材料を研究するための統一的な方法を提供したことにあります。以前は、研究者は近似法に頼るか、あるいは非常に小さな規模のシステムしか研究できず、それによって重要な大規模挙動を見逃してしまうことがありました。この新手法を用いることで、以前は不可能だった指数関数的に大きなスケールで、これらの複雑な格子の電子特性を探求できるようになりました。この手法が特に強力なのは、それが「厳密（exact）」であるためです。情報を捨てたり、システムの挙動について推測したりすることはありません。単に、ルールを記述するためのより効率的な方法を見つけ出したのです。これにより、光、回路、原子配列を用いた新しいタイプの量子シミュレータの構築においてますます重要となっている、非ユークリッド幾何学における電子の挙動を探求する道が開かれます。

研究者たちはまた、彼らの手法は研究した形状については完璧に機能するものの、それらの形状が単純で繰り返されるルールから構築されているという事実に依存していることも指摘しています。もし材料がそのような明確なパターンに従わない構造を持っている場合、この特定のアプローチはそれほどうまく機能しない可能性があります。しかし、再帰的な構造を持つ広大な家族に対しては、この手法は強力な新しいレンズを提供します。これにより、科学者たちは、数十億の原子を持つ材料の電子特性について、実世界のデバイスの理解や、特性をカスタマイズした新材料の設計に関連するスケールで問いを投げかけることができるようになります。これらのシステムを厳密にシミュレートできるということは、コンピュータによる予測が物理法則そのものと同様に信頼できることを意味し、量子物質の領域における将来の実験と発見のための強固な基礎を提供するのです。