[論文レビュー] Efficient and scalable Path Integral Monte Carlo Simulations with worm-type updates for Bose-Hubbard and XXZ models

Nicolas Sadoune, Lode Pollet|arXiv (Cornell University)|Apr 26, 2022

Physics of Superconductivity and Magnetism被引用数 4

ひとこと要約

本論文は、任意の格子上でボーズ＝ハッバード模型およびXXZ模型の経路積分モンテカルロシミュレーションを効率的かつスケーラブルに行うための、オープンソースでC++14準拠のワームアルゴリズムの実装を提示する。系のサイズおよび逆温度に対してほぼ線形のスケーリングを達成し、スーパーフルイド相では自己相関時間が1に近い。これにより、強い相関を持つ量子系における臨界現象の正確なシミュレーションが可能になる。

ABSTRACT

We present a novel and open-source implementation of the worm algorithm, which is an algorithm to simulate Bose-Hubbard and sign-positive spin models using a path integral representation of the partition function. The code can deal with arbitrary lattice structures and assumes spin-exchange terms, or bosonic hopping amplitudes, between nearest-neighbor sites, and local or nearest-neighbor interactions of the density-density type. We explicitly demonstrate the near-linear scaling of the algorithm with respect to the system volume and the inverse temperature and analyze the autocorrelation times in the vicinity of a U(1) second order phase transition. The code is written in such a way that extensions to other lattice models as well as closely-related sign-positive models can be done straightforwardly on top of the provided framework.

研究の動機と目的

• 符号が正の量子モデルのスケーラブルでオープンソースのワームアルゴリズムの実装を開発すること。
• 最近接相互作用と密度-密度相互作用を伴う任意の格子構造でのシミュレーションを可能にすること。
• ボーズ＝ハッバード模型およびXXZ模型において、系のサイズおよび逆温度に対してほぼ線形のスケーリングを達成し、臨界領域の研究を効率的に行えるようにすること。
• 将来の他の格子模型や関連する符号が正の系への拡張を容易にするモジュラーで拡張可能なフレームワークを提供すること。

提案手法

• アルゴリズムはフォック状態を用いた経路積分表現と、運動エネルギーにおける強い結合項展開を用いる。
• 単一粒子グリーン関数をサンプリングするために、INSERTWORM-GLUEWORMおよびMOVEWORM更新を含むワーム型更新スキームを採用する。
• データ構造は、多数体場理論における時間順序を保つために、時間順に並べたホッピングイベントを格納するC++のvector of listsを用いる。
• 実装はALPSCoreライブラリに依存しており、静的格子定義とXMLベースの格子入力を両方サポートすることで柔軟性を確保する。
• 粒子数、エネルギー、密度-密度相関関数、ウィンディング数、虚時間グリーン関数といった主要な物理量は、それぞれ専用の推定器によって計算される。
• コードは高性能を実現しており、メモリ使用量はホッピングイベントの数に比例し、単一ノードおよびスーパーコンピュータでも効率的にスケーリングする。

実験結果

リサーチクエスチョン

• RQ1 ボーズ＝ハッバード模型およびXXZ模型に対して、系のサイズおよび逆温度に対してほぼ線形のスケーリングを達成するワームアルゴリズムの実装方法は何か？
• RQ2 U(1)の第二種相転移付近における自己相関時間の挙動はどのように変化するか？また、系のサイズおよび逆温度にどのように依存するか？
• RQ3 特に臨界領域において、小規模から大規模な系へのスケーリングにおいても、性能とメモリ効率を維持できるか？
• RQ4 本実装は、以前のワームアルゴリズムコードと比較して、保守性、拡張性、パフォーマンスの観点でどのように優れているか？

主な発見

• アルゴリズムは系の体積および逆温度に対してほぼ線形のスケーリングを示し、大規模なシミュレーションを可能にする。
• スーパーフルイド相では自己相関時間が1に近いまま維持され、ワームが長距離相関を効率的に探索できることを示している。
• メモリ使用量はホッピングイベントの数に比例し、中程度のβにおける96×96系では約70 MBで、フラクチュエーションの分で倍増する。
• コードは高いパフォーマンスを維持しており、3.1 GHzのi5プロセッサを搭載した単一iMacノードで1秒間に2×10^7回以上の更新が可能である。
• 実装は非常に拡張可能であり、静的定義またはXMLベースの定義により任意の格子をサポートし、新しいモデルへの容易な適応が可能である。
• ボーズ＝ハッバード模型およびXXZ模型の自動テストにより、正確性が確認され、正確な対角化結果と一致している。