[論文レビュー] BIGSTICK: A flexible configuration-interaction shell-model code (updated)

Calvin W. Johnson, W. E. Ormand+3 more|arXiv (Cornell University)|Jan 24, 2018

Parallel Computing and Optimization Techniques参考文献 28被引用数 41

ひとこと要約

BIGSTICKは、オンザフライのマトリクス要素計算、MPI/OpenMP並列処理、経験的およびab initio（ノーコア）空間をサポートする柔軟なオープンソースの配置相互作用シェルモデルコードです。

ABSTRACT

We present BIGSTICK, a flexible configuration-interaction open-source shell-model code for the many-fermion problem. Written mostly in Fortran 90 with some later extensions, BIGSTICK utilizes a factorized on-the-fly algorithm for computing many-body matrix elements, and has both MPI (distributed memory) and OpenMP (shared memory) parallelization, and can run on platforms ranging from laptops to the largest parallel supercomputers. It uses a flexible yet efficient many-body truncation scheme, and reads input files in multiple formats, allowing one to tackle both phenomenological (major valence shell space) and ab initio (the so-called no-core shell model) calculations. BIGSTICK can generate energy spectra, static and transition one-body densities, and expectation values of scalar operators. Using the built-in Lanczos algorithm one can compute transition probability distributions and decompose wave functions into components defined by group theory. This manual provides a general guide to compiling and running BIGSTICK, which comes with numerous sample input files, as well as some of the basic theory underlying the code. Updated November 2025 to version 8.0.0

研究の動機と目的

• Phenomenologicalおよびab initioスペースの両方に対応する柔軟なCIシェルモデルツールの開発を動機づける。
• 非常に大きな多体空間を扱うための高速でメモリ効率の良いアルゴリズムを提供する。
• ノート PC からスーパーコンピュータまで、エネルギースペクトル、密度、演算子期待値の計算を実現する。
• アクセス性と信頼性を広げるための入力フォーマット、使用指針、並列化戦略を紹介する。

提案手法

• 設定されたシェルモデル基底における多体問題をマトリクス固有値問題として定式化する。
• 全要素を貯蔵せずにハミルトニアンのマトリクス要素をオンザフライ計算するアルゴリズムを用い、二体（および任意で三体）相互作用を活用する。
• データを削減し演算を加速するため、基底を種別（例：陽子と中性子）と量子数で因子分解する。
• 低エネルギー固有値と固有ベクトルを効率的に得るためにLanczos対角化を採用する。
• 価電子空間およびノーコア計算の双方に対して、複数の基底スキーム（M-scheme、J-scheme、SU(3)-scheme）および縮約戦略（例：Nmax）をサポートする。
• MPIとOpenMPの並列化を提供し、ノートPCからリーダーシップ級スーパーコンピュータまでスケールする。

実験結果

リサーチクエスチョン

• RQ1 ハミルトニアン全体のストレージを持たずに、非常に大きな基底空間へ配置相互作用計算をどうスケールさせるか？
• RQ2 異なる物理問題に対して、どの基底表現（M-scheme、J-scheme、SU(3)-scheme）がメモリ使用と計算効率を最適化するか？
• RQ3 因子化されたオンザフライのマトリクス要素を用いて、Phenomenologicalおよびab initioスペースを横断して低エネルギースペクトルと遷移密度を信頼性高く計算できるか？
• RQ4 縮約スキーム（例：Nmax、n-particle n-hole）が大規模ノーコアおよび価電子空間計算の精度と実現性にどう影響するか？
• RQ5 ノートPCからスーパーコンピュータまで portabilityを確保する性能特性と並列化戦略はどうあるべきか？

主な発見

• BIGSTICKはエネルギー固有値スペクトル、静的および遷移一体密度、スカラー演算子の期待値を生成できる。
• ハミルトニアンを適用するために因子化されたオンザフライアルゴリズムを使用し、すべての非零要素の明示的な貯蔵を回避する。
• Phenomenologicalな主要価電子空間とab initioノーコアシェルモデル計算の両方をサポートし、三体力を含む場合はコストが高くなる。
• Lanczos対角化により、非常に大きな空間で低エネルギー固有対を効率的に抽出できる（プラットフォームに応じて数千万〜十億級の基底状態まで）。
• MPIとOpenMPによる並列化と基底セクター整理および量子数の因子化により、ノートPCからリーダーシップ級機械までスケーラブルな性能を発揮する。
• 基底は複数のスキーム（M-scheme、J-scheme、SU(3)-scheme）で表現でき、それぞれサイズと複雑さのトレードオフがあり、さまざまな核構造問題に合わせたアプローチを可能にする。