Quantum hierarchical Fokker–Planck equations with U(1) gauge fields [U(1)-QHFPE]: A computational framework for Aharonov–Bohm effects

We present a software package that solves the quantum Fokker-Planck equation with gauge fields, formulated within the hierarchical equations of motion framework [U(1)-QHFPE]. The framework rigorously preserves gauge invariance and rotational symmetry under non-Markovian and non-perturbative system-bath (S-B) interactions, enabling accurate simulations of transport phenomena such as the Aharonov-Bohm effect under thermal environments. In a strong S-B coupling regime, quantum S-B entanglement emerges naturally. The demonstration programs perform calculations of response functions in Aharonov-Bohm ring geometries with a mechanical potential that may induce quantum tunneling, thereby illustrating the software's capability to resolve topological quantum interference in dissipative open systems. Written in C++, the code includes a central processing unit version with highly readable OpenMP directives and a graphics processing unit version tailored for high-performance computing.