Bose-Hubbard模型时间演化模拟异常求助

• 基于同一表象完成了Bose-Hubbard模型的哈密顿量、初始态、粒子数算符矩阵构建，理论逻辑自洽。
• 5格点晶格系统中，2个玻色子初始分布为01000和00010，按理论预期应呈现对称演化，但实际模拟出现单侧粒子数增多的不对称异常。
• 自行排查数日未定位问题，怀疑是矩阵与向量相乘环节存在代数错误，相关代码来自Bose-Hubbard项目下的2 particles new.py文件。

内容的提问来源于stack exchange，提问作者Mikel Puy