ScatterWorks: A Python package for building and solving scattering network models

技術概要：ScatterWorks

問題提起
ネットワークモデルは、量子ホール遷移のためのチャーカー・コドントン（Chalker–Coddington）モデル、モアレ系、およびフロケ問題といった、無秩序およびトポロジカルな系における非相互作用の輸送および波の伝搬現象を研究するための効率的な枠組みを提供します。その汎用性にもかかわらず、これらのモデルは、凝縮系物理学の一部において、連続体近似やタイトバインディング的手法よりも広く採用されているとは言えません。この利用が進まない背景には、汎用ツールの欠如があります。研究者はしばしば、再利用が困難な、個別の問題に特化したプライベートなコードに依存しています。さらに、標準的な転送行列法は、指数関数的に増大または減衰するモードによって数値的不安定性を生じる可能性があります。幾何学的形状、境界条件、または無秩序プロトコルの変更を行うたびに、単にパラメータを調整するだけでなく、コアとなるセットアップ論理を書き直す必要が生じます。したがって、記号的（シンボリック）および数値的なネットワークモデルの両方を、柔軟に構築し、堅牢に解くことができる、統一されたオープンソースの環境が求められています。特に、高密度な線形代数が計算上の限界となるような大規模なシステムにおいて、その必要性は顕著です。

手法
本論文では、散乱ネットワークモデルを構築し解くために設計された、オープンソースのPythonパッケージであるScatterWorksを紹介します。その手法は、コンパクトなネットワーク表現とモジュール化されたワークフローを中心としています。

1. ネットワーク表現: 本パッケージは、標準的なリンク（canonical_links）のリストによって定義されるNetworkオブジェクトを使用します。各リンクは、ソースノードからシンクノードへと接続され、周期境界条件やタイリングを処理するためのオプションのシフトベクトルを含みます。この表現により、ネットワークのトポロジカルな定義と、散乱行列の割り当てが分離されます。
2. 構成可能な操作: パッケージは、散乱行列の基底を保持したまま、標準的なリンクの順序を変更せずにネットワークの幾何学的形状を操作するための一連の操作を提供します。主な操作は以下の通りです：
   • tile(): ユニットセルを複製して、より大きな周期的または有限のネットワークを作成します。
   • union(): 離れたネットワークを結合します。
   • cut(): 特定のリンクまたはノードのサブセットを選択します。
   • relink(): link_map配列を修正することでリンクの端点を再接続し、標準的な基底を再定義することなく複雑な幾何学的形状の変更を可能にします。
3. 散乱方程式のコンパイル: ローカルな散乱行列は、ラベル（例：ノードファミリー、サブ格子）または明示的なインデックスに基づいてノードに割り当てられます。これらのローカル行列は、グローバルな散乱演算子 $S$ へと組み立てられます。パッケージは、解析的な導出のための記号的行列（SymPyを使用）と、シミュレーションのための数値的行列（密または疎）の両方をサポートしています。
4. ソルバーと観測量: 輸送観測量を計算するために、パッケージはバルクリンクとリードリンク（入射および放出）を区別します。バルクリンクを除去することで、リード間の有効散乱行列（$S_{eff}$）を導出します。
   • 一般的なケースでは、solve_scattering_equationsが有効行列を計算します。
   • 大規模な疎行列系の場合、schur_solveはSchur補完法（MUMPSライブラリ経由）を利用して、バルク部分ブロック $(I - S_{bb})^{-1}$ の逆行列計算を効率的に処理します。このアプローチは、転送行列の数値的不安定性を回避し、大規模なネットワークに対して優れたスケーラビリティを実現します。
   • 輸送観測量、例えば無次元コンダクタンス $g$ は、$S_{eff}$ の透過ブロックに適用されたランダウアー公式を用いて計算されます。
5. 可視化: パッケージには、ラベルを認識したホバー機能および静的ラベルを備えたネットワーク幾何学をレンダリングするプロットユーティリティが含まれており、複雑なネットワーク構築のデバッグを支援します。

主な貢献

• 統一されたフレームワーク: ScatterWorksは、記号的および数値的なネットワークモデルの両方を構築するための単一のインターフェースを提供し、解析的な導出と大規模な数値シミュレーションの間の溝を埋めます。
• 柔軟な幾何学構築: 構成可能な操作（tile, cut, relink, union）を用いることで、ユーザーはコアとなる論理を書き直すことなく、複雑なユニットセル、リボン、および閉じた幾何学的形状を定義できます。
• 疎行列ソルバー: 疎な散乱行列のためのSchurベースのソルバーを実装することで、パッケージは密な手法では実用的ではない大規模なネットワークサイズでの研究を可能にします。これは、ネットワークモデルのシミュレーションにおける主要なボトルネックに対処するものです。
• 記号的機能: 記号的な散乱行列を割り当てる能力により、ネットワークモデルのコードから直接、解析的な式を復元することが可能です。
• ラベルベースのワークフロー: ノードおよびリンクのラベルにpandas DataFrameを使用することで、散乱行列の割り当てやリードの選択が容易になり、カスタム幾何学の設定が簡素化されます。

結果
本論文では、以下の2つの主要な例を通じてパッケージの能力を実証しています。

1. 記号的ファブリ・ペロー干渉計: 本パッケージを使用して、ミラーと位相ノードを表す記号的な散乱行列を持つ小さなネットワークを組み立てました。得られた有効散乱行列は解析的に導出され、往復の共鳴分母を含む、標準的なファブリ・ペローの透過および反射構造を再現しました。
2. 数値的チャーカー・コドントン・モデル: 量子ホール相転移を研究するために、Chalker–Coddingtonネットワークの数値的実装を構築しました。
   • パッケージは準エネルギーバンド構造を正常に再現し、固有位相が交差する臨界点 $\theta = \pi/4$ を特定しました。
   • 臨界点付近の様々なシステムサイズ ($L$) におけるコンダクタンス計算は、期待されるピーク値である約 $0.5 e^2/h$ を示しました。
   • 実行時間の分析により、Schurベースのソルバーがシステムサイズに対して期待通りの $O(L^3)$ スケーリングを示すことが確認され、密な手法と比較した際の効率性が検証されました。

意義
本論文は、ScatterWorksをネットワークモデル研究への参入障壁を下げるための実用的なツールとして位置付けています。その意義は、疎な定式化を通じてパフォーマンスを維持しながら、柔軟なオープンソース環境を提供することにあります。ネットワーク幾何学の構築と散乱行列の割り当て、および解法を切り離すことにより、本パッケージを使用することで、研究者はコードを再設計することなく、複雑な幾何学的形状や無秩序プロトコルを探索できます。プロットユーティリティとラベル認識レンダリングの導入は、大規模でカスタムなネットワーク設定のデバッグと再現性をさらにサポートします。著者らは、本パッケージが拡張可能であり、今後の課題として、ユーザーインターフェースの改善や、フロケ駆動物理学やランダムウォーク研究などの他の領域への機能拡張を挙げていると述べています。