[论文解读] 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并行，以及对现象学与从头计算(no-core)空间的支持。

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

研究动机与目标

• 推动开发一个灵活的 CI 壳模型工具，涵盖现象学与从头计算（ab initio）空间。
• 提供一个快速、节省内存的算法，以处理极大的多体空间。
• 让能量谱、密度及算符本征值的期望值计算在从笔记本电脑到超级计算机的多平台上可行。
• 引入输入格式、使用指南和并行化策略，以提升可访问性与可靠性。

提出的方法

• 将多体问题在配置的壳模型基中表述为矩阵特征值问题。
• 使用即时算法计算哈密顿量矩阵元素，避免完全存储，利用二体（以及可选的三体）相互作用。
• 按粒子种类（如质子与中子）和量子数对基底进行分解，以减少数据并加速运算。
• 采用 Lanczos 对角化高效获得低激发的特征值和特征向量。
• 支持多种基底方案（M-scheme、J-scheme、SU(3)-scheme）和截断（如 Nmax），用于价空间与无核计算。
• 提供 MPI 和 OpenMP 并行化，使规模从笔记本电脑到领军级超级计算机可扩展。

实验结果

研究问题

• RQ1 如何在不完全存储哈密顿量的情况下，将配置相互作用壳模型计算扩展到极其大的基空间？
• RQ2 不同物理问题下，哪种基表示（M-scheme、J-scheme、SU(3)-scheme）能优化内存使用和计算效率？
• RQ3 是否能使用因式分解的即时矩阵元素，可靠地在现象学与从头计算空间中计算低-激发谱和跃迁密度？
• RQ4 截断方案（如 Nmax、n-particle n-hole）对大规模无核与价空间计算的准确性和可行性有何影响？
• RQ5 哪些性能特征和并行化策略能够确保从笔记本电脑到领军级机器的可移植性？

主要发现

• BIGSTICK 能生成能量谱、静态与跃迁一体密度，以及标量算符的期望值。
• 该代码使用分解的即时算法来作用哈密顿量，避免显式存储所有非零元素。
• 它同时支持现象学的主价空间和从头无核壳模型计算，包括三体力（成本较高）。
• Lanczos 对角化使在极大空间中高效提取低-激发特征对成为可能（具体平台下可达数千万到数十亿基态状态）。
• 通过 MPI 和 OpenMP 的并行化，加上基 Sector 组织和量子数分解，获得从笔记本到领军级机器的可扩展性能。
• 基底可以用多种方案表示（M-scheme、J-scheme、SU(3)-scheme），每种在规模和复杂性上有权衡，便于针对不同核结构问题采取定制化方法。